武大李星星团队开源的GNSS UPD估计模块（GREAT-UPD）window 下调试

目录

【GREAT-UPD简介】

【在window下调试】

1 环境要求和许可证

2 下载

3 编译（Windows下）

【数据的处理】

1 数据准备

2 数据预处理

2.1 单个项目的预处理

2.2 批量预处理

3 估计UPD

3.1 单个项目的估计UPD

3.2 批量估计UPD

4 结果绘图和分析

5 补充说明

5.1 ambupd文件格式说明

【参考】

---

【GREAT-UPD简介】

武汉大学李星星团队在【GPS Solution】开源了GNSS UPD估计模块（GREAT-UPD）。可为 PPP 模糊度固定解计算小数周偏差改正数。根据介绍，该模块可以估计目前四大主流的卫星导航系统（GPS、GLONASS、Galileo 和 BDS）的超宽巷、宽巷和窄巷 UPD。软件处理的流程图如下：

注意：所有的路径不能含有中文

【在window下调试】

1 环境要求和许可证

环境为：Microsoft Visual Studio 2017 on Windows 10（64bit）

许可证：GNU通用公共许可证

2 下载

下载链接：https://geodesy.noaa.gov/gps-toolbox/

下图为GREAT-UPD的目录

3 编译（Windows下）

示意图如上所示，具体步骤说明见下文

3.1 下载并安装CMake，打开CMake界面

3.2 点击【Browse Source】打开GREAT-UPD_1.0文件夹；在GREAT-UPD_1.0下新建build文件夹，点击【Browse Build】打开build文件夹；

3.3 接着点击左下角的【Configure】出现IDE选择界面；

3.4 接着点击主界面的【Generate】将build files写入3.2中新建的build文件夹；接着点击【Open Project】（路径必须英文，否则没反应）；最后在build文件夹下面打开项目即可调试

 

【数据的处理】

约定：YYYY: 4-digit year; YY: 2-digit year; MM: 2-digit month; DD: 2-digit day; DOY: 3-digit DOY (Day of Year). 

1 数据准备

估计超宽巷EWL UPD和宽巷WL UPD需要GNSS观测文件（O文件）、DCB改正文件、

有的观测文件缺少GLONASS的频率号和BDS的SICB（satellite-induced code bias），所以估计GLO和BDS的EWL和WL时还需要广播星历文件

估计GPS Block IIF 卫星EWL的UPD还需要IFCB文件（需要自己事先估计，格式参考GREAT-UPD官方手册A.6）

估计窄巷NL UPD时需要ambupd文件和WL UPD文件（格式参考A.3和A.5）

NL UPD通过宽巷模糊度和无电离层模糊度求得，宽巷模糊度通过MW组合观测值求得，无电离层组合模糊度通过浮点PPP求得

GREAT-UPD提供python脚本下载GNSS观测文件、导航文件、DCB文件。（脚本在GREAT-UPD_1.0/util/batch_process）

（注意：python 脚本需要Python 3.* 的环境）

2 数据预处理

GREAT-UPD提供数据预处理工具GREAT-PreEdit，能探测GF和MW周跳、提出异常观测值、eliminate short arcs等；，输出结果为ambflag文件（参考A.4）

准备预处理配置文件，XML格式（在GREAT-UPD_1.0\util\PreEdit\PreEdit_config目录下，配置文件的设置参考A.2）

2.1和2.2互斥，选一个执行即可

2.1 单个项目的预处理

a）新建项目project_test（例如E:\project_test），将观测文件（.**o）和导航文件（.**n）、预处理配置文件（.xml）放在目录下，如下图所示

b）修改预处理配置文件（.xml）

c）打开cmd命令行，输入下图所示命令

这样就得到了ambflag文件

2.2 批量预处理

GREAT-UPD提供了一个可以逐站逐天处理的python脚本"PreEdit.py"（GREAT-UPD_1.0/util/batch_process）。

a）准备预处理初始化文件（.ini）、参考XML文件、“sitelist”文件（参考GREAT-UPD_1.0/util/PreEdit/PreEdit_sample目录下）。这里先把PreEdit_Windows.ini和site_list两个文件拷贝到E:\project_test目录下，并做以下修改

b）定位到预处理配置文件（.ini）所在目录下，执行PreEdit.py文件

首先cd到PreEdit_Windows.ini所在目录

执行：python ..\..\PreEdit.py -c PreEdit_Windows.ini

或者：python ..\..\PreEdit.py --config=PreEdit_Windows.ini

注意：所有的路径不能含有中文

 

3 估计UPD

3.1和3.2互斥，选一个执行即可

3.1 单个项目的估计UPD

准备：配置文件UPD XML（具体格式详见A.1，文件在GREAT-UPD_1.0/sample data/UPD_2020001/XML路径下）

注意：处理宽巷UPD和窄巷UPD时的配置文件是不同的，处理窄巷UPD时需要先算出WL UPD

修改配置文件如下图所示

打开命令行，定位到GREAT-UPD.exe的目录下，执行配置文件

当命令行出现“Normal End!”字样代表计算完成，结果默认在GREAT-UPD.exe的目录下

3.2 批量估计UPD

像预处理一样，GREAT-UPD提供python脚本“upd.py”（在GREAT-UPD_1.0/util/batch_process目录下）

准备一个初始化文件（.ini），做如下图修改

定位到配置文件（.ini）目录下，使用python脚本处理配置文件

注意：所有的路径不能含有中文

 

4 结果绘图和分析

参考手册

5 补充说明

5.1 ambupd文件格式说明

IF模糊度可以取稳定后的平均值，比如100个历元之后的平均值。若连续历元都存在x卫星，则文件中的这些历元的IF值、WL值、标准差都是一样的。

 

【参考】

1. 公众号【不迷途导航程序员】
2. GREAT-UPD官方手册（GREAT-UPD_1.0/doc/GREAT-UPD_1.0.pdf）
3.