手抄报 安全手抄报 手抄报内容 手抄报图片 英语手抄报 清明节手抄报 节约用水手抄报

MATLAB使用sinusoidal和mollweide投影绘制地图

时间:2024-10-12 04:19:41

1、第一,首先使用sinusoidal投影方式绘制分插的世界地图。启动MATALB,新建脚本(觊皱筠桡Ctrl+N),输入如下代码:close al造婷用痃l; clear all; clcslongs=[-100 0;-75 25;-5 45;25 145;45 100;145 295;100 290];slats= [ 8 80;-80 8; 8 80;-80 8; 8 80;-80 0; 0 80];for i=1:7m_proj('sinusoidal','long',slongs(i,:),'lat',slats(i,:));m_grid('fontsize',6,'linest','-','color',[.7,.7,.7],...'xtick',[-180:30:360],'xticklabels',[],...'ytick',[-80:20:80],'yticklabels',[]);m_coast('patch','g');endxlabel('Interrupted sinusoidal projection of world oceans');set(gca,'xlimmode','auto','ylimmode','auto');其中通过for-end循环设置7组不同经度longitude和纬度latitude的绘制区域,所使用的投影方式为sinusoidal。另外,m_proj设置投影方式和经纬度区域。m_grid设置网格线属性。m_coast以填充(补丁patch)的方式设置海岸线。需要注意的时,set(gca,'xlimmode','auto','ylimmode','auto')需要加上,不然不能同时显示7组绘制区域。

MATLAB使用sinusoidal和mollweide投影绘制地图

2、第二,保存和运行上述脚本,得到如下图形:sinusoidal投影方式的分插的(上下共7组)世界地图。

MATLAB使用sinusoidal和mollweide投影绘制地图

3、第三,下面采用mollweide投影方式绘制分插的世界地图。新建脚本(Ctrl+N),输入如下代码:close all; clear all; clcslongs=[-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];for i=1:6m_proj('mollweide','long',slongs(i,:),'lat',slats(i,:));m_grid('fontsize',6,'linest','-','color','k',...'xtick',[-180:30:360],'xticklabels',[],...'ytick',[-80:20:80],'yticklabels',[]);m_coast('patch',[.6,.6,.6]);endxlabel('Interrupted mollweide projection of world oceans');set(gca,'xlimmode','auto','ylimmode','auto');其中通过for-end循环设置6组不同经度longitude和纬度latitude的绘制区域,所使用的投影方式为mollweide。关于代码的其他解释参照第一步。

MATLAB使用sinusoidal和mollweide投影绘制地图

4、第四,保存和运行第三步中的脚本,得到如下图形:mollweide投影方式的分插的(上下共6组)世界地图。

MATLAB使用sinusoidal和mollweide投影绘制地图

5、第五,把第一步和第三步的脚本合并在一起,从而将第二步和第四步中的图形绘制在一个绘图区。合并后的脚本如下,主要是增加了subpl泠贾高框ot( )命令用于在一个绘图区绘制子图。close all; clear all; clcsubplot(2,1,1)slongs=[-100 0;-75 25;-5 45;25 145;45 100;145 295;100 290];slats= [ 8 80;-80 8; 8 80;-80 8; 8 80;-80 0; 0 80];for i=1:7m_proj('sinusoidal','long',slongs(i,:),'lat',slats(i,:));m_grid('fontsize',6,'linest','-','color',[.7,.7,.7],...'xtick',[-180:30:360],'xticklabels',[],...'ytick',[-80:20:80],'yticklabels',[]);m_coast('patch','g');endxlabel('Interrupted sinusoidal projection of world oceans');set(gca,'xlimmode','auto','ylimmode','auto');subplot(2,1,2)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];for i=1:6m_proj('mollweide','long',slongs(i,:),'lat',slats(i,:));m_grid('fontsize',6,'linest','-','color','k',...'xtick',[-180:30:360],'xticklabels',[],...'ytick',[-80:20:80],'yticklabels',[]);m_coast('patch',[.6,.6,.6]);endxlabel('Interrupted mollweide projection of world oceans');set(gca,'xlimmode','auto','ylimmode','auto');

MATLAB使用sinusoidal和mollweide投影绘制地图

6、第六,保存和运行上述合并的脚本,最终得到如下图形:上方为sinusoidal投影方式绘制的分插世界地图,下方为mollweide投影方式绘制的分插世界地图。

MATLAB使用sinusoidal和mollweide投影绘制地图
© 手抄报圈