[发明专利]一种卫星重力观测数据异常值提取方法、装置及电子设备

说明书

本发明实施例涉及测绘工程技术领域，公开了一种卫星重力观测数据异常值提取方法、装置及电子设备。本发明中的卫星重力观测数据异常值提取方法包括：获取第一卫星重力观测数据，采用小波连续函数降噪方法对第一卫星重力观测数据进行异常值探测，根据第一卫星重力观测数据中位于异常值位置的前后数值计算一个新的数值，并将计算出的新的数值替换异常值，获得第二卫星重力观测数据，采用狄克逊检验法对第二卫星重力观测数据进行异常值探测检验。通过小波连续函数降噪方法与狄克逊检验方法的结合，能够实现对卫星重力观测数据中异常值的准确探测。

技术领域

本发明实施例涉及测绘工程技术领域，特别涉及一种卫星重力观测数据异常值提取方法、装置及电子设备。

背景技术

数据是进行科学研究的基础，因此人们运用各种手段获得测量或实验数据。由于环境或人为因素的干扰，测量数据即卫星重力观测数据与被测原型的真值不完全一致，即存在误差。测量学领域将测量数据中的误差按照大小、特征及产生的原因分为系统误差、偶然误差和粗差。在应用测量数据时，要求数据仅含有偶然误差，系统误差可以采用一定的方法发现和修正，粗差即异常值则是由于实验工作者操作时疏忽大意或测量仪器本身有缺陷而测量者又未能及时发现，而导致的严重歪曲测量结果的误差。因此，如果卫星重力观测数据中存在粗差，不论是在科学实验还是航天应用等实际生产领域，都会产生严重的影响。因此在用于科学处理之前，需要对卫星重力观测数据进行异常值探测，将卫星重力观测数据中的异常值标定，剔除或者用合理的数值代替。

发明人发现现有技术中至少存在如下问题：在对卫星重力观测数据进行粗差探测时，目前应用比较广泛的探测方法有小波降噪方法，如小波收缩阈值去噪法。在数据处理过程中，小波收缩阈值去噪法虽然对离散态粗差和区域态粗差比较敏感，但是在粗差探测过程中仍然不能实现准确探测。虽然有人针对小波降噪方法中的小波收缩阈值去噪法进行了改进，但改进后的小波收缩阈值去噪方法的粗差探测效果仍不理想。并且，单独使用现有技术中的任意一类单一粗差探测方法都不能实现对粗差的准确探测。

发明内容

本发明实施方式的目的在于提供一种卫星重力观测数据异常值提取方法、装置及电子设备，通过小波连续函数降噪方法与狄克逊检验方法的结合，实现对卫星重力观测数据中异常值的准确探测。

为解决上述技术问题，本发明的实施方式提供了一种卫星重力观测数据异常值提取方法，包括：获取第一卫星重力观测数据；采用小波连续函数降噪方法对第一卫星重力观测数据进行异常值探测；根据第一卫星重力观测数据中位于异常值位置的前后数值计算一个新的数值，并将计算出的新的数值替换异常值，获得第二卫星重力观测数据；采用狄克逊检验法对第二卫星重力观测数据进行异常值探测检验。

本发明的实施方式还提供了一种卫星重力观测数据异常值提取装置，包括：获取模块，用于获取第一卫星重力观测数据；第一处理模块，用于采用小波连续函数降噪方法对第一卫星重力观测数据进行异常值探测；第二处理模块，用于根据所述第一卫星重力观测数据中位于异常值位置的前后数值计算一个新的数值，并将计算出的新的数值替换异常值，获得第二卫星重力观测数据；第三处理模块，用于采用狄克逊检验法对第二卫星重力观测数据进行异常值探测检验。

本发明的实施方式还提供了一种电子设备，包括至少一个处理器；以及，与至少一个处理器通信连接的存储器；其中，存储器存储有可被至少一个处理器执行的指令，指令被至少一个处理器执行，以使至少一个处理器能够执行上述实施方式中的卫星重力观测数据异常值提取方法。

本发明实施方式相对于现有技术而言，先通过小波连续函数降噪方法对卫星重力观测数据中的异常值进行探测，根据探测出的异常值前后数值的特征计算一个新的数值来替换探测出的异常值，再使用狄克逊检验法对异常值替换后的卫星重力观测数据进行进一步探测。通过将小波连续函数降噪方法与狄克逊检验法进行结合，实现了对卫星重力观测数据中异常值的准确探测。