一种基于点特征的单目视觉位姿测量方法

文档序号:1626270 发布日期:2020-01-14 浏览:9次 >En<

阅读说明:本技术 一种基于点特征的单目视觉位姿测量方法 (Monocular vision pose measuring method based on point characteristics ) 是由 王中宇 李亚茹 程银宝 于 2019-10-29 设计创作,主要内容包括:发明涉及一种基于点特征的单目视觉位姿测量方法,属于单目视觉位姿测量领域,其特征在于:(1)根据共面4个特征点之间的位置关系,分别对平行和相交两种情况进行分析;(2)根据特征点的空间坐标、图像坐标和空间位置关系,推导出世界坐标系中3个坐标轴上的向量变换到摄像机坐标系中的单位向量,进而求解出物体相对于摄像机的初始位姿;(3)用LM算法对初始位姿进行优化,得到最终位姿。上述方法仅需给出任意4个特征点的空间相对位置,其中任意3个特征点不共线,不需要给定其他限制条件,就可求解出物体相对于摄像机的位姿。解决单目视觉位姿求解过程中存在的解不唯一、选解难和解的精度不高等问题。(The invention relates to a monocular vision pose measurement method based on point characteristics, which belongs to the field of monocular vision pose measurement and is characterized in that: (1) according to the position relation among the coplanar 4 characteristic points, analyzing the two conditions of parallel and intersection respectively; (2) deducing vectors on 3 coordinate axes in a world coordinate system to be converted into unit vectors in a camera coordinate system according to the spatial coordinates, the image coordinates and the spatial position relation of the feature points, and further solving the initial pose of the object relative to the camera; (3) and optimizing the initial pose by using an LM algorithm to obtain a final pose. The method can solve the pose of the object relative to the camera only by giving the spatial relative positions of any 4 characteristic points, wherein any 3 characteristic points are not collinear, and other limiting conditions are not required to be given. The method solves the problems that the solution is not unique, the solution is difficult to select, the precision is not high and the like in the monocular vision pose solving process.)

一种基于点特征的单目视觉位姿测量方法

技术领域

本发明属于单目视觉位姿测量领域,涉及一种基于点特征的单目视觉位姿测量方法,解决单目视觉位姿求解过程中存在的解不唯一、选解难和解的精度不高等问题。

技术背景

视觉测量技术是基于计算机视觉将三维空间中的物体通过相机投射到二维图像中,然后借助于数学模型从二维图像还原为三维信息的技术,能够实现物体的几何尺寸以及物体在空间中的位置、姿态等测量。位置和姿态简称为位姿,通常用6个自由度来表示,即两个坐标系之间的3个旋转量和3个平移量。位姿测量目前已经广泛地应用于航空航天、目标跟踪、增强现实、自动导航、视觉伺服控制和工业制造等领域。求解位姿时首先需要选取特征,常用的有点特征、直线特征和二次曲线特征等。基于直线和曲线特征的算法一般是由点特征算法演变过来的,因此研究点特征具有重要意义。

基于点特征的位姿估计问题又称为PnP(Perspective-n-Point)问题。国内外学者对此进行了大量研究,主要集中于点的个数n及其算法,初步结论为当n<3时,PnP问题有无穷多个解;当n=3时,PnP问题最多有4个解,如果3个点组成的三角形为等腰三角形,并且摄像机在特定区域内时,PnP问题有唯一解;当n=4时,如果4点共面,则PnP问题有唯一解,如果4点不共面,则PnP问题最多有4个解;当n=5时,PnP问题最多有2个解;当n>5时,PnP问题有唯一解,但解的形式不封闭,因此求解的算法具有局限性。除此之外,依据呈四边形分布的共面4特征点,用相应的几何关系求解位姿等方法也有应用,但是存在特征点提取困难、解不唯一、鲁棒性差和迭代结果不收敛等缺陷,并且在实际的测量场景中,特征点往往是未知的,此时特征点提取精度对位姿解算的精度有极大的影响。

