物探论坛

 找回密码
 立即注册

QQ登录

只需一步,快速开始

查看: 2004|回复: 2

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

[复制链接]
发表于 2013-4-13 22:16:01 | 显示全部楼层 |阅读模式
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*

2835860390376893869.jpg

困惑之处:

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 是什么格式的数据呀,为什么要做成这种形式呀,不懂
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

QQ|Archiver|手机版|小黑屋|物探论坛 ( 鄂ICP备12002012号 微信号:iwutan )

GMT+8, 2024-3-29 05:08 , Processed in 0.136525 second(s), 18 queries .

Powered by Discuz! X3.4

© 2001-2017 Comsenz Inc.

快速回复 返回顶部 返回列表