GMT绘制三维地形图

huolongguo

bat代码：

gmtset ANOT_FONT_SIZE 8 LABEL_FONT_SIZE 12 ANOT_FONT 1 LABEL_FONT 1
set psf=topo_cdm.ps
set inf=topo.dat

psbasemap -R0/16/0/24 -Jx1 -B0 -P  -K -U"xiaoqiang" > %psf%
echo 5 22.5 30 0 33 1 Qaidam Basin Topo | pstext -R -J -O -K -V>> %psf%

surface -I2m -V %inf% -Gtmp.grd -R90/100/35/40
grd2cpt tmp.grd -Crainbow -Z -V > tmp.cpt
grdgradient tmp.grd -A270 -Gtmp.grad
grdhisteq tmp.grad -Gtmp.his -N
grdmath tmp.his 0.4 x = tmp.int

grdview -K -O -S -P tmp.grd -JM10 -JZ2.5 -Ba1/a1/2000:"Topo(m)":SEwnZ -Qs -N-1/235 -Itmp.int -X2 -Ctmp.cpt -R -E150/45 -Y14 >> %psf%

grdcontour tmp.grd -J -P  -O -K -E220/50 -R -C100 -A1000+p+g255+s5 -Wa3/0  -X-1 -Y-12 >> %psf%

grdview -K -O -S -P tmp.grd -J -JZ2.5i -Ba1/a1/2000SeWnz  -Qs -N-100 -Itmp.int -Ctmp.cpt -R -K -E220/50 >> %psf%

echo 88 36 0 0 1.1 | psxyz -J -E220/50 -R -SV0.2/0.5/0.4 -W1p -O -K -N -V>> %psf%
echo 88.25 37.3 20 180 1 5 N | pstext -J -E220/50 -R -O -K -N >>%psf%

psscale -E -D13/5/7/0.3 -Ctmp.cpt -Ba1000g500/:"(m)": -O -P -V -K>> %psf%

del .gmt* tmp.*

效果图：