针对这些问题,本发明提出了一种基于点特征的单目视觉位姿测量方法,根据任意不共线的4个特征点的空间相对位置,不需要给定特征点与光心连线形成的夹角、不需要给定光心到特征点的距离以及限制4个特征点呈矩形或特殊形状的分布等条件,利用给定4个不共线特征点相对位置已知的特性,求解出物体相对于摄像机的位姿。

发明内容

本发明的目的在于:提供一种基于点特征的单目视觉位姿测量方法,针对单目视觉位姿求解过程中存在的解不唯一、选解难和解的精度不高等问题。

本发明属于单目视觉位姿测量领域,涉及一种基于点特征的单目视觉位姿测量方法,解决单目视觉位姿求解过程中存在的解不唯一、选解难和解的精度不高等问题。

本发明的技术解决方案:一种基于点特征的单目视觉位姿测量方法,由以下步骤实现:

步骤一:获取摄像机拍摄的图像,用张正友标定法对摄像机进行标定,获取摄像机内参数矩阵;

步骤二:提取图像中任意4个不共面的特征点Pi(i=1,2,3,4),计算特征点在成像平面坐标系中的坐标{qi(xi,yi,f),i=1,2,3,4},根据两点间距离公式计算出任意两个特征点之间的距离:d12=|P1P2|、d23=|P2P3|、d13=|P1P3|、d14=|P1P4|;

步骤三:设P1为世界坐标系的原点,Oc为光心点,分析给定的任意不共线的4个特征点之间的相对位置关系:当4个特征点所在的两条直线平行时,定义特征点相关距离因子{ki,i=1,2,3,4};当4个特征点所在的两条直线相交时,分别计算两条直线交点在世界坐标系中的坐标及成像坐标系中的坐标;

步骤四:分别求解出世界坐标系中Xw、Yw、Zw轴上的向量变换到摄像机坐标系中的单位向量(r11,r21,r31);

步骤五:根据世界坐标系的原点在摄像机中的坐标求解出物体相对于摄像机的初始位置,并通过LM算法对初始位姿进行优化,得到物体相对于摄像机的初始姿态。

其中,步骤三和步骤四所述的分析给定的4个特征点的相对位置关系,求解出世界坐标系中Xw、Yw、Zw轴上的向量变换到摄像机坐标系中的单位向量(r11,r21,r31)的具体实现过程如下:

根据任意3点不共线的共面4个特征点Pi(i=1,2,3,4)的位置关系,对P1、P2两点所在的直线P1P2与P3、P4两点所在的直线P3P4之间平行和相交两种情况分别分析:

(1)当P1P2与P3P4平行时

设P1为世界坐标系的原点,Oc为光心点,点Oc、q3、q4确定的平面为π1,点Oc、q1、q2确定的平面为π2,点Oc到靶标特征点Pi(i=1,2,3,4)的距离与到像点qi(i=1,2,3,4)的距离之比分别为ki=|OcPi|/|Ocqi|(i=1,2,3,4)。

根据向量

Figure BDA0002251021100000021

和相机内参数矩阵计算得到平面π1的标准法向量

Figure BDA0002251021100000022

由于直线P1P2与直线P3P4平行,可得方程:

由|P1P2|=d12可得方程:

Figure BDA0002251021100000031

联立方程(1)和方程(2)解得:

同理,根据|P3P4|=d34解得:

所以,世界坐标系中Xw轴上的向量变换到摄像机坐标系中的单位向量(r11,r21,r31)为

Figure BDA0002251021100000034

Zw轴上的向量变换到摄像机坐标系中的单位向量(r13,r23,r33)为Yw轴上的向量变换到摄像机坐标系中的单位向量(r12,r22,r32)为(r13,r23,r33)×(r11,r21,r31)。

从而得到世界坐标系到摄像机坐标系的旋转矩阵R为:

Figure BDA0002251021100000036

由P1为世界坐标系的原点可知,世界坐标系到摄像机坐标系的平移向量为

Figure BDA0002251021100000037

(2)当P1P2与P3P4相交时

