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');
成果图
❤️欢迎点赞收藏❤️
)
