Matlab M_map工具箱绘制Interrupted Mollweide Projection

GMT自带了许多的地图投影，但是对于Interrupted Mollweide投影效果却不好。

作为平替的m_map工具箱中带有的投影类型可完美解决这一问题。

Interrupted Mollweide Projection长这样

全球陆地

全球海洋

使用Matlab工具箱m_map展示全球海平面变化的空间分布

addpath('/Users/chistrongwen/Desktop/m_map')

Slongs=[-100 43;-75 20; 20 145;43 100;145 295;100 295];

Slats= [  0  90;-90  0;-90   0; 0  90;-90   0;  0  90];

x = -179.875:0.25:179.875;

y = -89.875:0.25:89.875;

[x,y] = meshgrid(x,y);

for l=1:6

    m_proj('mollweide','long',Slongs(l,:),'lat',Slats(l,:));

    m_grid('fontsize',6,'xticklabels',[],'xtick',[-180:30:360],...

                'ytick',[-80:20:80],'yticklabels',[],'linest','-','color','k')

    m_coast('patch',[.9 .9 .9],'edgecolor','black');

    hold on

    m_pcolor(aviso_msla_gia.lon,aviso_msla_gia.lat,aviso_msla_gia.rg(:,:,1));

    m_pcolor(x,y,[aviso_msla_gia.rg(:,721:1440,1),aviso_msla_gia.rg(:,1:720,1)]);

    colormap jet

end

xlabel('Interrupted Mollweide Projection of World Oceans');

set(gca,'xlimmode','auto','ylimmode','auto');

成果图

️欢迎点赞收藏️