【GMT】用GMT绘制测高卫星Topex_Poseidon_Jason的地面轨迹

运行环境:Ubuntu16.04+GMT5.2.1

1.思路

AVISO网站提供了多个卫星的kml文件,这是现成的轨迹数据,想办法用GMT读取出来即可,正好GMT的kml2gmt模块可以将Google Earth的KML文件转换为GMT表数据。如果是kmz文件,GMT似乎不能转换,在Google Earth中将其另存为kml即可。

注:Topex、Jason1、Jason2的轨迹是重合的。

2.代码

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
# 变量
ps=global.ps
J=M25c
R=-180/180/-70/70
# ps文件头
gmt psxy -J$J -R$R -T -K > $ps
# 将Google Earth的KML文件转换为GMT表数据
gmt kml2gmt Visu_J2J1TP_Tracks_GE_V3.kml > TPTrack.txt
# 海岸线
gmt pscoast -R$R -J$J -Ba60g30 -I1/1p,blue -I2/0.25p,blue -N1/0.25p,- -W0.25p,white -Ggreen -Sblue -K -O >> $ps
# 卫星轨迹
gmt psxy TPTrack.txt -R -J -W1p,lightred -K -O>> $ps
# ps文件尾
gmt psxy -R$R -J$J -T -O >> $ps
# 导出jpg
gmt ps2raster $ps -Tj -E300 -A -P
# 用gs打开ps文件
gs $ps

3.结果

global