GMT绘制不同投影方式的地形图

小二

以ETOPO5为底图数据，使用GMT绘制多种投影的全球地形图。

ETOPO有多种类型的数据，之前有ETOPO5和ETOPO2版本，最新的是ETOPO1版本，其分辨率是1′，虽然不如SRTM的分辨率，但是ETOPO数据集成了海底地形和极地冰盖高程数据，对于绘制一些海洋和极地区域的地形图很有价值。

(下载：http://www.ngdc.noaa.gov/mgg/global/global.html )

（1）Robinson投影

首先以Robinson投影为例子，绘制全球陆地地形和海底地形图。

下面是bat代码：

xyz2grd ETOPO5.DAT -Getopo5.grd -I5m -Rg -ZTLxyhw

makecpt -Cglobe -T-10000/10000/1000 -Z > colors.cpt

grdimage etopo5.grd -Ba60g30 -Rg -JN180/25c -Ccolors.cpt -K > world.ps

psscale -Ba2500f500::/:"m": -Ccolors.cpt -D12.5c/0.3c/15c/.35ch -Y-0.75i -O -U"Author:YangLei">> world.ps

（2） Eckert IV projection[equal-area]

该投影也专门用于绘制全球地图，属于伪圆柱等积投影。其特征是中央子午线为直线，最外侧子午线为半圆形，Scale is true along the parallel at 40:30 North and South.

需要修改grdimage 程序中的投影为-JKf180/25c，GMT绘制的全球地图如下：

（3）Sinusoidal projection

修改投影为 -JI0/25c，结果如下：

（4）Interrupted Sinusoidal projection

其绘制方法有点区别，需要分次绘制

grdimage etopo5.grd -Ba60g30 -R-160/-20/-90/90 -Ji-90/0.02i -Ccolors.cpt -K > world.ps

grdimage etopo5.grd -Ba60g30 -R-20/60/-90/90 -Ji20/0.02i -X2.8i -Ccolors.cpt -K -O>> world.ps

grdimage etopo5.grd -Ba60g30 -R60/200/-90/90 -Ji130/0.02i -X1.6i -Ccolors.cpt -K -O>> world.ps

（5）墨卡托投影

投影方式设置为-JM25c，还要注意-R的设置，墨卡托投影不能设置为-Rg。