CN111104890B - 一种空间飞行器广义模型的识别和重构方法及装置 - Google Patents

一种空间飞行器广义模型的识别和重构方法及装置 Download PDF

Info

Publication number
CN111104890B
CN111104890B CN201911286093.8A CN201911286093A CN111104890B CN 111104890 B CN111104890 B CN 111104890B CN 201911286093 A CN201911286093 A CN 201911286093A CN 111104890 B CN111104890 B CN 111104890B
Authority
CN
China
Prior art keywords
ellipses
ellipse
constraint condition
circle center
target
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.)
Active
Application number
CN201911286093.8A
Other languages
English (en)
Other versions
CN111104890A (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.)
Shanghai Aerospace System Engineering Institute
Original Assignee
Shanghai Aerospace System Engineering Institute
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 Shanghai Aerospace System Engineering Institute filed Critical Shanghai Aerospace System Engineering Institute
Priority to CN201911286093.8A priority Critical patent/CN111104890B/zh
Publication of CN111104890A publication Critical patent/CN111104890A/zh
Application granted granted Critical
Publication of CN111104890B publication Critical patent/CN111104890B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00Scenes; Scene-specific elements
    • G06V20/10Terrestrial scenes
    • G06V20/13Satellite images
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00Scenes; Scene-specific elements
    • G06V20/10Terrestrial scenes

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Multimedia (AREA)
  • Theoretical Computer Science (AREA)
  • Astronomy & Astrophysics (AREA)
  • Remote Sensing (AREA)
  • Image Analysis (AREA)

Abstract

本发明实施例提供了一种空间飞行器广义模型的识别和重构方法及装置。所述方法包括:获取所述空间飞行器对应的目标图像;检测所述目标图像中的椭圆,生成椭圆集;遍历所述椭圆集,从所述椭圆集中获取满足轮廓约束条件和空间圆心约束条件的目标椭圆;基于所述目标椭圆,重构所述空间飞行器的广义模型。本发明实施例通过将空间飞行器简化为一种广义模型,该模型将空间飞行器的主体简化为满足轴对称结构的截面圆特征的集合,通过识别轴对称截面圆特征重构空间飞行器的结构。对于一幅空间飞行器图像,在完全未知目标信息的条件下,能够重构空间飞行器主体的归一化模型,既可用于后续位姿测量也可完成帧间跟踪。

Description

