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

目标

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

代码

准备文件:

点文件:example_point.txt

1
2
3
-53.854 -1.999 0.2c
116.594 40.065 0.2c
-96.111 39.672 0.2c

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

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

代码:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
# 在全球地图上标一些点
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

结果