【GMT5】在全球地图上绘制点和文字

目标

在全球背景地图上绘制点和文字。

代码

准备文件：

点文件：example_point.txt

-53.854 -1.999  0.2c
116.594 40.065  0.2c
-96.111 39.672  0.2c

文字文件：example_text.txt，文字的坐标由点坐标的纬度加2度得到。

-53.854 0.001   10p 0   BC  Amazon
116.594 42.065  10p 0   BC  BeiJing
-96.111 41.672  10p 0   BC  America

代码：

# 在全球地图上标一些点
ps=example.ps
R=-180/180/-80/80
J=X16d/9d #d : degree of arc
point=example_point.txt
text=example_text.txt

# ps文件头
gmt psxy -J$J -R$R -T -K > $ps

# 海岸线
gmt pscoast -R$R -J$J -Ba60  -K -Ggrey -Swhite -O >> $ps

# 文字
gmt pstext $text -R$R -J$J -F+f+a+j -K -O >> $ps
# 画点
gmt psxy $point -R$R -J$J -Sa -Gblack -K -O >> $ps

# ps文件尾
gmt psxy -R$R -J$J -T -O >> $ps

# 保存
gmt psconvert $ps -A -Tg -P

结果