一种空间飞行器广义模型的识别和重构方法及装置
技术领域
本发明涉及航天技术领域,特别是一种空间飞行器广义模型的识别和重构方法及装置。
背景技术
空间飞行器的识别与重构是重要的在轨航天任务,是航天领域非常具有难度的任务之一,目前已有的方法需要已知的先验知识较多,应用范围较窄,能够识别的目标空间飞行器非常有限。
在交会对接、敌我目标识别、在轨服务等航天任务中,都需要识别、重构目标空间飞行器。通常的技术手段主要有以下几种:(1)合作目标识别法。该方法在空间飞行器上安装合作靶标,通过识别这些靶标完成空间飞行器的识别任务;(2)模板匹配法。该方法不需要在空间飞行器上安装靶标,先验地对目标的各种姿态进行模型训练并存储在硬件***中,通过匹配已知模板和得到的图像完成识别任务;(3)目标特征识别法。该方法识别目标固有的特征,利用这些特征之间可能存在的内在联系完成识别任务。综上,第一种方法需要目标上具有合作靶标,第二种方法只能识别训练过的目标并且消耗的硬件***的存储资源较大,这两种方法的应用限制条件较多。第三种方法具更广的适用范围,通过改进识别算法可以达到非常好的性能。
综上,现有的目标特征识别法具有较大限制性,大部分着眼于局部特征的检测,因此在先验知识有限的条件下无法对目标空间飞行器有着整体的认识,也无法将这些局部特征和干扰特征区分。
发明内容
本发明解决的技术问题是:克服现有技术中椭圆特征检测着眼于局部特征的检测,因此在先验知识有限的条件下无法对目标空间飞行器有着整体的认识,也无法将这些局部特征和干扰特征区分的不足,提供了一种空间飞行器广义模型的识别和重构方法及装置。
本发明的技术解决方案是:
第一方面,本发明实施例提供了一种空间飞行器广义模型的识别和重构方法,包括:
获取所述空间飞行器对应的目标图像;
检测所述目标图像中的椭圆,生成椭圆集;
遍历所述椭圆集,从所述椭圆集中获取满足轮廓约束条件和空间圆心约束条件的目标椭圆;
基于所述目标椭圆,重构所述空间飞行器的广义模型。
可选地,所述遍历所述椭圆集,从所述椭圆集中获取满足轮廓约束条件和空间圆心约束条件的目标椭圆,包括:
从所述椭圆集中随机获取任意两个椭圆;
计算所述任意两个椭圆对应的空间圆所在平面的法向量;
根据所述法向量,确定所述任意两个椭圆是否满足所述轮廓约束条件;
在所述任意两个椭圆满足所述轮廓约束条件的情况下,确定所述任意两个椭圆是否满足所述空间圆心约束条件;
在所述任意两个椭圆满足所述空间圆心约束条件的情况下,将所述任意两个椭圆作为所述目标椭圆;
重复执行所述从所述椭圆集中随机获取任意两个椭圆,至所述在所述任意两个椭圆满足所述空间圆心约束条件的情况下,将所述任意两个椭圆作为所述目标椭圆的步骤;
在不同的目标椭圆对中含有相同的椭圆时,合并所述不同的目标椭圆对。
可选地,所述计算所述任意两个椭圆对应的空间圆所在平面的法向量,包括:
设定所述任意两个椭圆分别对应的矩阵形式;
根据两个所述矩阵形式,确定所述任意两个椭圆形成的椭圆束的椭圆束矩阵形式;
根据所述椭圆束矩阵形式,确定无穷远直线的投影;
获取所述无穷远直线的投影对应的投影坐标;
基于所述投影坐标,计算得到所述法向量。
可选地,所述根据所述法向量,确定所述任意两个椭圆是否满足所述轮廓约束条件,包括:
计算以所述目标图像对应的获取设备为圆心,所述任意两个椭圆映射至空间中形成的锥面方程;
利用坐标系变换原则,将所述锥面方程转换为标准方程;
根据所述标准方程,计算截面轮廓方程;
基于所述截面轮廓方程,确定所述任意两个椭圆是否满足所述轮廓约束条件。
可选地,所述在所述任意两个椭圆满足所述轮廓约束条件的情况下,确定所述任意两个椭圆是否满足所述空间圆心约束条件的步骤,包括:
计算所述空间圆圆心投影到图像上的投影点的投影齐次坐标;
计算所述空间圆圆心的齐次坐标;
基于所述投影齐次坐标和所述齐次坐标,确定所述任意两个椭圆是否满足所述圆心约束条件。
可选地,所述基于所述目标椭圆,重构所述空间飞行器的广义模型,包括:
计算所述圆心到所述任意两个椭圆所形成的空间截面之间的第一相对距离;
计算所述任意两个椭圆所形成的空间截面圆的相对半径;
计算所述任意两个椭圆所形成的空间截面之间的第二相对距离;
基于所述第一相对距离、所述相对半径和所述第二相对距离,重构所述广义模型。
第二方面,本发明实施例提供了一种空间飞行器广义模型的识别和重构装置,包括:
目标图像获取模块,用于获取所述空间飞行器对应的目标图像;
椭圆集生成模块,用于检测所述目标图像中的椭圆,生成椭圆集;
目标椭圆获取模块,用于遍历所述椭圆集,从所述椭圆集中获取满足轮廓约束条件和空间圆心约束条件的目标椭圆;
广义模型重构模块,用于基于所述目标椭圆,重构所述空间飞行器的广义模型。
可选地,所述目标椭圆获取模块包括:
任意椭圆获取子模块,用于从所述椭圆集中随机获取任意两个椭圆;
法向量计算子模块,用于计算所述任意两个椭圆对应的空间圆所在平面的法向量;
轮廓约束确定子模块,用于根据所述法向量,确定所述任意两个椭圆是否满足所述轮廓约束条件;
圆心约束确定子模块,用于在所述任意两个椭圆满足所述轮廓约束条件的情况下,确定所述任意两个椭圆是否满足所述空间圆心约束条件;
目标椭圆获取子模块,用于在所述任意两个椭圆满足所述空间圆心约束条件的情况下,将所述任意两个椭圆作为所述目标椭圆;
重复执行子模块,用于重复执行所述任意椭圆获取子模块,至所述目标椭圆获取子模块;
椭圆对合并子模块,用于在不同的目标椭圆对中含有相同的椭圆时,合并所述不同的目标椭圆对。
可选地,所述法向量计算子模块包括:
矩阵形式设定子模块,用于设定所述任意两个椭圆分别对应的矩阵形式;
椭圆束矩阵形式确定子模块,用于根据两个所述矩阵形式,确定所述任意两个椭圆形成的椭圆束的椭圆束矩阵形式;
直线投影确定子模块,用于根据所述椭圆束矩阵形式,确定无穷远直线的投影;
投影坐标获取子模块,用于获取所述无穷远直线的投影对应的投影坐标;
法向量获取子模块,用于基于所述投影坐标,计算得到所述法向量。
可选地,所述轮廓约束确定子模块包括:
锥面方程计算子模块,用于计算以所述目标图像对应的获取设备为圆心,所述任意两个椭圆映射至空间中形成的锥面方程;
标准方程转换子模块,用于利用坐标系变换原则,将所述锥面方程转换为标准方程;
轮廓方程计算子模块,用于根据所述标准方程,计算截面轮廓方程;
轮廓约束条件确定子模块,用于基于所述截面轮廓方程,确定所述任意两个椭圆是否满足所述轮廓约束条件。
可选地,所述圆心约束确定子模块包括:
圆心投影坐标计算子模块,用于计算所述空间圆圆心投影到图像上的投影点的投影齐次坐标;
圆心坐标计算子模块,用于计算所述空间圆圆心的齐次坐标;
圆心约束条件确定子模块,用于基于所述投影齐次坐标和所述齐次坐标,确定所述任意两个椭圆是否满足所述圆心约束条件。
可选地,所述广义模型重构模块包括:
第一相对距离计算子模块,用于计算所述圆心到所述任意两个椭圆所形成的空间截面之间的第一相对距离;
相对半径计算子模块,用于计算所述任意两个椭圆所形成的空间截面圆的相对半径;
第二相对距离计算子模块,用于计算所述任意两个椭圆所形成的空间截面之间的第二相对距离;
广义模型重构子模块,用于基于所述第一相对距离、所述相对半径和所述第二相对距离,重构所述广义模型。
本发明与现有技术相比的优点在于:
本发明实施例提供的椭圆特征检测方案,通过获取空间飞行器对应的目标图像,检测目标图像中的椭圆,生成椭圆集,遍历椭圆集,从椭圆集中获取满足轮廓约束条件和空间圆心约束条件的目标椭圆,基于目标椭圆重构空间飞行器的广义模型。本发明分析了现有空间飞行器的整体结构特征,将空间飞行器简化为一种广义模型,该模型将空间飞行器的主体简化为满足轴对称结构的截面圆特征的集合,通过识别轴对称截面圆特征重构空间飞行器的结构。对于一幅空间飞行器图像,在完全未知目标信息的条件下,能够重构空间飞行器主体的归一化模型,既可用于后续位姿测量也可完成帧间跟踪。
附图说明
图1为本发明实施例提供的一种空间飞行器广义模型的识别和重构方法的步骤流程图;
图2为本发明实施例提供的一种轮廓约束的示意图;
图3为本发明实施例提供的一种空间圆心约束的示意图;
图4为本发明实施例提供的一种重建空间飞行器的广义模型的示意图;
图5为本发明实施例提供的一种空间飞行器广义模型的识别和重构装置的结构示意图。
具体实施方式
下面将结合本发明的实施例中的附图,对本发明的实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明的实施例一部分实施例,而不是全部的实施例。基于本发明的实施例中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明的实施例保护的范围。
实施例一
参照图1,示出了本发明实施例提供的一种空间飞行器广义模型的识别和重构方法的步骤流程图,如图1所示,该空间飞行器广义模型的识别和重构方法具体可以包括如下步骤:
步骤101:获取所述空间飞行器对应的目标图像。
本发明实施例可以应用于对空间飞行器上的椭圆特征进行检测的场景中。
空间飞行器是指在地球大气层以外的宇宙空间,基本上按照天体力学的规律运行的各类飞行器。
本发明中,空间飞行器可以为人造地球卫星、空间探测器、载人航天器等飞行器中的任一种,具体地,可以根据实际情况而定,本发明实施例对此不加以限制。
目标图像是指空间飞行器的表面图像。本发明中可以采用相机采集空间飞行器的目标图像,如采用飞行器航拍设备拍摄得到空间飞行器的目标图像等。
在获取空间飞行器对应的目标图像之后,执行步骤102。
步骤102:检测所述目标图像中的椭圆,生成椭圆集。
椭圆集是指根据空间飞行器上的椭圆形成的集合。
在获取空间飞行器对应的目标图像之后,可以检测目标图像中的椭圆图像特征,并结合椭圆图像特征可以生成多个空间飞行器的表面上的椭圆,结合这些椭圆可以得到椭圆集。
在生成椭圆集之后,执行步骤103。
步骤103:遍历所述椭圆集,从所述椭圆集中获取满足轮廓约束条件和空间圆心约束条件的目标椭圆。
轮廓约束条件是指椭圆集中任意两个椭圆的轮廓所需满足的约束条件。
空间圆心约束条件是指椭圆集中任意两个椭圆在空间所形成的椭圆平面的圆心所需满足的约束条件。
目标椭圆是指椭圆集中满足轮廓约束条件和空间圆心约束条件的椭圆。
在生成椭圆集之后,可以遍历椭圆集,从椭圆集中获取满足轮廓约束条件和空间圆心约束条件的目标椭圆,具体地,可以结合下述优选实施例进行详细描述。
在本发明的一种优选实施例中,上述步骤103可以包括:
子步骤A1:从所述椭圆集中随机获取任意两个椭圆;
子步骤A2:计算所述任意两个椭圆对应的空间圆所在平面的法向量。
在本发明实施例中,首先可以从椭圆集中随机获取任意两个椭圆,并计算任意两个椭圆对应的空间圆所在平面的法向量,具体地:
若两个椭圆由空间飞行器广义模型截面圆投影形成,则它们形成的椭圆束必然退化为截面无穷远直线的投影因此法向量计算可以如下述优选实施例的描述。
在本发明的另一种优选实施例中,上述子步骤A2可以包括:
子步骤B1:设定所述任意两个椭圆分别对应的矩阵形式;
子步骤B2:根据两个所述矩阵形式,确定所述任意两个椭圆形成的椭圆束的椭圆束矩阵形式;
子步骤B3:根据所述椭圆束矩阵形式,确定无穷远直线的投影;
子步骤B4:获取所述无穷远直线的投影对应的投影坐标;
子步骤B5:基于所述投影坐标,计算得到所述法向量。
设两个椭圆参数的矩阵形式为和/>由这两个椭圆形成的椭圆束的矩阵形式可依据公式(1)计算
上述公式(1)中,D为任意两个椭圆形成的椭圆束的矩阵。
(2)对矩阵D进行特征值分解得到椭圆束的三种退化形式Di(按照小数公式(2)计算),三组解中包含了待求的无穷远直线的投影:
(3)计算退化二次曲线D的相对符号(公式(3)),并依据相对符号和退化直线的关系(公式(4))得到含有无穷远直线的投影的解:
∑(D)=|η-ξ| (3)
上述公式(3)和(4)中,∑(D)为二次曲线D的相对符号。
(4)求无穷远直线的投影坐标不妨令z3是解/>对应的特征向量,则/>可由以下原则从退化直线的参数矢量{L,m}中选出:特征向量z1和z2表示的点分布在/>的一侧,分布在/>的两侧,即:
(5)求法向量:
上述公式中,K为摄像机标定矩阵。
在计算任意两个椭圆对应的空间圆所在平面的法向量之后,执行子步骤A3。
子步骤A3:根据所述法向量,确定所述任意两个椭圆是否满足所述轮廓约束条件。
在计算任意两个椭圆对应的空间圆所在平面的法向量之后,可以根据任意两个椭圆的法向量确定任意两个椭圆是否满足轮廓约束条件,具体地,结合下述优选实施例进行如下描述。
在本发明的一种优选实施例中,上述子步骤A3可以包括:
子步骤C1:计算以所述目标图像对应的获取设备为圆心,所述任意两个椭圆映射至空间中形成的锥面方程;
子步骤C2:利用坐标系变换原则,将所述锥面方程转换为标准方程;
子步骤C3:根据所述标准方程,计算截面轮廓方程;
子步骤C4:基于所述截面轮廓方程,确定所述任意两个椭圆是否满足所述轮廓约束条件。
在本发明实施例中,如图2所示,以目标图像对应的获取设备(如相机等)的中心为顶点,将图像中的椭圆映射到空间中形成圆锥,和法向量垂直的平面截取该圆锥形成的轮廓必然为圆。具体算法流程如下:
(1)计算以相机中心为顶点,将图像中的椭圆映射到空间中形成锥面方程如下述公式(8):
上述公式(8)中,为空间锥面的矩阵形式,可由下述公式(9)计算得到:
(2)利用坐标系变换,将空间锥面方程转化为标准形式方程,使得新坐标系z'轴垂直于空间中截面圆所在的平面。坐标系变换矩阵如公式(10)所示,在三维摄影空间IP3中由摄像机坐标系O-xyz变换到标准坐标系O-x'y'z'的射影变换见公式(11),在标准坐标系O-x'y'z'下空间锥面的齐次坐标可由公式(12)计算。
(3)计算截面轮廓方程。如图2锥面和任意平面z'=k(k>0)相交形成的二次曲线方程为:
(4)该二次曲线必然为圆,给出轮廓约束,优选地,阈值Th1和Th2取1%。
(5)若满足轮廓约束,则确定任意两个椭圆满足轮廓约束条件;否则返回,执行重新从椭圆集中随机获取任意两个椭圆的步骤。
在确定任意两个椭圆满足轮廓约束条件之后,执行子步骤A4。
子步骤A4:在所述任意两个椭圆满足所述轮廓约束条件的情况下,确定所述任意两个椭圆是否满足所述空间圆心约束条件。
在确定任意两个椭圆满足轮廓约束条件之后,则可以判定任意两个椭圆是否满足空间圆心约束条件,具体地,可以结合下述优选实施例进行详细描述。
在本发明的一种优选实施例中,上述子步骤A5可以包括:
子步骤D1:计算所述空间圆圆心投影到图像上的投影点的投影齐次坐标;
子步骤D2:计算所述空间圆圆心的齐次坐标;
子步骤D3:基于所述投影齐次坐标和所述齐次坐标,确定所述任意两个椭圆是否满足所述圆心约束条件。
在本发明实施例中,投影齐次坐标是指空间圆圆心投影到图像上的投影点的齐次坐标。
计算任意两个椭圆是否满足空间圆心约束条件的依据如下:如图3所示,在标准坐标系O-x'y'z'下,空间圆圆心坐标相等。
具体算法流程如下:
(1)计算空间圆圆心在图像上投影的齐次坐标。设空间中截面圆C1、C2的圆心为p1、p2,对应在摄像机像平面上的投影点为点/>在摄像机像面坐标系o-uv下的齐次坐标可直接由对极关系得到:
上述公式(15)和(16)中,su1、su2为齐次坐标的比例因子。
这两点在标准坐标系O-x'y'z'下的坐标可表示为:
上述公式(17)和(18)中,分别表示任意两个椭圆的空间圆心在标准坐标系中的坐标。
(2)计算空间圆圆心的相对坐标。如图3,在三角形OP1P2中,向量和/>的夹角α12可由式(19)计算得到:
向量和平面法向量以及/>和平面法向量构成的角度可能为:
三角形op1p2内的角度{α121212}可由式(22)确定:
β12∈{β′,β″},γ12∈{γ′,γ″};subject to α121212=π (22)
定义:线段op1的长度d1为标准长度。线段op2的长度可由正弦定理得出,即线段op2的相对长度为:
因此,在坐标系O-x'y'z'下,空间中圆C1、C2圆心p1、p2的坐标可表示为:
(3)圆心坐标必然相等,给出约束,优选地,阈值Th3和Th4取1%。
子步骤A5:在所述任意两个椭圆满足所述空间圆心约束条件的情况下,将所述任意两个椭圆作为所述目标椭圆。
在任意两个椭圆同时满足轮廓约束条件和空间圆心约束条件时,则将任意两个椭圆确定为目标椭圆。
子步骤A6:重复执行所述子步骤A1,至所述子步骤A5。
然后,重复执行子步骤A1至子步骤A5,即从椭圆集中再随机获取任意两个椭圆,执行轮廓约束条件和空间圆心约束条件的判定步骤,并确定后获取的任意两个椭圆是否为目标椭圆;
子步骤A7:在不同的目标椭圆对中含有相同的椭圆时,合并所述不同的目标椭圆对。
在检测得到的多个目标椭圆对中可能存在不同的目标椭圆对中含有相同的椭圆,此时,可以合并含有相同的椭圆的目标椭圆对,从而能够得到形成广义模型的所有椭圆。
在确定出椭圆集中的目标椭圆之后,执行步骤104。
步骤104:基于所述目标椭圆,重构所述空间飞行器的广义模型。
在得到椭圆集中的目标椭圆之后,可以基于目标椭圆重构空间飞行器的广义模型,具体地,可以结合下述优选实施例进行详细描述。
在本发明的一种优选实施例中,上述步骤104可以包括:
子步骤E1:计算所述圆心到所述任意两个椭圆所形成的空间截面之间的第一相对距离;
子步骤E2:计算所述任意两个椭圆所形成的空间截面圆的相对半径;
子步骤E3:计算所述任意两个椭圆所形成的空间截面之间的第二相对距离;
子步骤E4:基于所述第一相对距离、所述相对半径和所述第二相对距离,重构所述广义模型。
在本发明实施例中,在获取目标椭圆之后,可以结合目标椭圆对空间飞行器的广义模型进行重构,得到空间飞行器的归一化模型。具体算法流程如下:
(1)计算原点到空间截面的相对距离k1和k2(即第一相对距离):
(2)计算空间截面圆的相对半径R1和R2
(3)计算空间截面间的第二相对距离
重构空间飞行器的广义模型。约去各式中的非零标准值:坐标系原点O到C1圆心的距离d1,可计算空间轴对称结构中圆特征和矩形特征的归一化半径和距离(公式(32)和公式(33)),沿法线n方向对所有特征重新编号,可实现以空间中n个截面圆重建广义模型,如图4所示。
RR=R1:R2:...:Rn (32)
RL=L12:L13:...:L1n (33)
本发明不需要已知目标飞行器上安装有合作标志、也不需要提前对飞行器训练建模,仅要求空间飞行器具有圆截面的轴对称结构即可完成识别和重构问题,这种截面是圆的轴对称结构是绝大部分空间飞行器的主体,因而本发明可以识别绝大部分空间飞行器。
本发明创建的广义模型是一种定量的计算模型,因而非常适合进行拓展,如:在本发明的基础上建立目标坐标系,可完成后续的位姿测量,是相对导航、目标动态跟踪等任务的基础。
本发明实施例提供的空间飞行器广义模型的识别和重构方法,通过获取空间飞行器对应的目标图像,检测目标图像中的椭圆,生成椭圆集,遍历椭圆集,从椭圆集中获取满足轮廓约束条件和空间圆心约束条件的目标椭圆,基于目标椭圆重构空间飞行器的广义模型。本发明分析了现有空间飞行器的整体结构特征,将空间飞行器简化为一种广义模型,该模型将空间飞行器的主体简化为满足轴对称结构的截面圆特征的集合,通过识别轴对称截面圆特征重构空间飞行器的结构。对于一幅空间飞行器图像,在完全未知目标信息的条件下,能够重构空间飞行器主体的归一化模型,既可用于后续位姿测量也可完成帧间跟踪。
实施例二
参照图5,示出了本发明实施例提供的一种空间飞行器广义模型的识别和重构装置的结构示意图,如图5所示,该空间飞行器广义模型的识别和重构装置具体可以包括如下模块:
目标图像获取模块510,用于获取所述空间飞行器对应的目标图像;
椭圆集生成模块520,用于检测所述目标图像中的椭圆,生成椭圆集;
目标椭圆获取模块530,用于遍历所述椭圆集,从所述椭圆集中获取满足轮廓约束条件和空间圆心约束条件的目标椭圆;
广义模型重构模块540,用于基于所述目标椭圆,重构所述空间飞行器的广义模型。
可选地,所述目标椭圆获取模块530包括:
任意椭圆获取子模块,用于从所述椭圆集中随机获取任意两个椭圆;
法向量计算子模块,用于计算所述任意两个椭圆对应的空间圆所在平面的法向量;
轮廓约束确定子模块,用于根据所述法向量,确定所述任意两个椭圆是否满足所述轮廓约束条件;
圆心约束确定子模块,用于在所述任意两个椭圆满足所述轮廓约束条件的情况下,确定所述任意两个椭圆是否满足所述空间圆心约束条件;
目标椭圆获取子模块,用于在所述任意两个椭圆满足所述空间圆心约束条件的情况下,将所述任意两个椭圆作为所述目标椭圆;
重复执行子模块,用于重复执行所述任意椭圆获取子模块,至所述目标椭圆获取子模块。
可选地,所述法向量计算子模块包括:
矩阵形式设定子模块,用于设定所述任意两个椭圆分别对应的矩阵形式;
椭圆束矩阵形式确定子模块,用于根据两个所述矩阵形式,确定所述任意两个椭圆形成的椭圆束的椭圆束矩阵形式;
直线投影确定子模块,用于根据所述椭圆束矩阵形式,确定无穷远直线的投影;
投影坐标获取子模块,用于获取所述无穷远直线的投影对应的投影坐标;
法向量获取子模块,用于基于所述投影坐标,计算得到所述法向量。
可选地,所述轮廓约束确定子模块包括:
锥面方程计算子模块,用于计算以所述目标图像对应的获取设备为圆心,所述任意两个椭圆映射至空间中形成的锥面方程;
标准方程转换子模块,用于利用坐标系变换原则,将所述锥面方程转换为标准方程;
轮廓方程计算子模块,用于根据所述标准方程,计算截面轮廓方程;
轮廓约束条件确定子模块,用于基于所述截面轮廓方程,确定所述任意两个椭圆是否满足所述轮廓约束条件。
可选地,所述圆心约束确定子模块包括:
圆心投影坐标计算子模块,用于计算所述空间圆圆心投影到图像上的投影点的投影齐次坐标;
圆心坐标计算子模块,用于计算所述空间圆圆心的齐次坐标;
圆心约束条件确定子模块,用于基于所述投影齐次坐标和所述齐次坐标,确定所述任意两个椭圆是否满足所述圆心约束条件。
可选地,所述广义模型重构模块540包括:
第一相对距离计算子模块,用于计算所述圆心到所述任意两个椭圆所形成的空间截面之间的第一相对距离;
相对半径计算子模块,用于计算所述任意两个椭圆所形成的空间截面圆的相对半径;
第二相对距离计算子模块,用于计算所述任意两个椭圆所形成的空间截面之间的第二相对距离;
广义模型重构子模块,用于基于所述第一相对距离、所述相对半径和所述第二相对距离,重构所述广义模型。
本发明实施例提供的空间飞行器广义模型的识别和重构装置,通过获取空间飞行器对应的目标图像,检测目标图像中的椭圆,生成椭圆集,遍历椭圆集,从椭圆集中获取满足轮廓约束条件和空间圆心约束条件的目标椭圆,基于目标椭圆重构空间飞行器的广义模型。本发明分析了现有空间飞行器的整体结构特征,将空间飞行器简化为一种广义模型,该模型将空间飞行器的主体简化为满足轴对称结构的截面圆特征的集合,通过识别轴对称截面圆特征重构空间飞行器的结构。对于一幅空间飞行器图像,在完全未知目标信息的条件下,能够重构空间飞行器主体的归一化模型,既可用于后续位姿测量也可完成帧间跟踪。
以上所述仅为本发明的实施例的较佳实施例而已,并不用以限制本发明的实施例,凡在本发明的实施例的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的实施例的保护范围之内。

