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

权利要求书

1.一种卫星重力观测数据异常值提取方法，其特征在于，包括：

获取第一卫星重力观测数据；

采用小波连续函数降噪方法对所述第一卫星重力观测数据进行异常值探测；

根据所述第一卫星重力观测数据中位于所述异常值位置的前后数值计算一个新的数值，并将计算出的所述新的数值替换所述异常值，获得第二卫星重力观测数据；

采用狄克逊检验法对所述第二卫星重力观测数据进行异常值探测检验。

2.根据权利要求1所述的卫星重力观测数据异常值提取方法，其特征在于，所述采用小波连续函数降噪方法对所述第一卫星重力观测数据进行异常值探测，具体包括：

对所述第一卫星重力观测数据进行小波分解，得到所述第一卫星重力观测数据在各尺度上的小波系数；

对所述小波系数进行修改，其中，将大于预设最大阈值的所述小波系数进行保留，将小于预设最小阈值的所述小波系数进行无限接近于零值处理，将大于所述预设最小阈值同时小于所述预设最大阈值之间的所述小波系数，加上修正量之后进行保留，其中，所述修正量为：所述预设最小阈值和所述预设最大阈值差值的绝对值乘以修改之前的所述小波系数所得的结果；

根据修改后的所述小波系数进行小波系数合成，得到第三卫星重力观测数据；

通过计算所述第三卫星重力观测数据与所述第一卫星重力观测数据的差值得到残差信号；

根据所述残差信号对所述第一卫星重力观测数据中的所述异常值进行定位。

3.根据权利要求2所述的卫星重力观测数据异常值提取方法，其特征在于，所述对所述小波系数进行修改，具体包括：

利用阈值函数公式：计算得到修改后的所述小波系数；

其中，η表示修改后的所述小波系数，x表示修改之前的所述小波系数，λ_Max表示所述预设最大阈值：λ_Min表示所述预设最小阈值：其中，σ表示所述第一卫星重力观测数据的信号标准差，N表示所述第一卫星重力观测数据的信号长度，j表示所述第一卫星重力观测数据的分解尺度，参数a、b、m表示所述阈值函数的调节因子，m的取值在区间(0,1)之间。

4.根据权利要求1所述的卫星重力观测数据异常值提取方法，其特征在于，所述根据所述第一卫星重力观测数据中位于所述异常值位置的前后数值计算一个新的数值，具体包括：

利用差值公式：计算所述异常值对应的新的数值，

其中，i表示所述第一卫星重力观测数据中数值的序列号，y'_i表示所述第一卫星重力观测数据中序列号为i的数值对应的新的数值，y_i-1表示位于所述异常值y_i之前的第一个数值，y_i-2表示位于所述异常值y_i之前的第二个数值，y_i+1表示位于所述异常值y_i之后的第一个数值，y_i+2表示位于所述异常值y_i之后的第二个数值。

5.根据权利要求4所述的卫星重力观测数据异常值提取方法，其特征在于，所述将计算出的所述新的数值替换所述异常值，获得第二卫星重力观测数据具体包括：

将计算出的所述新的数值y'_i代替探测出的所述异常值y_i；

通过对所述第一卫星重力观测数据进行所述异常值y_i的替换后，获得所述第二卫星重力观测数据。