CN108871349A - 一种深空探测器光学导航位姿加权确定方法 - Google Patents
一种深空探测器光学导航位姿加权确定方法 Download PDFInfo
- Publication number
- CN108871349A CN108871349A CN201810766656.2A CN201810766656A CN108871349A CN 108871349 A CN108871349 A CN 108871349A CN 201810766656 A CN201810766656 A CN 201810766656A CN 108871349 A CN108871349 A CN 108871349A
- Authority
- CN
- China
- Prior art keywords
- crater
- pose
- deep space
- matrix
- equation
- 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
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/24—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for cosmonautical navigation
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- General Physics & Mathematics (AREA)
- Astronomy & Astrophysics (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开的一种深空探测器光学导航位姿加权确定方法,涉及深空探测器位姿确定方法,属于深空探测技术领域。本发明实现方法如下:在导航相机图像中利用图像处理算法检测和提取行星表面陨石坑信息;对检测提取的陨石坑信息进行椭圆拟合与陨石坑中心定位;利用拟合的椭圆方程,根据不确定性传播特性,求取测量误差协方差矩阵来描述陨石坑中心的定位不确定性;利用测量不确定性对陨石坑中心进行加权,使得不同陨石坑中心在位姿估计算法中的权值大小不同,进而提高位姿估计的精确性,即基于测量不确定性加权方法实现深空探测器位姿的确定。本发明能够提高位姿确定的精确性,具有实现方法简单的优点。本发明适用于利用陨石坑作为导航路标进行位姿确定的导航系统。
Description
技术领域
本发明涉及深空探测器位姿确定方法,具体涉及适用于利用陨石坑作为导航路标进行位姿确定的导航系统的深空探测器位姿确定方法,属于深空探测技术领域。
背景技术
深空探测航行距离远、时间长,采用传统的测控方式存在较大的测量延迟,技术操作难度大,深空探测通信成本高。传统的惯性导航虽不需要地面深空测控网的支持,但是星载惯性器件在使用过程中存在不同程度的零漂误差、随机误差以及计算误差积累,无法满足深空探测的精确导航任务要求。随着计算机硬件技术的突破和光学敏感器件的发展,使得基于星载计算机和光学导航相机的自主光学导航方法成为研究热点。其中,行星表面陨石坑作为天然地理地形特征路标,无需探测器额外携带路标载荷,有效减少任务复杂度,应用前景广泛。以陨石坑光学信息为主的位姿估计方法已经成为深空探测自主导航的重要方式,如何利用陨石坑光学信息来提高位姿估计精度是当前各国航天科研人员重点关注的研究方向之一。
在已发展的深空探测器位姿估计方法中,在先技术[1](Dementhon D F,Davis LS.Model-based object pose in 25lines of code[J].International Journal ofComputer Vision,1995,15(1-2):123-141.),针对单张图像,利用四个或更多非共面特征点,在迭代过程中反复地把透视投影模型近似为2次投影的合成,先通过正交投影把所有的三维物体点平行于光轴投影到一个与光轴垂直的合适的空间平面上,然后将该空间平面上的点再透视投影到像平面上,不断找到近似的姿态,直至收敛到精确的估计值。但是该算法没有考虑到用于估计的物体观测误差对迭代的影响,随着迭代次数的增加会产生误差积累的问题。
在先技术[2](参见W.Shao,et al.,Attitude and position determinationbased on craters for precision landing,Proc IMechE Part G:J AerospaceEngineering.2016,Vol.230(10)1934-1942.),针对光学自主导航问题,提出了一种实时的利用至少三条陨石坑提取边缘的椭圆曲线的位姿估计方法,利用椭圆曲线和二维图像的几何与代数约束,并利用克罗内克积建立了着陆器位置姿态的线性方程。但该方法并没有考虑陨石坑提取时边缘点误差对位姿估计的影响,因此该方法只适用于边缘提取足够精确且没有误匹配的情况。
在先技术[3](参见周润,张征宇,黄叙辉.相机位姿估计的加权正交迭代算法[J].光学学报,2018,38(05):193-199.),针对测量误差,提出了一种相机位姿估计的加权正交迭代算法,以加权共线误差为目标函数,根据像面重投影误差自动调节权重系数,优化相机位姿估计结果。但是该方法是以共线误差方程为基础加入的权重,并不适用于深空探测以陨石坑边缘特征点为观测量,最终要得到陨石坑中心不确定性的位姿确定。
利用陨石坑光学信息进行深空探测器位姿确定过程中,陨石坑边缘检测不可避免的存在观测误差,使得陨石坑中心存在定位不确定性,现有的深空探测器位姿确定方法未考虑陨石坑中心定位不确定性的影响,进而导致探测器位姿解算精度低。
发明内容
针对现有的深空探测器位姿确定方法未考虑陨石坑中心定位不确定性的影响,进而导致探测器位姿解算精度低的问题。本发明公开的一种深空探测器光学导航位姿加权确定方法要解决的技术问题是:基于测量不确定性加权方法实现深空探测器位姿的确定,通过考虑陨石坑中心定位的不确定性,将原始的陨石坑中心加入测量不确定性的权重,使得不同陨石坑中心在位姿估计算法中的权值大小不同,进而提高位姿确定的精确性,具有实现方法简单的优点。
本发明是通过下述技术方案实现的。
本发明公开的一种深空探测器光学导航位姿加权确定方法,在导航相机图像中利用图像处理算法检测和提取行星表面陨石坑信息。对检测提取的陨石坑信息进行椭圆拟合与陨石坑中心定位。利用拟合的椭圆方程,根据不确定性传播特性,求取测量误差协方差矩阵来描述陨石坑中心的定位不确定性。利用测量不确定性对陨石坑中心进行加权,使得不同陨石坑中心在位姿估计算法中的权值大小不同,进而提高位姿估计的精确性。即基于测量不确定性加权方法实现深空探测器位姿的确定。
本发明公开的一种深空探测器光学导航位姿加权确定方法,包括如下步骤:
步骤1:在导航相机图像中利用图像处理算法检测和提取行星表面陨石坑信息。
读取光学相机拍摄到的目标天体表面地形图像之后,基于图像分割和形态学处理算法等计算机图形处理技术对图像进行陨石坑边缘的检测和提取。
步骤1行星表面陨石坑信息包括图像处理得到的陨石坑边缘点的像素值。
步骤2:对步骤1中检测提取的陨石坑信息进行椭圆拟合与陨石坑中心定位。
当观测椭圆边界点不少于五个时,通过最小二乘法确定陨石坑拟合椭圆方程的系数B,C,D,E,F,进而确定椭圆定位定向参数,得到陨石坑边缘的拟合椭圆方程为x2+2Bxy+Cy2+2Dx+2Ey+F=0。拟合椭圆的中心(x0,y0)为
其中,B,C,D,E,F为椭圆方程的系数。
步骤3:利用步骤2中拟合的椭圆方程,根据不确定性传播特性,求取测量误差协方差矩阵来描述陨石坑中心的定位不确定性。
对n个边缘点写出n个椭圆方程,并用矩阵表示n个误差方程,则求得各特征点的协方差矩阵P,将各特征点视为具有相同方差σ2来处理。最后,根据协方差传播率,由P产生每个陨石坑中心的误差协方差矩阵来描述陨石坑中心的定位不确定性。
其中,K是由步骤1中的形心公式求偏导所得的多项式参数的系数B,C,D,E,F构成的矩阵。
作为优选,步骤3实现方法为:
对n个边缘点写出n个椭圆一般方程,并用矩阵表示误差方程为
令
则误差方程写为V=AX+B。
各特征点的协方差矩阵P为
其中,为V的测量噪声方差阵,I是n×n的单位矩阵,σ2为n个特征点的方差,将各特征点视为具有相同方差σ2来处理。最后,根据协方差传播率,由P产生各陨石坑中心的误差协方差矩阵来描述陨石坑中心的定位不确定性。
其中,K是由步骤2中的形心公式求偏导所得的多项式参数的系数B,C,D,E,F构成的矩阵。
将每个陨石坑中心计算得到的误差椭圆画在拍摄的导航图像上,能够直观地显示出拟合陨石坑中心的方向和大小不确定性。
步骤4:利用测量不确定性对陨石坑中心进行加权并解算位置姿态,即基于测量不确定性加权方法实现深空探测器位姿的确定。
对步骤3中求取的测量误差协方差矩阵进行奇异值分解,
利用协方差矩阵定义仿射变换矩阵F
F=diag(1/σ1,1/σ2,1)UT
(u′,v′,1)T=F(u,v,1)T
经过仿射变换后,得到加权后的像素坐标点(u′,v′)T,结合已建立的地形数据库中上述像素点对应的陨石坑三维坐标点,将这些数据输入到改进的比例正交投影迭代变换POSIT算法中,即得到更为精确的位姿估计解,即基于测量不确定性加权方法实现深空探测器位姿的确定。
步骤4所述的输入到改进的比例正交投影迭代变换POSIT算法中即得到更为精确的位姿估计解,具体实现方法如下:
待求的探测器的旋转矩阵R和平移矩阵t分别为
定义检测得到的陨石坑中心在导航相机图像平面像素坐标系下的齐次坐标为Pi=(ui,vi,1)T,非齐次坐标为像点对应的陨石坑中心在目标天体固连坐标系下的匹配得到的齐次坐标为非齐次坐标为其中,i=1,2,…,n,n为陨石坑匹配数。图像点Pi和数据库匹配点Fi w的透视投影关系为
λiPi=K[R t]Pi w
式中,为摄像机内参矩阵,f为导航相机焦距,λi为数据库匹配点Pi w的深度值。
定义(x,y)是特征点在图像平面物理坐标系的像点坐标,(x′,y′)是特征点经过仿射变换后在图像平面物理坐标系的像点坐标,(XC,YC,ZC)是特征点在光学导航相机坐标系的坐标,特征点和像点的关系为
则透视投影关系表示为
上式两边同时除以tz得到
式中,λ=ZC/tz,s=f/tz
加入不确定性权阵,将等式左右两边分别左乘仿射变换矩阵得
展开得
由于在探测器任务期间,行星表面各特征点在导航相机坐标系中的z轴变化范围,所述z轴变化范围远远小于行星表面特征点到相机的平均深度。因此认为λ=ZC/tz≈1。
λ=ZC/tz≈1是比例正交投影迭代变换POSIT算法的迭代初值。在λ=ZC/tz≈1初始状态下,物体的所有点在同一个深度上时,透视变换退化为一个比例正交投影。由于给定λ一个估计值,因此先将其看做已知量,删掉第三行方程中就少了4个未知量,更方便求解,即得到
展开得迭代方程为
求解8个未知量[sR11 sR12 sR13 sR21 sR22 sR23 stx sty]T利用最小二乘法求解如下方程
由于λ被视为已知量,一个特征点及对应图像坐标得到2个独立方程,方程中总共有8个独立变量,拆成3个向量分别为
由此计算得到向量sR1,sR2,同时R1,R2均为单位向量。由于存在计算误差,向量sR1,sR2的模长不一定相等,取因此计算得到
tz=f/s
tx=(stx)/s
ty=(sty)/s
R1=(sR1)/s
R2=(sR2)/s
R3=R1×R2
由于迭代初始时,λ≈1,通过上式得到初始的旋转矩阵和平移矩阵,代入下式
λ=[R3T/tz 1][XW YW ZW 1]T
通过迭代更新使每个特征点都有更加精确的迭代系数λ,将更新值代入迭代方程中,循环往复,逐步逼近位姿估计的精确解,直至旋转矩阵R和平移向量t稳定时停止迭代。
还包括步骤5:将步骤4基于测量不确定性加权方法确定的深空探测器位姿,用于目标天体表面存在陨石坑导航路标的深空探测器位姿解算,进而实现用于目标天体表面存在陨石坑导航路标的深空探测器的导航。
有益效果:
1、本发明公开的一种深空探测器光学导航位姿加权确定方法,由于考虑陨石坑中心定位的不确定性,将原始的陨石坑中心加入测量不确定性的权重,不确定性大的权重小,不确定性小的权重大,使得不同陨石坑中心在位姿估计算法中的权值大小不同,进而提高位姿估计的精确性。
2、本发明公开的一种深空探测器光学导航位姿加权确定方法,能自主地确定探测器的位置和姿态,能够用于目标天体表面存在陨石坑导航路标的深空探测器,实现用于目标天体表面存在陨石坑导航路标的深空探测器的导航任务。
附图说明
图1是本发明的一种深空探测器光学导航位姿加权确定方法流程示意图;
图2是本发明位姿确定过程中导航相机与像平面和行星表面关系示意图;
图3是本发明步骤1图像处理效果示意图,其中图3a)是导航相机拍摄的原始图像,图3b)是原始图像经过图像处理之后的效果图;
图4是本发明步骤2椭圆拟合定位效果示意图;
图5是本发明步骤3陨石坑边缘拟合椭圆与陨石坑中心误差椭圆示意图;
图6是本发明实例中步骤3模拟的陨石坑边缘的拟合椭圆与陨石坑中心误差椭圆的像平面图。
具体实施方式
为了更好的说明本发明的目的和优点,下面结合附图和实例对发明内容做进一步说明。
为了验证本发明的可行性,利用拍摄到的小行星Eros 433表面陨石坑图像信息,将陨石坑经图像处理之后的拟合椭圆及中心作为数学仿真的真实陨石坑边缘和中心,如图5所示,并定义各陨石坑中心在小天体固连坐标系下的深度,得到像点对应的陨石坑中心在世界坐标系下的三维坐标取噪声值为5来模拟图像处理过程中的误差,并假设世界坐标系原点为[500,-300,2000]T,探测器在世界坐标系下的初始位置为[0,0,2000]T,初始姿态为[0°,180°,0°]。视场角大小为30°,导航相机焦距为f=8mm,进行数学仿真验证。
本实施例公开的一种深空探测器光学导航位姿加权确定方法,包括如下步骤:
步骤1:在导航相机图像中利用图像处理算法检测和提取行星表面陨石坑信息。
读取光学相机拍摄到的Eros 433表面地形图像之后,对图像进行高斯噪声滤除,利用LOG算子自动选择阈值进行边缘检测,之后基于八连通法保留连通区域,剔除非边缘点,将光照方向作为先验信息,最后利用光照方向与梯度方向的内积去除伪边缘,完成对24个陨石坑边缘的检测和提取。
步骤2:对步骤1中检测提取的陨石坑信息进行椭圆拟合与陨石坑中心定位。
通过最小二乘法确定24组陨石坑的拟合椭圆方程的系数Bi,Ci,Di,Ei,Fi,i=1,2,3...24,进而确定椭圆定位定向参数,得到陨石坑边缘的拟合椭圆方程为x2+2Bxy+Cy2+2Dx+2Ey+F=0。拟合椭圆的中心(x0,y0)为
步骤3:利用步骤2中拟合的椭圆方程,根据不确定性传播特性,求取测量误差协方差矩阵来描述陨石坑中心的定位不确定性。
对每个陨石坑,n个边缘点写出n个椭圆方程,并用矩阵表示n个误差方程,则求得各特征点的协方差矩阵P,将各特征点视为具有相同方差σ2=5·In×n来处理。最后,根据协方差传播率,由P产生每个陨石坑中心的误差协方差矩阵来描述陨石坑中心的定位不确定性。
其中,K是由步骤1中的形心公式求偏导所得的多项式参数的系数B,C,D,E,F构成的矩阵。
步骤3具体实现方法为:
对n个边缘点写出n个椭圆一般方程,并用矩阵表示误差方程为
令
则误差方程写为V=AX+B。
各特征点的协方差矩阵P为
其中,为V的测量噪声方差阵,I是n×n的单位矩阵,σ2为n个特征点的方差。最后,根据协方差传播率,由P产生各陨石坑中心的误差协方差矩阵
其中,K是由步骤2中的形心公式求偏导所得的多项式参数的系数B,C,D,E,F构成的矩阵。
如图6所示,将拟合椭圆和中心误差椭圆画在模拟的真实陨石坑图像上,能够直观地显示出拟合陨石坑中心的方向和大小不确定性。
步骤4:利用测量不确定性对陨石坑中心进行加权并解算位置姿态,即基于测量不确定性加权方法实现深空探测器位姿的确定。
对24个陨石坑,利用步骤3中求取的测量误差协方差矩阵进行奇异值分解,
利用协方差矩阵求仿射变换矩阵F
(u′,v′,1)T=F(u,v,1)T
经过仿射变换后,得到加权后的像素坐标点(u′,v′)T,结合已定义的像素点对应的陨石坑三维坐标点,将这些数据输入到改进的比例正交投影迭代变换POSIT算法中,即得到更为精确的位姿估计解,即基于测量不确定性加权方法实现深空探测器位姿的确定。
步骤4所述的输入到改进的比例正交投影迭代变换POSIT算法中即得到更为精确的位姿估计解,具体实现方法如下:
设待求的探测器的旋转矩阵R和平移矩阵t分别为
定义检测得到的陨石坑中心在导航相机图像平面像素坐标系下的齐次坐标为Pi=(ui,vi,1)T,非齐次坐标为像点对应的陨石坑中心在目标天体固连坐标系下的齐次坐标为非齐次坐标为其中,i=1,2,…,24。图像点Pi和数据库匹配点Fi w的透视投影关系为
λiPi=K[R t]Pi w
式中,为摄像机内参矩阵,f为导航相机焦距,λi为Pi w的深度值。
定义(x,y)是特征点在图像平面物理坐标系的像点坐标,(x′,y′)是特征点经过仿射变换后在图像平面物理坐标系的像点坐标,(XC,YC,ZC)是特征点在光学导航相机坐标系的坐标,特征点和像点的关系为
则透视投影关系表示为
上式两边同时除以tz得到
式中,λ=ZC/tz,s=f/tz
加入不确定性权阵,将等式左右两边分别左乘仿射变换矩阵得
展开得
由于在探测器任务期间,行星表面各特征点在导航相机坐标系中的z轴变化范围,所述z轴变化范围远远小于行星表面特征点到相机的平均深度。因此认为λ=ZC/tz≈1。
λ=ZC/tz≈1是比例正交投影迭代变换POSIT算法的迭代初值。在λ=ZC/tz≈1初始状态下,物体的所有点在同一个深度上时,透视变换退化为一个比例正交投影。由于给定λ一个估计值,因此先将其看做已知量,删掉第三行方程中就少了4个未知量,更方便求解,即得到
展开得迭代方程为
求解8个未知量[sR11 sR12 sR13 sR21 sR22 sR23 stx sty]T,利用最小二乘法求解如下方程
由于λ被视为已知量,一个特征点及对应图像坐标得到2个独立方程,方程中总共有8个独立变量,拆成3个向量分别为
由此计算得到向量sR1,sR2,同时R1,R2均为单位向量。由于存在计算误差,向量sR1,sR2的模长不一定相等,取因此计算得到
tz=f/s
tx=(stx)/s
ty=(sty)/s
R1=(sR1)/s
R2=(sR2)/s
R3=R1×R2
由于迭代初始时,λ≈1,通过上式得到初始的旋转矩阵和平移矩阵,代入下式
通过迭代更新使每个特征点都有更加精确的迭代系数λ,将更新值代入迭代方程中,循环往复,逐步逼近位姿估计的精确解,直至旋转矩阵R和平移向量t稳定时停止迭代。
通过1000次的蒙特卡罗仿真得到两种方法的位姿估计精度如下表所示
表1探测器位姿估计结果
姿态误差 | 位置误差 | |
经典算法 | [-1.0519,-0.5774,0.0672]T | [0.8291,-0.7906,2.1456]T |
加权算法 | [-0.0198,0.6649,0.0525]T | [0.5275,0.1482,1.6599]T |
由计算结果,本发明提出的方法在相同条件下比传统方法姿态精度能够提高44.48%,位置精度提高28.13%,仿真结果验证了本发明的有效性。由此可知,不确定性越大,陨石坑中心误差就越大,该点在位姿算法中的权重就越小,其误差对整体估算的影响就越小,深空探测器位姿估计的精度就越高。
包括步骤5:将步骤4基于测量不确定性加权方法确定的深空探测器位姿,用于目标天体表面存在陨石坑导航路标的深空探测器位姿解算,进而实现用于目标天体表面存在陨石坑导航路标的深空探测器的导航。
以上所述的具体描述,对发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (8)
1.一种深空探测器光学导航位姿加权确定方法,其特征在于:包括如下步骤,
步骤1:在导航相机图像中利用图像处理算法检测和提取行星表面陨石坑信息;
步骤2:对步骤1中检测提取的陨石坑信息进行椭圆拟合与陨石坑中心定位;
步骤3:利用步骤2中拟合的椭圆方程,根据不确定性传播特性,求取测量误差协方差矩阵来描述陨石坑中心的定位不确定性;
步骤4:利用测量不确定性对陨石坑中心进行加权并解算位置姿态,即基于测量不确定性加权方法实现深空探测器位姿的确定。
2.如权利要求1所述的一种深空探测器光学导航位姿加权确定方法,其特征在于:还包括步骤5:将步骤4基于测量不确定性加权方法确定的深空探测器位姿,用于目标天体表面存在陨石坑导航路标的深空探测器位姿解算,进而实现用于目标天体表面存在陨石坑导航路标的深空探测器的导航。
3.如权利要求1或2所述的一种深空探测器光学导航位姿加权确定方法,其特征在于:步骤1具体实现方法为,
读取光学相机拍摄到的目标天体表面地形图像之后,基于图像分割和形态学处理算法对图像进行陨石坑边缘的检测和提取;
步骤1行星表面陨石坑信息包括图像处理得到的陨石坑边缘点的像素值。
4.如权利要求3所述的一种深空探测器光学导航位姿加权确定方法,其特征在于:步骤2具体实现方法为,
当观测椭圆边界点不少于五个时,通过最小二乘法确定陨石坑拟合椭圆方程的系数B,C,D,E,F,进而确定椭圆定位定向参数,得到陨石坑边缘的拟合椭圆方程为x2+2Bxy+Cy2+2Dx+2Ey+F=0;拟合椭圆的中心(x0,y0)为
其中,B,C,D,E,F为椭圆方程的系数。
5.如权利要求4所述的一种深空探测器光学导航位姿加权确定方法,其特征在于:步骤3具体实现方法为,
对n个边缘点写出n个椭圆方程,并用矩阵表示n个误差方程,则求得各特征点的协方差矩阵P,将各特征点视为具有相同方差σ2来处理;最后,根据协方差传播率,由P产生每个陨石坑中心的误差协方差矩阵来描述陨石坑中心的定位不确定性;
其中,K是由步骤1中的形心公式求偏导所得的多项式参数的系数B,C,D,E,F构成的矩阵。
6.如权利要求5所述的一种深空探测器光学导航位姿加权确定方法,其特征在于:步骤4具体实现方法为,
对步骤3中求取的测量误差协方差矩阵进行奇异值分解,
利用协方差矩阵定义仿射变换矩阵F
F=diag(1/σ1,1/σ2,1)UT
(u′,v′,1)T=F(u,v,1)T
经过仿射变换后,得到加权后的像素坐标点(u′,v′)T,结合已建立的地形数据库中上述像素点对应的陨石坑三维坐标点,将这些数据输入到改进的比例正交投影迭代变换POSIT算法中,即得到更为精确的位姿估计解,即基于测量不确定性加权方法实现深空探测器位姿的确定。
7.如权利要求6所述的一种深空探测器光学导航位姿加权确定方法,其特征在于:步骤3优选实现方法为,
对n个边缘点写出n个椭圆一般方程,并用矩阵表示误差方程为
令
则误差方程写为V=AX+B;
各特征点的协方差矩阵P为
其中,为V的测量噪声方差阵,I是n×n的单位矩阵,σ2为n个特征点的方差,将各特征点视为具有相同方差σ2来处理;最后,根据协方差传播率,由P产生各陨石坑中心的误差协方差矩阵来描述陨石坑中心的定位不确定性;
其中,K是由步骤2中的形心公式求偏导所得的多项式参数的系数B,C,D,E,F构成的矩阵;
将每个陨石坑中心计算得到的误差椭圆画在拍摄的导航图像上,能够直观地显示出拟合陨石坑中心的方向和大小不确定性。
8.如权利要求7所述的一种深空探测器光学导航位姿加权确定方法,其特征在于:步骤4所述的输入到改进的比例正交投影迭代变换POSIT算法中即得到更为精确的位姿估计解,具体实现方法如下,
待求的探测器的旋转矩阵R和平移矩阵t分别为
定义检测得到的陨石坑中心在导航相机图像平面像素坐标系下的齐次坐标为Pi=(ui,vi,1)T,非齐次坐标为像点对应的陨石坑中心在目标天体固连坐标系下的匹配得到的齐次坐标为非齐次坐标为其中,i=1,2,…,n,n为陨石坑匹配数;图像点Pi和数据库匹配点Fi w的透视投影关系为
λiPi=K[R t]Pi w
式中,为摄像机内参矩阵,f为导航相机焦距,λi为数据库匹配点Pi w的深度值;
定义(x,y)是特征点在图像平面物理坐标系的像点坐标,(x′,y′)是特征点经过仿射变换后在图像平面物理坐标系的像点坐标,(XC,YC,ZC)是特征点在光学导航相机坐标系的坐标,特征点和像点的关系为
则透视投影关系表示为
上式两边同时除以tz得到
式中,λ=ZC/tz,s=f/tz
加入不确定性权阵,将等式左右两边分别左乘仿射变换矩阵得
展开得
由于在探测器任务期间,行星表面各特征点在导航相机坐标系中的z轴变化范围,所述z轴变化范围远远小于行星表面特征点到相机的平均深度;因此认为λ=ZC/tz≈1;
λ=ZC/tz≈1是比例正交投影迭代变换POSIT算法的迭代初值;在λ=ZC/tz≈1初始状态下,物体的所有点在同一个深度上时,透视变换退化为一个比例正交投影;由于给定λ一个估计值,因此先将其看做已知量,删掉第三行方程中就少了4个未知量,更方便求解,即得到
展开得迭代方程为
求解8个未知量[sR11 sR12 sR13 sR21 sR22 sR23 stx sty]T利用最小二乘法求解如下方程
由于λ被视为已知量,一个特征点及对应图像坐标得到2个独立方程,方程中总共有8个独立变量,拆成3个向量分别为
由此计算得到向量sR1,sR2,同时R1,R2均为单位向量;由于存在计算误差,向量sR1,sR2的模长不一定相等,取因此计算得到
tz=f/s
tx=(stx)/s
ty=(sty)/s
R1=(sR1)/s
R2=(sR2)/s
R3=R1×R2
由于迭代初始时,λ≈1,通过上式得到初始的旋转矩阵和平移矩阵,代入下式
通过迭代更新使每个特征点都有更加精确的迭代系数λ,将更新值代入迭代方程中,循环往复,逐步逼近位姿估计的精确解,直至旋转矩阵R和平移向量t稳定时停止迭代。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810766656.2A CN108871349B (zh) | 2018-07-13 | 2018-07-13 | 一种深空探测器光学导航位姿加权确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810766656.2A CN108871349B (zh) | 2018-07-13 | 2018-07-13 | 一种深空探测器光学导航位姿加权确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108871349A true CN108871349A (zh) | 2018-11-23 |
CN108871349B CN108871349B (zh) | 2021-06-15 |
Family
ID=64301357
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810766656.2A Active CN108871349B (zh) | 2018-07-13 | 2018-07-13 | 一种深空探测器光学导航位姿加权确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108871349B (zh) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109697734A (zh) * | 2018-12-25 | 2019-04-30 | 浙江商汤科技开发有限公司 | 位姿估计方法及装置、电子设备和存储介质 |
CN109708643A (zh) * | 2019-01-14 | 2019-05-03 | 北京理工大学 | 小行星表面光学导航路标评价选取方法 |
CN110132281A (zh) * | 2019-05-21 | 2019-08-16 | 哈尔滨工程大学 | 一种基于询问应答模式的水下高速目标高精度自主声学导航方法 |
CN110146093A (zh) * | 2019-06-19 | 2019-08-20 | 北京理工大学 | 双体小行星探测自主协同光学导航方法 |
CN110619368A (zh) * | 2019-09-23 | 2019-12-27 | 北京理工大学 | 行星表面导航特征成像匹配检测方法 |
CN111721302A (zh) * | 2020-06-24 | 2020-09-29 | 北京理工大学 | 不规则小行星表面复杂地形特征识别与感知方法 |
CN111735460A (zh) * | 2020-08-05 | 2020-10-02 | 北京控制与电子技术研究所 | 基于小天体中心提取的航天器导航方法、系统及装置 |
CN111861873A (zh) * | 2020-07-20 | 2020-10-30 | 北京航天飞行控制中心 | 仿真图像的生成方法和装置 |
CN112365541A (zh) * | 2020-11-24 | 2021-02-12 | 北京航空航天大学青岛研究院 | 一种基于相似变换的大场景相机姿态注册方法 |
CN112819794A (zh) * | 2021-02-04 | 2021-05-18 | 青岛科技大学 | 一种基于深度学习的小天体陨石坑检测方法 |
CN112906573A (zh) * | 2021-02-10 | 2021-06-04 | 北京理工大学 | 基于轮廓点集的行星表面导航路标匹配方法 |
CN114683279A (zh) * | 2022-02-17 | 2022-07-01 | 广东博智林机器人有限公司 | 导航精度的确定方法、装置和电子设备 |
CN116306185A (zh) * | 2023-02-14 | 2023-06-23 | 北京理工大学 | 基于不规则离散元的陨星进入地球大气层后运动仿真方法及系统 |
CN116524031A (zh) * | 2023-07-03 | 2023-08-01 | 盐城数智科技有限公司 | 一种基于yolov8的大范围月球车定位建图方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1847792A (zh) * | 2006-05-12 | 2006-10-18 | 哈尔滨工业大学 | 深空探测器软着陆视觉导航算法模拟测试装置 |
CN102999915A (zh) * | 2012-12-03 | 2013-03-27 | 哈尔滨工业大学 | 一种基于面积比的陨石坑匹配方法 |
CN107144278A (zh) * | 2017-04-24 | 2017-09-08 | 北京理工大学 | 一种基于多源特征的着陆器视觉导航方法 |
-
2018
- 2018-07-13 CN CN201810766656.2A patent/CN108871349B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1847792A (zh) * | 2006-05-12 | 2006-10-18 | 哈尔滨工业大学 | 深空探测器软着陆视觉导航算法模拟测试装置 |
CN102999915A (zh) * | 2012-12-03 | 2013-03-27 | 哈尔滨工业大学 | 一种基于面积比的陨石坑匹配方法 |
CN107144278A (zh) * | 2017-04-24 | 2017-09-08 | 北京理工大学 | 一种基于多源特征的着陆器视觉导航方法 |
Non-Patent Citations (6)
Title |
---|
TAKEHIRO KOYAGUCHI 等: "Effects of the crater on eruption column dynamics", 《JGR SOLID EARTH》 * |
ZHENGSHI YU 等: "Design and optimization of navigation and guidance techniques for Mars pinpoint landing: review and prospect", 《PROGRESS IN AEROSPACE SCIENCE》 * |
冯军华: "月球探测器软着陆自主光学导航方法研究", 《中国博士学位论文全文数据库》 * |
邵巍 等: "一种基于星表特征直线匹配的着陆器位姿估计方法", 《深空探测学报》 * |
陈建清: "应用灰度特征的行星陨石坑自主检测方法与着陆导航研究", 《宇航学报》 * |
高锡珍: "基于陨石坑拟合椭圆的着陆器位姿估计方法研究", 《中国优秀硕士学位论文全文数据库》 * |
Cited By (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109697734A (zh) * | 2018-12-25 | 2019-04-30 | 浙江商汤科技开发有限公司 | 位姿估计方法及装置、电子设备和存储介质 |
CN109708643A (zh) * | 2019-01-14 | 2019-05-03 | 北京理工大学 | 小行星表面光学导航路标评价选取方法 |
CN109708643B (zh) * | 2019-01-14 | 2020-07-07 | 北京理工大学 | 小行星表面光学导航路标评价选取方法 |
CN110132281A (zh) * | 2019-05-21 | 2019-08-16 | 哈尔滨工程大学 | 一种基于询问应答模式的水下高速目标高精度自主声学导航方法 |
CN110132281B (zh) * | 2019-05-21 | 2023-10-20 | 哈尔滨工程大学 | 一种基于询问应答模式的水下高速目标高精度自主声学导航方法 |
CN110146093A (zh) * | 2019-06-19 | 2019-08-20 | 北京理工大学 | 双体小行星探测自主协同光学导航方法 |
CN110619368A (zh) * | 2019-09-23 | 2019-12-27 | 北京理工大学 | 行星表面导航特征成像匹配检测方法 |
CN111721302B (zh) * | 2020-06-24 | 2021-11-09 | 北京理工大学 | 不规则小行星表面复杂地形特征识别与感知方法 |
CN111721302A (zh) * | 2020-06-24 | 2020-09-29 | 北京理工大学 | 不规则小行星表面复杂地形特征识别与感知方法 |
CN111861873B (zh) * | 2020-07-20 | 2023-12-29 | 北京航天飞行控制中心 | 仿真图像的生成方法和装置 |
CN111861873A (zh) * | 2020-07-20 | 2020-10-30 | 北京航天飞行控制中心 | 仿真图像的生成方法和装置 |
CN111735460A (zh) * | 2020-08-05 | 2020-10-02 | 北京控制与电子技术研究所 | 基于小天体中心提取的航天器导航方法、系统及装置 |
CN112365541B (zh) * | 2020-11-24 | 2022-09-02 | 北京航空航天大学青岛研究院 | 一种基于相似变换的大场景相机姿态注册方法 |
CN112365541A (zh) * | 2020-11-24 | 2021-02-12 | 北京航空航天大学青岛研究院 | 一种基于相似变换的大场景相机姿态注册方法 |
CN112819794B (zh) * | 2021-02-04 | 2022-12-30 | 青岛科技大学 | 一种基于深度学习的小天体陨石坑检测方法 |
CN112819794A (zh) * | 2021-02-04 | 2021-05-18 | 青岛科技大学 | 一种基于深度学习的小天体陨石坑检测方法 |
CN112906573B (zh) * | 2021-02-10 | 2022-06-28 | 北京理工大学 | 基于轮廓点集的行星表面导航路标匹配方法 |
CN112906573A (zh) * | 2021-02-10 | 2021-06-04 | 北京理工大学 | 基于轮廓点集的行星表面导航路标匹配方法 |
CN114683279A (zh) * | 2022-02-17 | 2022-07-01 | 广东博智林机器人有限公司 | 导航精度的确定方法、装置和电子设备 |
CN114683279B (zh) * | 2022-02-17 | 2023-05-23 | 广东博智林机器人有限公司 | 导航精度的确定方法、装置和电子设备 |
CN116306185A (zh) * | 2023-02-14 | 2023-06-23 | 北京理工大学 | 基于不规则离散元的陨星进入地球大气层后运动仿真方法及系统 |
CN116306185B (zh) * | 2023-02-14 | 2023-11-03 | 北京理工大学 | 基于不规则离散元的陨星进入地球大气层后运动仿真方法及系统 |
CN116524031A (zh) * | 2023-07-03 | 2023-08-01 | 盐城数智科技有限公司 | 一种基于yolov8的大范围月球车定位建图方法 |
CN116524031B (zh) * | 2023-07-03 | 2023-09-22 | 盐城数智科技有限公司 | 一种基于yolov8的大范围月球车定位建图方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108871349B (zh) | 2021-06-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108871349A (zh) | 一种深空探测器光学导航位姿加权确定方法 | |
CN107749079B (zh) | 一种面向无人机扫描重建的点云质量评价和轨迹规划方法 | |
Pizarro et al. | Large area 3-D reconstructions from underwater optical surveys | |
CN104748750B (zh) | 一种模型约束下的在轨三维空间目标姿态估计方法及系统 | |
CN106023298B (zh) | 基于局部泊松曲面重建的点云刚性配准方法 | |
CN110866969B (zh) | 基于神经网络与点云配准的发动机叶片重构方法 | |
CN109712172A (zh) | 一种初始位姿测量结合目标追踪的位姿测量方法 | |
CN108230247B (zh) | 基于云端的三维地图的生成方法、装置、设备及计算机可读的存储介质 | |
CN106485690A (zh) | 基于点特征的点云数据与光学影像的自动配准融合方法 | |
CN105913489A (zh) | 一种利用平面特征的室内三维场景重构方法 | |
CN103988226A (zh) | 用于估计摄像机运动以及用于确定实境三维模型的方法 | |
CN111429494B (zh) | 一种基于生物视觉的点云高精度自动配准方法 | |
CN110490917A (zh) | 三维重建方法及装置 | |
CN107492107B (zh) | 基于平面与空间信息融合的物体识别与重建方法 | |
CN109708643B (zh) | 小行星表面光学导航路标评价选取方法 | |
CN110095123B (zh) | 不规则小天体表面路标观测信息评价优化方法 | |
CN113393577B (zh) | 一种倾斜摄影地形重建方法 | |
Ren et al. | A unified energy minimization framework for model fitting in depth | |
CN110243390A (zh) | 位姿的确定方法、装置及里程计 | |
Boulekchour et al. | Robust L∞ convex optimisation for monocular visual odometry trajectory estimation | |
CN112906573B (zh) | 基于轮廓点集的行星表面导航路标匹配方法 | |
CN109003295A (zh) | 一种无人机航空影像快速匹配方法 | |
Takimoto et al. | Shape reconstruction from multiple RGB-D point cloud registration | |
Baker et al. | Limb-based shape modeling and localization for autonomous navigation around small bodies | |
Oumer | Visual tracking and motion estimation for an on-orbit servicing of a satellite |
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 |