[发明专利]一种航空重力测量数据向下延拓方法及系统

权利要求书

1.一种航空重力测量数据向下延拓方法，其特征在于，所述方法包括：

获取各空中点的原始重力观测数据；

对所述原始重力观测数据进行预处理，得到格网化的空中点重力异常数据；

根据将格网化的空中点重力异常数据向地面延拓，得到地面重力异常数据，所述地面重力异常数据为所述空中点重力异常数据向地面的延拓，其中，为所述地面重力异常数据的第n+1次迭代，为所述地面重力异常数据的第n次迭代，为设定格网面元内的地面重力异常数据的第n次迭代，A为系数矩阵，所述系数矩阵的对角线元素为非对角线元素为ρ为空中点和地面点之间的距离，r为空中点的地心向径，θ为空中点地心余纬，λ为空中点的地心经度；R表示地面点的地心向径，θ'表示地面点的地心余纬，λ'表示地面点的地心经度；ψ表示地面点和空中点之间的球面角距，cosψ＝cosθ'cosθ+sinθ'sinθcos(λ′-λ)，ds_j为空中点对应的地面格网面元，ψ₀为设定积分半径，N_c为所述设定积分半径ψ₀中包含的地面点的个数。

2.根据权利要求1所述的航空重力测量数据向下延拓方法，其特征在于，所述对所述原始重力观测数据进行预处理，得到格网化的空中点重力异常数据，具体包括：

对所述原始重力观测数据进行当量转换、误差补偿和数据格式转换，得到第一预处理数据；

对所述第一预处理数据进行惯性导航解算、GNSS解算以及惯性/GNSS组合导航解算，得到第二预处理数据；

对所述第二预处理数据进行测线划分、采样率统一、测线扰动重力解算，得到第三预处理数据；

对所述第三预处理数据进行低通滤波、测线网平差和格网化处理，得到格网化的空中点重力异常数据。

3.根据权利要求1所述的航空重力测量数据向下延拓方法，其特征在于，所述根据将格网化的空中点重力异常数据向地面延拓，具体包括：

设定迭代初始条件

根据对所述地面重力异常数据进行迭代，直至相邻两次迭代得到的所述地面重力异常数据之间的差值小于设定阈值时，停止迭代。

4.根据权利要求1所述的航空重力测量数据向下延拓方法，其特征在于，所述方法还包括：

对所述地面重力异常数据进行精度评估。

5.根据权利要求4所述的航空重力测量数据向下延拓方法，其特征在于，所述对所述地面重力异常数据进行精度评估，具体包括：

获取地面重力异常数据；

获取所述地面重力异常数据对应的延拓面的已知重力数据；

根据所述已知重力数据对所述地面重力异常数据进行精度评估。

6.一种航空重力测量数据向下延拓系统，其特征在于，所述系统包括：

空中点原始重力观测数据获取单元，用于获取各空中点的空中重力异常观测数据；

预处理单元，用于对所述空中重力异常观测数据进行预处理，得到格网化的空中点重力异常数据；

延拓单元，用于根据将格网化的空中点重力异常数据向地面延拓，得到地面重力异常数据，所述地面重力异常数据为所述空中点重力异常数据向地面的延拓，其中，为所述地面重力异常数据的第n+1次迭代，为所述地面重力异常数据的第n次迭代，为设定格网面元内的地面重力异常数据的第n次迭代，A为系数矩阵，所述系数矩阵的对角线元素为非对角线元素为ρ为空中点和地面点之间的距离，r为空中点的地心向径，θ为空中点地心余纬，λ为空中点的地心经度；R表示地面点的地心向径，θ'表示地面点的地心余纬，λ'表示地面点的地心经度；ψ表示地面点和空中点之间的球面角距，cosψ＝cosθ'cosθ+sinθ'sinθcos(λ′-λ)，ds_j为空中点对应的地面格网面元，ψ₀为设定积分半径，N_c为所述设定积分半径ψ₀中包含的地面点的个数。