matlab 填充陆地的m_map工具包的调用代码（例子）

Hua_Zhang

matlab 填充陆地的m_map工具包的调用代码（例子）m_map 工具包（太大了，无法上传），只能自行下载了，数据为处理后PDO-SOI数据，本贴主要提供调用m_map工具包代码

m_map下载地址http://www.eos.ubc.ca/~rich/map.html

1. %% soi
   
2. clear all; close all; clc;
   
3. load('D:\file_data\课件_秋季\气候统计方法和应用\气候统计方法与应用_homework\气候统计2\PDO_SOI.mat');
   
4. load('b2r.txt');b2r=b2r/256;
   
5. soi_r(soi_r==0)=nan;
   
6. soi_r(soi_r>=5 | soi_r<=-5)=nan;
   
7. lon_w = lon(1:180);lon_e = lon(181:360);
   
8. lon_shape = [lon_e;lon_w+360];
   
9. soi_w = soi_r(1:180,:);soi_e = soi_r(181:360,:);
   
10. soi_shape = [soi_e;soi_w];
    
11. [x,y]=meshgrid(lon_shape,lat);
    
12. 
    
13. figure
    
14. m_proj('Equidistant Cylindrical','lat',[-30 60],'long',[60 300]);
    
15. m_coast('patch',[.7 .7 .7],'edgecolor','none');%f\h\l\iq
    
16. m_grid('linest','none','linewidth',2,'tickdir','in');
    
17. hold on
    
18. %caxis([-.01 0.45]);
    
19. m_contourf(x,y,pdo_shape',66,'linecolor','none');
    
20. % m_pcolor(x,y,soi_shape');
    
21. % shading interp
    
22. colormap((b2r));
    
23. colorbar
    
24. caxis([-1 1]);
    
25. m_coast('patch',[.7 .7 .7],'edgecolor','none');%f\h\l\iq
    
26. xlabel('Longitude(^OE)','fontsize',12);
    
27. ylabel('Latitude(^ON)','fontsize',12);
    
28. title('SOI-SSTA','fontsize',15);

复制代码

si18

请问 边界处的空白怎样处理？