GMT范例之三维地形图上加文字

葡萄

gmtset PLOT_DEGREE_FORMAT ddd:mm:ssF
set filename=gmt2011_03a_text.ps
grd2cpt tw.grd -Crainbow -Z > tmp.cpt
grdgradient tw.grd -A45 -Gtw.grd.int -Ne1 -M
grdview tw.grd -Jm1.8 -Jz0.0003 -R117/125/19/27/-10000/4000 -Ctmp.cpt -Itw.grd.int  -E150/40 -Qi -N-10001/200/200/200 -B1/1/2500:"meter":SEwnz -K -V -X3 -Y4 >> %filename%
pstext whiteword.dat -Jm -R -G255/255/255 -N -O -K -V >> %filename%
pstext bk.txt -Jm -R -G255/0/255 -K -V -O>>%filename%
psscale -D10/-1/15/0.3h -Ctmp.cpt -Bf200a1000/:meter: -I -O -V >> %filename%
del .gmt*
del *.int
del *.cpt
附：pstext所用数据
==========================================================
whiteword.dat文件内容
# x,y,size,angle,fontno,justify,text
119.40 22.50 15 0      5 6 South China Sea
119.40 22.10 15 0      5 6 Basin
123.50 20.05 13 50    5 1 West Philippine Basin
119.90 21.50 15 35    5 1 Manlila Trench
122.95 21.90 10 50    5 1 Huatang Basin
125.17 22.13 10 335  5 6 Ryukyu Trench
126.50 24.00 15 0      5 6 Okinawa
126.80 23.75 15 0      5 6 Trouth
121.60 24.90 17 15    5 6 Taiwan Strait
=============================================================
bk.txt文件内容
# x,y,size,angle,fontno,justify,text
122.5 26 30 0 33 CT http://xxqhome.blog.163.com

文件中最好加上一句注释，否则容易搞混各列是啥数据。