Claims (4)

1.一种空间飞行器广义模型的识别和重构方法,其特征在于,包括:
获取所述空间飞行器对应的目标图像;
检测所述目标图像中的椭圆,生成椭圆集;
遍历所述椭圆集,从所述椭圆集中获取满足轮廓约束条件和空间圆心约束条件的目标椭圆;
基于所述目标椭圆,重构所述空间飞行器的广义模型;
所述遍历所述椭圆集,从所述椭圆集中获取满足轮廓约束条件和空间圆心约束条件的目标椭圆,包括:
从所述椭圆集中随机获取任意两个椭圆;
计算所述任意两个椭圆对应的空间圆所在平面的法向量;
根据所述法向量,确定所述任意两个椭圆是否满足所述轮廓约束条件,包括:计算以所述目标图像对应的获取设备为圆心,所述任意两个椭圆映射至空间中形成的锥面方程;利用坐标系变换原则,将所述锥面方程转换为标准方程;根据所述标准方程,计算截面轮廓方程;基于所述截面轮廓方程,确定所述任意两个椭圆的截面是否为圆,如果是圆则满足所述轮廓约束条件,如果不是圆则不满足所述轮廓约束条件;
在所述任意两个椭圆满足所述轮廓约束条件的情况下,确定所述任意两个椭圆是否满足所述空间圆心约束条件,包括:计算所述空间圆圆心投影到图像上的投影点的投影齐次坐标;计算所述空间圆圆心在标准坐标系O-x'y'z'下的相对坐标;确定所述任意两个椭圆的圆心相对坐标是否满足所述任意两个椭圆的圆心在x’轴的坐标值分别为/>和/>所述任意两个椭圆的圆心在y’轴的坐标值分别为/>和/>Th3和Th4为阈值,如果是则满足所述圆心约束条件,如果不是则不满足所述圆心约束条件;
在所述任意两个椭圆满足所述空间圆心约束条件的情况下,将所述任意两个椭圆作为所述目标椭圆;
重复执行所述从所述椭圆集中随机获取任意两个椭圆,至所述在所述任意两个椭圆满足所述空间圆心约束条件的情况下,将所述任意两个椭圆作为所述目标椭圆的步骤;
在不同的目标椭圆对中含有相同的椭圆时,合并所述不同的目标椭圆对;
所述基于所述目标椭圆,重构所述空间飞行器的广义模型,包括:
计算所述圆心到所述任意两个椭圆所形成的空间截面之间的第一相对距离;
计算所述任意两个椭圆所形成的空间截面圆的相对半径;
计算所述任意两个椭圆所形成的空间截面之间的第二相对距离;
基于所述第一相对距离、所述相对半径和所述第二相对距离,重构所述广义模型,包括:
根据空间轴对称结构中相对半径和第二相对距离,沿法线n方向对所有特征重新编号,得到广义模型:RR=R1:R2:...:Rn,RL=L12:L13:...:L1n,下标n为沿法线n方向所有特征的数量,RR=R1:R2:...:Rn表示各个特征的归一化半径,RL=L12:L13:...:L1n表示各个特征的距离。
2.根据权利要求1所述的方法,其特征在于,所述计算所述任意两个椭圆对应的空间圆所在平面的法向量,包括:
设定所述任意两个椭圆分别对应的矩阵形式;
根据两个所述矩阵形式,确定所述任意两个椭圆形成的椭圆束的椭圆束矩阵形式;
根据所述椭圆束矩阵形式,确定无穷远直线的投影;
获取所述无穷远直线的投影对应的投影坐标;
基于所述投影坐标,计算得到所述法向量。
3.一种空间飞行器广义模型的识别和重构装置,其特征在于,包括:
目标图像获取模块,用于获取所述空间飞行器对应的目标图像;
椭圆集生成模块,用于检测所述目标图像中的椭圆,生成椭圆集;
目标椭圆获取模块,用于遍历所述椭圆集,从所述椭圆集中获取满足轮廓约束条件和空间圆心约束条件的目标椭圆;
广义模型重构模块,用于基于所述目标椭圆,重构所述空间飞行器的广义模型;所述目标椭圆获取模块包括:
任意椭圆获取子模块,用于从所述椭圆集中随机获取任意两个椭圆;
法向量计算子模块,用于计算所述任意两个椭圆对应的空间圆所在平面的法向量;
轮廓约束确定子模块,用于根据所述法向量,确定所述任意两个椭圆是否满足所述轮廓约束条件;
圆心约束确定子模块,用于在所述任意两个椭圆满足所述轮廓约束条件的情况下,确定所述任意两个椭圆是否满足所述空间圆心约束条件;
目标椭圆获取子模块,用于在所述任意两个椭圆满足所述空间圆心约束条件的情况下,将所述任意两个椭圆作为所述目标椭圆;
重复执行子模块,用于重复执行所述任意椭圆获取子模块,至所述目标椭圆获取子模块;
椭圆对合并子模块,用于在不同的目标椭圆对中含有相同的椭圆时,合并所述不同的目标椭圆对;
所述轮廓约束确定子模块包括:
锥面方程计算子模块,用于计算以所述目标图像对应的获取设备为圆心,所述任意两个椭圆映射至空间中形成的锥面方程;
标准方程转换子模块,用于利用坐标系变换原则,将所述锥面方程转换为标准方程;
轮廓方程计算子模块,用于根据所述标准方程,计算截面轮廓方程;
轮廓约束条件确定子模块,用于基于所述截面轮廓方程,确定所述任意两个椭圆的截面是否为圆,如果是圆则满足所述轮廓约束条件,如果不是圆则不满足所述轮廓约束条件;
所述圆心约束确定子模块包括:
圆心投影坐标计算子模块,用于计算所述空间圆圆心投影到图像上的投影点的投影齐次坐标;
圆心坐标计算子模块,用于计算所述空间圆圆心在标准坐标系O-x'y'z'下的相对坐标;
圆心约束条件确定子模块,用于确定所述任意两个椭圆的圆心相对坐标是否满足所述任意两个椭圆的圆心在x’轴的坐标值分别为/>和/>所述任意两个椭圆的圆心在y’轴的坐标值分别为/>和/>Th3和Th4为阈值,如果是则满足所述圆心约束条件,如果不是则不满足所述圆心约束条件;
所述广义模型重构模块包括:
第一相对距离计算子模块,用于计算所述圆心到所述任意两个椭圆所形成的空间截面之间的第一相对距离;
相对半径计算子模块,用于计算所述任意两个椭圆所形成的空间截面圆的相对半径;
第二相对距离计算子模块,用于计算所述任意两个椭圆所形成的空间截面之间的第二相对距离;
广义模型重构子模块,用于基于所述第一相对距离、所述相对半径和所述第二相对距离,重构所述广义模型,包括:根据空间轴对称结构中相对半径和第二相对距离,沿法线n方向对所有特征重新编号,得到广义模型:RR=R1:R2:...:Rn,RL=L12:L13:...:L1n,下标n为沿法线n方向所有特征的数量,RR=R1:R2:...:Rn表示各个特征的归一化半径,RL=L12:L13:...:L1n表示各个特征的距离。
4.根据权利要求3所述的装置,其特征在于,所述法向量计算子模块包括:
矩阵形式设定子模块,用于设定所述任意两个椭圆分别对应的矩阵形式;
椭圆束矩阵形式确定子模块,用于根据两个所述矩阵形式,确定所述任意两个椭圆形成的椭圆束的椭圆束矩阵形式;
直线投影确定子模块,用于根据所述椭圆束矩阵形式,确定无穷远直线的投影;
投影坐标获取子模块,用于获取所述无穷远直线的投影对应的投影坐标;
法向量获取子模块,用于基于所述投影坐标,计算得到所述法向量。
CN201911286093.8A 2019-12-13 2019-12-13 一种空间飞行器广义模型的识别和重构方法及装置 Active CN111104890B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911286093.8A CN111104890B (zh) 2019-12-13 2019-12-13 一种空间飞行器广义模型的识别和重构方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911286093.8A CN111104890B (zh) 2019-12-13 2019-12-13 一种空间飞行器广义模型的识别和重构方法及装置

Publications (2)

Publication Number Publication Date
CN111104890A CN111104890A (zh) 2020-05-05
CN111104890B true CN111104890B (zh) 2023-09-29

Family

ID=70423014

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911286093.8A Active CN111104890B (zh) 2019-12-13 2019-12-13 一种空间飞行器广义模型的识别和重构方法及装置

Country Status (1)

Country Link
CN (1) CN111104890B (zh)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4692878A (en) * 1985-03-29 1987-09-08 Ampower Technologies, Inc. Three-dimensional spatial image system
JPH06231249A (ja) * 1993-02-05 1994-08-19 Fujitsu Ltd 物体認識方式
CA2241383A1 (en) * 1997-06-20 1998-12-20 Ronen Ezra Basri A system for visual homing and navigation
RU69674U1 (ru) * 2007-09-19 2007-12-27 Мурадин Абубекирович Кумахов Устройство для формирования направленного пучка рентгеновского излучения
CN104517291A (zh) * 2014-12-15 2015-04-15 大连理工大学 基于目标同轴圆特征的位姿测量方法
CN108492333A (zh) * 2018-03-30 2018-09-04 哈尔滨工业大学 基于星箭对接环图像信息的航天器姿态估计方法
CN109102567A (zh) * 2018-10-11 2018-12-28 北京理工大学 一种基于重建误差最小化的位姿参数高精度求解方法
CN109405835A (zh) * 2017-08-31 2019-03-01 北京航空航天大学 基于非合作目标直线与圆单目图像的相对位姿测量方法
CN109471108A (zh) * 2018-12-17 2019-03-15 中国人民解放军空军工程大学 进动椎体目标三维成像方法、装置及电子设备

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8619074B2 (en) * 2010-12-10 2013-12-31 Xerox Corporation Rendering personalized text on curved image surfaces

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4692878A (en) * 1985-03-29 1987-09-08 Ampower Technologies, Inc. Three-dimensional spatial image system
JPH06231249A (ja) * 1993-02-05 1994-08-19 Fujitsu Ltd 物体認識方式
CA2241383A1 (en) * 1997-06-20 1998-12-20 Ronen Ezra Basri A system for visual homing and navigation
RU69674U1 (ru) * 2007-09-19 2007-12-27 Мурадин Абубекирович Кумахов Устройство для формирования направленного пучка рентгеновского излучения
CN104517291A (zh) * 2014-12-15 2015-04-15 大连理工大学 基于目标同轴圆特征的位姿测量方法
CN109405835A (zh) * 2017-08-31 2019-03-01 北京航空航天大学 基于非合作目标直线与圆单目图像的相对位姿测量方法
CN108492333A (zh) * 2018-03-30 2018-09-04 哈尔滨工业大学 基于星箭对接环图像信息的航天器姿态估计方法
CN109102567A (zh) * 2018-10-11 2018-12-28 北京理工大学 一种基于重建误差最小化的位姿参数高精度求解方法
CN109471108A (zh) * 2018-12-17 2019-03-15 中国人民解放军空军工程大学 进动椎体目标三维成像方法、装置及电子设备

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Multiple circle recognition and pose estimation for aerospace application;Wu, B;《OPTIK》;第145卷;第148-157页 *
POSITION AND ORIENTATION MEASUREMENT FOR AUTONOMOUS AERIAL REFUELING BASED ON MONOCULAR VISION;Wu, WQ等;《INTERNATIONAL JOURNAL OF ROBOTICS & AUTOMATION》;第32卷(第1期);第13-21页 *
吴斌等.基于航天应用的多圆特征识别和姿态估计.《光学学报》.2017,第37卷(第09期),第0915004-1~0915004-11页. *
郑顺义等.基于模型和广义点摄影测量的圆柱体自动三维重建与检测.《测绘学报》.2011,第40卷(第04期),第477-482页. *

