发明名称 基于目标旋转对称特征的位姿测量方法
摘要 基于目标旋转对称特征的位姿测量方法,属于计算机视觉测量技术领域。其特征是本发明所要测量的目标上的图案具有旋转对称性,拍摄一幅包含实际图案投影的目标图片,根据实际图案的旋转对称性通过径向展开变换和低秩特征恢复的方法由目标图片恢复出对正图片,根据实际图案上的点在目标图片和对正图片上投影点之间的关系解出目标的位姿参数。本发明的效果和益处是,利用一幅具有旋转对称特征的图案即可完成目标位姿的测量,将视觉测量问题与图像本身的对称性建立了联系,不需要提取图像的底层特征,减少了人工干预,易于实现自动化,能够保证测量精度。
申请公布号 CN104596486A 申请公布日期 2015.05.06
申请号 CN201410787089.0 申请日期 2014.12.15
申请人 大连理工大学 发明人 孙怡;刘睿智;卢健
分类号 G01C11/08(2006.01)I 主分类号 G01C11/08(2006.01)I
代理机构 大连理工大学专利中心 21200 代理人 侯明远
主权项 基于目标旋转对称特征的位姿测量方法,采用视觉传感器采集测量目标的图像,通过图像处理系统对采集的图像进行分析处理来完成目标在空间的位置、姿态信息的测量,其特征在于,步骤如下:S1获取目标上的旋转对称图案的灰度图片,选择图片中具有旋转对称图案投影的矩形区域,该矩形区域构成一个矩阵,记作I;获取相机的焦距f和像元尺寸d<sub>x</sub>,则f<sub>d</sub>=f/d<sub>x</sub>;获取目标上的旋转对称图案的直径d<sub>real</sub>;S2将参数始初化,角度参数初始值θ<sub>x</sub>=0,θ<sub>y</sub>=0,θ<sub>z</sub>=0,归一化位移参数的初始值<img file="FDA0000633844900000011.GIF" wi="481" he="101" />参数用向量的形式表示为<img file="FDA0000633844900000012.GIF" wi="575" he="101" />初始化的参数用向量的形式表示为τ<sub>init</sub>=[0,0,0,0,0,1];S3将步骤S1中选择的矩形区域I和初化的τ<sub>init</sub>作为条件,求解式(11)的优化问题:<img file="FDA0000633844900000013.GIF" wi="1496" he="134" />其中<img file="FDA0000633844900000014.GIF" wi="380" he="96" />E是约束方程,Iοτ表示根据向量τ中的参数对矩阵I进行变换,设矩阵I上的一点为I(u<sub>2</sub>,v<sub>2</sub>),经过变换后该点的坐标为(u<sub>1</sub>,v<sub>1</sub>),则(u<sub>1</sub>,v<sub>1</sub>)与(u<sub>2</sub>,v<sub>2</sub>)的关系用式(6)来表示:<img file="FDA0000633844900000015.GIF" wi="1780" he="331" /><img file="FDA0000633844900000016.GIF" wi="1513" he="322" />在式(6)中R<sub>11</sub>=C<sub>y</sub>C<sub>z</sub>,R<sub>21</sub>=S<sub>x</sub>S<sub>y</sub>C<sub>z</sub>‑C<sub>x</sub>S<sub>z</sub>,R<sub>31</sub>=C<sub>x</sub>S<sub>y</sub>C<sub>z</sub>+S<sub>x</sub>S<sub>z</sub>,R<sub>12</sub>=C<sub>y</sub>S<sub>z</sub>,R<sub>22</sub>=S<sub>x</sub>S<sub>y</sub>S<sub>z</sub>+C<sub>x</sub>C<sub>z</sub>,R<sub>32</sub>=C<sub>x</sub>S<sub>y</sub>S<sub>z</sub>‑S<sub>x</sub>C<sub>z</sub>,其中,C<sub>x</sub>=cosθ<sub>x</sub>,S<sub>x</sub>=sinθ<sub>x</sub>,C<sub>y</sub>=cosθ<sub>y</sub>,S<sub>y</sub>=sinθ<sub>y</sub>,C<sub>z</sub>=cosθ<sub>z</sub>,S<sub>z</sub>=sinθ<sub>z</sub>;F{Iοτ}表示 对Iοτ进行径向展开变换,设变换后得到图像I<sub>p</sub>,I<sub>p</sub>上任意一点为I<sub>p</sub>(i,j),在进行FE变换时,并不是对整幅图片进行变换,而是根据需要选取一对同心圆,对内圆和外圆之间的环面做FE变换;根据式(7)写出环面的FE变换的关系如下:<img file="FDA0000633844900000021.GIF" wi="1905" he="178" />其中,m,n分别表示图像I<sub>p</sub>的行数和列数,(i,j)表示图像I<sub>p</sub>中任意一点,i=1,2,3,…,m,j=1,2,3,…,n,r<sub>1</sub>、r<sub>2</sub>分别是内外环的半径;对Iοτ而言,m表示径向采样点的个数,则径向采样间隔为<img file="FDA0000633844900000022.GIF" wi="181" he="147" />n为旋转方向的采样个数,则角度采样间隔为<img file="FDA0000633844900000023.GIF" wi="125" he="158" />式(11)中的<img file="FDA0000633844900000024.GIF" wi="438" he="101" />是目标函数,<img file="FDA0000633844900000025.GIF" wi="232" he="99" />表示求矩阵<img file="FDA0000633844900000026.GIF" wi="63" he="100" />的秩,λ是一个常数,E为误差矩阵,||E||<sub>0</sub>表示求误差矩阵E的零范数;<img file="FDA0000633844900000027.GIF" wi="66" he="102" />为I<sub>p</sub>减去误差矩阵E之后得到的矩阵;式(11)求解的具体步骤如下:S3.1对式(11)约束方程中的F{Iοτ}做归一化处理:<img file="FDA0000633844900000028.GIF" wi="989" he="197" />表示求F{Iоτ}的F范数,即<img file="FDA0000633844900000029.GIF" wi="966" he="207" />S3.2对式(11)的目标函数进行凸松弛并对其约束方程进行线性化,得到式(13)如下:<img file="FDA00006338449000000210.GIF" wi="1684" he="138" />其中<img file="FDA0000633844900000031.GIF" wi="124" he="113" />表示求矩阵<img file="FDA0000633844900000032.GIF" wi="68" he="101" />的核范数,||E||<sub>1</sub>表示求矩阵E的1范数,J表示||F{Iοτ}||<sub>norm</sub>的雅克比矩阵,即<img file="FDA0000633844900000033.GIF" wi="644" he="210" />Δτ表示τ附近的增量,<img file="FDA0000633844900000034.GIF" wi="876" he="97" />S3.3利用增广拉格朗日乘子法来求解式(13)的优化问题<img file="FDA0000633844900000035.GIF" wi="1383" he="142" />其对应的拉格朗日函数是:<img file="FDA0000633844900000036.GIF" wi="1741" he="145" />其中,<img file="FDA0000633844900000037.GIF" wi="1696" he="117" />Y是一个拉格朗日乘子阵,μ>0,用来折中约束部分与目标函数部分,&lt;·,·&gt;表示矩阵内积算子,||·||<sub>F</sub>表示矩阵的F范数;则式(13)中的问题转化为同解的无约束的凸优化问题:<img file="FDA0000633844900000038.GIF" wi="1533" he="157" />利用交替方向迭代法求解式(14),分别迭代求解各个变量的最优解,步骤如下:S3.3.1参数的初始值设置为:Y<sub>0</sub>=0,E<sub>0</sub>=0,Δτ<sub>0</sub>=0,μ<sub>0</sub>>0,ρ>1,k=0,t<sub>inner</sub>>0,t<sub>outer</sub>>0,f<sub>pre</sub>=0;S3.3.2<img file="FDA0000633844900000039.GIF" wi="1180" he="111" />其中<img file="FDA00006338449000000310.GIF" wi="864" he="109" />表示对<img file="FDA00006338449000000311.GIF" wi="727" he="101" />做SVD分解;S3.3.3<img file="FDA00006338449000000312.GIF" wi="637" he="111" />其中<img file="FDA00006338449000000313.GIF" wi="240" he="110" />表示对Σ<sub>k</sub>进行奇异值收缩,收缩阈值为<img file="FDA00006338449000000314.GIF" wi="123" he="93" />收缩算子D<sub>ξ</sub>(·)也称作软阈值操作算子,若定义在标量上,其表示形式写作 D<sub>ξ</sub>(x)=max(sign(x)·(|x|‑ξ),0);S3.3.4<img file="FDA0000633844900000041.GIF" wi="1286" he="109" />其中<img file="FDA0000633844900000042.GIF" wi="1106" he="121" />表示对矩阵的每个元素做收缩,收缩阈值为<img file="FDA0000633844900000043.GIF" wi="151" he="92" />S3.3.5<img file="FDA0000633844900000044.GIF" wi="1356" he="104" />S3.3.6<img file="FDA0000633844900000045.GIF" wi="1168" he="117" />S3.3.7μ<sub>k+1</sub>=ρμ<sub>k</sub>;S3.3.8<img file="FDA0000633844900000046.GIF" wi="902" he="113" />判断s<t<sub>inner</sub>是否为真,若为假,返回步骤S3.3.2继续执行,若为真,继续执行下一步;S3.3.9更新变换参数τ=τ+Δτ<sub>k+1</sub>,<img file="FDA0000633844900000047.GIF" wi="664" he="134" />判断 |f<sub>p</sub>‑f<sub>pre</sub>|<t<sub>outer</sub>是否为真,若为假,f<sub>pre</sub>=f<sub>p</sub>并返回步骤S3.1继续执行,若为真,继续执行下一步;S3.3.10输出<img file="FDA0000633844900000048.GIF" wi="321" he="104" />τ,E=E<sub>k+1</sub>作为最终解;S4根据步骤S3的计算结果求解目标的位置和姿态参数;其中姿态参数θ<sub>x</sub>,θ<sub>y</sub>,θ<sub>z</sub>直接从τ中得到;位置参数<img file="FDA0000633844900000049.GIF" wi="276" he="93" />是经过归一化的,<img file="FDA00006338449000000410.GIF" wi="244" he="138" /><img file="FDA00006338449000000411.GIF" wi="492" he="145" />用霍夫变换获取对正图片Iοτ中旋转对称图形的直径d<sub>image</sub>,则d=f·d<sub>real</sub>/(d<sub>image</sub>·d<sub>x</sub>),其中f为焦距,d<sub>x</sub>为像元尺寸;根据求出的d得到:<img file="FDA00006338449000000412.GIF" wi="847" he="100" />
地址 116024 辽宁省大连市甘井子区凌工路2号