- 积分
- 237
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2018-10-24
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
matlab 填充陆地的m_map工具包的调用代码(例子)m_map 工具包(太大了,无法上传),只能自行下载了,数据为处理后PDO-SOI数据,本贴主要提供调用m_map工具包代码
m_map下载地址http://www.eos.ubc.ca/~rich/map.html
- %% soi
- clear all; close all; clc;
- load('D:\file_data\课件_秋季\气候统计方法和应用\气候统计方法与应用_homework\气候统计2\PDO_SOI.mat');
- load('b2r.txt');b2r=b2r/256;
- soi_r(soi_r==0)=nan;
- soi_r(soi_r>=5 | soi_r<=-5)=nan;
- lon_w = lon(1:180);lon_e = lon(181:360);
- lon_shape = [lon_e;lon_w+360];
- soi_w = soi_r(1:180,:);soi_e = soi_r(181:360,:);
- soi_shape = [soi_e;soi_w];
- [x,y]=meshgrid(lon_shape,lat);
-
- figure
- m_proj('Equidistant Cylindrical','lat',[-30 60],'long',[60 300]);
- m_coast('patch',[.7 .7 .7],'edgecolor','none');%f\h\l\iq
- m_grid('linest','none','linewidth',2,'tickdir','in');
- hold on
- %caxis([-.01 0.45]);
- m_contourf(x,y,pdo_shape',66,'linecolor','none');
- % m_pcolor(x,y,soi_shape');
- % shading interp
- colormap((b2r));
- colorbar
- caxis([-1 1]);
- m_coast('patch',[.7 .7 .7],'edgecolor','none');%f\h\l\iq
- xlabel('Longitude(^OE)','fontsize',12);
- ylabel('Latitude(^ON)','fontsize',12);
- title('SOI-SSTA','fontsize',15);
复制代码
|
-
|