Also Published As

Publication number Publication date
CN111104890A (zh) 2020-05-05

Similar Documents

Publication Publication Date Title
Sharma et al. Pose estimation for non-cooperative rendezvous using neural networks
Sharma et al. Neural network-based pose estimation for noncooperative spacecraft rendezvous
CN111080627B (zh) 一种基于深度学习的2d+3d大飞机外形缺陷检测与分析方法
CN108562274B (zh) 一种基于标志器的非合作目标位姿测量方法
CN107063228B (zh) 基于双目视觉的目标姿态解算方法
RU2609434C2 (ru) Выявление расположения и местоположения объектов
Zhang et al. Vision-based pose estimation for textureless space objects by contour points matching
CN108596980B (zh) 圆形目标视觉定位精度评估方法、装置、存储介质和处理设备
CN110866969B (zh) 基于神经网络与点云配准的发动机叶片重构方法
JP2016161569A (ja) オブジェクトの3d姿勢およびオブジェクトのランドマーク点の3dロケーションを求める方法、およびオブジェクトの3d姿勢およびオブジェクトのランドマークの3dロケーションを求めるシステム
CN103649674A (zh) 测量设备以及信息处理设备
CN112347550A (zh) 耦合式室内三维语义建图及建模方法
Sell et al. Pose performance of LIDAR-based navigation for satellite servicing
CN107481281B (zh) 相对位姿计算方法和装置及航天飞行器交会对接***
He et al. Pose tracking of spacecraft based on point cloud DCA features
Takeishi et al. Evaluation of interest-region detectors and descriptors for automatic landmark tracking on asteroids
Vaillant Using occluding contours for 3D object modeling
Hewitt et al. LiDAR-Based map relative localization performance analysis for landing on Europa
Chen et al. A new pose estimation method for non-cooperative spacecraft based on point cloud
CN111104890B (zh) 一种空间飞行器广义模型的识别和重构方法及装置
Koizumi et al. Development of attitude sensor using deep learning
Troiani et al. 1-point-based monocular motion estimation for computationally-limited micro aerial vehicles
CN115205462A (zh) 一种针对光切式的基于半有序点云的三维重建方法
Yan et al. Horizontal velocity estimation via downward looking descent images for lunar landing
Wang et al. Full Period Three-dimensional (3-D) Reconstruction Method for a Low Cost Singlelayer Lidar.

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