不知道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 $ps3.结果
