matlab中地图边界与掩膜(去掉边界外区域)的实现（基于shape文件）

147

找了好久的方法，顶一个

小寶

在直白一点就好了，我还是看不明白撒

斥鷃

小寶 发表于 2013-10-24 18:18

                               
登录/注册后可看大图

在直白一点就好了，我还是看不明白撒

哪里不明白可以留言~

小寶

斥鷃 发表于 2013-10-24 23:19

                               
登录/注册后可看大图

哪里不明白可以留言~

谢谢可能问题比较多，我以前用IDL来画气溶胶数据，不过发现不能只画中国的一些标注也不是很好字体等等都不能改，现在要用MATLAB，问题是你写的可能还是有一点不明白如何可以帮我回到下，谢谢了
[lon lat]=meshgrid([97:0.1:107],[21:0.1:30]); 没什么说的
Z=griddata(gy_locat(:,2),gy_locat(:,3),EOF_used(:,1),lon,lat,'v4');
这个里面的V4是？你这个EOF_used(:,1)算是一维的，如果我有一个二维的如何写？EOF[360,180]
yunnan=shaperead('D:\myitem\arcgis\yunnan.shp');
这个的问题，我在其他文件下了中国的shp应该就可以了吧
R=makerefmat('RasterSize',size(Z),'Lonlim',[97 107],'Latlim',[21 30]);
MASK=vec2mtx(yunnan.Y,yunnan.X,Z,R,'filled');
MASK(find(MASK>1))=nan;
MASK(find(MASK==1))=0;
contourf(lon,lat,Z+MASK,30);
shading flat
colorbar
hold on
plot(yunnan.X,yunnan.Y,'-k','linewidth',3)
hold off
如上这些和你写的一样是不是就可以了出图了？
PS如果用M_MAP可能会比较简单么

小寶

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. 
   
8. l1=double(l333(72:137,16:55))';
   
9. 
   
10. [lon lat]=meshgrid([72:1:137],[16:1:55]);
    
11. %Z=griddata(l1);
    
12. yunnan=shaperead('bou1_4l.shp');
    
13. 
    
14. R=makerefmat('RasterSize',size(l1),'Lonlim',[72 137],'Latlim',[16 55]);
    
15. MASK=vec2mtx(yunnan(:).Y,yunnan(:).X,l1,R,'filled');
    
16. MASK(find(MASK>1))=nan;
    
17. MASK(find(MASK==1))=0;
    
18. 
    
19. contourf(lon,lat,l1+MASK,30);
    
20. shading flat
    
21. colorbar
    
22. hold on
    
23. plot(yunnan.X,yunnan.Y,'-k','linewidth',3)
    
24. hold off

复制代码

Error using vec2mtx (line 89)
Too many input arguments.

Error in ff (line 17)
MASK=vec2mtx(yunnan(:).Y,yunnan(:).X,l1,R,'filled');这个错误哎，不明白了

斥鷃

小寶 发表于 2013-10-25 08:53

                               
登录/注册后可看大图

谢谢可能问题比较多，我以前用IDL来画气溶胶数据，不过发现不能只画中国的一些标注也不是很好字体等等都不 ...

griddata是一个二维插值函数，'v4'是插值算法名，因为我的数据是格点数据，所以需要先插值到平面以后才能使用，如果你的直接是网格数据就不影响了。MASK=vec2mtx(yunnan(:).Y,yunnan(:).X,l1,R,'filled');中yunnan后面不需要加(:)表示引用，这里yunnan是个结构数组，.x是它的属性，去掉这部分你试试吧。另外问一下你的matlab是哪个版本的，版本问题可能也有影响。还有就是m_map没有怎么用过，但是工具箱挺强大的样子，所以也没有什么发言权，会不会简单些就不清楚了。希望能帮到你。

小寶

斥鷃 发表于 2013-10-25 23:24

                               
登录/注册后可看大图

griddata是一个二维插值函数，'v4'是插值算法名，因为我的数据是格点数据，所以需要先插值到平面以后才能 ...

我的是R2011b

是不是上面我写的代码没什么问题？

小寶

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. 
   
8. l1=double(l333(72:137,16:55))';
   
9. 
   
10. [lon lat]=meshgrid([72:1:137],[16:1:55]);
    
11. %Z=griddata(l1);
    
12. yunnan=shaperead('bou1_4l.shp');
    
13. 
    
14. R=makerefmat('RasterSize',size(l1),'Lonlim',[72 137],'Latlim',[16 55]);
    
15. MASK=vec2mtx(yunnan.Y,yunnan.X,l1,R,'filled');
    
16. MASK(find(MASK>1))=nan;
    
17. MASK(find(MASK==1))=0;
    
18. inpolygon
    
19. contourf(lon,lat,l1+MASK,30);
    
20. shading flat
    
21. colorbar
    
22. hold on
    
23. plot(yunnan.X,yunnan.Y,'-k','linewidth',3)
    
24. hold off

复制代码

Error using vec2mtx (line 89)
Too many input arguments.

Error in ff (line 15)
MASK=vec2mtx(yunnan.Y,yunnan.X,l1,R,'filled');

xiaowan

学好matlab出图是不是就可以不用学习surfer出图了呢~~

onion5376

太棒了，一直苦于纠结于用matlab画shp图