本发明属于gnss观测数据处理的技术领域,涉及一种gnss精密单点定位预处理技术,具体涉及一种奇异点保留滤波降噪的周跳探测方法。
背景技术:
目前,随着gnss(globalnavigationsatellitesystem,全球导航卫星系统)在各种领域中的应用和发展,不同的用户需要不同精度的定位服务。在测绘、大气、地震监测等领域要求的服务精度越来越高。相比于伪距观测量,载波相位观测量具有精度高、噪声低等优点,被广泛的应用于高精度定位定轨服务中。但由于载波相位本身存在周跳、整周模糊度不固定等问题,使得在高精度定位定轨服务中,必须首先实现周跳的精准探测和模糊度的解算。
周跳探测是模糊度解算的前提,只有做到精准的周跳探测才能正确的参数化模糊度,构造正确的法方程求解模糊度,提高gnss精密定位的精度和可靠性。当gnss接收机锁相环失锁就会产生周跳,在低信噪比以及卫星仰角过低时都有可能发生周跳。特别的,对于我国北斗卫星导航系统,geo(geosynchronousorbit,地球同步轨道)卫星有着卫星轨道高、卫星几何差、需要频繁进行机动调整等原因导致周跳偏多。特别的对于一些低仰角(15°-20°)的geo卫星,更是存在频繁的一周小周跳,经典的turboedit方法不能准确可靠地探测该周跳。
mw组合作为周跳探测经典组合量,有着不含几何量和电离层残差误差影响的优点,但噪声偏大不利于小周跳的探测。
turboedit方法采用mw组合作为主要探测量,对于某一历元
其中
技术实现要素:
本发明的目的在于避免现有技术中的不足而提供一种奇异点保留滤波降噪的周跳探测方法,其可以实现对mw组合降噪的同时保持周跳信号不被滤除,从而达到提高探测精度的效果,解决mw组合由于噪声过大导致的一到二周小周跳无法精准探测的问题。
本发明的目的通过以下技术方案实现:
提供一种奇异点保留滤波降噪的周跳探测方法,包括如下步骤:
step1、建立mw组合观测序列
step2、初始化,选取合适的滑动平均窗长
step3、通过中值滤波器剔除mw组合观测序列中的野值;
step4:对序列:
进行一阶差分得到:
其中
step5:对step4中得到的两个一阶差分序列分别进行窗长为
其中
step6:计算窗长
其中,
得到历元
step7:若
其中
并且,若
step8:若
作为进一步的改进,在步骤step1中,由rinex观测文件读取获得伪距与载波相位数据,做线性组合得到mw组合观测序列
作为进一步的改进,在步骤step3中,通过一个窗长为3的中值滤波器剔除mw组合观测序列中的野值。
作为进一步的改进,在步骤step5中,所述极值中值滤波过程如下:
若
step5.1:搜索出集合
中的最大值
其中
step5.2:按照公式:
更新
step5.3:若
作为进一步的改进,由于mw组合的一阶差分序列几乎是零均值的,故在step5.2中,更新
通过本发明的方法,可以得到保留周跳的降噪序列
附图说明
利用附图对本发明作进一步说明,但附图中的实施例不构成对本发明的任何限制,对于本领域的普通技术人员,在不付出创造性劳动的前提下,还可以根据以下附图获得其它的附图。
图1是本发明奇异点保留滤波降噪的周跳探测方法的流程图。
图2是本发明的效果图。
具体实施方式
为了使本领域的技术人员更好地理解本发明的技术方案,下面结合附图和具体实施例对本发明作进一步详细的描述,需要说明的是,在不冲突的情况下,本申请的实施例及实施例中的特征可以相互组合。
如图1所示,本发明实施例提供的一种奇异点保留滤波降噪的周跳探测方法,包括如下步骤:
step1、由rinex(receiverindependentexchangeformat,即与接收机无关的交换格式)观测文件读取获得伪距与载波相位数据,做线性组合得到mw组合观测序列
step2、选取合适的滑动平均窗长
step3、初始化,对滑窗l内的数据,通过一个窗长为3的中值滤波器剔除mw组合观测序列中的野值。按照下述step4到step7计算第
step4、对序列
进行一阶差分得到:
其中
step5、进行上述两个一阶差分序列分别进行窗长为
其中
step6、计算窗
得到历元
step7、若
其中
并且若
step8、若
通过上述步骤,可以得到保留周跳的降噪序列
作为进一步优选的实施方式,在步骤step5中,所述极值中值滤波过程如下:
若
step5.1:搜索出集合
中的最大值
其中
step5.2:按照公式:
更新
step5.3:若
作为进一步优选的实施方式,由于mw组合的一阶差分序列几乎是零均值的,故在step5.2中,更新
如图2所示,图2中的上图为jfng站c01卫星mw组合周跳探测仿真效果图,图中黄色点图为mw组合原始数据点,蓝色实线为新算法滤波降噪后的mw组合平滑曲线,红色虚线为
上面的描述中阐述了很多具体细节以便于充分理解本发明,但是,本发明还可以采用其他不同于在此描述的其他方式来实施,因此,不能理解为对本发明保护范围的限制。
总之,本发明虽然列举了上述优选实施方式,但是应该说明,虽然本领域的技术人员可以进行各种变化和改型,除非这样的变化和改型偏离了本发明的范围,否则都应该包括在本发明的保护范围内。