CN109283591B - 以地面点作为控制的航空重力数据向下延拓方法及系统 - Google Patents

以地面点作为控制的航空重力数据向下延拓方法及系统 Download PDF

Info

Publication number
CN109283591B
CN109283591B CN201811228354.6A CN201811228354A CN109283591B CN 109283591 B CN109283591 B CN 109283591B CN 201811228354 A CN201811228354 A CN 201811228354A CN 109283591 B CN109283591 B CN 109283591B
Authority
CN
China
Prior art keywords
gravity
data
ground
aerial
point
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201811228354.6A
Other languages
English (en)
Other versions
CN109283591A (zh
Inventor
刘晓刚
翟振和
秦显平
管斌
马健
段渭超
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chinese People's Liberation Army 61540
Original Assignee
Chinese People's Liberation Army 61540
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Chinese People's Liberation Army 61540 filed Critical Chinese People's Liberation Army 61540
Priority to CN201811228354.6A priority Critical patent/CN109283591B/zh
Publication of CN109283591A publication Critical patent/CN109283591A/zh
Application granted granted Critical
Publication of CN109283591B publication Critical patent/CN109283591B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V7/00—Measuring gravitational fields or waves; Gravimetric prospecting or detecting
    • G01V7/16—Measuring gravitational fields or waves; Gravimetric prospecting or detecting specially adapted for use on moving platforms, e.g. ship, aircraft

Landscapes

  • Engineering & Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Ocean & Marine Engineering (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种以地面点作为控制的航空重力数据向下延拓方法及系统。所述评估方法包括：根据航空重力测量系统获取空中格网重力测量数据，并根据地面重力测量系统获取地面重力点测量数据；根据空中格网重力测量数据以及地面重力点测量数据，采用泊松Poisson积分方程构建联合地面重力测量点的航空重力数据向下延拓的泊松Poisson积分法观测模型；判断泊松Poisson积分法观测模型是否存在病态性，若是，采用正则化技术解算观测方程，确定第一延拓结果；根据第一延拓结果对延拓面地面重力数据进行评估；若否，对矩阵进行求逆计算，确定第二延拓结果；根据第二延拓结果对延拓面地面重力数据进行评估。采用本发明的评估方法及系统能够提高最终延拓结果精度，有效控制延拓结果的稳定性。

Description

以地面点作为控制的航空重力数据向下延拓方法及系统
技术领域
本发明涉及航空重力测量领域，特别是涉及一种以地面点作为控制的航空重力数据向下延拓方法及系统。
背景技术
航空重力测量技术，因其可以在沙漠、沼泽、冰川、丘陵、原始森林、陆海交界等一些难以开展地面重力测量的区域进行作业，快速经济地获取高精度高分辨率地球重力场中高频信息，从而成为获取地球重力场测量数据最有效的技术手段之一。
航空重力测量数据在应用时，一般需要延拓到地球表面或大地水准面上，主要用于不同类型重力测量数据的融合、全球或区域地球重力场模型的构建、全球或区域(似)大地水准面的精化、水下重力匹配辅助导航中重力基准图的生成等；以地面的重力测量点作为控制，对航空重力测量数据进行向下延拓，可以减弱向下延拓的不稳定性，抑制边界效应对延拓结果的影响，提高延拓结果精度；目前常用的联合地面重力测量点的航空重力数据向下延拓的方法主要是最小二乘配置法。
最小二乘配置法的基本原理是利用少量的高质量地面重力测量点联合航空重力数据来推估未测地面重力点，从而实现将航空重力测量数据延拓到地面的目的；优点是相比单一的航空重力测量数据直接延拓来说，可以提高延拓结果精度，并且能给出延拓结果的精度估计信息；缺点是地面重力测量点和航空重力数据的协方差函数构建难度较大，并且构建效果将对最终延拓结果(未测地面重力点的地面重力数据)的精度造成直接影响，无法有效控制延拓结果的稳定性。
发明内容
本发明的目的是提供一种以地面点作为控制的航空重力数据向下延拓方法及系统，以解决现有航空重力数据向下延拓的最小二乘配置法评估未测地面重力点的地面重力数据时的评估精度低，无法有效控制延拓结果稳定性的问题。
为实现上述目的，本发明提供了如下方案：
一种以地面点作为控制的航空重力数据向下延拓方法，所述以地面点作为控制的航空重力数据向下延拓方法基于航空重力测量系统以及地面重力测量系统，包括：
根据所述航空重力测量系统获取空中格网重力测量数据，并根据所述地面重力测量系统获取地面重力点测量数据；
根据所述空中格网重力测量数据以及所述地面重力点测量数据，采用泊松Poisson积分方程构建联合地面重力测量点的航空重力数据向下延拓的泊松Poisson积分法观测模型；
判断所述泊松Poisson积分法观测模型是否存在病态性，得到第一判断结果；
若所述第一判断结果表示为所述泊松Poisson积分法观测模型存在病态性，采用正则化技术解算观测方程，确定第一延拓结果；
根据所述第一延拓结果对延拓面地面重力数据进行评估；
若所述第一判断结果表示为所述泊松Poisson积分法观测模型不存在病态性，对矩阵进行求逆计算，确定第二延拓结果；
根据所述第二延拓结果对延拓面地面重力数据进行评估。
可选的，所述根据所述空中格网重力测量数据以及所述地面重力点测量数据，采用泊松Poisson积分方程构建联合地面重力测量点的航空重力数据向下延拓的泊松Poisson积分法观测模型，具体包括：
利用重力场模型确定所述重力场模型之外的地面区域重力数据对所述空中格网重力测量数据的第一空中重力测量影响数据；
滤除所述空中格网重力测量数据内与所述第一空中重力测量影响数据相同的空中网格重力测量数据，确定空中残差格网重力数据；
根据所述地面重力点测量数据确定空中重力测量数据；
根据所述空中重力测量数据，采用泊松Poisson积分方程确定第二空中重力测量影响数据；
滤除所述空中残差格网重力数据内与所述第二空中重力测量影响数据相同的空中网格重力测量数据，确定滤除后的空中重力测量数据；
根据所述滤除后的空中重力测量数据构建联合地面重力测量点的航空重力数据向下延拓的泊松Poisson积分法观测模型。
可选的，所述根据所述滤除后的空中重力测量数据构建联合地面重力测量点的航空重力数据向下延拓的泊松Poisson积分法观测模型，具体包括：
利用公式
Figure BDA0001836527800000031
构建联合地面重力测量点的航空重力数据向下延拓的泊松Poisson积分法观测模型；其中，
Figure BDA0001836527800000032
为空中重力测量数据；R1为地面平均高程面的球面半径；r为空中点的地心向径；N为地面重力点的个数；Bj为地面点j的纬度；Lj为地面点j的经度；
Figure BDA0001836527800000033
ΔBj为纬度的格网分辨率，ΔLj为经度的格网分辨率；
Figure BDA0001836527800000034
为地面重力点测量数据；
Figure BDA0001836527800000035
为地面点j和空中点i之间的距离。
可选的，所述采用正则化技术解算观测方程，确定第一延拓结果，具体包括：
根据所述泊松Poisson积分法观测模型建立正则化法观测方程；
对所述正则化法观测方程进行谱分解，确定矩阵谱分解方程；
利用曲线法，确定所述矩阵谱分解方程的最优正则化参数；所述曲线法包括L曲线法以及U曲线法；
根据所述最优正则化参数以及所述正则化法观测方程确定第一延拓结果。
一种以地面点作为控制的航空重力数据向下延拓系统，所述以地面点作为控制的航空重力数据向下延拓方法基于航空重力测量系统以及地面重力测量系统，包括：
测量数据获取模块，用于根据所述航空重力测量系统获取空中格网重力测量数据，并根据所述地面重力测量系统获取地面重力点测量数据；
泊松Poisson积分法观测模型构建模块，用于根据所述空中格网重力测量数据以及所述地面重力点测量数据，采用泊松Poisson积分方程构建联合地面重力测量点的航空重力数据向下延拓的泊松Poisson积分法观测模型；
第一判断模块，用于判断所述泊松Poisson积分法观测模型是否存在病态性，得到第一判断结果；
第一延拓结果确定模块，用于若所述第一判断结果表示为所述泊松Poisson积分法观测模型存在病态性，采用正则化技术解算观测方程，确定第一延拓结果；
第一评估模块，用于根据所述第一延拓结果对延拓面地面重力数据进行评估；
第二延拓结果确定模块，用于若所述第一判断结果表示为所述泊松Poisson积分法观测模型不存在病态性，对矩阵进行求逆计算，确定第二延拓结果；
第二评估模块，用于根据所述第二延拓结果对延拓面地面重力数据进行评估。
可选的，所述泊松Poisson积分法观测模型确定模块具体包括：
第一空中重力测量影响数据确定单元，用于利用重力场模型确定所述重力场模型之外的地面区域重力数据对所述空中格网重力测量数据的第一空中重力测量影响数据；
空中残差格网重力数据确定单元，用于滤除所述空中格网重力测量数据内与所述第一空中重力测量影响数据相同的空中网格重力测量数据，确定空中残差格网重力数据；
空中重力测量数据确定单元，用于根据所述地面重力点测量数据确定空中重力测量数据；
第二空中重力测量影响数据确定单元，用于根据所述空中重力测量数据，采用泊松Poisson积分方程确定第二空中重力测量影响数据；
滤除后的空中重力测量数据确定单元，用于滤除所述空中残差格网重力数据内与所述第二空中重力测量影响数据相同的空中网格重力测量数据，确定滤除后的空中重力测量数据；
泊松Poisson积分法观测模型确定单元，用于根据所述滤除后的空中重力测量数据构建联合地面重力测量点的航空重力数据向下延拓的泊松Poisson积分法观测模型。
可选的，所述泊松Poisson积分法观测模型确定单元具体包括：
泊松Poisson积分法观测模型确定子单元，用于利用公式
Figure BDA0001836527800000051
构建联合地面重力测量点的航空重力数据向下延拓的泊松Poisson积分法观测模型；其中，
Figure BDA0001836527800000052
为空中重力测量数据；R1为地面平均高程面的球面半径；r为空中点的地心向径；N为地面重力点的个数；Bj为地面点j的纬度；Lj为地面点j的经度；
Figure BDA0001836527800000053
ΔBj为纬度的格网分辨率，ΔLj为经度的格网分辨率；
Figure BDA0001836527800000054
为地面重力点测量数据；
Figure BDA0001836527800000055
为地面点j和空中点i之间的距离。
可选的，所述第一延拓结果确定模块具体包括：
正则化法观测方程建立单元，用于根据所述泊松Poisson积分法观测模型建立正则化法观测方程；
矩阵谱分解方程确定单元，用于对所述正则化法观测方程进行谱分解，确定矩阵谱分解方程；
最优正则化参数确定单元，用于利用曲线法，确定所述矩阵谱分解方程的最优正则化参数；所述曲线法包括L曲线法以及U曲线法；
第一延拓结果确定单元，用于根据所述最优正则化参数以及所述正则化法观测方程确定第一延拓结果。
根据本发明提供的具体实施例，本发明公开了以下技术效果：采用本发明所提供的以地面点作为控制的航空重力数据向下延拓方法及系统不涉及到航空重力测量数据和地面重力测量点的协方差函数构建，大大节约了工作量；以地面的重力测量点作为辅助控制，能够抑制边界效应对延拓结果的影响，减弱向下延拓的不稳定性，提高延拓结果精度；同时能够削弱航空重力测量数据中的系统误差的影响。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明所提供的联合地面重力测量点的以地面点作为控制的航空重力数据向下延拓系统组成图；
图2为本发明所提供的以地面点作为控制的航空重力数据向下延拓方法流程图；
图3为本发明所提供的另一种以地面点作为控制的航空重力数据向下延拓方法流程图；
图4为本发明所提供的正则化技术解算观测方程流程图；
图5为本发明所提供的以地面点作为控制的航空重力数据向下延拓系统结构图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明的目的是提供一种以地面点作为控制的航空重力数据向下延拓方法及系统，能够提高最终延拓结果(未测地面重力点的地面重力数据)的精度，有效控制延拓结果的稳定性。
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1为本发明所提供的联合地面重力测量点的以地面点作为控制的航空重力数据向下延拓系统组成图，如图1所示，航空重力测量系统由重力传感器分系统、全球导航定位系统(Global Navigation Satellite System，GNSS)定位分系统、姿态分系统、激光测高分系统、数据采集分系统等组成；本发明利用航空重力测量系统获得的空中格网点数据，进而采用联合地面重力测量点的航空重力数据向下延拓的点质量法，将空中格网点数据向下延拓到地球表面或大地水准面上，并将该数据作下一步的应用。
图2为本发明所提供的以地面点作为控制的航空重力数据向下延拓方法流程图，如图2所示，一种以地面点作为控制的航空重力数据向下延拓方法，所述以地面点作为控制的航空重力数据向下延拓方法基于航空重力测量系统以及地面重力测量系统，包括：
步骤201：根据所述航空重力测量系统获取空中格网重力测量数据，并根据所述地面重力测量系统获取地面重力点测量数据。
步骤202：根据所述空中格网重力测量数据以及所述地面重力点测量数据，采用泊松Poisson积分方程构建联合地面重力测量点的航空重力数据向下延拓的泊松Poisson积分法观测模型。
在球近似下，扰动重力与扰动位有以下关系：
Figure BDA0001836527800000071
其中，r是地心向径。
假设航空重力测量是在空中飞行平面σ上，地面Σ上进行了部分重力测量，通常在σ覆盖的区域，是重力空白区，或是测点较少精度较弱的地区。我们希望将空中扰动引力延拓到该区域的地面平均高程面Ω上，这样更有利于以后的应用。可以将地面平均高程面Ω视为半径为R1的球面，在这个面上的扰动重力为
Figure BDA0001836527800000072
由式关系知，
Figure BDA0001836527800000073
为调和函数，因此可以采用Poisson积分由该球面上的扰动重力求得球外一点的值：
Figure BDA0001836527800000074
其中l是空中测点到Ω上积分面元的空间距离，
Figure BDA0001836527800000075
首先设测量区域σ覆盖的地面无重力点，在其外围有重力数据(通常为平均扰动重力)。式为一积分方程，待求的是积分内对应于σ的Ω'上的扰动重力。因为积分是全球的，需要做一改化。对于扰动重力，令
Figure BDA0001836527800000081
其中，
Figure BDA0001836527800000082
为已知球谐重力位模型计算部分：
Figure BDA0001836527800000083
Figure BDA0001836527800000084
为实际空中测量与地面已知扰动重力扣除了模型扰动重力之后的剩余扰动重力应用式，将
Figure BDA0001836527800000085
分为两部分：
其中：
Figure BDA0001836527800000086
Figure BDA0001836527800000087
角标i表示在空中观测点处，
Figure BDA0001836527800000088
为测量区以外区域的积分，由于该区域有重力资料，可以用数值积分方法计算得到，此时积分只需选择适当范围则可。虽然已知的平均扰动重力是在地面上，但由于距离较远，可以视为平均高程面上的值。将式的积分离散化，从而得到观测方程：
Figure BDA0001836527800000089
其中，
Figure BDA00018365278000000810
可以写成矩阵形式：G＝Al。
其中G是M维空中扰动重力观测量向量，l是N维待求的平均高程面上扰动重力向量，A是M×N矩阵。当M＞N时，可采用最小二乘方法解得l：l＝(ATA)-1ATG。
当测区范围较大时，式存在着病态性，需要采用正则化处理技术来削弱延拓模型病态性的影响，提高延拓结果的精度。
所述步骤202具体包括：利用重力场模型确定所述重力场模型之外的地面区域重力数据对所述空中格网重力测量数据的第一空中重力测量影响数据；滤除所述空中格网重力测量数据内与所述第一空中重力测量影响数据相同的空中网格重力测量数据，确定空中残差格网重力数据；根据所述地面重力点测量数据确定空中重力测量数据；根据所述空中重力测量数据，采用泊松Poisson积分方程确定第二空中重力测量影响数据；滤除所述空中残差格网重力数据内与所述第二空中重力测量影响数据相同的空中网格重力测量数据，确定滤除后的空中重力测量数据；根据所述滤除后的空中重力测量数据构建联合地面重力测量点的航空重力数据向下延拓的泊松Poisson积分法观测模型。
所述根据所述滤除后的空中重力测量数据构建联合地面重力测量点的航空重力数据向下延拓的泊松Poisson积分法观测模型，具体包括：利用公式
Figure BDA0001836527800000091
构建联合地面重力测量点的航空重力数据向下延拓的泊松Poisson积分法观测模型；其中，
Figure BDA0001836527800000092
为空中重力测量数据；R1为地面平均高程面的球面半径；r为空中点的地心向径；N为地面重力点的个数；Bj为地面点j的纬度；Lj为地面点j的经度；
Figure BDA0001836527800000093
ΔBj为纬度的格网分辨率，ΔLj为经度的格网分辨率；
Figure BDA0001836527800000094
为地面重力点测量数据；
Figure BDA0001836527800000095
为地面点j和空中点i之间的距离。
步骤203：判断所述泊松Poisson积分法观测模型是否存在病态性，若是，执行步骤204；若否，执行步骤206。
步骤204：采用正则化技术解算观测方程，确定第一延拓结果。
所述采用正则化技术解算观测方程，确定第一延拓结果，具体包括：根据所述泊松Poisson积分法观测模型建立正则化法观测方程；对所述正则化法观测方程进行谱分解，确定矩阵谱分解方程；利用曲线法，确定所述矩阵谱分解方程的最优正则化参数；所述曲线法包括L曲线法以及U曲线法；根据所述最优正则化参数以及所述正则化法观测方程确定第一延拓结果。
步骤205：根据所述第一延拓结果对延拓面地面重力数据进行评估。
步骤206：对矩阵进行求逆计算，确定第二延拓结果。
步骤207：根据所述第二延拓结果对延拓面地面重力数据进行评估。
根据本发明所提供的以地面点作为控制的航空重力数据向下延拓方法还可以延伸出以下评估方法，可分为以下步骤：
第一步：通过航空重力测量系统和地面重力测量仪分别获得空中格网重力测量数据和地面重力点测量数据。
第二步利用已知重力场模型计算搜索半径以外的地面区域重力数据对航空重力测量数据的影响，并从航空重力测量数据中扣除，获得空中残差格网重力数据。
第三步：利用地面重力测量点，采用Poisson积分方程计算空中重力测量数据，获得地面重力点对空中重力测量数据的影响。
第四步：从空中残差格网重力数据中减去地面重力点对空中重力测量数据的影响值。
第五步：构建联合地面重力测量点的航空重力数据向下延拓的Poisson积分法观测方程。
第六步：对经过一系列处理后的空中重力测量数据进行向下延拓。
第七步：判断观测方程是否存在病态性，如果是，则采用正则化技术解算观测方程，获得延拓结果，如果不存在病态性，则直接进行矩阵求逆获得延拓结果。
第八步：对延拓面重力数据进行精度评估，获得精度评估结果；
具体评估流程如图3所示，其中，第七步中的正则化技术解算观测方程的具体实现过程如图4所示：
根据联合地面重力测量点的航空重力数据向下延拓的Poisson积分法观测方程，建立正则化法观测方程；对正则化法观测方程进行谱分解，得到矩阵的谱分解方程，可以提高计算速度；根据L曲线法、U曲线法等方法，确定正则化方程的最优正则化参数；将最优正则化参数代入正则化法观测方程，得到矩阵计算结果。
图5为本发明所提供的以地面点作为控制的航空重力数据向下延拓系统结构图，一种以地面点作为控制的航空重力数据向下延拓系统，所述以地面点作为控制的航空重力数据向下延拓方法基于航空重力测量系统以及地面重力测量系统，包括：
测量数据获取模块501，用于根据所述航空重力测量系统获取空中格网重力测量数据，并根据所述地面重力测量系统获取地面重力点测量数据。
泊松Poisson积分法观测模型构建模块502，用于根据所述空中格网重力测量数据以及所述地面重力点测量数据，采用泊松Poisson积分方程构建联合地面重力测量点的航空重力数据向下延拓的泊松Poisson积分法观测模型。
所述泊松Poisson积分法观测模型确定模块502具体包括：第一空中重力测量影响数据确定单元，用于利用重力场模型确定所述重力场模型之外的地面区域重力数据对所述空中格网重力测量数据的第一空中重力测量影响数据；空中残差格网重力数据确定单元，用于滤除所述空中格网重力测量数据内与所述第一空中重力测量影响数据相同的空中网格重力测量数据，确定空中残差格网重力数据；空中重力测量数据确定单元，用于根据所述地面重力点测量数据确定空中重力测量数据；第二空中重力测量影响数据确定单元，用于根据所述空中重力测量数据，采用泊松Poisson积分方程确定第二空中重力测量影响数据；滤除后的空中重力测量数据确定单元，用于滤除所述空中残差格网重力数据内与所述第二空中重力测量影响数据相同的空中网格重力测量数据，确定滤除后的空中重力测量数据；泊松Poisson积分法观测模型确定单元，用于根据所述滤除后的空中重力测量数据构建联合地面重力测量点的航空重力数据向下延拓的泊松Poisson积分法观测模型。
所述泊松Poisson积分法观测模型确定单元具体包括：泊松Poisson积分法观测模型确定子单元，用于利用公式
Figure BDA0001836527800000111
构建联合地面重力测量点的航空重力数据向下延拓的泊松Poisson积分法观测模型；其中，
Figure BDA0001836527800000112
为空中重力测量数据；R1为地面平均高程面的球面半径；r为空中点的地心向径；N为地面重力点的个数；Bj为地面点j的纬度；Lj为地面点j的经度；
Figure BDA0001836527800000113
ΔBj为纬度的格网分辨率，ΔLj为经度的格网分辨率；
Figure BDA0001836527800000114
为地面重力点测量数据；
Figure BDA0001836527800000115
为地面点j和空中点i之间的距离。
第一判断模块503，用于判断所述泊松Poisson积分法观测模型是否存在病态性，得到第一判断结果。
第一延拓结果确定模块504，用于若所述第一判断结果表示为所述泊松Poisson积分法观测模型存在病态性，采用正则化技术解算观测方程，确定第一延拓结果。
所述第一延拓结果确定模块504具体包括：正则化法观测方程建立单元，用于根据所述泊松Poisson积分法观测模型建立正则化法观测方程；矩阵谱分解方程确定单元，用于对所述正则化法观测方程进行谱分解，确定矩阵谱分解方程；最优正则化参数确定单元，用于利用曲线法，确定所述矩阵谱分解方程的最优正则化参数；所述曲线法包括L曲线法以及U曲线法；第一延拓结果确定单元，用于根据所述最优正则化参数以及所述正则化法观测方程确定第一延拓结果。
第一评估模块505，用于根据所述第一延拓结果对延拓面地面重力数据进行评估。
第二延拓结果确定模块506，用于若所述第一判断结果表示为所述泊松Poisson积分法观测模型不存在病态性，对矩阵进行求逆计算，确定第二延拓结果。
第二评估模块507，用于根据所述第二延拓结果对延拓面地面重力数据进行评估。
由于本发明所提供的以地面点作为控制的航空重力数据向下延拓方法及系统不涉及到协方差函数的构建，可以对延拓结果起到有效的控制作用，抑制边界效应对延拓结果的影响，改善向下延拓的效果。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。

Claims (6)

1.一种以地面点作为控制的航空重力数据向下延拓方法，其特征在于，所述以地面点作为控制的航空重力数据向下延拓方法基于航空重力测量系统以及地面重力测量系统，包括：
根据所述航空重力测量系统获取空中格网重力测量数据，并根据所述地面重力测量系统获取地面重力点测量数据；
根据所述空中格网重力测量数据以及所述地面重力点测量数据，采用泊松Poisson积分方程构建联合地面重力测量点的航空重力数据向下延拓的泊松Poisson积分法观测模型；所述根据所述空中格网重力测量数据以及所述地面重力点测量数据，采用泊松Poisson积分方程构建联合地面重力测量点的航空重力数据向下延拓的泊松Poisson积分法观测模型，具体包括：利用重力场模型确定搜索半径之外的地面区域重力数据对所述空中格网重力测量数据的第一空中重力测量影响数据；滤除所述空中格网重力测量数据内与所述第一空中重力测量影响数据相同的空中网格重力测量数据，确定空中残差格网重力数据；根据所述地面重力点测量数据确定空中重力测量数据；根据所述空中重力测量数据，采用泊松Poisson积分方程确定第二空中重力测量影响数据；滤除所述空中残差格网重力数据内与所述第二空中重力测量影响数据相同的空中网格重力测量数据，确定滤除后的空中重力测量数据；根据所述滤除后的空中重力测量数据构建联合地面重力测量点的航空重力数据向下延拓的泊松Poisson积分法观测模型；
判断所述泊松Poisson积分法观测模型是否存在病态性，得到第一判断结果；
若所述第一判断结果表示为所述泊松Poisson积分法观测模型存在病态性，采用正则化技术解算观测方程，确定第一延拓结果；
根据所述第一延拓结果对延拓面地面重力数据进行评估；
若所述第一判断结果表示为所述泊松Poisson积分法观测模型不存在病态性，对矩阵进行求逆计算，确定第二延拓结果；
根据所述第二延拓结果对延拓面地面重力数据进行评估。
2.根据权利要求1所述的以地面点作为控制的航空重力数据向下延拓方法，其特征在于，所述根据所述滤除后的空中重力测量数据构建联合地面重力测量点的航空重力数据向下延拓的泊松Poisson积分法观测模型，具体包括：
利用公式
Figure FDA0002284371910000021
构建联合地面重力测量点的航空重力数据向下延拓的泊松Poisson积分法观测模型；其中，
Figure FDA0002284371910000022
为空中重力测量数据；R1为地面平均高程面的球面半径；r为空中点的地心向径；N为地面重力点的个数；Bj为地面点j的纬度；Lj为地面点j的经度；
Figure FDA0002284371910000023
ΔBj为纬度的格网分辨率，ΔLj为经度的格网分辨率；
Figure FDA0002284371910000024
为地面重力点测量数据；
Figure FDA0002284371910000025
为地面点j和空中点i之间的距离。
3.根据权利要求2所述的以地面点作为控制的航空重力数据向下延拓方法，其特征在于，所述采用正则化技术解算观测方程，确定第一延拓结果，具体包括：
根据所述泊松Poisson积分法观测模型建立正则化法观测方程；
对所述正则化法观测方程进行谱分解，确定矩阵谱分解方程；
利用曲线法，确定所述矩阵谱分解方程的最优正则化参数；所述曲线法包括L曲线法以及U曲线法；
根据所述最优正则化参数以及所述正则化法观测方程确定第一延拓结果。
4.一种以地面点作为控制的航空重力数据向下延拓系统，其特征在于，所述以地面点作为控制的航空重力数据向下延拓系统基于航空重力测量系统以及地面重力测量系统，包括：
测量数据获取模块，用于根据所述航空重力测量系统获取空中格网重力测量数据，并根据所述地面重力测量系统获取地面重力点测量数据；
泊松Poisson积分法观测模型构建模块，用于根据所述空中格网重力测量数据以及所述地面重力点测量数据，采用泊松Poisson积分方程构建联合地面重力测量点的航空重力数据向下延拓的泊松Poisson积分法观测模型；所述泊松Poisson积分法观测模型确定模块具体包括：第一空中重力测量影响数据确定单元，用于利用重力场模型确定搜索半径之外的地面区域重力数据对所述空中格网重力测量数据的第一空中重力测量影响数据；空中残差格网重力数据确定单元，用于滤除所述空中格网重力测量数据内与所述第一空中重力测量影响数据相同的空中网格重力测量数据，确定空中残差格网重力数据；空中重力测量数据确定单元，用于根据所述地面重力点测量数据确定空中重力测量数据；第二空中重力测量影响数据确定单元，用于根据所述空中重力测量数据，采用泊松Poisson积分方程确定第二空中重力测量影响数据；滤除后的空中重力测量数据确定单元，用于滤除所述空中残差格网重力数据内与所述第二空中重力测量影响数据相同的空中网格重力测量数据，确定滤除后的空中重力测量数据；泊松Poisson积分法观测模型确定单元，用于根据所述滤除后的空中重力测量数据构建联合地面重力测量点的航空重力数据向下延拓的泊松Poisson积分法观测模型；
第一判断模块，用于判断所述泊松Poisson积分法观测模型是否存在病态性，得到第一判断结果；
第一延拓结果确定模块，用于若所述第一判断结果表示为所述泊松Poisson积分法观测模型存在病态性，采用正则化技术解算观测方程，确定第一延拓结果；
第一评估模块，用于根据所述第一延拓结果对延拓面地面重力数据进行评估；
第二延拓结果确定模块，用于若所述第一判断结果表示为所述泊松Poisson积分法观测模型不存在病态性，对矩阵进行求逆计算，确定第二延拓结果；
第二评估模块，用于根据所述第二延拓结果对延拓面地面重力数据进行评估。
5.根据权利要求4所述的以地面点作为控制的航空重力数据向下延拓系统，其特征在于，所述泊松Poisson积分法观测模型确定单元具体包括：
泊松Poisson积分法观测模型确定子单元，用于利用公式
Figure FDA0002284371910000041
构建联合地面重力测量点的航空重力数据向下延拓的泊松Poisson积分法观测模型；其中，
Figure FDA0002284371910000042
为空中重力测量数据；R1为地面平均高程面的球面半径；r为空中点的地心向径；N为地面重力点的个数；Bj为地面点j的纬度；Lj为地面点j的经度；
Figure FDA0002284371910000043
ΔBj为纬度的格网分辨率，ΔLj为经度的格网分辨率；
Figure FDA0002284371910000044
为地面重力点测量数据；
Figure FDA0002284371910000045
为地面点j和空中点i之间的距离。
6.根据权利要求5所述的以地面点作为控制的航空重力数据向下延拓系统，其特征在于，所述第一延拓结果确定模块具体包括：
正则化法观测方程建立单元，用于根据所述泊松Poisson积分法观测模型建立正则化法观测方程；
矩阵谱分解方程确定单元，用于对所述正则化法观测方程进行谱分解，确定矩阵谱分解方程；
最优正则化参数确定单元，用于利用曲线法，确定所述矩阵谱分解方程的最优正则化参数；所述曲线法包括L曲线法以及U曲线法；
第一延拓结果确定单元，用于根据所述最优正则化参数以及所述正则化法观测方程确定第一延拓结果。
CN201811228354.6A 2018-10-22 2018-10-22 以地面点作为控制的航空重力数据向下延拓方法及系统 Active CN109283591B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811228354.6A CN109283591B (zh) 2018-10-22 2018-10-22 以地面点作为控制的航空重力数据向下延拓方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811228354.6A CN109283591B (zh) 2018-10-22 2018-10-22 以地面点作为控制的航空重力数据向下延拓方法及系统

Publications (2)

Publication Number Publication Date
CN109283591A CN109283591A (zh) 2019-01-29
CN109283591B true CN109283591B (zh) 2020-03-13

Family

ID=65178241

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811228354.6A Active CN109283591B (zh) 2018-10-22 2018-10-22 以地面点作为控制的航空重力数据向下延拓方法及系统

Country Status (1)

Country Link
CN (1) CN109283591B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109856691B (zh) * 2019-03-14 2020-08-14 中国人民解放军61540部队 一种基于梯度法的航空重力矢量向下延拓方法及系统
CN110161582B (zh) * 2019-05-24 2020-05-12 中国地质科学院 空中与地面数据结合的重力换算方法及系统
CN111650659A (zh) * 2020-06-17 2020-09-11 宁夏大学 一种基于空中重力异常的垂线偏差估计方法及系统
CN112881743B (zh) * 2021-01-13 2023-05-02 徐培亮 一种基于正则化算法的速度和加速度计算方法及测量装置

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2009137931A1 (en) * 2008-05-14 2009-11-19 Geotech Airborne Limited Airborne geophysical survey using airship
WO2017135183A1 (ja) * 2016-02-01 2017-08-10 滋樹 水谷 重力偏差データの表層密度値推定方法
CN108415879A (zh) * 2018-01-19 2018-08-17 中国人民解放军92859部队 基于向上延拓的航空重力最小二乘向下延拓解析算法
CN108594319A (zh) * 2018-05-11 2018-09-28 中国人民解放军61540部队 一种航空重力测量数据向下延拓方法及系统

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2009137931A1 (en) * 2008-05-14 2009-11-19 Geotech Airborne Limited Airborne geophysical survey using airship
WO2017135183A1 (ja) * 2016-02-01 2017-08-10 滋樹 水谷 重力偏差データの表層密度値推定方法
CN108415879A (zh) * 2018-01-19 2018-08-17 中国人民解放军92859部队 基于向上延拓的航空重力最小二乘向下延拓解析算法
CN108594319A (zh) * 2018-05-11 2018-09-28 中国人民解放军61540部队 一种航空重力测量数据向下延拓方法及系统

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
航空重力测量数据向下延拓方法比较;王兴涛 等;《地球物理学报》;20041130;第47卷(第6期);第1017-1022页 *
航空重力测量数据向下延拓模型的病态性分析;刘晓刚 等;《第四届高分辨率对地观测学术年会论文集》;20170917;第1358-1373页 *
航空重力测量数据向下延拓的改进Poisson积分迭代法;刘晓刚 等;《测绘学报》;20180930;第47卷(第9期);第1185-1195页 *

Also Published As

Publication number Publication date
CN109283591A (zh) 2019-01-29

Similar Documents

Publication Publication Date Title
CN109283591B (zh) 以地面点作为控制的航空重力数据向下延拓方法及系统
CN103900565B (zh) 一种基于差分gps的惯导系统姿态获取方法
Cai et al. The first airborne scalar gravimetry system based on SINS/DGPS in China
CN109507706B (zh) 一种gps信号丢失的预测定位方法
CN104535080B (zh) 大方位失准角下基于误差四元数的传递对准方法
CN107728182A (zh) 基于相机辅助的柔性多基线测量方法和装置
CN104075699A (zh) 三维固态电子罗盘及其传感器的零点和比例系数核正方法
CN113916222B (zh) 基于卡尔曼滤波估计方差约束的组合导航方法
CN109856691B (zh) 一种基于梯度法的航空重力矢量向下延拓方法及系统
CN111595345A (zh) 一种基于多维重力梯度灯塔的潜艇导航方法及系统
CN110017812A (zh) 一种垂线偏差、重力加速度和重力梯度的测量装置及方法
CN115439469B (zh) 基于无人机的建筑缺陷检测方法、装置及电子设备
CN118067157B (zh) 一种惯性测量单元性能评估方法、装置、设备及介质
CN112577518A (zh) 一种惯性测量单元标定方法及装置
CN103983274B (zh) 一种适用于低精度无方位基准双轴转位设备的惯性测量单元标定方法
CN109581523B (zh) 一种采用卫星跟踪卫星装置对加速度计校准的方法和系统
Lu et al. Improved calibration of IMU biases in analytic coarse alignment for AHRS
CN113447993A (zh) 磁力矢量测量的补偿飞行方法、系统及磁补偿方法、系统
CN113670266A (zh) 利用无人机倾斜摄影开展不动产权籍测量技术
Wen et al. Study on the key technology and application of UAV surveying and mapping data processing
CN103064128B (zh) 基于星间距离误差模型的地球重力场恢复方法
CN107702718B (zh) 一种基于瞬间可观测度模型的机载pos机动优化方法与装置
CN111174974B (zh) 一种车辆悬架升沉测量方法及系统
CN104655132A (zh) 一种基于加速度计的机体弹性变形角估计方法
CN108919371B (zh) 一种联合地面重力点的航空重力数据向下延拓方法及系统

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant