使用GMT绘制GPS站点分布图
本文讲解GMT在大地测量学中一个重要的用途，绘制点位分布图。以全球IGS站和中国某地区的GPS连续运行站位例子，最终分布图如下。图（1）：带有省界的中国某地区XXX分布。
图（2）全球部分IGS站点分布
文件格式：带有坐标和GPS测站点名的文本文件，前两列为经纬度，第三列为GPS站点的缩写名称。http://hiphotos.baidu.com/yangleir/pic/item/8b6ef7ff5266d0169f86fe8a972bd40734fa3528.jpg
代码如下：Rem windows 7，GMT4.5.5Rem 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.psrem psbasemap -BNEWSf1a2g1/f1a2g1 -R114/120/36/42 -JM5i -P -K >site.psRem 下面一行添加指定形状的点位，颜色和形状通过-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.psdel *.ds4REM Edit by YANG Lei//**************************************************************//如图形和坐标涉及保密信息，请告之作者。
页:
[1]