- 积分
- 8029
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2017-9-18
- 最后登录
- 1970-1-1
|
发表于 2018-5-11 11:27:52
|
显示全部楼层
本帖最后由 callmezoro 于 2018-5-11 11:29 编辑
figure
% m_proj('stereographic','lat',-90,'radius',70,'rotangle',360)
% m_proj('Equidistant','lat',[min(lat),max(lat)],'lon',[min(lon),max(lon)])
m_proj('Lambert','lat',[min(lat),max(lat)],'lon',[min(lon),max(lon)])
[lon,lat]=meshgrid(lon,lat);
m_pcolor(lon,lat,new)
shading interp
hold on
[C,H]=m_contour(lon,lat,new,'k');
clabel(C,H,'LabelSpacing',200,'Rotation',0,'FontSize',8);
clear C H
% coastline
m_coast('patch',[.9 .9 .9]);
m_grid('box','fancy','xaxislocation','top');
% xlabel('longitude');ylabel('latitude');
title(['EOF',num2str(num)])
caxis([-0.6 0.61]);
colorbar;
% load('lightRWB','lightRWB')
% colormap(gca,lightRWB)
colormap(jet)
m_map工具箱了解一下,这是我用来画图的程序,new就是你的网格点,lon,lat不用说了吧,把colorbar范围改一下就行,如果你觉得我填充的陆地颜色不好看也可以随意改,[.9 .9 .9]那个地方最上面提供了三种投影方式,你可以试试哪种比较好看,不过第一种我没写成普适性的你自己看看m_map的介绍改一改就行
|
|