CN102706348A - 一种基于三角形的重力图快速匹配方法 - Google Patents
一种基于三角形的重力图快速匹配方法 Download PDFInfo
- Publication number
- CN102706348A CN102706348A CN2012101534863A CN201210153486A CN102706348A CN 102706348 A CN102706348 A CN 102706348A CN 2012101534863 A CN2012101534863 A CN 2012101534863A CN 201210153486 A CN201210153486 A CN 201210153486A CN 102706348 A CN102706348 A CN 102706348A
- Authority
- CN
- China
- Prior art keywords
- triangle
- coupling
- point
- search
- standard
- 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
Images
Landscapes
- Navigation (AREA)
Abstract
重力辅助惯性导航系统中一种基于三角形的重力图快速匹配方法,主要有粗匹配和精匹配两个过程。根据惯导系统短时精度高的特点,首先利用载体三个连续位置点间的相对距离信息建立标准三角形,然后结合载体上重力仪实时测量的重力值,通过粗匹配和精匹配两个过程,在已知重力图中建立匹配三角形;最后将匹配三角形与标准三角形进行匹配解算获得惯导系统长时间累积的位置误差和航向角误差。本发明是一种快速的重力图匹配方法,具有实时性好、精度高、鲁棒性好的特点,可应用于重力辅助惯性导航等无源导航系统。
Description
技术领域
本发明涉及了重力辅助惯性导航系统中一种基于三角形的重力图快速匹配方法,主要针对重力辅助惯性导航中重力图的快速匹配问题。
背景技术
就地球物理场辅助惯性导航的无源导航技术来说,并不是一个新鲜的课题,最早开始的是陆地上的地形辅助惯性导航研究,并于20世纪70年代末80年代初取得明显的成果。近年来,借助大地测量技术成果,惯性/重力、惯性/地磁无源组合导航的理论和方法同样引起了国内外研究机构和学者的注意,相继开展了重力匹配导航、地磁匹配导航等技术的研究。目前重力辅助惯性导航的技术常采用的方法大多是沿用地形匹配的相关技术,主要有SITAN算法、TERCOM算法等。无论是SITAN算法,还是TERCOM算法,其基本流程都是通过运载体上承载的传感器实时测量运载体运动轨迹上的地球物理场数据;同时,根据惯导系统提供的运载体位置信息,从导航用参考重力图中获取地球物理场数据;然后将这两组数据集送给匹配解算计算机,利用匹配算法确定运载体的最佳匹配位置。
数据集匹配问题最早可以追溯到二十世纪五十年代末的遥感图像分析,那时Hobrough首次对遥感模拟图像进行了相似性分析,提出了图像匹配的概念,至六七十年代,数据集匹配的研究方法大多基于对图像的特征提取,如,Anuta的快速Fourier变换(FFT)方法,但是,由于提取特征方法的限制,使得在特征提取时数据集的信息不能完全得到反映。直到八十年代初,Lucas针对视频图像匹配给出了一种迭代算法,由于该方法直接针对图像数据集进行搜索,因此Lucas的方法完全保留了图像的性质,从而确保了匹配的精确性。到了九十年代初,Chen和Medioni,Besl和McKay以及Zhang等人分别在Lucas的工作基础上将各种具体匹配问题上升到数据集之间的匹配问题,并独立提出了求解数据集匹配问题的ICP(Iterated Closest Point,简称:ICP)算法,尤其是Besl和McKey的工作已经成为数据集匹配问题的一个基础,ICP方法也已经成为求解数据集匹配问题的一般性方法而被广泛采用,而且已被广泛应用于地球物理场辅助惯性导航的重力图(地形图/地磁图)匹配。
ICP算法是一个关于全局收敛性的研究,在理论上可以在任何初始位置误差情况都能工作,而在实际应用中由于目标函数、对应关系选择、变换类型以及初始位置的选择等因素的影响,几乎所有基于ICP的现有方法都只能得到局部最优匹配,是否能取得全局匹配最优的关键在于初始匹配参数的选择。针对初始匹配的选择问题,Hugli等人提出了基于初始参数划分的SIC-range方法,Chetverikov等人通过对每一步迭代中对应点的选择范围进行“裁剪”,并随着迭代的进行逐步缩小和修正对应点的选择范围,以上方法主要源于对静态遥感图像的相似性匹配,它不能满足重力辅助惯性导航系统中对实时性的要求。
发明内容
本发明解决的技术问题:克服现有技术的不足,提供一种重力辅助惯性导航系统中重力图匹配的快速匹配方法,可有效地提高重力辅助惯性导航系统的可靠性和实时性,从而可消除和克服在惯导系统输出的位置信息存在大误差条件下重力图匹配算法失效的缺陷。
本发明采用的技术方案步骤如下:
第一步,根据导航信息得出标准三角形。
从惯导系统的导航解算结果中提取经度、纬度作为位置点的信息,进而连续获取载体的三个位置点信息,以其两两间的航行距离作为标准三角形的三条边。
第二步,结合实时重力量测值,在重力图中进行大范围搜索,获取粗匹配点集。
在已知重力图中建立块阵列进行大范围搜索,块阵列和其中的每个块都是正方形,块阵列按每边上块的个数定义其级数,块阵列的最大级数根据惯导系统的估计位置误差确定。
以第一步中的三个位置点分别为中心在重力图中建立三个块阵列,然后结合第一步三个位置点对应的实时重力值,分别在三个块阵列中以块为最小单元搜索等值点,其搜索过程按照块阵列级数从小到大的顺序进行逐级搜索,由这些搜索的等值点构成粗匹配的三组点集。
第三步,通过标准三角形的约束,建立匹配三角形集,并对其进行筛选。
(1)以第一步标准三角形的三边距离为约束条件,将第二步粗匹配的三组点集按照排列组合的原则建立一组匹配三角形集,其约束条件中三边距离的误差限由第二步块阵列的单个块大小来确定。
(2)计算第三步(1)中各个匹配三角形和标准三角形对应三边的误差之和,并将所有匹配三角形按误差之和从小到大排序;
(3)筛选误差之和最小的匹配三角形用于后续的精匹配过程。
第四步,基于粗匹配筛选结果,在重力图中进行小范围搜索,获取精匹配点集。
以第三步(3)中筛选的匹配三角形的三个点分别为搜索中心,在重力图中建立单个块进行小范围搜索,其块的大小与第二步的块大小相等,然后仍结合第一步三个位置点对应的实时重力值,分别在三个块中搜索等值点,由这些等值点直接构成精匹配的三组点集。
第五步,通过标准三角形的约束,建立新的匹配三角形集,并进行筛选。
(1)以第一步标准三角形的三边距离为约束条件,利用第四步精匹配的三组点集建立一组新的匹配三角形集,其约束条件中三边距离的误差限由重力测量误差和重力图网格精度来确定;
(2)计算第五步(1)中各个新的匹配三角形和标准三角形对应三边的误差之和,并将所有新的匹配三角形按误差之和从小到大排序;
(3)筛选误差之和最小的匹配三角形用于后续的匹配解算。
第六步,利用精匹配筛选的匹配三角形进行匹配解算。
根据两个点集序列间误差平方和最小的原理建立匹配解算的价值函数,将第一步中标准三角形对应的点序列与第五步(3)中精匹配过程筛选的匹配三角形点序列进行匹配解算,以上两个三角形间存在平移和旋转变换,可以通过匹配解算可求得以上两个三角形间的平移向量和旋转矩阵,该平移向量和旋转矩阵中包含惯导系统长时间累积的位置误差和航向角误差信息。
本发明原理说明如下:
无源重力导航就是利用局部重力异常分布的任意性,对惯导系统的累积误差进行校正。虽然惯导系统给出的导航信息存在随时间而不断增大的累积误差,但是在短时间内惯导能够提供高精度的相对距离信息,所以本发明充分利用惯导短时高精度的特点,选取惯导提供的三个连续位置点建立标准三角形,该三角形与载体对应的三个真实点构成的三角形是全等关系,所以他们之间可以通过旋转和平移进行相互转换;然后结合重力仪在三个对应位置点量测的实时重力值,在重力图中获取满足标准三角形约束的匹配三角形集,并通过粗匹配和精匹配过程筛选对应距离误差之和最小的匹配三角形进行最终的匹配解算;最后基于两组点集序列间误差平方和最小的原理建立价值函数,并通过微分求解获得惯导系统长时间累积的位置误差和航向角误差。
附图说明
图1为本发明一种基于三角形的重力图匹配方法流程图;
图2为本发明中标准三角形与载体真实位置点的关系图;
图3为本发明粗匹配过程中逐级搜索的四级块阵列;
图4为本发明粗匹配过程建立的三角形集;
图5为本发明匹配三角形集中存在的奇异三角形;
图6为本发明三角形的旋转和平移原理图;
具体实施方式
本发明将按照以下步骤实现重力图的匹配过程:
第一步,根据导航信息得出标准三角形。
从惯导系统的导航解算结果中提取经度、纬度作为位置点的信息,进而连续获取载体的三个位置点信息,以其两两间的航行距离作为标准三角形的三条边。如图2所示:运载体在三个连续时刻tk、tk+1和tk+2对应的真实位置点分别为A'、B'、C',而由惯导解算出的对应位置点分别为A(λA,LA)、B(λB,LB)、C(λC,LC),其中λA、λB、λC分别表示A、B、C三点的经度,LA、LB、LC分别表示A、B、C三点的纬度,图中细实线表示载体的真实轨迹,点划线表示惯导指示的轨迹,则载体两两位置点间的航行距离分别为AB、BC、AC。根据球面三角形的公式:
cosS1=cosS2cosS3+sinS2sinS3cosA1
S1=arccos(sinLAsinLB+cosLAcosLBcos(λB-λA))
于是得到边AB的长度为:
AB=Re·S1=Re·arccos(sinLAsinLB+cosLAcosLBcos(λB-λA))
其中Re为地球半径,同理得到BC的长度和AC的长度计算公式为:
BC=Re·arccos(sinLBsinLC+cosLBcosLCcos(λC-λB))
AC=Re·arccos(sinLAsinLC+cosLAcosLCcos(λC-λA))
这样根据惯导系统提供的位置信息,计算得到了标准三角形ΔABC。
第二步,结合实时测量重力值,在重力图中进行大范围搜索,获取粗匹配点集。
在已知重力图中建立块阵列进行大范围搜索,块阵列和其中的每个块都是正方形,本发明选取10×10的重力图网格方阵做为基本块(即e=10),块阵列按每边上块的个数定义其级数n,块阵列的最大级数nmax根据惯导系统的估计位置误差确定。
如图3所示:(a)、(b)、(c)、(d)分别为一级、二级、三级、四级的块阵列,块阵列中央空心点表示惯导提供的位置点(如A、B、C),图中曲线表示重力图中的重力等值线,粗实曲线表示与中心位置点的测量重力值相等的等值线,粗实线上的黑实点表示对应块内的粗匹配点,每一级阵列的搜索范围都只限于阵列边缘的块(阴影部分块)。
(1)根据惯导的误差特性确定最大的块阵列级数:
其中,tk为位置点A'对应的时刻,δ为惯导常值漂移误差,单位是°/h,σ为重力图分辨率,单位是°/格,e为块边长,[·]为取整运算。
确定了块阵列的最大级数后,在重力图中分别以位置点A、B、C为中心建立三个块阵列。
(2)结合载体在A'、B'、C'量测的重力值gA、gB、gC,以块阵列级数从小到大的顺序进行逐级搜索,如果在某一级阵列中搜索到重力值的等值点,则停止下一级的搜索。
(3)分别求取单个块中所有等值点的平均值作为该块的代表点,没有等值点的块没有代表点,由所有块的代表点组成粗匹配的点集。这样可以获取三个位置点A、B、C对应的粗匹配点集P1,P2,P3分别为:
其中n1,n2,n3分别表示搜索A、B、C对应等值点的块阵列级数,imax,jmax,kmax分别表示三组点集P1,P2,P3中点的个数。
第三步,通过标准三角形的约束,建立匹配三角形集,并对其进行筛选,其具体步骤为:
(1)以第一步标准三角形的三边距离为约束条件,利用第二步粗匹配的三组点集建立一组匹配三角形集。如图4所示,假设位置点A在块阵列第二级搜到两个待匹配点,位置点B也在块阵列第二级搜到两个待匹配点,位置点C在块阵列第二级搜到一个待匹配点,则以位置点A、B、C(空心点)对应的三个点集(实心点)排列组合可生成粗匹配的匹配三角形集{ΔP1,iP2,jP3,k}。
三角形三边距离的约束条件为:
其中,R为距离误差限。
距离误差的来源包括:块单元引入的误差、重力观测值的误差和重力图网格精度引入的误差。粗匹配过程中的主要误差是由块单元引入的误差,所以误差限R的确定如下:
其中,π=3.1415926,Re=6378165m为地球半径。
(2)计算第三步(1)中各个匹配三角形和标准三角形ΔABC对应三边距离的误差之和,并将所有匹配三角形按误差之和从小到大排序。其误差之和Ml的计算公式为:
M1=|P1,iP2,j-AB|+|P2,jP3,k-BC+|+|P1,iP3,k-AC|
其中,l表示匹配三角形点集中的第l个三角形。
(3)剔除奇异三角形,并筛选误差之和最小的匹配三角形用于精匹配过程。
针对匹配三角形集中存在的翻转奇异三角形,其剔除条件为:
如图5所示,ΔP1,1P2,1P3,1和ΔP1,2P2,2P3,2都满足第三步(1)中的约束条件,但是ΔP1,2P2,2P3,2与标准三角形ΔABC存在翻转的对应关系,所以是一个奇异三角形。
通过上式剔除了奇异三角形后,选择距离误差之和最小的三角形(定义为ΔP1,AP2,BP3,C)用于后续精匹配过程。
第四步,基于粗匹配筛选结果,在重力图中进行小范围搜索,获取精匹配点集。
以第三步(3)中筛选的匹配三角形的三个点P1,A、P2,B、P3,C分别为搜索中心,在重力图中建立三个独立的块进行小范围搜索,其块的大小与第二步的块大小相等,然后仍结合第一步载体在A'、B'、C'三个位置点的重力值gA、gB、gC,分别在三个块中搜索等值点,由这些等值点直接构成精匹配的三组点集P1',P2',P3'分别为:
其中,n1',n2',n3'分别表示搜索P1,A、P2,B、P3,C对应等值点的块阵列级数。
第五步,通过标准三角形的约束,建立新的匹配三角形集,并进行筛选。
(1)以第一步标准三角形的三边距离为约束条件,利用第四步精匹配的三组点集建立一组新的匹配三角形集{ΔP1,i′P2,j′P3,k′}。其约束条件为:
其中R'为精匹配三边距离的误差限。
误差限R'由重力测量误差和重力图网格精度来确定:
其中,ξ为重力仪的测量误差(单位mGal),τ为匹配点附近重力图单位格对应重力观测值的变化量(单位mGal/格)。
(2)计算第五步(1)中各新匹配三角形和标准三角形ΔABC对应三边的误差之和,并将所有新的匹配三角形按误差之和从小到大排序。其误差之和Ml′的公式为:
M1'=|P1,i'P2,j′-AB|+|P2,j'P3,k′-BC|+|P1,i′P3,k-AC|
(3)筛选误差之和最小的匹配三角形(定义为ΔP1,A'P2,B'P3,C')作为精匹配的结果,用于后续的匹配解算。
第六步,利用精匹配筛选的匹配三角形进行匹配解算。
根据两组点集序列间误差平方和最小的原理建立匹配解算的价值函数,将第一步中惯导提供的标准三角形对应点序列与精匹配筛选的匹配三角形对应点序列进行匹配解算,求得两个三角形间的位置平移向量和旋转矩阵。以下给出平移向量和旋转矩阵的求解过程:
(2)定义三角形匹配序列和价值函数。
如图6所示,为方便表示,将标准三角形定义为ΔQAQBQC(即上述的ΔABC),将精匹配给出的最终匹配三角形定义为ΔPAPBPC(即第五步(3)中ΔP1,A'P2,B'P3,C'),而标准三角形ΔQAQBQC通过旋转和平移后对应的三角形定义为ΔQA'QB'QC'。设三维矩阵R为ΔQAQBQC的旋转矩阵(绕ΔQAQBQC的重心旋转,逆时针旋转为正方向),同时定义d=[ΔλΔLΔh]T为ΔQAQBQC旋转后到ΔQA'QB'QC'的平移量,其中Δλ表示惯导的经度误差,ΔL表示惯导的纬度误差,Δh表示惯导的高度误差。
最佳匹配的准则是:两个三角形对应点的距离差的平方和最小,即标准三角形ΔQAQBQC经旋转平移后的三角形ΔQA'QB'QC'与ΔPAPBPC对应点的距离差的平方和达到最小值。其中,定义ΔPAPBPC各点的坐标向量为pi,ΔQAQBQC各点的坐标向量为qi,ΔQA'QB'QC'各点的坐标向量为qi'。则价值函数T(d,R)为:
同时,以点q0为旋转中心的旋转平移公式可以为:
qi'=d+q0+R(qi-q0)
于是,最终的价值函数为:
(2)求解平移向量d。
为了让目标函数取最小值,平移向量d应满足:
即
所以求解平移向量d得
d=p0-q0
(3)求解旋转矩阵R。
定义 则
同时分别定义 由于MV为常量,要使目标函数取最小值,即求取MV的最大值。根据矩阵分析知识有:
将S奇异值分解为S=UWVT,其中U,V都是正交矩阵,则有
Trace(SRT)=Trace(UWVTRT)
=Trace(UWVTRTUUT)
=Trace(WVTRTU)
令Z=VTRTU,由于U,V,R都是正交矩阵,所以矩阵Z也是正交矩阵则
Trace(SRT)=Trace(WVTRTU)
=Trace(WZ)
≤Trace(W)
当且仅当Z=VTRTU=I时,Trace(SRT)取得最大值,即MV取得最大值,目标函数T(d,R)取得最小值。于是可求得旋转矩阵R的计算公式为:
R=UVT
这样就解算出了从标准三角形ΔQAQBQC到匹配三角形ΔPAPBPC对应的平移向量d和旋转矩阵R。而平移向量d和旋转矩阵R与惯导系统长时间的位置误差和航向角误差Δθ存在如下关系式:
本发明说明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。
Claims (2)
1.一种基于三角形的重力图快速匹配方法,其特征步骤如下:
第一步,根据导航信息得出标准三角形;
从惯导系统的导航解算结果中提取经度、纬度作为位置点的信息,进而连续获取载体的三个位置点信息,以其两两间的航行距离作为标准三角形的三条边;
第二步,结合实时重力量测值,在重力图中进行大范围搜索,获取粗匹配点集;
在已知重力图中建立块阵列进行大范围搜索,块阵列和其中的每个块都是正方形,块阵列按每边上块的个数定义其级数,块阵列的最大级数根据惯导系统的估计位置误差确定;
以第一步中的三个位置点分别为中心在重力图中建立三个块阵列,然后结合第一步三个位置点对应的实时重力值,分别在三个块阵列中以块为最小单元搜索等值点,其搜索过程按照块阵列级数从小到大的顺序进行逐级搜索,由这些搜索的等值点构成粗匹配的三组点集;
第三步,通过标准三角形的约束,建立匹配三角形集,并对其进行筛选,其具体步骤为:
(1)以第一步标准三角形的三边距离为约束条件,将第二步粗匹配的三组点集按照排列组合的原则建立一组匹配三角形集,其约束条件中三边距离的误差限由第二步块阵列的单个块大小来确定;
(2)计算第三步(1)中各个匹配三角形和标准三角形对应三边的误差之和,并将所有匹配三角形按误差之和从小到大排序;
(3)筛选误差之和最小的匹配三角形用于后续的精匹配过程;
第四步,基于粗匹配筛选结果,在重力图中进行小范围搜索,获取精匹配点集;
以第三步(3)中筛选的匹配三角形的三个点分别为搜索中心,在重力图中建立单个块进行小范围搜索,其块的大小与第二步的块大小相等,然后仍结合第一步三个位置点对应的实时重力值,分别在三个块中搜索等值点,由这些等值点直接构成精匹配的三组点集;
第五步,通过标准三角形的约束,建立新的匹配三角形集,并进行筛选,其具体步骤为:
(1)以第一步标准三角形的三边距离为约束条件,利用第四步精匹配的三组点集建立一组新的匹配三角形集,其约束条件中三边距离的误差限由重力测量误差和重力图网格精度来确定;
(2)计算第五步(1)中各个新的匹配三角形和标准三角形对应三边的误差之和,并将所有新的匹配三角形按误差之和从小到大排序;
(3)筛选误差之和最小的匹配三角形用于后续的匹配解算;
第六步,利用第五步(3)中精匹配筛选的匹配三角形进行匹配解算;
根据两个点集序列间误差平方和最小的原理建立匹配解算的价值函数,将第一步中标准三角形对应的点序列与第五步精匹配过程筛选的匹配三角形点序列进行匹配解算,可求得两个三角形间的平移向量和旋转矩阵,该平移向量和旋转矩阵中包含惯导系统长时间累积的位置误差和航向角误差信息。
2.根据权利要求1所述的一种基于三角形的重力图快速匹配方法,其特征在于:所述第二步采用块阵列逐级搜索的方法获取粗匹配的点集,其步骤为:
(1)根据惯导的误差特性估计出惯导解算位置误差,进而确定出块阵列的最大级数,以第一步中惯导给出的位置点为中心,在重力图中建立块阵列;
(2)以块阵列的级数从小到大的顺序进行逐级搜索,如果在某一级阵列中搜索到实时重力值的等值点,则停止下一级的搜索;
(3)分别求取单个块中所有等值点的平均值作为该块的代表点,没有等值点的块没有代表点,由所有块的代表点组成粗匹配的点集。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210153486.3A CN102706348B (zh) | 2012-05-16 | 2012-05-16 | 一种基于三角形的重力图快速匹配方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210153486.3A CN102706348B (zh) | 2012-05-16 | 2012-05-16 | 一种基于三角形的重力图快速匹配方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102706348A true CN102706348A (zh) | 2012-10-03 |
CN102706348B CN102706348B (zh) | 2015-05-20 |
Family
ID=46899300
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210153486.3A Expired - Fee Related CN102706348B (zh) | 2012-05-16 | 2012-05-16 | 一种基于三角形的重力图快速匹配方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102706348B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104697523A (zh) * | 2015-03-31 | 2015-06-10 | 哈尔滨工业大学 | 基于迭代计算的惯性/地磁匹配定位方法 |
CN108225322A (zh) * | 2017-12-25 | 2018-06-29 | 北京理工大学 | 一种基于重力异常基准图的重力异常值连续化方法 |
CN108896040A (zh) * | 2018-03-30 | 2018-11-27 | 中国空间技术研究院 | 天空海一体化水下潜器惯性/重力组合导航方法和系统 |
CN110068829A (zh) * | 2019-04-11 | 2019-07-30 | 上海快仓智能科技有限公司 | 简化的基于图形快速匹配的激光定位方法 |
CN111044041A (zh) * | 2019-12-31 | 2020-04-21 | 北京理工大学 | 基于重力场三维特征的重力辅助惯性导航适配区选取方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102128625A (zh) * | 2010-12-08 | 2011-07-20 | 北京航空航天大学 | 重力辅助惯性导航系统中重力图匹配的初始匹配方法 |
CN102168979A (zh) * | 2010-12-08 | 2011-08-31 | 北京航空航天大学 | 一种基于三角形约束模型的无源导航的等值线匹配方法 |
-
2012
- 2012-05-16 CN CN201210153486.3A patent/CN102706348B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102128625A (zh) * | 2010-12-08 | 2011-07-20 | 北京航空航天大学 | 重力辅助惯性导航系统中重力图匹配的初始匹配方法 |
CN102168979A (zh) * | 2010-12-08 | 2011-08-31 | 北京航空航天大学 | 一种基于三角形约束模型的无源导航的等值线匹配方法 |
Non-Patent Citations (3)
Title |
---|
ZHAOXIA LIU ET AL.: "A new algorithm of global feature matching based on triangle regions for image registration", 《IEEE:ICSP2010 PROCEEDINGS》, 31 December 2010 (2010-12-31) * |
ZHENLI YANG ET AL.: "A Triangle Matching Algorithm for Gravity-aided Navigation for Underwater Vehicles", 《THE JOURNAL OF NAVIGATION》, vol. 67, 17 October 2013 (2013-10-17) * |
朱庄生等: "无源重力导航的三角形匹配算法及仿真", 《仪器仪表学报》, vol. 33, no. 10, 31 October 2012 (2012-10-31) * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104697523A (zh) * | 2015-03-31 | 2015-06-10 | 哈尔滨工业大学 | 基于迭代计算的惯性/地磁匹配定位方法 |
CN104697523B (zh) * | 2015-03-31 | 2017-05-31 | 哈尔滨工业大学 | 基于迭代计算的惯性/地磁匹配定位方法 |
CN108225322A (zh) * | 2017-12-25 | 2018-06-29 | 北京理工大学 | 一种基于重力异常基准图的重力异常值连续化方法 |
CN108225322B (zh) * | 2017-12-25 | 2021-04-27 | 北京理工大学 | 一种基于重力异常基准图的重力异常值连续化方法 |
CN108896040A (zh) * | 2018-03-30 | 2018-11-27 | 中国空间技术研究院 | 天空海一体化水下潜器惯性/重力组合导航方法和系统 |
CN108896040B (zh) * | 2018-03-30 | 2021-03-26 | 中国空间技术研究院 | 天空海一体化水下潜器惯性/重力组合导航方法和系统 |
CN110068829A (zh) * | 2019-04-11 | 2019-07-30 | 上海快仓智能科技有限公司 | 简化的基于图形快速匹配的激光定位方法 |
CN111044041A (zh) * | 2019-12-31 | 2020-04-21 | 北京理工大学 | 基于重力场三维特征的重力辅助惯性导航适配区选取方法 |
CN111044041B (zh) * | 2019-12-31 | 2021-04-27 | 北京理工大学 | 基于重力场三维特征的重力辅助惯性导航适配区选取方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102706348B (zh) | 2015-05-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102128625B (zh) | 重力辅助惯性导航系统中重力图匹配的初始匹配方法 | |
CN103047985B (zh) | 一种空间目标的快速定位方法 | |
CN106526590A (zh) | 一种融合多源sar影像工矿区三维地表形变监测及解算方法 | |
CN108061889A (zh) | Ais与雷达角度系统偏差的关联方法 | |
CN102508245B (zh) | 一种星载多频率与多基线InSAR高程估计精度等效性确定方法 | |
CN102706348B (zh) | 一种基于三角形的重力图快速匹配方法 | |
CN107065025A (zh) | 一种基于重力梯度不变量的轨道要素估计方法 | |
CN105093299A (zh) | 一种基于炮检距向量片技术优化观测系统的方法及装置 | |
CN104236546A (zh) | 一种卫星星光折射导航误差确定与补偿方法 | |
CN104061932A (zh) | 一种利用引力矢量和梯度张量进行导航定位的方法 | |
CN102096944A (zh) | 地质体结构面三维激光扫描点云识别方法 | |
CN107709926A (zh) | 自动化的移动岩土测绘 | |
CN102446354A (zh) | 一种高精度多源地面激光点云的整体配准方法 | |
CN106546235A (zh) | 一种基于载体补偿的磁性目标定位方法 | |
CN202209953U (zh) | 用于水下载体的地磁辅助惯性导航系统 | |
CN105740505A (zh) | 一种基于gps-rtk技术的道路空间线形恢复方法 | |
CN103954996B (zh) | 一种基于旅行时法确定地层裂隙裂缝走向的方法及装置 | |
CN106168682A (zh) | 一种基于旋转重力场的运动目标体监测方法 | |
CN101266153B (zh) | 测绘工程类陀螺全站仪精度评定方法 | |
CN103353612A (zh) | 一种地下目标物体的测量定位设备及测量定位方法 | |
CN107703527A (zh) | 一种基于北斗三频单历元宽巷/超宽巷的组合定位方法 | |
CN107340529A (zh) | 一种星载测频定位方法、装置和系统 | |
CN102236101A (zh) | 多层裂缝预测方法和装置 | |
CN112378376B (zh) | 一种基于传感阵列和测斜仪的海底变形联合监测方法 | |
CN106323271A (zh) | 基于特征奇异值的航天器相对姿态测量矢量选取方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150520 Termination date: 20210516 |