|||
GMT博大精深,功能强大,我也算是一个刚入门的初学者,希望能将我的学习GMT的经历与大家分享,希望能帮助到没有头绪的初学者,当然对于GMT的高手,我也希望GMTer们能多多交流,共同提高。
NOTE: 如果有更好的建议或者改进请及时通知原作者 zhouforme@gmail.com。
主要针对中国地区IGS站分布图
GMT脚本(igs_cn.gmt):
#!/bin/bash
#
# Purpose:Plot IGS sites in China
# Author : Feng Zhou
# Date : 2012-12-10
#
ps=igs_cn.eps
gmtset HEADER_FONT_SIZE 10 X_ORIGIN 4c Y_ORIGIN 8c
ANNOT_FONT_SIZE_PRIMARY 6 FRAME_WIDTH 0.05c
pscoast -R70/140/15/55 -JM5i -Glightyellow -Slightblue -U"Made by Feng Zhou" -A1000
-Ba10f5/5:."IGS Stations in China": -K -Dh -N1/1p -Wthinnest -P > $ps
cat $1 | awk '{print $1, $2}' | psxy -R -J -O -Sc0.15 -Gred -Wthinnest >> $ps
# Clean up after ourselves:
rm -f .gmt*
NOTE:$1是可执行脚本后面带的参数,比如本例子./igs_cn.gmt igs_cn.list这样才可以正确运行。
数据文件(igs_cn.list):
115.89 39.61 BJFS LM
121.20 31.10 SHAO LM
114.36 30.53 WUHN LM
87.60 43.81 URUM LM
125.44 43.79 CHAN LM
109.22 34.37 XIAN RM
102.80 25.03 KUNM LM
91.10 29.66 LHAZ LM
87.18 43.47 GUAO RM
NOTE:格式为经度,纬度(北纬为正,南纬为负),站点名,站点名与本例子的圆圈表示的相对位置 LM-left middle
效果图如下:
上面的例子画出来的没有站点名(针对很多个站点的时候比较好),如果站点数不多,想加上站点名,数据文件不变,脚本需要改变一下。
改后脚本:
#!/bin/bash
#
# Purpose:Plot IGS sites in China
# Author : Feng Zhou
# Date : 2012-12-10
#
ps=igsname_cn.eps
gmtset HEADER_FONT_SIZE 10 X_ORIGIN 4c Y_ORIGIN 8c \
ANNOT_FONT_SIZE_PRIMARY 6 FRAME_WIDTH 0.05c
pscoast -R70/140/15/55 -JM5i -Glightyellow -Slightblue -U"Made by Feng Zhou" -A1000 \
-Ba10f5/5:."IGS Stations in China": -K -Dh -N1/1p -Wthinnest -P > $ps
cat $1 | awk '{print $1, $2}' | psxy -R -J -O -K -Sc0.15 -Gred -Wthinnest >> $ps
cat $1 | awk '{print $1, $2, 5, 0, 5, $4, $3}' | pstext -R -J -O -Dj0.15/0 -Gred -N >> $ps
# Clean up after ourselves:
rm -f .gmt*
效果图:
注明:由于上述脚本出现了问题,因此上传脚本和数据文件,供大家下载。
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-9-28 01:16
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社