如何将中国区域外面的值去掉

秋风@秋雨

我利用matlab将中国区域的年气温平均图画出来，但是由于插值的原因，在中国以外区域也出现值，不知如何去掉，菜鸟求帮忙

whmooncool

下载中国的basemap

whmooncool

貌似家园有个帖子说了，你去查查

Aires

楼上说的是，楼主搜一下帖子

金陵蔡肉包

有了shapefile，然后有个inpolygon函数可以干这个。

royalosyin

做一个cookie文件，就是0，1的那种，1代表中国。
然后Data（cookie==0）= nan 就ok了。

huzengyun

看来matlab还是很强大的，还是需要好好学习！

小寶

金陵蔡肉包 发表于 2013-10-22 00:25

                               
登录/注册后可看大图

有了shapefile，然后有个inpolygon函数可以干这个。

你好我想问下这个inpolygon怎么用呢

1. %中国地图附加南海地图
   
2. clc;clear
   
3. 
   
4. s='F:\studyData\MODIS\Terra\L3MonthyAerosol\2004\001\MOD08_M3.A2004001.051.2010284170915.hdf';
   
5. l333=hdfread(s,'Optical_Depth_Land_And_Ocean_Mean_Mean');
   
6.             [alon,alat]=meshgrid(72:1:137,16:1:55);
   
7. guojie=shaperead('bou1_4l.shp');%读取国界shp文件的内容
   
8. bou1_4lx=[guojie(:).X];%提取经度信息
   
9. bou1_4ly=[guojie(:).Y];%提取纬度信息
   
10. l1=double(l333(72:137,16:55))';
    
11. %%
    
12. m_proj('miller','lon',[72,137],'lat',[16,55])%选择投影方式
    
13. m_pcolor(alon,alat,l1);
    
14. hold on;
    
15. m_plot(bou1_4lx,bou1_4ly)%绘国界
    
16. 
    
17. inpolygon 这个地方要怎么加 在哪里加？
    
18. 
    
19. m_grid('linestyle','none','box','fancy','tickdir','out','fontsize',16)%添加坐标
    
20. h1=axes('Position',[0.14 0.155 0.09 0.2]);%创建坐标系时返回它的句柄
    
21. axes(h1);%将坐标系h1置为当前坐标系
    
22. %set(gcf,'PaperPositionMode','auto')
    
23. %m_proj('miller','lon',[106,122],'lat',[2,26]) %设置南海区域
    
24. m_plot(bou1_4lx,bou1_4ly)%绘国界
    
25. m_grid('XTick',[],'YTick',[])%添加坐标
    

复制代码

金陵蔡肉包

小寶 发表于 2013-10-25 10:26

                               
登录/注册后可看大图

你好我想问下这个inpolygon怎么用呢

你找个中国国境线的shapefile，就是不要有省区划的。然后用shaperead函数读进去。再
map_path = shaperead('')   %引号中间写shapefile的路径。
map_X = [map_path(:).X];  
map_Y = [map_path(:).Y];  
IN(:,:)=inpolygon(lon(:,:),lat(:,:),map_X,map_Y);

小寶

金陵蔡肉包 发表于 2013-10-26 00:54

                               
登录/注册后可看大图

你找个中国国境线的shapefile，就是不要有省区划的。然后用shaperead函数读进去。再
map_path = shapere ...

谢谢你在论坛其他帖子下到中国的shp文件了，我试试看看能不能弄明白，不会在问你哦好心人谢啦