最近我们在船舶上单独安装了一套双频GNSS,通过PPP解算用来验证GPS验潮的可能性。
GNSS观测文件是长格式:IGS000USA_R_20231920000_01D_01S_MO.rnx
观测时间为2023年7月11日(GPS时间)。
从ftp://igs.ign.fr/pub/igs/products/下载精密星历、钟差和地球自转参数:
ESA0OPSFIN_20231920000_01D_05M_ORB.SP3
ESA0OPSFIN_20231920000_01D_30S_CLK.CLK
ESA0OPSFIN_20231900000_07D_01D_ERP.ERP
1、使用TerraPos(GPS+GLONASS)
解算的结果如下图,即使选择GPS+GLONASS+BeiDou,解算结果与GPS+GLONASS一摸一样。
下载IGS精密星历、钟差和地球自转参数:
IGS0OPSFIN_20231920000_01D_15M_ORB.SP3
IGS0OPSFIN_20231920000_01D_30S_CLK.CLK
IGS0OPSFIN_20231900000_07D_01D_ERP.ERP
处理过程显示很多数据超限,结果很差。
下载GFZ精密星历、钟差和地球自转参数:
GFZ0OPSFIN_20231920000_01D_05M_ORB.SP3
GFZ0OPSFIN_20231920000_01D_30S_CLK.CLK
GFZ0OPSFIN_20231900000_07D_01D_ERP.ERP
结果见下图,比ESA稍差。
将计算结果进行EGM2008改正到85面,与实测潮位的对比如下图,有点惨不忍睹。
2、使用PRIDE-PPPAR处理
报一个错误:
prepare products:no IGS Antex file:/home/Administrator/PRIDE-PPPAR-master/table/igs20_2247.atx。
更新PRIDE-PPPAR,首先删除PRIDE-PPPAR的安装和执行文件,从https://github.com/PrideLab/PRIDE-PPPAR下载源文件并在msys2环境下安装。
报编译错误:
make[1]: Leaving directory '/home/Administrator/PRIDE-PPPAR-master/src/lib'
make -j4 -C orbit
make[1]: Entering directory '/home/Administrator/PRIDE-PPPAR-master/src/orbit'
gfortran -g -O3 sp3orb.f90 get_sp3orb_args.f90 ../lib/libcom.a -o sp3orb
collect2: fatal error: cannot find ‘ld’
compilation terminated.
make[1]: *** [Makefile:17: sp3orb] Error 1
make[1]: Leaving directory '/home/Administrator/PRIDE-PPPAR-master/src/orbit'
make: *** [Makefile:14: orbit] Error 2
cp: cannot stat './scripts/pdp3.sh': No such file or directory
cp: target '/home/Administrator/.PRIDE_PPPAR_BIN' is not a directory
errror: PRIDE-PPPAR installation failed!
解决办法 :将usr/bin/ld.bfd.exe改名为ld.exe,源文件也保留,不要删除,再次安装成功。
执行命令:
# pdp3 IGS000USA_R_20231920000_01D_01S_MO.rnx
处理过程报错,见下图,似乎是GLONASS卫星R22有问题。
剔除GLONASS卫星,重新计算。
# pdp3 -sys GECJ IGS000USA_R_20231920000_01D_01S_MO.rnx
解算成功,结果经EGM2008改正以后,与实测潮位的对比见下图。
3、小结
从上面的结果来看,即使在船舶上单独安装了GNSS双频观测站,无论是使用TerraPos还是PRIDE-PPPAR经行PPP解算,计算结果都不理想,两者相比较,矮子里拔将军,PRIDE-PPPAR的结果与实测潮位形态更为贴合些。
参考文献
https://blog.csdn.net/qq_41896822/article/details/123739888
https://blog.csdn.net/weibuu/article/details/109328428