基于偶然误差特性的时间序列粗差探测方法及装置与流程

技术特征：
1.一种基于偶然误差特性的时间序列粗差探测方法，其特征在于，包括：根据时间序列拟合模型对原始数据时间序列进行拟合，并得到时间时序的残差序列r；计算残差序列r的绝对值，并将结果进行排序；选取rbin值对残差序列r进行分段；根据残差序列r<rbin分别进行正态分布拟合，并分别求出不同rbin分段的σ
i
；判断是否存在粗差；若不存在粗差，则σ
i
服从同一正态分布，并随着rbin范围所包含的数据增多，σ
i
逐渐趋于固定值；若σ
i
显现出不同趋势，则存在粗差，并对相邻σ
i
进行差分计算；对σ
i
进行插值，并对其去趋势，得到不含趋势项的插值σ
i
；通过局部峰值查找函数，确定查找差分结果表现出不同运动趋势时刻的σ
x
；根据偶然误差分布特性进行粗差去除。2.根据权利要求1所述的基于偶然误差特性的时间序列粗差探测方法，其特征在于，所述时间序列拟合模型公式为：设任意观测序列y(t
i
)，则：其中，i为坐标历元时刻标识；y(t
i
)为测站某一分量t
i
时刻的坐标；b为截距；v0为线性速度；m0为谐波个数；a
m
和b
m
是频率为f
m
时周期项的振幅；为随机噪声，即3.根据权利要求1所述的基于偶然误差特性的时间序列粗差探测方法，其特征在于，所述计算残差序列r的绝对值，并将结果进行排序时是按照从小到大的顺序进行排序的。4.根据权利要求1所述的基于偶然误差特性的时间序列粗差探测方法，其特征在于，所述选取rbin值对残差序列r进行分段时，rbin值在gnss高精度解算中，水平方向可取至1mm。5.根据权利要求1所述的基于偶然误差特性的时间序列粗差探测方法，其特征在于，所述对σ
i
进行插值时可进行不同类型插值。6.一种基于偶然误差特性的时间序列粗差探测装置，其特征在于，包括：拟合模块，用于根据时间序列模型对原始数据时间序列进行拟合，并得到时间时序的残差序列r；排序模块，用于计算残差序列r的绝对值，并将结果进行排序；分段模块，用于选取rbin值对残差序列r进行分段；计算模块，用于根据残差序列r<rbin分别进行正态分布拟合，并分别求出不同rbin分段的σ
i
；判断模块，用于判断是否存在粗差；若σ
i
服从同一正态分布，并随着rbin范围所包含的数据增多，逐渐趋于固定值，则不存在粗差；若σ
i
显现出不同趋势，则存在粗差并对相邻σ
i
进行差分计算；去趋势模块，用于对σ
i
进行插值，并对其去趋势，得到不含趋势项的插值σ
i
；查找模块，用于通过局部峰值查找函数，确定查找差分结果表现出不同运动趋势时刻
的σ
x
；粗差去除模块，用于根据偶然误差分布特性进行粗差去除。7.根据权利要求6所述的基于偶然误差特性的时间序列粗差探测装置，其特征在于，所述时间序列拟合模型公式为：设任意观测序列y(t
i
)，则：其中，i为坐标历元时刻标识；y(t
i
)为测站某一分量t
i
时刻的坐标；b为截距；v0为线性速度；m0为谐波个数；a
m
和b
m
是频率为f
m
时周期项的振幅；为随机噪声，即8.根据权利要求6所述的基于偶然误差特性的时间序列粗差探测装置，其特征在于，所述计算残差序列r的绝对值，并将结果进行排序时是按照从小到大的顺序进行排序的。9.一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至5中任一项所述的方法步骤。10.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现权利要求1至5中任一项所述的方法步骤。

技术总结
本发明公开了一种基于偶然误差特性的时间序列粗差探测方法及装置，包括：根据时间序列拟合模型对原始数据时间序列进行拟合，得到R；计算R的绝对值，并将结果进行排序；选取rbin值对R进行分段；根据R<rbin分别进行正态分布拟合并求出不同rbin分段的σ


技术研发人员：张明 高涵 邱云峰 黄彦锋 陈海林
受保护的技术使用者：云南省测绘产品检测站
技术研发日：2022.02.21
技术公布日：2022/5/30