使用GMT绘制GPS站点分布图

小二

本文讲解GMT在大地测量学中一个重要的用途，绘制点位分布图。

以全球IGS站和中国某地区的GPS连续运行站位例子，最终分布图如下。

图（1）：带有省界的中国某地区XXX分布。

图（2）全球部分IGS站点分布

文件格式：

带有坐标和GPS测站点名的文本文件，前两列为经纬度，第三列为GPS站点的缩写名称。

代码如下：

Rem windows 7，GMT4.5.5

Rem creat fig(1):XXX site distribution of XXX area in China.

pscoast -R114/120/36/42 -JM5i -K -Ba2g1/a2g1  -Di -A5000 -G221/206/167 -S135/240/56 -Na -Ia>site.ps

rem psbasemap -BNEWSf1a2g1/f1a2g1 -R114/120/36/42 -JM5i -P -K >site.ps

Rem 下面一行添加指定形状的点位，颜色和形状通过-G和-S设置。                                                         

psxy HB.txt  -R  -JM  -K -O  -Gred -St0.08i  -W0.25p >>site.ps

Rem 下面一行绘制中国的省界，还需另备省界文件。                     

psxy prov.dat -JM -R -W1.0p/red -m/nan -B -O -K -P>>site.ps

Rem 下面一行添加站点名称，需要linux gawk的支持，新版本GMT内置该小程序。老版Rem 本GMT需安装cygwin（一种windows 下的linux虚拟环境）

gawk "{print $1, $2, 12, 0, 1, 1, $3}" HB.txt | pstext -R -JM -D0.1i/-0.115i -O>>site.ps

del *.ds4

REM Edit by YANG Lei

//**************************************************************//

如图形和坐标涉及保密信息，请告之作者。