MATLAB学习与使用：8种地图投影方式

MATLAB学习与使用：8种地图投影方式

MATLAB（MATrix &LABoratory）是一款强大的科学软件，可用于数据处理、科学计算、图形绘制、模拟仿真等等。今天，简单展示MATLAB包含的8种地图投影方式：圆柱投影（Cylindrical）、伪圆柱投影（Pseudocylindrical）、圆锥投影（Conic）、多圆锥投影（PolyConic）、伪圆锥投影（PseudoConic）、方位角投影（Azimuthal）、伪方位角投影（Pseudoazimuthal）和改进的方位角投影（Modified Azimuthal）。

圆柱投影（Cylindrical）。在命令行窗口输入maps后回车，可以查看MATLAB提供得8大种地图投影方式，共计71小种方案。下面挑一些方案展示其投影效果。运行下面脚本，依次看到balthsrt,behrmann,lambcyln,mercator四种方案的圆柱投影。

close all,clear all,clc

load coast %加载海岸线

%圆柱形投影，有balthsrt,behrmann,lambcyln,mercator等方案

subplot(2,2,1);axesm balthsrt;patchm(lat,long,'g')

subplot(2,2,2);axesm behrmann;patchm(lat,long,'g')

subplot(2,2,3);axesm lambcyln;patchm(lat,long,'g')

subplot(2,2,4);axesm mercator;patchm(lat,long,'g')

伪圆柱投影（Pseudocylindrical）。运行下面脚本，依次看到apianus，collig，fournier，robinson四种方案的伪圆柱投影。

close all,clear all,clc

load coast %加载海岸线

%伪圆柱投影，有apianus，collig，fournier，robinson等方案

subplot(2,2,1);axesm apianus;patchm(lat,long,'g')

subplot(2,2,2);axesm collig ;patchm(lat,long,'g')

subplot(2,2,3);axesm fournier;patchm(lat,long,'g')

subplot(2,2,4);axesm robinson;patchm(lat,long,'g')

圆锥投影（Conic）。运行下面脚本，依次看到eqaconicstd，eqaconic，eqdconic，lambert，murdoch1，murdoch3六种方案的圆锥投影。

close all,clear all,clc

load coast %加载海岸线

%圆锥投影，有eqaconicstd，eqaconic，eqdconic，lambert，murdoch1，murdoch3等方案

subplot(2,3,1);axesm eqaconicstd;patchm(lat,long,'g')

subplot(2,3,2);axesm eqaconic ;patchm(lat,long,'g')

subplot(2,3,3);axesm eqdconic ;patchm(lat,long,'g')

subplot(2,3,4);axesm lambert;patchm(lat,long,'g')

subplot(2,3,5);axesm murdoch1;patchm(lat,long,'g')

subplot(2,3,6);axesm murdoch3;patchm(lat,long,'g')

多圆锥投影（PolyConic）和伪圆锥投影（PseudoConic）。运行下面脚本，依次看到polyconstd，polycon，vgrint1三种方案的多圆锥投影和bonne，werner两种方案的伪圆锥投影。

close all,clear all,clc

load coast %加载海岸线

%多圆锥投影，有polyconstd，polycon，vgrint1三种方案

subplot(2,3,1);axesm polyconstd;patchm(lat,long,'g')

subplot(2,3,2);axesm polycon ;patchm(lat,long,'g')

subplot(2,3,3);axesm vgrint1 ;patchm(lat,long,'g')

%伪圆锥投影，有bonne，werner两种方案

subplot(2,3,4);axesm bonne;patchm(lat,long,'g')

subplot(2,3,5);axesm werner;patchm(lat,long,'g')

方位角投影（Azimuthal）。运行下面脚本，依次看到breusing，eqaazim，eqdazim，globe，gnomonic，ortho六种方案的方位角投影。

close all,clear all,clc

load coast %加载海岸线

%方位角投影，有breusing，eqaazim，eqdazim，globe，gnomonic，ortho等方案

subplot(2,3,1);axesm breusing;patchm(lat,long,'g')

subplot(2,3,2);axesm eqaazim ;patchm(lat,long,'g')

subplot(2,3,3);axesm eqdazim ;patchm(lat,long,'g')

subplot(2,3,4);axesm globe;patchm(lat,long,'g')

subplot(2,3,5);axesm gnomonic;patchm(lat,long,'g')

subplot(2,3,6);axesm ortho ;patchm(lat,long,'g')

伪方位角投影（Pseudoazimuthal）和改进的方位角投影（Modified Azimuthal）。运行下面脚本，依次看到wiechel一种方案的伪方位角投影和aitoff，bries，hammer三种方案的改进的方位角投影。

close all,clear all,clc

load coast %加载海岸线

%伪方位角投影，有wiechel一种方案

subplot(2,2,1);axesm wiechel;patchm(lat,long,'g')

%改进的方位角投影，有aitoff，bries，hammer三种方案

subplot(2,2,2);axesm aitoff ;patchm(lat,long,'g')

subplot(2,2,3);axesm bries ;patchm(lat,long,'g')

subplot(2,2,4);axesm hammer;patchm(lat,long,'g')