本文给出了使用GMT绘制GPS站点分布的方法和代码
安装GMT的方法请参考https://docs.gmt-china.org/latest/install/。其中 grdPath 文件夹下的 earth_relief_15s.grd 文件需要自行到http://topex.ucsd.edu/WWW_html/mar_topo.html下载。
1.数据准备
- 准备站点数据文件 site.txt
为了简洁，这里仅给出3个站点以作示例
数据文件的格式为
#站名 经度 纬度
8118 120.5 23.4
aknd 120.3 22.8
banc 121.4 24.9
2.编写代码
本脚本文件需要在Linux系统下运行
#!/bin/bash
#plot sites distribution
set -o errexit
set -o nounset
gmt gmtset FONT_ANNOT_PRIMARY=12p,Times-Bold # 坐标轴字体，小五号字体，Times-Roman
gmt gmtset FONT_LABEL=9p,25,black
# gmt gmtset MAP_FRAME_TYPE plain # 地理地图边框设置
gmt gmtset FORMAT_GEO_MAP ddd:mm:ssF # 地图的单位标注格式
gmt gmtset MAP_TITLE_OFFSET 0p # 标题离坐标轴的距离
gmt gmtset PS_MEDIA=a5
gmt gmtset PS_PAGE_ORIENTATION=portrait # portrait
gmt gmtset MAP_FRAME_WIDTH=0.15c # set width for axies
R=119.90/122.10/21.80/25.40
J=M10c
B=a1f0.5
## setting the path of topographic data
grdPath=./
## To change the sta file for plotting specific station distribution
InputDataPath=./
staDistribution=${InputDataPath}/site.txt
##setting the path of results
results=./result
OutputPS=${results}/sites.ps
##preparing the prerequisite file for basemap
gmt grdcut ${grdPath}/earth_relief_15s.grd -R$R -GcutTopo.grd
gmt grdgradient cutTopo.grd -Ne0.7 -A50 -GcutTopo_i.grd
gmt grd2cpt cutTopo.grd -Cglobe -S-10000/10000/200 -Z -D > colorTopo.cpt
## plotting the basemap
gmt psbasemap -R$R -J$J -B$B -BWS -X1.5c -Y1.5c -K > $OutputPS
gmt psbasemap -R$R -J$J -B$B -Ben -K -O >> $OutputPS
gmt grdimage cutTopo.grd -R$R -J$J -IcutTopo_i.grd -CcolorTopo.cpt -Q -O -K >> $OutputPS
gmt pscoast -R$R -J$J -Da -W1/0.35p -I1/0.25p -O -K >> $OutputPS
## plotting station Distribution
awk '{print $2,$3}' $staDistribution | gmt psxy -R -Gred -Sc0.2c -J$J -B+n -O -K >> $OutputPS
gmt psxy -R -J -T -O >> $OutputPS
# gmt psxy -R -J -T -O -K >> $OutputPS
## plotting name of station (optional)
# awk '{print $2,$3,$1}' $staDistribution | gmt pstext -R -J -F+f4p,1,black+jTL -O >> $OutputPS
## convert PostScript file to a specific format image file
gmt psconvert -A0.5c -Tf -P $OutputPS
rm gmt.* cutTopo*.grd colorTopo.cpt
echo "Done!"
你可以根据自己的需要进行更改。
3 效果图
站点分布图
