CN107504974A - 地形分块与地形测点加权的地形匹配定位方法 - Google Patents

地形分块与地形测点加权的地形匹配定位方法 Download PDF

Info

Publication number
CN107504974A
CN107504974A CN201710832755.1A CN201710832755A CN107504974A CN 107504974 A CN107504974 A CN 107504974A CN 201710832755 A CN201710832755 A CN 201710832755A CN 107504974 A CN107504974 A CN 107504974A
Authority
CN
China
Prior art keywords
mrow
msub
landform
msup
node
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.)
Granted
Application number
CN201710832755.1A
Other languages
English (en)
Other versions
CN107504974B (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.)
Harbin Engineering University
Original Assignee
Harbin Engineering University
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 Harbin Engineering University filed Critical Harbin Engineering University
Priority to CN201710832755.1A priority Critical patent/CN107504974B/zh
Publication of CN107504974A publication Critical patent/CN107504974A/zh
Application granted granted Critical
Publication of CN107504974B publication Critical patent/CN107504974B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/20Instruments for performing navigational calculations
    • G01C21/203Specially adapted for sailing ships

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Automation & Control Theory (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Image Processing (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)
  • Navigation (AREA)

Abstract

本发明的目的在于提供地形分块与地形测点加权的地形匹配定位方法,影响地形匹配定位精度的有两方面的因素:地形的适配性和地形测量误差。对于适配性较大的局部地形区域可以提供较多的地形信息,而测量误差会使局部地形产生畸变,给地形匹配定位带来负面影响。由于适配性和地形的测量误差都是局部地形的特征,为了体现适配性和局部地形测量误差对匹配定位的影响,将地形进行分块并利用适配性对分块子地形图中的节点进行加权,同时在匹配过程中利用分快子地形的残差统计方差估计子地形图的测量误差,并对子地形图中的节点进行加权处理,利用适配性和测量误差得到的权值同时对子地形图中的节点进行加权,通过迭代过程获得最优的地形匹配定位结果。

Description

地形分块与地形测点加权的地形匹配定位方法
技术领域
本发明涉及的是一种地形匹配定位方法,具体地说是深海地形匹配定位方法。
背景技术
深海地形的测绘主要依靠AUV来完成,由于AUV作业能力和导航精度的限制,在大范围的绘图过程中需要先得到一些局部小区域的地形图最后进行拼接处理得到大范围的水下地形图。此外,在测量过程中需要利用已获得的地形图进行地形匹配定位以此修正导航误差。地形拼接和地形匹配定位过程都需要高精度的地形匹配定位方法,由于地形匹配定位的精度主要有两个方面的影响:地形的适配性和地形测量误差及地形插值重构过程中导致的地形畸变误差。现有的地形匹配定位方法并没有考虑地形的适配性,而且对于地形的测量误差也只是考虑成高斯噪声,这些简化处理使得定位结果的和似然函数的收敛性并不好,同时定位的精度也容易受到地形适配性和地形畸变误差的影响。
本发明主要是出于对匹配过程中地形适配性和地形畸变误差对地形匹配定位精度的考虑,通过对先验地形的分块对地形的测量点进行加权处理。
发明内容
本发明的目的在于提供考虑对匹配过程中地形适配性和地形畸变误差对地形匹配定位精度的地形分块与地形测点加权的地形匹配定位方法。
本发明的目的是这样实现的:
本发明地形分块与地形测点加权的地形匹配定位方法,其特征是:
(1)进行初匹配,并将测量地形的测点分块:
假设分块子地形图的边界点为k,设定先验地形图地形图的分块大小,将先验地形图分割成M×N块子地形图,每一个子地形图的边界上的地形节点数为k,计算每一个地形点的适配性量化参数,量化参数采用地形节点的8方向信噪比:
式中:i,j分别表示鲜艳地形图节点的行和列索引号,d表示先验地形的网格边长,k表示表示8个方向的索引号,表示地形节点i,j在k方向的梯度,σ表示地形测量误差;
计算每一个分块的适配性参数:
式中:p表示分块自地图的边界地形节点个数,I,J表示地形分块的行和列索引号;
利用最大化原则将8个方向上的信息量化成一个量
同时根据下式对测量地形进行初步的定位,得到初步定位偏差(dx1,dy1),然后利用定位偏差对测量地形进行位置修正,得到初步修正后的测量地形:
式中:Xp表示地形匹配定位的位置,i,j表示搜索区域内搜索点的索引号, zk表示MTM地形图中的点,表示在搜索点(i,j)处测量地形序列zk在DEM中的差值点高度;
(2)获取测量地形的适配性权值和测量误差权值:
根据修正后的测量地形和先验地形,得到先验地形与测量地形的重叠区域的节点序列,Zn和Zd分别表示重叠区域中位于先验地形和测量地形图中的节点;
根据步骤(1)中得到的先验地形的分块信息得到Zn中的地形节点所在的地形分块索引,根据节点所在的地形分块的适配性参数对先验地形和测量地形重叠区域的每一个节点序列Zn进行加权,假设节点位于分块 (I,J)中,则节点的权值为即位于(I,J)分块内的地形节点取值均为获取所有的Zn序列点的权值λi后对权值进行归一化处理,该权值为地形适配性权值,
根据步骤(1)中的匹配定位计算位置修正后的先验地形和测量地形的高度偏差序列残差:
Δh=[Z-h(Xp)]
残差的均值和方差估计:
根据计算得到的每一个地形分块中的残差方差对Zn地形节点进行加权,与权值λi的确定方法一致,同一个地形块中的地行节点的测量误差权值大小相同,用1/σi表示Zn中的节点地形测量误差权值,对地形测量误差权值进行归一化处理:
(3)根据步骤(2)中得到的适配性权重和测量误差权值计算归一化的匹配权值:
(4)根据步骤(3)中得到的先验地形和测量地形交叠区域内节点的权值重新计算定位点;
其中:qi表示最终得到的节点权值;
(5)判断是否达到迭代终点,如果达到则返回定位结果Xp,若没有达到迭代终点则回到步骤(2)。
本发明的优势在于:由于适配性和地形的测量误差都是局部地形的特征,为了体现适配性和局部地形测量误差对匹配定位的影响,将地形进行分块并利用适配性对分块子地形图中的节点进行加权,同时在匹配过程中利用分快子地形的残差统计方差估计子地形图的测量误差,并对子地形图中的节点进行加权处理,这样利用适配性和测量误差得到的权值同时对子地形图中的节点进行加权,通过迭代过程获得最优的地形匹配定位结果。
附图说明
图1为本发明的流程图;
图2为地形节点适配性计算时的8个离散方向;
图3为地形分块结果;
图4为地形重叠结果;
图5为适配性权值的计算方法;
图6为地形测量误差权值的计算方法;
图7为地形分块与加权匹配定位拼接方法流程图。
具体实施方式
下面结合附图举例对本发明做更详细地描述:
结合图1-7,地形分块与地形测点加权的地形匹配定位方法主要步骤包括,将先验地形图001进行分块,计算分块后的每一个子地图的适配性参数002,对测量地形004进行初步匹配定位计算003根据匹配结果对测量地形进行修正,然后根据分块结果002和修正后的地形进行适配性权值008和测量误差权值的计算006的计算,然后对权值分别进行归一化处理,再利用权值进行加权匹配定位010,判断是否达到了迭代次数011如果没有则一直迭代直到迭代完成,虚线框内表示迭代过程。最后输出匹配定位修正后的测量地形012。
1、初匹配和测量地形的测点分块
设定先验地形图地形图201的分块大小(假设分块子地形图的边界点为k)。将先验地形图(DEM)201分割成M×N块子地形图,每一个子地形图的边界上的地形节点数为k,如图1所示表示先验地形的分块情况。计算每一个地形点 301、202的适配性量化参数,量化参数采用地形节点的8方向302信噪比:
式中:
i,j分别表示鲜艳地形图节点的行和列索引号;
d表示先验地形的网格边长;
k表示表示8个方向的索引号;
表示地形节点i,j在k方向的梯度;
σ表示地形测量误差;
计算每一个分块的适配性参数,计算公式如下:
式中:
p表示分块自地图的边界地形节点个数;
I,J表示地形分块的行和列索引号;
利用最大化原则将8个方向上的信息量化成一个量
同时根据下面的公式对测量地形进行初步的定位,得到初步定位偏差 (dx1,dy1),然后利用定位偏差对测量地形进行位置修正,得到初步修正后的测量地形。
式中:
Xp表示地形匹配定位的位置;
i,j表示搜索区域内搜索点的索引号;
zk表示MTM地形图中的点;
表示在搜索点(i,j)处测量地形序列zk在DEM中的差值点高度;
2、测量地形的适配性权值和测量误差权值计算
根据修正后的测量地形(MTM)301和先验地形(DEM),得到DEM与MTM的重叠区域的节点序列302,Zn和Zd分别表示重叠区域中位于DEM地形和MTM地形图中的节点。
根据1中得到的DEM的分块303信息得到Zn中的地形节点所在的地形分块索引,根据节点所在的地形分块的适配性参数对DEM201和MTM301重叠区域302的每一个节点序列Zn进行加权,假设节点位于分块(I,J)中,则节点的权值为也就是说位于(I,J)分块内的地形节点取值均为获取所有的Zn序列点的权值λi后对权值进行归一化处理,该权值称为地形适配性权值。
然后,根据1中的匹配定位计算位置修正后的DEM201和MTM301的高度偏差序列残差:
Δh=[Z-h(Xp)]
残差的均值和方差估计:
根据计算得到的每一个地形分块中的残差方差对Zn地形节点进行加权,与权值λi的确定方法一致,同一个地形块中的地行节点的测量误差权值大小相同,用1/σi表示Zn中的节点地形测量误差权值,对地形测量误差权值进行归一化处理:
3、根据2中得到的适配性权重和测量误差权值计算归一化的匹配权值:
4、根据3中得到的DEM201和MTM301交叠区域302内节点的权值重新计算定位点;
其中:qi表示最终得到的节点权值;
5、判断是否达到迭代终点,如果达到则返回定位结果Xp,若没有达到迭代终点则回到步骤2。

Claims (1)

1.地形分块与地形测点加权的地形匹配定位方法,其特征是:
(1)进行初匹配,并将测量地形的测点分块:
假设分块子地形图的边界点为k,设定先验地形图地形图的分块大小,将先验地形图分割成M×N块子地形图,每一个子地形图的边界上的地形节点数为k,计算每一个地形点的适配性量化参数,量化参数采用地形节点的8方向信噪比:
<mrow> <msubsup> <mi>SSNR</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mi>k</mi> </msubsup> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <msubsup> <mi>&amp;Delta;z</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mi>k</mi> </msubsup> <mrow> <mo>(</mo> <mo>&amp;CenterDot;</mo> <mo>)</mo> </mrow> </mrow> <mrow> <mi>d</mi> <mi>&amp;sigma;</mi> </mrow> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>,</mo> <mi>k</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>3</mn> <mo>,</mo> <mn>5</mn> <mo>,</mo> <mn>7</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <msubsup> <mi>&amp;Delta;z</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mi>k</mi> </msubsup> <mrow> <mo>(</mo> <mo>&amp;CenterDot;</mo> <mo>)</mo> </mrow> </mrow> <mrow> <msqrt> <mn>2</mn> </msqrt> <mi>d</mi> <mi>&amp;sigma;</mi> </mrow> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>,</mo> <mi>k</mi> <mo>=</mo> <mn>2</mn> <mo>,</mo> <mn>4</mn> <mo>,</mo> <mn>6</mn> <mo>,</mo> <mn>8</mn> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>
式中:i,j分别表示鲜艳地形图节点的行和列索引号,d表示先验地形的网格边长,k表示表示8个方向的索引号,表示地形节点i,j在k方向的梯度,σ表示地形测量误差;
计算每一个分块的适配性参数:
<mrow> <msubsup> <mover> <mrow> <mi>S</mi> <mi>N</mi> <mi>R</mi> </mrow> <mo>&amp;OverBar;</mo> </mover> <mrow> <mi>I</mi> <mi>J</mi> </mrow> <mi>k</mi> </msubsup> <mo>=</mo> <mfrac> <mn>1</mn> <msup> <mi>p</mi> <mn>2</mn> </msup> </mfrac> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>p</mi> </munderover> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>p</mi> </munderover> <mrow> <mo>(</mo> <msubsup> <mi>SSNR</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mi>k</mi> </msubsup> <mo>)</mo> </mrow> </mrow>
式中:p表示分块自地图的边界地形节点个数,I,J表示地形分块的行和列索引号;
利用最大化原则将8个方向上的信息量化成一个量
<mrow> <mover> <mrow> <msub> <mi>SSNR</mi> <mrow> <mi>I</mi> <mi>J</mi> </mrow> </msub> </mrow> <mo>&amp;OverBar;</mo> </mover> <mo>=</mo> <mi>m</mi> <mi>a</mi> <mi>x</mi> <mrow> <mo>(</mo> <msub> <mi>SSNR</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow>
同时根据下式对测量地形进行初步的定位,得到初步定位偏差(dx1,dy1),然后利用定位偏差对测量地形进行位置修正,得到初步修正后的测量地形:
<mrow> <msub> <mi>L</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mn>2</mn> <msup> <mi>&amp;pi;</mi> <mrow> <mi>N</mi> <mo>/</mo> <mn>2</mn> </mrow> </msup> <msup> <mi>&amp;sigma;</mi> <mi>N</mi> </msup> </mrow> </mfrac> <mo>&amp;CenterDot;</mo> <mi>exp</mi> <mo>&amp;lsqb;</mo> <mo>-</mo> <mfrac> <mn>1</mn> <mrow> <mn>2</mn> <msup> <mi>&amp;sigma;</mi> <mn>2</mn> </msup> </mrow> </mfrac> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>K</mi> </munderover> <msup> <mrow> <mo>(</mo> <msub> <mi>z</mi> <mi>k</mi> </msub> <mo>-</mo> <msubsup> <mi>h</mi> <mi>k</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msubsup> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>&amp;rsqb;</mo> </mrow>
<mrow> <msup> <mi>X</mi> <mi>p</mi> </msup> <mo>=</mo> <munder> <mi>argmax</mi> <mrow> <msup> <mi>X</mi> <mi>p</mi> </msup> <mo>&amp;Element;</mo> <msub> <mi>X</mi> <mi>s</mi> </msub> </mrow> </munder> <mrow> <mo>(</mo> <msub> <mi>L</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow>
式中:Xp表示地形匹配定位的位置,i,j表示搜索区域内搜索点的索引号,zk表示MTM地形图中的点,表示在搜索点(i,j)处测量地形序列zk在DEM中的差值点高度;
(2)获取测量地形的适配性权值和测量误差权值:
根据修正后的测量地形和先验地形,得到先验地形与测量地形的重叠区域的节点序列,Zn和Zd分别表示重叠区域中位于先验地形和测量地形图中的节点;
根据步骤(1)中得到的先验地形的分块信息得到Zn中的地形节点所在的地形分块索引,根据节点所在的地形分块的适配性参数对先验地形和测量地形重叠区域的每一个节点序列Zn进行加权,假设节点位于分块(I,J)中,则节点的权值为即位于(I,J)分块内的地形节点取值均为获取所有的Zn序列点的权值λi后对权值进行归一化处理,该权值为地形适配性权值,
根据步骤(1)中的匹配定位计算位置修正后的先验地形和测量地形的高度偏差序列残差:
Δh=[Z-h(Xp)]
残差的均值和方差估计:
<mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mover> <mrow> <mi>&amp;Delta;</mi> <mi>h</mi> </mrow> <mo>&amp;OverBar;</mo> </mover> <mo>=</mo> <mi>m</mi> <mi>e</mi> <mi>a</mi> <mi>n</mi> <mo>(</mo> <mi>&amp;Delta;</mi> <mi>h</mi> <mo>)</mo> </mtd> </mtr> <mtr> <mtd> <mi>&amp;sigma;</mi> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mi>M</mi> <mi>N</mi> </mrow> </mfrac> <mi>s</mi> <mi>u</mi> <mi>m</mi> <mo>(</mo> <mi>&amp;Delta;</mi> <mi>h</mi> <mo>-</mo> <mover> <mrow> <mi>&amp;Delta;</mi> <mi>h</mi> </mrow> <mo>&amp;OverBar;</mo> </mover> <mo>)</mo> </mtd> </mtr> </mtable> </mfenced>
根据计算得到的每一个地形分块中的残差方差对Zn地形节点进行加权,与权值λi的确定方法一致,同一个地形块中的地行节点的测量误差权值大小相同,用1/σi表示Zn中的节点地形测量误差权值,对地形测量误差权值进行归一化处理:
<mrow> <mover> <msub> <mi>&amp;sigma;</mi> <mi>i</mi> </msub> <mo>&amp;OverBar;</mo> </mover> <mo>=</mo> <mfrac> <mrow> <mn>1</mn> <mo>/</mo> <msub> <mi>&amp;sigma;</mi> <mi>i</mi> </msub> </mrow> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>k</mi> </munderover> <mrow> <mo>(</mo> <mn>1</mn> <mo>/</mo> <msub> <mi>&amp;sigma;</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>;</mo> </mrow>
(3)根据步骤(2)中得到的适配性权重和测量误差权值计算归一化的匹配权值:
<mrow> <msub> <mi>q</mi> <mi>i</mi> </msub> <mo>=</mo> <mfrac> <mrow> <mover> <msub> <mi>&amp;lambda;</mi> <mi>i</mi> </msub> <mo>&amp;OverBar;</mo> </mover> <mo>&amp;CenterDot;</mo> <mover> <msub> <mi>&amp;sigma;</mi> <mi>i</mi> </msub> <mo>&amp;OverBar;</mo> </mover> </mrow> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>k</mi> </munderover> <mrow> <mo>(</mo> <mover> <msub> <mi>&amp;lambda;</mi> <mi>i</mi> </msub> <mo>&amp;OverBar;</mo> </mover> <mo>&amp;CenterDot;</mo> <mover> <msub> <mi>&amp;sigma;</mi> <mi>i</mi> </msub> <mo>&amp;OverBar;</mo> </mover> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>;</mo> </mrow>
(4)根据步骤(3)中得到的先验地形和测量地形交叠区域内节点的权值重新计算定位点;
<mrow> <msub> <mi>L</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mn>2</mn> <msup> <mi>&amp;pi;</mi> <mrow> <mi>N</mi> <mo>/</mo> <mn>2</mn> </mrow> </msup> <msup> <mi>&amp;sigma;</mi> <mi>N</mi> </msup> </mrow> </mfrac> <mo>&amp;CenterDot;</mo> <mi>exp</mi> <mo>&amp;lsqb;</mo> <mo>-</mo> <mfrac> <msub> <mi>q</mi> <mi>i</mi> </msub> <mrow> <mn>2</mn> <msup> <mi>&amp;sigma;</mi> <mn>2</mn> </msup> </mrow> </mfrac> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>k</mi> </munderover> <msup> <mrow> <mo>(</mo> <msub> <mi>z</mi> <mi>i</mi> </msub> <mo>-</mo> <msubsup> <mi>h</mi> <mi>i</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msubsup> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>&amp;rsqb;</mo> </mrow>
<mrow> <msup> <mi>X</mi> <mi>p</mi> </msup> <mo>=</mo> <munder> <mi>argmax</mi> <mrow> <msup> <mi>X</mi> <mi>p</mi> </msup> <mo>&amp;Element;</mo> <msub> <mi>X</mi> <mi>s</mi> </msub> </mrow> </munder> <mrow> <mo>(</mo> <msub> <mi>L</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow>
其中:qi表示最终得到的节点权值;
(5)判断是否达到迭代终点,如果达到则返回定位结果Xp,若没有达到迭代终点则回到步骤(2)。
CN201710832755.1A 2017-09-15 2017-09-15 地形分块与地形测点加权的地形匹配定位方法 Active CN107504974B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710832755.1A CN107504974B (zh) 2017-09-15 2017-09-15 地形分块与地形测点加权的地形匹配定位方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710832755.1A CN107504974B (zh) 2017-09-15 2017-09-15 地形分块与地形测点加权的地形匹配定位方法

Publications (2)

Publication Number Publication Date
CN107504974A true CN107504974A (zh) 2017-12-22
CN107504974B CN107504974B (zh) 2020-09-25

Family

ID=60696981

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710832755.1A Active CN107504974B (zh) 2017-09-15 2017-09-15 地形分块与地形测点加权的地形匹配定位方法

Country Status (1)

Country Link
CN (1) CN107504974B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110207721A (zh) * 2019-06-06 2019-09-06 哈尔滨工程大学 一种考虑残差分布的无效地形匹配结果识别方法
CN110441760A (zh) * 2019-09-16 2019-11-12 河海大学 一种基于先验地形图的大范围海底地形图拓展构图方法
CN111784717A (zh) * 2020-06-28 2020-10-16 中国计量大学 一种智慧台灯及基于该台灯的自动阅卷方法
CN113297344A (zh) * 2021-06-16 2021-08-24 重庆富沛和科技有限公司 基于三维遥感图像的地性线匹配方法、装置及地物目标位置定位方法
WO2022107126A1 (en) * 2020-11-23 2022-05-27 Israel Aerospace Industries Ltd. Location identification based on terrain model with distance measurement

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105352496A (zh) * 2015-11-17 2016-02-24 中国海洋大学 基于声呐辅助自主导航的auv导航方法
CN106875486A (zh) * 2017-02-22 2017-06-20 哈尔滨工程大学 一种基于节点信息量统计的多波束地形分块方法
CN106885576A (zh) * 2017-02-22 2017-06-23 哈尔滨工程大学 一种基于多点地形匹配定位的auv航迹偏差估计方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105352496A (zh) * 2015-11-17 2016-02-24 中国海洋大学 基于声呐辅助自主导航的auv导航方法
CN106875486A (zh) * 2017-02-22 2017-06-20 哈尔滨工程大学 一种基于节点信息量统计的多波束地形分块方法
CN106885576A (zh) * 2017-02-22 2017-06-23 哈尔滨工程大学 一种基于多点地形匹配定位的auv航迹偏差估计方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
YE LI ET AL.: "Terrain Matching Positioning Method Based on Node Multi-information Fusion", 《THE JOURNAL OF NAVIGATION》 *
陈小龙: "AUV水下地形匹配辅助导航技术研究", 《中国博士学位论文全文数据库 工程科技Ⅱ辑》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110207721A (zh) * 2019-06-06 2019-09-06 哈尔滨工程大学 一种考虑残差分布的无效地形匹配结果识别方法
CN110207721B (zh) * 2019-06-06 2022-06-21 哈尔滨工程大学 一种考虑残差分布的无效地形匹配结果识别方法
CN110441760A (zh) * 2019-09-16 2019-11-12 河海大学 一种基于先验地形图的大范围海底地形图拓展构图方法
CN110441760B (zh) * 2019-09-16 2023-04-04 河海大学 一种基于先验地形图的大范围海底地形图拓展构图方法
CN111784717A (zh) * 2020-06-28 2020-10-16 中国计量大学 一种智慧台灯及基于该台灯的自动阅卷方法
CN111784717B (zh) * 2020-06-28 2024-01-16 中国计量大学 一种智慧台灯及基于该台灯的自动阅卷方法
WO2022107126A1 (en) * 2020-11-23 2022-05-27 Israel Aerospace Industries Ltd. Location identification based on terrain model with distance measurement
CN113297344A (zh) * 2021-06-16 2021-08-24 重庆富沛和科技有限公司 基于三维遥感图像的地性线匹配方法、装置及地物目标位置定位方法
CN113297344B (zh) * 2021-06-16 2021-11-30 重庆富沛和科技有限公司 基于三维遥感图像的地性线匹配方法、装置及地物目标位置定位方法

Also Published As

Publication number Publication date
CN107504974B (zh) 2020-09-25

Similar Documents

Publication Publication Date Title
CN107504974A (zh) 地形分块与地形测点加权的地形匹配定位方法
CN106885576B (zh) 一种基于多点地形匹配定位的auv航迹偏差估计方法
CN104061932B (zh) 一种利用引力矢量和梯度张量进行导航定位的方法
CN106052688B (zh) 基于地形轮廓匹配的惯性导航系统速度累积误差修正方法
CN106871880B (zh) 采用三维坐标计算肋位、半宽和高度进行船舶合拢的方法
CN104075715B (zh) 一种结合地形和环境特征的水下导航定位方法
CN108896040B (zh) 天空海一体化水下潜器惯性/重力组合导航方法和系统
CN103217688B (zh) 一种基于不规则三角网机载激光雷达点云平差计算方法
CN111982532B (zh) 一种自动泊车性能测试方法
CN105091907B (zh) Sins/dvl组合中dvl方位安装误差估计方法
CN103369466B (zh) 一种地图匹配辅助室内定位方法
CN109507706B (zh) 一种gps信号丢失的预测定位方法
CN109708629B (zh) 一种用于差异定位性能条件的飞行器集群协同导航方法
CN103379619A (zh) 一种定位方法和系统
CN111694040B (zh) 一种卫星/惯性组合导航系统定位方法与装置
CN104075711B (zh) 一种基于CKF的IMU/Wi‑Fi信号紧组合室内导航方法
CN114061591B (zh) 一种基于滑动窗数据回溯的等值线匹配方法
CN107966145B (zh) 一种基于稀疏长基线紧组合的auv水下导航方法
CN104697523A (zh) 基于迭代计算的惯性/地磁匹配定位方法
CN109059964A (zh) 一种基于重力峰的惯性导航与重力测量双校准方法
CN107643082A (zh) 基于多波束的多路径并行iccp水下地形匹配方法
CN103033822B (zh) 移动信息确定装置、方法以及接收机
CN109855623A (zh) 基于Legendre多项式和BP神经网络的地磁模型在线逼近方法
CN102661742B (zh) 基于曲率特征加权质心点约束的自适应标志点布局方法
CN105865444A (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
GR01 Patent grant