CN107289943A - 构建水下重力辅助导航背景场的Shepard插值算法 - Google Patents
构建水下重力辅助导航背景场的Shepard插值算法 Download PDFInfo
- Publication number
- CN107289943A CN107289943A CN201710475567.8A CN201710475567A CN107289943A CN 107289943 A CN107289943 A CN 107289943A CN 201710475567 A CN201710475567 A CN 201710475567A CN 107289943 A CN107289943 A CN 107289943A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- gravity
- mfrac
- gravity anomaly
- 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.)
- Withdrawn
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/20—Instruments for performing navigational calculations
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Automation & Control Theory (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明涉及一种构建水下重力辅助导航背景场的Shepard插值算法，属于顾及移去恢复技术构建水下重力辅助导航背景场的电子信息技术领域。包括以海洋重力观测值为基础数据，联合全球重力场模型，采用顾及移去恢复技术的Shepard插值算法，构建满足水下重力辅助导航需求的海洋重力背景场。本发明包括以下技术优点，消除边缘效应的影响，提高了十倍的算法精度。
Description
技术领域
本发明涉及一种构建水下重力辅助导航背景场的Shepard插值算法，属于顾及移去恢复技术构建水下重力辅助导航背景场的电子信息技术领域。
背景技术
随着国家深海战略的推进，水下导航技术受到了越来越多的关注。水下导航技术中应用最广泛的方法是惯性导航，该方法具有精度高、可靠性好等优点，但是在进行轨迹解算和定位解算时系统误差随着时间积累，如不定期进行重调，定位误差会无限制增长，不适合水下潜器长时间自主航行。重力属于地球固有物理特征，利用海洋重力信息开展水下重力辅助导航，是克服惯性导航系统定位误差随时间积累的缺陷的有效途径。高精度高分辨率的海洋重力背景场是进行水下重力辅助导航的前提和必要条件，海洋重力测量由于经济成本和时间成本等方面的原因无法直接获取达到水下重力辅助导航分辨率要求的重力信息，寻求构建海洋重力背景场的高精度插值算法，是实现水下重力辅助导航的必要步骤。
Shepard插值法是以重力异常在空间位置上的相对距离为基础的有理函数逼近方法，具有整体空间插值和局部空间插值的双重特性，且该方法不需要求解线性方程组，时间效率高，在诸多工程领域受到重视。地球重力信息有其固有的物理特性，以拉普拉斯方程为基础，可以表达成由不同阶次组成的球谐形式，以EGM2008为代表的全球重力场模型提供2160阶次的中长波重力信息，为建立海洋重力背景场提供了物理特性信息。Shepard插值算法仅考虑了待估点与观测点之间的空间位置信息，为提高算法的插值精度，有必要研究顾及海洋重力物理特性构建水下重力辅助导航背景场的Shepard插值算法。
发明内容
本发明的目的在于提供一种顾及移去恢复技术，构建水下重力辅助导航背景场的Shepard插值算法，克服了现有技术的不足。
本发明的原理在于，针对Shepard插值算法仅考虑了待估点与观测点之间的空间位置信息，没有顾及地球重力场固有物理特性的问题，本发明提出一种顾及移去恢复技术构建水下重力辅助导航背景场的Shepard插值新算法，即利用全球重力场模型，首先移去重力异常观测值中的重力异常模型值，得到观测点处的残差重力异常，然后以观测点处的残差重力异常为基础，应用Shepard插值算法，得到待估点处的残差重力异常，最后恢复待估点处残差重力异常的重力异常模型值，得到待估点处的重力异常值，构建满足水下重力辅助导航需求的海洋重力背景场。
本发明通过以下技术方案实现，包括以海洋重力观测值为基础数据，联合全球重力场模型，采用顾及移去恢复技术的Shepard插值算法，构建满足水下重力辅助导航需求的海洋重力背景场。
构建水下重力辅助导航背景场的Shepard插值算法包括以下步骤：
步骤1、移去重力异常观测值中的模型重力异常，得到观测点处的残差重力异常，
δΔgi＝Δgi-ΔgiM
式中，Δgi是观测点i处重力异常实际测量值，ΔgiM是由全球重力场模型计算的观测点i处的模型重力异常，δΔgi是观测点i处的残差重力异常；
模型重力异常ΔgM的计算公式为：
式中，(r,θ,λ)分别是该点处的地心向径、余纬和经度，a为参考椭球长半径，GM为万有引力常数和地球总质量的乘积，和是n阶m次完全规格化位系数，为缔合Legendre函数，Nmax是最高阶数；
步骤2、应用Shepard插值算法，得到待估点处的残差重力异常，
式中，为待估点X0处的残差重力异常估计值，δΔg(Xi)为观测点Xi处的残差重力异常，n为观测点个数，ρ(di)为观测处重力异常的权函数，di为待估点X0和观测点Xi之间的距离；
Shepard插值算法的权函数ρ(di)的计算公式为：
式中，R为以待估点为圆心的插值半径，本发明R取0.5度；
步骤3、恢复待估点残差重力异常的模型重力异常，得到待估点处的重力异常值，
Δg0＝δΔg0+Δg0M
式中，Δg0是待估点处重力异常计算值，Δg0M是由全球重力场模型计算的待估点0处的模型重力异常，δΔg0是待估点处的残差重力异常。
本发明包括以下技术优点，消除边缘效应的影响，提高了十倍的算法精度；极大地提高了算法的精确性，具有突出的实质性特点和显著的进步。
附图说明
图1为重力异常数据的特征统计图。
图2为重力异常标准值的比对图。
图3为传统Shepard插值算法与DTU10重力异常标准值的差异值示意图。
图4为本发明算法与DTU10重力异常标准值的差异值示意图。
具体实施方式
本发明是一种顾及移去恢复技术构建水下重力辅助导航背景场的Shepard插值新方法，包括以海洋重力观测值为基础数据，联合全球重力场模型，采用顾及移去恢复技术的Shepard插值算法，构建满足水下重力辅助导航需求的海洋重力背景场。
构建水下重力辅助导航背景场的Shepard插值算法包括以下步骤：
步骤1、移去重力异常观测值中的模型重力异常，得到观测点处的残差重力异常，
δΔgi＝Δgi-ΔgiM
式中，Δgi是观测点i处重力异常实际测量值，ΔgiM是由全球重力场模型计算的观测点i处的模型重力异常，δΔgi是观测点i处的残差重力异常；
模型重力异常ΔgM的计算公式为：
式中，(r,θ,λ)分别是该点处的地心向径、余纬和经度，a为参考椭球长半径，GM为万有引力常数和地球总质量的乘积，和是n阶m次完全规格化位系数，为缔合Legendre函数，Nmax是最高阶数；
步骤2、应用Shepard插值算法，得到待估点处的残差重力异常，
式中，为待估点X0处的残差重力异常估计值，δΔg(Xi)为观测点Xi处的残差重力异常，n为观测点个数，ρ(di)为观测处重力异常的权函数，di为待估点X0和观测点Xi之间的距离；
Shepard插值算法的权函数ρ(di)的计算公式为：
式中，R为以待估点为圆心的插值半径，本发明R取0.5度；
步骤3、恢复待估点残差重力异常的模型重力异常，得到待估点处的重力异常值，
Δg0＝δΔg0+Δg0M
式中，Δg0是待估点处重力异常计算值，Δg0M是由全球重力场模型计算的待估点0处的模型重力异常，δΔg0是待估点处的残差重力异常。
为使本发明的目的、技术方案及优点更加清楚明白，以下结合附图举例对本发明内容作详细说明。以DTU10重力异常数值模型作为基础数据设计试验，该模型是由丹麦科技大学联合多代卫星测高数据计算的分辨率为1′×1′的全球海域重力场数值模型。采用2160阶次的EGM2008重力场模型作为移去恢复技术的参考场。选择了一个3°×3°区块作为主要试验区进行数值计算，该海域的DTU10重力异常值、EGM2008重力异常模型值及残差重力异常值的特征统计见图1：试验区重力异常数据的特征统计表/mGal。
试验方案为格网化加密应用，即设原有的格网化数据分辨率不满足水下重力辅助导航应用的需求，需要进行格网化加密以得到更高分辨率的海域重力异常背景场。将DTU10重力异常数值模型抽稀成5'×5'格网化数据，然后利用格网化方法将其加密成分辨率为1'×1'的格网化数据。利用未参与格网化计算的DTU10重力异常数据作为参考对本发明方法得到的计算结果进行精度评价，为了消除边缘效应的影响，参与精度评价的数据范围中心2°×2°的区域。为比较分析本发明算法的有效性，引入没有应用移去恢复技术的传统Shepard插值算法进行格网化计算。图2给出了两种插值方法的计算结果与DTU10重力异常标准值的比对结果，结合图3、4可看出，本发明算法构建的水下重力辅助导航背景场精度为0.29mGal，优于传统Shepard插值算法精度的2.739mGal，提高了十倍的算法精度，验证了本发明算法的先进性，基于本发明算法构建的海域背景场能有效满足水下重力辅助导航的需求。
Claims (2)
1.构建水下重力辅助导航背景场的Shepard插值算法，包括以海洋重力观测值为基础数据，联合全球重力场模型，采用顾及移去恢复技术的Shepard插值算法，构建满足水下重力辅助导航需求的海洋重力背景场。
2.根据权利要求1所述的构建水下重力辅助导航背景场的Shepard插值算法，包括以下步骤：
步骤1、移去重力异常观测值中的模型重力异常，得到观测点处的残差重力异常，
δΔgi＝Δgi-ΔgiM
式中，Δgi是观测点i处重力异常实际测量值，ΔgiM是由全球重力场模型计算的观测点i处的模型重力异常，δΔgi是观测点i处的残差重力异常。
模型重力异常ΔgM的计算公式为：
<mrow>
<msub>
<mi>&Delta;g</mi>
<mi>M</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>G</mi>
<mi>M</mi>
</mrow>
<msup>
<mi>r</mi>
<mn>2</mn>
</msup>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>2</mn>
</mrow>
<msub>
<mi>N</mi>
<mrow>
<mi>M</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
</munderover>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mi>a</mi>
<mi>r</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mi>n</mi>
</msup>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>n</mi>
</munderover>
<mrow>
<mo>(</mo>
<msubsup>
<mover>
<mi>C</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>n</mi>
<mi>m</mi>
</mrow>
<mo>*</mo>
</msubsup>
<mi>cos</mi>
<mi> </mi>
<mi>m</mi>
<mi>&lambda;</mi>
<mo>+</mo>
<msub>
<mover>
<mi>S</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>n</mi>
<mi>m</mi>
</mrow>
</msub>
<mi>sin</mi>
<mi> </mi>
<mi>m</mi>
<mi>&lambda;</mi>
<mo>)</mo>
</mrow>
<msub>
<mover>
<mi>P</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>n</mi>
<mi>m</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
</mrow>
式中，(r,θ,λ)分别是该点处的地心向径、余纬和经度，a为参考椭球长半径，GM为万有引力常数和地球总质量的乘积，和是n阶m次完全规格化位系数，为缔合Legendre函数，Nmax是最高阶数；
步骤2、应用Shepard插值算法，得到待估点处的残差重力异常，
<mrow>
<mi>&delta;</mi>
<mi>&Delta;</mi>
<mover>
<mi>g</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<msub>
<mi>X</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<mi>&rho;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>d</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mi>&delta;</mi>
<mi>&Delta;</mi>
<mi>g</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>X</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
式中，为待估点X0处的残差重力异常估计值，δΔg(Xi)为观测点Xi处的残差重力异常，n为观测点个数，ρ(di)为观测处重力异常的权函数，di为待估点X0和观测点Xi之间的距离；
Shepard插值算法的权函数ρ(di)的计算公式为：
<mrow>
<mi>&rho;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>d</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mfrac>
<mn>1</mn>
<msub>
<mi>d</mi>
<mi>i</mi>
</msub>
</mfrac>
</mtd>
<mtd>
<mrow>
<mn>0</mn>
<mo><</mo>
<msub>
<mi>d</mi>
<mi>i</mi>
</msub>
<mo>&le;</mo>
<mfrac>
<mi>R</mi>
<mn>3</mn>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mn>27</mn>
<mrow>
<mn>4</mn>
<mi>R</mi>
</mrow>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>d</mi>
<mi>i</mi>
</msub>
<mi>R</mi>
</mfrac>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mtd>
<mtd>
<mrow>
<mfrac>
<mi>R</mi>
<mn>3</mn>
</mfrac>
<mo><</mo>
<msub>
<mi>d</mi>
<mi>i</mi>
</msub>
<mo>&le;</mo>
<mi>R</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<msub>
<mi>d</mi>
<mi>i</mi>
</msub>
<mo>></mo>
<mi>R</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
式中，R为以待估点为圆心的插值半径，本发明R取0.5度；
步骤3、恢复待估点残差重力异常的模型重力异常，得到待估点处的重力异常值，
Δg0＝δΔg0+Δg0M
式中，Δg0是待估点处重力异常计算值，Δg0M是由全球重力场模型计算的待估点0处的模型重力异常，δΔg0是待估点处的残差重力异常。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710475567.8A CN107289943A (zh) | 2017-06-21 | 2017-06-21 | 构建水下重力辅助导航背景场的Shepard插值算法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710475567.8A CN107289943A (zh) | 2017-06-21 | 2017-06-21 | 构建水下重力辅助导航背景场的Shepard插值算法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107289943A true CN107289943A (zh) | 2017-10-24 |
Family
ID=60097650
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710475567.8A Withdrawn CN107289943A (zh) | 2017-06-21 | 2017-06-21 | 构建水下重力辅助导航背景场的Shepard插值算法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107289943A (zh) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111428347A (zh) * | 2020-03-06 | 2020-07-17 | 中国人民解放军92859部队 | 基于频域特性构建海域垂线偏差模型的Shepard格网算法 |
| CN111460618A (zh) * | 2020-03-06 | 2020-07-28 | 中国人民解放军92859部队 | 基于频域特性构建海域垂线偏差模型的反距离加权格网法 |
| CN111830585A (zh) * | 2020-06-12 | 2020-10-27 | 国家海洋局南海调查技术中心(国家海洋局南海浮标中心) | 基于移动平均格网算法的海域垂线偏差获取方法及其装置 |
| CN111830586A (zh) * | 2020-06-12 | 2020-10-27 | 国家海洋局南海调查技术中心(国家海洋局南海浮标中心) | 基于克里金格网算法的海域垂线偏差获取方法及其装置 |
| CN112949049A (zh) * | 2021-02-08 | 2021-06-11 | 中国人民解放军92859部队 | 一种利用带限思想计算重力异常低阶径向导数的方法 |
| CN112965124A (zh) * | 2021-02-08 | 2021-06-15 | 中国人民解放军92859部队 | 一种顾及局域保障条件计算外部重力异常垂直梯度的方法 |
-
2017
- 2017-06-21 CN CN201710475567.8A patent/CN107289943A/zh not_active Withdrawn
Non-Patent Citations (1)
| Title |
|---|
| 李姗姗: "水下重力辅助惯性导航的理论与方法研究", 《中国博士学位论文全文数据库基础科学辑》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111428347A (zh) * | 2020-03-06 | 2020-07-17 | 中国人民解放军92859部队 | 基于频域特性构建海域垂线偏差模型的Shepard格网算法 |
| CN111460618A (zh) * | 2020-03-06 | 2020-07-28 | 中国人民解放军92859部队 | 基于频域特性构建海域垂线偏差模型的反距离加权格网法 |
| CN111830585A (zh) * | 2020-06-12 | 2020-10-27 | 国家海洋局南海调查技术中心(国家海洋局南海浮标中心) | 基于移动平均格网算法的海域垂线偏差获取方法及其装置 |
| CN111830586A (zh) * | 2020-06-12 | 2020-10-27 | 国家海洋局南海调查技术中心(国家海洋局南海浮标中心) | 基于克里金格网算法的海域垂线偏差获取方法及其装置 |
| CN112949049A (zh) * | 2021-02-08 | 2021-06-11 | 中国人民解放军92859部队 | 一种利用带限思想计算重力异常低阶径向导数的方法 |
| CN112965124A (zh) * | 2021-02-08 | 2021-06-15 | 中国人民解放军92859部队 | 一种顾及局域保障条件计算外部重力异常垂直梯度的方法 |
| CN112949049B (zh) * | 2021-02-08 | 2021-11-30 | 中国人民解放军92859部队 | 一种利用带限思想计算重力异常低阶径向导数的方法 |
| CN112965124B (zh) * | 2021-02-08 | 2022-10-11 | 中国人民解放军92859部队 | 一种顾及局域保障条件计算外部重力异常垂直梯度的方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107289943A (zh) | 构建水下重力辅助导航背景场的Shepard插值算法 | |
| CN107063269A (zh) | 构建水下重力辅助导航背景场的反距离加权插值算法 | |
| CN107167142A (zh) | 构建水下重力辅助导航背景场的径向基函数插值算法 | |
| CN107300387A (zh) | 构建水下重力辅助导航背景场的Kriging插值算法 | |
| CN104035138B (zh) | 一种全球及局部海洋扰动重力的精确快速计算方法 | |
| CN104390646B (zh) | 水下潜器地形辅助惯性导航系统的位置匹配方法 | |
| CN103940446B (zh) | 一种基于横坐标系的极区航行船舶捷联惯性导航系统重调方法 | |
| CN103697854A (zh) | 一种非接触式结构面产状测量方法 | |
| CN104075715A (zh) | 一种结合地形和环境特征的水下导航定位方法 | |
| CN109345875A (zh) | 一种提高船舶自动识别系统测量精度的估计方法 | |
| CN103941290B (zh) | 一种海底电缆运动轨迹模拟方法及系统 | |
| CN107894235A (zh) | 一种超高速飞行器自主导航系统的模型误差补偿方法 | |
| CN103630139A (zh) | 一种基于地磁梯度张量测量的水下载体全姿态确定方法 | |
| CN106197409A (zh) | 一种海底管道三维地理坐标测量方法 | |
| CN104199115B (zh) | 一种基于线性模型的地磁矢量测量误差综合补偿方法 | |
| CN103900564A (zh) | 一种潜深辅助地磁异常反演测速/水下连续定位方法 | |
| CN103968838A (zh) | 一种基于极坐标系的auv曲线运动状态下的协同定位方法 | |
| Huang et al. | Weight self-adjustment Adams implicit filtering algorithm for attitude estimation applied to underwater gliders | |
| CN104040378B (zh) | 气象预测装置以及气象预测方法 | |
| CN103364842B (zh) | 一种捷联式航空重力仪误差分离方法 | |
| CN107270937A (zh) | 一种离线小波降噪快速初始对准方法 | |
| CN107014374B (zh) | 一种基于互补滤波的水下滑翔器节能算法 | |
| CN106705967A (zh) | 一种基于行人航位推算的精度改善的室内定位和方法 | |
| CN106441309A (zh) | 基于协方差交叉融合的火星进入段分布式自主导航方法 | |
| CN103344252B (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 | ||
| WW01 | Invention patent application withdrawn after publication |
Application publication date: 20171024 |
|
| WW01 | Invention patent application withdrawn after publication |