周跳检测
tedit基于Blewitt(1990)的预处理算法,用于识别RINEX观测文件中的不良测量值和周跳,生成RINEX健康诊断文件(rhd)。
使用方法:
tedit o_file -int -rnxn -rhd -xyz -short -lc_check yes/no/only -pc_check -elev -len -time
tedit.f90
此程序是周跳检测可执行程序的主程序。其处理流程如下:
- 调用
get_control_parameter()
函数将tedit后面的参数读入。主要包括RINEX观测文件、采样间隔、广播星历文件、开始结束的时间、编辑数据的长度、最小高度角、初始测站坐标等等。 - 如果使用广播星历,调用
rdrnxn
函数读取广播星历,并将星历内容保存到ephem
结构体中. - 读取rinex观测文件,调用
read_rinex_file
函数,在读取过程中,通过无几何距离组合给存在跳跃的历元给打上标记。检测并修复接收机钟跳。保存MW观测值并检查宽巷观测值是否为坏值;检查高度角是否过低。 - 调用
remove_short()
函数标记大的周跳。
# 函数目的:移除短弧段和标记大gap
subroutine remove_short(keep_end, nepo, ti, flg, len_short, len_gap, interval, flag_shrt, removed)
* 通过时间差判断是否存在大的gap,并标记为“gap”
* 循环查找每个模糊度对应的观测弧段,检测护短长度,并标记小于阈值的弧段,标记为移除。
- 如果需要检测单差观测值,则调用
check_sd()
函数检测单差观测值。 - 调用
edit_wideline()
函数,利用MW组合观测值检测周跳,调用check_ionosphere()
函数检查无电离层组合观测值。 - 生成rhd文件。