CN107504974A - 地形分块与地形测点加权的地形匹配定位方法 - Google Patents
地形分块与地形测点加权的地形匹配定位方法 Download PDFInfo
- 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
- terrain
- msub
- msup
- mfrac
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 23
- 238000005259 measurement Methods 0.000 claims abstract description 41
- 238000012876 topography Methods 0.000 claims abstract description 19
- 238000013139 quantization Methods 0.000 claims description 5
- 238000012804 iterative process Methods 0.000 abstract description 3
- 238000004364 calculation method Methods 0.000 description 6
- 230000006978 adaptation Effects 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 230000011218 segmentation Effects 0.000 description 2
- 238000011002 quantification Methods 0.000 description 1
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
- G01C21/203—Specially 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)
- Navigation (AREA)
- Position Fixing By Use Of Radio Waves (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>&Delta;z</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mi>k</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mo>&CenterDot;</mo>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>d</mi>
<mi>&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>&Delta;z</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mi>k</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mo>&CenterDot;</mo>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msqrt>
<mn>2</mn>
</msqrt>
<mi>d</mi>
<mi>&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>&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>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>p</mi>
</munderover>
<munderover>
<mo>&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>&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>&pi;</mi>
<mrow>
<mi>N</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
<msup>
<mi>&sigma;</mi>
<mi>N</mi>
</msup>
</mrow>
</mfrac>
<mo>&CenterDot;</mo>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<msup>
<mi>&sigma;</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<munderover>
<mo>&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>&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>&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>&Delta;</mi>
<mi>h</mi>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>=</mo>
<mi>m</mi>
<mi>e</mi>
<mi>a</mi>
<mi>n</mi>
<mo>(</mo>
<mi>&Delta;</mi>
<mi>h</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>&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>&Delta;</mi>
<mi>h</mi>
<mo>-</mo>
<mover>
<mrow>
<mi>&Delta;</mi>
<mi>h</mi>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
根据计算得到的每一个地形分块中的残差方差对Zn地形节点进行加权,与权值λi的确定方法一致,同一个地形块中的地行节点的测量误差权值大小相同,用1/σi表示Zn中的节点地形测量误差权值,对地形测量误差权值进行归一化处理:
<mrow>
<mover>
<msub>
<mi>&sigma;</mi>
<mi>i</mi>
</msub>
<mo>&OverBar;</mo>
</mover>
<mo>=</mo>
<mfrac>
<mrow>
<mn>1</mn>
<mo>/</mo>
<msub>
<mi>&sigma;</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<munderover>
<mo>&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>&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>&lambda;</mi>
<mi>i</mi>
</msub>
<mo>&OverBar;</mo>
</mover>
<mo>&CenterDot;</mo>
<mover>
<msub>
<mi>&sigma;</mi>
<mi>i</mi>
</msub>
<mo>&OverBar;</mo>
</mover>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<mrow>
<mo>(</mo>
<mover>
<msub>
<mi>&lambda;</mi>
<mi>i</mi>
</msub>
<mo>&OverBar;</mo>
</mover>
<mo>&CenterDot;</mo>
<mover>
<msub>
<mi>&sigma;</mi>
<mi>i</mi>
</msub>
<mo>&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>&pi;</mi>
<mrow>
<mi>N</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
<msup>
<mi>&sigma;</mi>
<mi>N</mi>
</msup>
</mrow>
</mfrac>
<mo>&CenterDot;</mo>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mo>-</mo>
<mfrac>
<msub>
<mi>q</mi>
<mi>i</mi>
</msub>
<mrow>
<mn>2</mn>
<msup>
<mi>&sigma;</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<munderover>
<mo>&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>&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>&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)。
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)
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)
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航迹偏差估计方法 |
-
2017
- 2017-09-15 CN CN201710832755.1A patent/CN107504974B/zh active Active
Patent Citations (3)
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)
Title |
---|
YE LI ET AL.: "Terrain Matching Positioning Method Based on Node Multi-information Fusion", 《THE JOURNAL OF NAVIGATION》 * |
陈小龙: "AUV水下地形匹配辅助导航技术研究", 《中国博士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (9)
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) | 地形分块与地形测点加权的地形匹配定位方法 | |
CN106123850B (zh) | Auv配载多波束声呐水下地形测绘修正方法 | |
CN108802770B (zh) | 一种ins增强gnss的高精度动态定位检定基准 | |
CN106885576B (zh) | 一种基于多点地形匹配定位的auv航迹偏差估计方法 | |
CN108007401A (zh) | 一种基于船载InSAR平台的河湖库沿岸形变检测装置及方法 | |
CN106871880A (zh) | 采用三维坐标计算肋位、半宽和高度进行船舶合拢的方法 | |
CN107643082A (zh) | 基于多波束的多路径并行iccp水下地形匹配方法 | |
CN103822634B (zh) | 一种基于改进的iccp算法的重力匹配辅助惯性导航方法 | |
CN101713649A (zh) | 一种基于扰动重力的似大地水准面陆海无缝拼接方法 | |
CN103363953B (zh) | 一种提高卫星测高交叉点精度的方法 | |
CN105157703A (zh) | 一种重力辅助惯性导航适配区可导航性的评价方法 | |
CN103927442A (zh) | 一种基于测角变换的超短基线安装角度误差抗粗差校准方法 | |
CN109193665A (zh) | 一种基于scada量测的电网支路静态参数辨识方法 | |
CN111368406A (zh) | 连续深度基准面构建方法 | |
CN110782452A (zh) | 一种t2定量图像成像方法及系统 | |
CN105160156B (zh) | 一种融合数据预处理的垂测电离图反演方法 | |
CN111721272A (zh) | 一种基于椭球面计算的工程表面测量方法 | |
CN106767834A (zh) | 一种基于模糊熵值的auv水下地形匹配适配区划分方法 | |
CN104462788A (zh) | 一种基于逆向子结构的有限元模型修正方法 | |
CN107909606A (zh) | 一种sar图像配准联系点粗差剔除方法 | |
CN104407366A (zh) | 一种对伪距进行平滑处理的方法 | |
CN107193043A (zh) | 一种起伏地表的地下构造成像方法 | |
CN117169979A (zh) | 一种基于机器学习融合海底地形数据的重力异常反演方法 | |
CN103940781B (zh) | 基于走航数据的smos卫星反演表层海水盐度的校正方法 | |
KR101941132B1 (ko) | 국지적 전리층지도의 가용영역 확장 장치 및 방법 |
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 |