直接给出代码：（需要有 ETOPO1_Ice_g_gmt4.grd 数据）

gmt gmtset FORMAT_GEO_MAP=ddd:mm:ssF

set PS="topo.ps"

gmt grdcut ETOPO1_Ice_g_gmt4.grd -R70/135/15/55 -GcutTopo.grd

gmt grdgradient cutTopo.grd -Ne0.7 -A50 -GcutTopo_i.grd

gmt grd2cpt cutTopo.grd -Cglobe -S-10000/10000/20 -Z -D>colorTopo.cpt

gmt psbasemap -R70/135/15/55 -JM7i -Bf5a10 -BWesN -Xc -Yc -K > %PS%

gmt grdimage cutTopo.grd -IcutTopo_i.grd -R -J -CcolorTopo.cpt -Q -O -K >>%PS%

gmt pscoast -R -J -Dh -W1/0.2p -I1/0.25p -N1/0.5p -O -K >>%PS%

gmt psscale -DjCB+w7i/0.15i+o0/-0.5i+h -R -J -CcolorTopo.cpt -Bxa2000f400+l"Elevation/m" -G-8000/8000 -O -K>>%PS%

gmt psxy -R -J -T -O >>%PS%

gmt psconvert %PS% -A -Tj -P

del gmt.*

http://s7/mw690/002XWA0fzy76yLTcmBE26&690