设P3P4与P1P2相交于点P0,根据给定的4个特征点Pi(i=1,2,3,4)在世界坐标系中的坐标得到P0在世界坐标系中的坐标,将点P0设为世界坐标系的原点,P0在图像坐标系中的投影点为q0(x0,y0,f),分别令|P0P1|=d01、|P0P2|=d02、|P0P3|=d03、|P0P4|=d04、|P0P5|=d05。根据投影点qi(xi,yi,f)(i=0,1,2,3,4)得到∠P1OcP0、∠P2OcP1、∠P2OcP0、∠P4OcP0、∠P4OcP3和∠P3OcP0

在ΔOcP0P2和ΔOcP0P1中,由正弦定理可得:

Figure BDA0002251021100000039

解得点P0在摄像机坐标系中的坐标

Figure BDA00022510211000000310

为:

Figure BDA00022510211000000311

同理,在ΔOcP0P2和ΔOcP1P2中,由正弦定理求解出点P2在摄像机坐标系下的坐标为:

Figure BDA0002251021100000042

在ΔOcP0P3和ΔOcP3P4中,由正弦定理求解出点P3在摄像机坐标系中的坐标

Figure BDA0002251021100000043

为:

Figure BDA0002251021100000044

则有:

Figure BDA0002251021100000045

根据式(8)和式(9)解得世界坐标系到摄像机坐标系的旋转矩阵R和平移向量t。则当cosβ≠0时,已知旋转矩阵R,求解得到相应的欧拉角为:

Figure BDA0002251021100000047

综上所述,从不同分布的4个特征点出发,由距离和角度的约束以及相应的几何关系均可得到唯一的位姿解。

附图说明

图1为本发明一种基于点特征的单目视觉位姿测量方法的流程图;

图2为本发明一种基于点特征的单目视觉位姿测量方法中当P1P2与P3P4平行时的位姿测量示意图;

图3为本发明一种基于点特征的单目视觉位姿测量方法中当P1P2与P3P4相交时的位姿测量示意图;

具体实施方式

下面结合附图和实施例对本发明作进一步说明。

本发明的一种基于点特征的单目视觉位姿测量方法,针对位姿求解过程中存在的解不唯一、选解难和解的精度不高等问题,提出一种基于点特征的单目视觉位姿测量算法。

实施例1

如图1所示,本发明的一种基于点特征的单目视觉位姿测量方法,具体包括以下步骤:

步骤一:获取摄像机拍摄的图像,用张正友标定法对摄像机进行标定,获取摄像机内参数矩阵;

步骤二:如图2所示提取图像中任意4个不共面的特征点Pi(i=1,2,3,4),计算特征点在成像平面坐标系中的坐标{qi(xi,yi,f),i=1,2,3,4},根据两点间距离公式计算出任意两个特征点之间的距离:d12=|P1P2|、d23=|P2P3|、d13=|P1P3|、d14=|P1P4|;

步骤三:由于P1、P2两点所在的直线P1P2与P3、P4两点所在的直线P3P4相互平行,设P1为世界坐标系的原点,Oc为光心点,点Oc、q3、q4确定的平面为π1,点Oc、q1、q2确定的平面为π2,点Oc到靶标特征点Pi(i=1,2,3,4)的距离与到像点qi(i=1,2,3,4)的距离之比分别为ki=|OcPi|/|Ocqi|(i=1,2,3,4)。

根据向量

Figure BDA0002251021100000051

和相机内参数矩阵计算得到平面π1的标准法向量

Figure BDA0002251021100000052

由于直线P1P2与直线P3P4平行,即

Figure BDA0002251021100000053

可得方程:

Figure BDA0002251021100000054

由|P1P2|=d12可得方程:

Figure BDA0002251021100000055

联立方程(1)和方程(2)解得:

Figure BDA0002251021100000056

Figure BDA0002251021100000057

从而得到点P1和P2在摄像机坐标系中的坐标分别为:

Figure BDA0002251021100000059

则:

Figure BDA0002251021100000061

同理,根据|P3P4|=d34可求解出点P3、P4在摄像机坐标系下的坐标分别为:

Figure BDA0002251021100000063

则:

Figure BDA0002251021100000064

步骤四:分别求解出世界坐标系中Xw轴上的向量变换到摄像机坐标系中的单位向量(r11,r21,r31)为

Figure BDA0002251021100000065

Zw轴上的向量变换到摄像机坐标系中的单位向量(r13,r23,r33)为

Figure BDA0002251021100000066

Yw轴上的向量变换到摄像机坐标系中的单位向量(r12,r22,r32)为(r13,r23,r33)×(r11,r21,r31)。

步骤五:由P1为世界坐标系的原点解得世界坐标系到摄像机坐标系的旋转矩阵R和平移向量t为:

Figure BDA0002251021100000067

并通过LM算法对初始位姿进行优化,得到物体相对于摄像机的初始位姿。

实施例2

如图1和图3所示,本发明的一种基于点特征的单目视觉位姿测量方法,具体包括以下步骤:

步骤一:获取摄像机拍摄的图像,用张正友标定法对摄像机进行标定,获取摄像机内参数矩阵;

步骤二:提取图像中任意4个不共面的特征点Pi(i=1,2,3,4),计算特征点在成像平面坐标系中的坐标{qi(xi,yi,f),i=1,2,3,4},根据两点间距离公式计算出任意两个特征点之间的距离:d12=|P1P2|、d23=|P2P3|、d13=|P1P3|、d14=|P1P4|;

步骤三:设P3P4与P1P2相交于点P0,根据给定的4个特征点Pi(i=1,2,3,4)在世界坐标系中的坐标得到P0在世界坐标系中的坐标,将点P0设为世界坐标系的原点,P0在图像坐标系中的投影点为q0(x0,y0,f),分别令|P0P1|=d01、|P0P2|=d02、|P0P3|=d03、|P0P4|=d04、|P0P5|=d05。根据投影点qi(xi,yi,f)(i=0,1,2,3,4)得到∠P1OcP0、∠P2OcP1、∠P2OcP0、∠P4OcP0、∠P4OcP3和∠P3OcP0

在ΔOcP0P2和ΔOcP0P1中,由正弦定理可得:

Figure BDA0002251021100000071

Figure BDA0002251021100000072

解得:

则点P0在摄像机坐标系中的坐标为:

同理,在ΔOcP0P2和ΔOcP1P2中,由正弦定理求解出点P2在摄像机坐标系下的坐标

Figure BDA0002251021100000076

为:

Figure BDA0002251021100000077

在ΔOcP0P3和ΔOcP3P4中,由正弦定理求解出点P3在摄像机坐标系中的坐标

Figure BDA0002251021100000078

为:

Figure BDA0002251021100000079

则有:

Figure BDA00022510211000000710

步骤四:根据式(8)和式(9)分别求解出世界坐标系中Xw轴上的向量变换到摄像机坐标系中的单位向量(r11,r21,r31)为

Figure BDA00022510211000000712

Zw轴上的向量变换到摄像机坐标系中的单位向量(r13,r23,r33)为

Figure BDA00022510211000000713

Yw轴上的向量变换到摄像机坐标系中的单位向量(r12,r22,r32)为(r13,r23,r33)×(r11,r21,r31)。

步骤五:当cos β≠0时,已知旋转矩阵R,根据世界坐标系的原点在摄像机中的坐标求解出物体相对于摄像机的初始位置,并通过LM算法对初始位姿进行优化,得到世界坐标系相对于摄像机的初始姿态的欧拉角为:

Figure BDA0002251021100000081

以上结合附图对本发明的具体实施方式作了说明,但这些说明不能被理解为限制了本发明的范围,本发明的保护范围由随附的权利要求书限定,任何在本发明权利要求基础上进行的改动都是本发明的保护范围。

12页详细技术资料下载
上一篇:一种医用注射器针头装配设备
下一篇:基于三维激光扫描的建筑物外立面测量方法及系统

网友询问留言

已有0条留言

还没有人留言评论。精彩留言会获得点赞!

精彩留言,会给你点赞!