GMT（Generic Mapping Tools）计算交叉点 --- 以Jason-2测高卫星数据为例

软件：GMT
1.将数据放在E\Crossovers_test\文件夹下的ssh.dat文件中，ssh.dat是由python或Matlab从Jason-2GDR的沿轨测量数据中读取的数据。数据格式为三列：lon lat ssh

数据格式示例：
   2.4270415e+02  -5.2266151e+01  -1.0009100e+01
   2.4275844e+02  -5.2225375e+01  -9.9849000e+00
   2.4281261e+02  -5.2184571e+01  -1.0037100e+01
   2.4286669e+02  -5.2143737e+01  -1.0052400e+01
   2.4292065e+02  -5.2102875e+01  -1.0036600e+01
   2.4297451e+02  -5.2061983e+01  -1.0029000e+01
   2.4302826e+02  -5.2021063e+01  -9.9870000e+00
   2.4308191e+02  -5.1980113e+01  -1.0042900e+01
   2.4313545e+02  -5.1939135e+01  -1.0017100e+01
   2.4318889e+02  -5.1898128e+01  -9.9669000e+00 

2.将以下命令写入E\Crossovers_test\文件夹下的Crossovers.bat文件。

gmt = set X2SYS_HOME=E\Crossovers_test\
gmt x2sys_init Jason -Edat   -Gg -Wd100 -F -Dnew.fmt 
gmt x2sys_cross ssh.dat -THY2A -Ia  -Rg -Z >out1.dat
gmt x2sys_list out1.dat  -THY2A -Cssh -Fxyz > ssh_result.dat

3.将以下命令写入E\Crossovers_test\文件夹下的new.fmt文件。该文件内容的功能参考：x2sys_init

Format define file for Jason project
---------------------------------------------------------------------
ASCII
SKIP 0
GEO
---------------------------------------------------------------------
name	intype	NaN-proxy?	NaN-proxy	scale	offset	oformat
lon 		a 		N 			0 		1 			0 		-
lat 		a 		N 			0 		1 			0 		-
ssh 		a 		N 			0 		1 			0 		-

4.数据（ssh.dat）程序（new.fmt、Crossovers.bat）放在同一级文件夹里（E\Crossovers_test\），双击Crossovers.bat运行即可得到计算结果（ssh_result.dat）。前两列为交叉点的经纬度，后两列分别为交叉点处两条轨迹的ssh。

部分计算结果： lon    lat    ssh_1     ssh_2
-86.440261876	-65.8088612287	-13.3419971252	-14.1282270191
-92.111080542	-65.5805235401	-15.9892816526	-18.1621567403
-90.695052881	-65.4434803959	-13.7933189859	-13.8282145292
-89.275942114	-65.2903534833	-15.4765926613	-15.1938802558
-92.112275256	-65.290609082	-18.3093075898	-16.4735820365
-87.860032778	-65.1211327228	-14.4035084107	-14.2521820198
-86.443450734	-65.2903294781	-13.1671800889	-13.1818843007
-89.277343444	-64.9350758252	-15.2219436863	-15.1567658511
-93.531234311	-65.1211360751	-17.5880653165	-16.2115077744
-94.947656369	-65.2903821444	-19.981725329	-20.0570627193
-86.444619349	-64.9352449113	-12.8342380364	-12.7591940506

参考：GMT x2sys模块
x2sys资料：Crossover point analysis using GMT x2sys
致谢：小李同学， yangleir师兄