使用GMT绘制GPS站点分布图_duweiji

个人资料

微博

正文字体大小：大中小

使用GMT绘制GPS站点分布图

(2012-05-08 18:05:28)

标签：

分类： GMT

使用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

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

阅读┊ 收藏 ┊ 喜欢 ▼ ┊打印┊举报/Report