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 -ZTLxyhwmakecpt -Cglobe -T-10000/10000/1000 -Z > colors.cptgrdimage etopo5.grd -Ba60g30 -Rg -JN180/25c -Ccolors.cpt -K > world.pspsscale -Ba2500f500::/:"m": -Ccolors.cpt -D12.5c/0.3c/15c/.35ch -Y-0.75i -O -U"Author:YangLei">> world.ps
（2） Eckert IV projection 该投影也专门用于绘制全球地图，属于伪圆柱等积投影。其特征是中央子午线为直线，最外侧子午线为半圆形，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.psgrdimage etopo5.grd -Ba60g30 -R-20/60/-90/90 -Ji20/0.02i -X2.8i -Ccolors.cpt -K -O>> world.psgrdimage etopo5.grd -Ba60g30 -R60/200/-90/90 -Ji130/0.02i -X1.6i -Ccolors.cpt -K -O>> world.ps
（5）墨卡托投影投影方式设置为-JM25c，还要注意-R的设置，墨卡托投影不能设置为-Rg。
页:
[1]