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

不知道CSDN现在接了多少广告，如果觉得界面有遮挡，请移步我的hexo博客:http://shanchuantian.com/: 【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.代码

# 变量
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.结果