红联Linux门户
Linux帮助

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

发布时间:2016-11-08 10:18:13来源:blog.csdn.net/shanchuan2012作者:贪玩的大川
运行环境:Ubuntu16.04+GMT5.2.1
 
1.思路
AVISO网站(http://www.aviso.altimetry.fr/en/data/tools/pass-locator.html)提供了多个卫星的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.结果
用GMT绘制测高卫星Topex_Poseidon_Jason的地面轨迹
 
本文永久更新地址:http://www.linuxdiyf.com/linux/25815.html