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*

效果图：

注明：由于上述脚本出现了问题，因此上传脚本和数据文件，供大家下载。

igs_cn.list

igsname_cn.gmt

igsnoname_cn.gmt