发明名称 一种基于四面体边长和体积加权约束的单目视觉位姿测量方法
摘要 本发明公开了一种基于四面体边长和体积加权约束的单目视觉位姿测量方法,该方法的步骤包括:步骤(1)、获取四面体目标图像;步骤(2)、提取图像中四面体目标特征点;步骤(3)、建立四面体体积和边长特征加权约束关系;步骤(4)、基于四面体体积边长加权约束关系解算位姿。该方法运用图像识别手段,以四面体体积和边长特征建立加权误差约束关系,提高了目标单目视觉位姿测量精度。
申请公布号 CN104880178A 申请公布日期 2015.09.02
申请号 CN201510290582.6 申请日期 2015.06.01
申请人 中国科学院光电技术研究所 发明人 赵汝进;王进;赵人杰;王明富;洪裕珍;颜坤;游迪
分类号 G01C11/04(2006.01)I 主分类号 G01C11/04(2006.01)I
代理机构 代理人
主权项 一种基于四面体边长和体积加权约束的单目视觉位姿测量方法,其特征在于该方法包括如下步骤:步骤(1)、对目标成像,并获取目标图像;首先设置目标由4个不共面特征点构成四面体,已知4个目标特征点在目标坐标系中3维坐标集{W<sub>1</sub>,W<sub>2</sub>,W<sub>3</sub>,W<sub>4</sub>},获取四面体边长参数L,体积参数V,<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>L</mi><mi>ij</mi></msub><mo>=</mo><mo>|</mo><mover><mrow><msub><mi>W</mi><mi>i</mi></msub><msub><mi>W</mi><mi>j</mi></msub></mrow><mo>&RightArrow;</mo></mover><mo>|</mo></mrow>]]></math><img file="FDA0000728784240000011.GIF" wi="238" he="113" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>V</mi><mo>=</mo><mrow><mo>(</mo><mover><mrow><msub><mi>W</mi><mi>i</mi></msub><msub><mi>W</mi><mi>j</mi></msub></mrow><mo>&RightArrow;</mo></mover><mo>&times;</mo><mover><mrow><msub><mi>W</mi><mi>i</mi></msub><msub><mi>W</mi><mi>g</mi></msub></mrow><mo>&RightArrow;</mo></mover><mo>)</mo></mrow><mo>&CenterDot;</mo><mover><mrow><msub><mi>W</mi><mi>i</mi></msub><msub><mi>W</mi><mi>h</mi></msub></mrow><mo>&RightArrow;</mo></mover></mrow>]]></math><img file="FDA0000728784240000012.GIF" wi="511" he="93" /></maths>其中i≠j≠g≠h,且i,j,g,h∈{1,2,3,4};步骤(2)、识别四面体目标特征点;将步骤(1)获取的图像进行识别并提取图像中目标特征点坐标{a<sub>i</sub>}(i=1…4);步骤(3)、建立四面体边长和体积加权约束关系;已知摄像机等效焦距为f,定义4个目标特征点在摄像机坐标系下坐标集为{A<sub>1</sub>,A<sub>2</sub>,A<sub>3</sub>,A<sub>4</sub>},定义特征点到摄像机光心p的距离为{l<sub>1</sub>,l<sub>2</sub>,l<sub>3</sub>,l<sub>4</sub>};由步骤(2)获取了目标特征点在图像平面成像点坐标a<sub>i</sub>=(u<sub>i</sub>,v<sub>i</sub>)(i=1…3),则成像点的3维坐标表示为(a<sub>i</sub>,f)=(u<sub>i</sub>,v<sub>i</sub>,f),则:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mover><mi>u</mi><mo>&RightArrow;</mo></mover><mi>i</mi></msub><mo>=</mo><mfrac><mover><mrow><mi>p</mi><mrow><mo>(</mo><msub><mi>a</mi><mi>i</mi></msub><mo>,</mo><mi>f</mi><mo>)</mo></mrow></mrow><mo>&RightArrow;</mo></mover><mrow><mo>|</mo><mover><mrow><mi>p</mi><mrow><mo>(</mo><msub><mi>a</mi><mi>i</mi></msub><mo>,</mo><mi>f</mi><mo>)</mo></mrow></mrow><mo>&RightArrow;</mo></mover><mo>|</mo></mrow></mfrac></mrow>]]></math><img file="FDA0000728784240000013.GIF" wi="303" he="194" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>A</mi><mi>i</mi></msub><mo>=</mo><msub><mi>l</mi><mi>i</mi></msub><mover><msub><mi>u</mi><mi>i</mi></msub><mo>&RightArrow;</mo></mover></mrow>]]></math><img file="FDA0000728784240000014.GIF" wi="180" he="87" /></maths>应用特征点构造边长L和体积V作为对刚体目标的约束,构成了基于边长L和体积V加权约束的误差项E<sub>LV</sub>,1个四面体有6条边,则共形成6组边长L误差约束项,对应形成6组约束加权误差约束项E<sub>LV</sub>;<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>e</mi><mi>L</mi></msub><mo>=</mo><msup><mrow><mo>(</mo><msub><mi>l</mi><mi>i</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>l</mi><mi>j</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><msub><mrow><mn>2</mn><mi>l</mi></mrow><mi>i</mi></msub><msub><mi>l</mi><mi>j</mi></msub><mrow><mo>(</mo><msub><mover><mi>u</mi><mo>&RightArrow;</mo></mover><mi>i</mi></msub><mo>&CenterDot;</mo><msub><mover><mi>u</mi><mo>&RightArrow;</mo></mover><mi>j</mi></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>L</mi><mi>ij</mi></msub></mrow>]]></math><img file="FDA0000728784240000015.GIF" wi="682" he="98" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>e</mi><mi>V</mi></msub><mo>=</mo><mo>[</mo><mrow><mo>(</mo><msub><mi>l</mi><mi>j</mi></msub><mover><msub><mi>u</mi><mi>j</mi></msub><mo>&RightArrow;</mo></mover><mo>-</mo><msub><mi>l</mi><mi>i</mi></msub><mover><msub><mi>u</mi><mi>i</mi></msub><mo>&RightArrow;</mo></mover><mo>)</mo></mrow><mo>&times;</mo><mrow><mo>(</mo><msub><mi>l</mi><mi>g</mi></msub><mover><msub><mi>u</mi><mi>g</mi></msub><mo>&RightArrow;</mo></mover><mo>-</mo><msub><mi>l</mi><mi>i</mi></msub><mover><msub><mi>u</mi><mi>i</mi></msub><mo>&RightArrow;</mo></mover><mo>)</mo></mrow><mo>]</mo><mo>&CenterDot;</mo><mrow><mo>(</mo><msub><mi>l</mi><mi>h</mi></msub><mover><msub><mi>u</mi><mi>h</mi></msub><mo>&RightArrow;</mo></mover><mo>-</mo><msub><mi>l</mi><mi>i</mi></msub><mover><msub><mi>u</mi><mi>i</mi></msub><mo>&RightArrow;</mo></mover><mo>)</mo></mrow><mo>-</mo><mi>V</mi></mrow>]]></math><img file="FDA0000728784240000016.GIF" wi="950" he="94" /></maths>E<sub>LV</sub>=αe<sub>L</sub>+βe<sub>V</sub>其中i≠j≠g≠h,且i,j,g,h∈{1,2,3,4};E为第k次迭代后由e<sub>L</sub>和e<sub>V</sub>形成的加权误差矩阵;步骤(4)、基于四面体边长和体积约束关系解算位姿;首次迭代解算出{l<sub>i</sub>}的修正量x:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>J</mi><mi>l</mi></msub><mo>=</mo><mfrac><msub><mrow><mo>&PartialD;</mo><mi>E</mi></mrow><mi>LV</mi></msub><mrow><mo>&PartialD;</mo><mi>L</mi></mrow></mfrac></mrow>]]></math><img file="FDA0000728784240000021.GIF" wi="221" he="129" /></maths>J<sub>l</sub>x=E<sub>LV</sub>x=(J<sub>l</sub><sup>T</sup>J<sub>l</sub>)<sup>‑1</sup>J<sub>l</sub><sup>‑1</sup>E<sub>LV</sub>其中J<sub>l</sub>为雅可比(Jacobian)矩阵,在迭代中解算出{l<sub>i</sub>}的修正量x,从而最终解算出准确的{l<sub>i</sub>};L<sup>(k+1)</sup>=L<sup>(k)</sup>‑x其中L<sup>(k)</sup>是第k次迭代后的{l<sub>i</sub>};当所有的边长和体积约束都能同时被满足,误差E<sub>LV</sub>在迭代过程中将会逐渐收敛到零,且最终解算得到最优的{l<sub>i</sub>};另一方面,此非线性最优化问题可以被认为是局部线性的,反复迭代可以被解算得到{l<sub>i</sub>};在{l<sub>i</sub>}得到的基础上,解算得到目标特征点在摄像机坐标系中坐标{A<sub>i</sub>},再由最小二乘平差方法,解算得到目标姿态角度θ=(A<sub>x</sub>,A<sub>y</sub>,A<sub>z</sub>)与目标平移向量T=(t<sub>x</sub>,t<sub>y</sub>,t<sub>z</sub>)。
地址 610209 四川省成都市双流350信箱