由于前一阵受一师兄拜托,需要解算5个测站站点时间序列,daily 解可以在NGL(内华达大地测量实验室)下载,同时NG了也提供了5min 解,但不完整,有的站不包括在内。很巧,这几个站就不在。解算时间为12-14年,早年间由于gamit起手,所以尝试用gamit 进行解算,然而由于数据量大,在跟新http://station.info 出各种问题,手动添加也不识别(推测是由于10.7版本的不稳定,10.5其实还好)大批量处理经常出现 too long 。如果分阶段处理又较为耗时。
于是放弃,后又尝试track模块bernese rtklib 等软件,最后选择武汉大学导航中心的pride——pppar,你问为什么?因为有人指导啊~其他的还得自己慢慢摸索。
我两天时间往导航中心跑了好几次,但是一直出各种问题,今天没睡午觉,跑到他们实验室跟软件的合作开发者也是我健身的师弟学习。最终解决所有问题,学成而归,用一下午的时间告诉了同师门一个师兄。这软件很强大,同时很智能,极易上手,但同时出现些小的问题,故记录在此。
(一)
PRIDE-PPPAR开源软件平台源自于葛茂荣教授最初的PPP-AR软件并经由耿江辉教授其后的改进发展而来。本软件包含了武汉大学GNSS中心众多科研人员集体的努力和付出,我们在此对他们表示由衷的感谢。我们开源本软件的目的是希望能够帮助到那些新进入GNSS领域的科研工作者以及推动地球物理和大地测量领域的研究工作。特别地,我们希望本软件能够贡献于地壳运动研究和大气变化监测等高精度地学应用中。
软件下载链接:
武汉大学PRIDE课题组pride.whu.edu.cn
注意,安装包里的table表不是最新,会出问题,需要在该页面下下载最新(15年更新)的table。也可以手动添加改正,如果你知道怎么改的话~
系统环境为Linux 乌班图 centos等较为常用的系统,软件构架和模块可以参考说明书。注:自带算例只是验证软件的安装成功与否,不能代表计算的结果。
(二)
安装:直接运行安装包里的install.sh(感谢潘元欣师弟的辛勤奉献,方便的大家),一步到位。
系统应预装awk(乌班图没有,很奇怪,一般都自带)还有gfortarn. 还需要预装一个32位库,不然后面rnx2rtkp 会运行出错。
OK ,开始运行:首先准备好自己的rinex文件和brdc文件,data/年份/年积日 放在一起,比如:
准备好table
注: 如果32位库装了还是出问题,说rnx2rtkp 没有权限,站点没有数据之类的,那么需要手动在table/sit_xyz 手动添加先验坐标。先验坐标可以在NGL 获得。
编辑脚本config_template,添加自己的数据路径(三行),添加自己的站点名称(最下面):
基本上不出问题这些就可以直接使用命令运行:
pride_pppar config_template YYYYMMDD(起始时间) YYYYMMDD(结束时间) N(要不要算模糊度,一般不用)
结果会生成以年命名的文件,pos 前缀就是单天解。
接下里需要计算小时解,这时候需要另一个config_template 以及一个脚本hpc.sh ,config_template也是上改路径下改站点。
另注意:table里的海潮文件ocean基本上包括全球的IGS点,但有的也是缺失,需要我们自己算一下然后添加进去,可以由下面的网址在线算:
Ocean tide loading providerholt.oso.chalmers.se
选择好模型,输入站点经纬度,然后输入自己的邮箱,就可以啦~
最后运行hpc.sh 脚本。就出来站点坐标的小时解啦。
不过是XYZ坐标,他自带了函数可以转换为ENU 。需要的自己再去看看。
OK ,这就是整个pride_PPPAR的心得。
分享一下,希望能帮助大家科研进步。