葡萄 发表于 2013-4-13 22:16:01

GMT范例之从地形图上切剖面

set filename=gmt2011_03c1.ps
set start=116/49
set end=135/40
project -C%start% -E%end% -G1 -Q > track.dat
grdcut etopo2_new.nc -Gtmp.grd -R110/140/35/50 -V
grd2cpt tmp.grd -Csealand -Z > tmp.cpt
grdgradient tmp.grd -A60 -Gtmp.grd.int -Ne1 -M -V
psbasemap -Jx1 -R-2/15/0/20 -B0 -V -K -P>>%filename%
pscoast -JM11 -R110/140/35/50 -Df -Ia/1p/100 -B5 -W3 -K -V -X2 -Y10 -O>>%filename%
grdimage tmp.grd -JM -R -Ctmp.cpt -Itmp.grd.int -B5 -K -V -O>>%filename%
psxy track.dat -JM -R -W5/red -O -V -K>>%filename%
psscale -D12.5/3/7/0.3 -Ctmp.cpt -Bf200a1000/:m: -I -V -O -K>> %filename%
grdtrack track.dat -Gtmp.grd | gawk "{print $3,$4}" > track_profile.dat
minmax -I5 track_profile.dat
psxy -R0/1170/0/1760 track_profile.dat -JX10/4 -Ba500f250:"Distance(km)":/a400f200:"H(m)": -W2/red-P -V -X1.5 -Y-7.5 -O>> %filename%
del .gmt*
del tmp.*
del track*

困惑之处: psxy -R0/1170/0/1760 track_profile.dat -JX10/4 -Ba500f250:"Distance(km)":/a400f200:"H(m)": -W2/red-P -V -X1.5 -Y-7.5 -O>> %filename%中Distance(km)及H(m)字体大小如何控制?不得其解。

假如 发表于 2014-1-18 18:58:58

楼主第五行的etopo2_new.nc 是什么格式的文件呀,解释一下呗,谢谢

假如 发表于 2014-1-18 21:36:45

第七行grdgradient tmp.grd -A60 -Gtmp.grd.int -Ne1 -M -V中的Gtmp.grd.int 是什么格式的数据呀,为什么要做成这种形式呀,不懂
页: [1]
查看完整版本: GMT范例之从地形图上切剖面