发明名称 显微视觉系统中运动误差矫正方法
摘要 本发明公开了一种显微视觉系统中运动误差矫正方法,应用于显微视觉系统,首先通过标定板对摄像机进行标定,建立摄像机模型,对标定板的角点进行跟踪,通过最小二乘法拟合出角点的跟踪轨迹和理想旋转轨迹,然后进行误差评估,建立误差矫正模型,从而利用误差矫正模型对显微图像进行矫正。本发明的方法只需通过一次误差评估就能建立起显微视觉系统的误差矫正模型,利用该模型能够矫正系统采集的任意显微运动视频,减小系统的运动误差,提高三维重建精度,具有较好的实用性。
申请公布号 CN102842117B 申请公布日期 2015.02.25
申请号 CN201210243071.5 申请日期 2012.07.13
申请人 浙江工业大学 发明人 刘盛;翟斌斌;金海强;陈胜勇
分类号 G06T5/00(2006.01)I;G06T7/00(2006.01)I 主分类号 G06T5/00(2006.01)I
代理机构 杭州天勤知识产权代理有限公司 33224 代理人 胡红娟
主权项 一种显微视觉系统中运动误差矫正方法,应用于显微视觉系统,所述显微视觉系统包括用于获取被观测物体的显微图像的单目光学显微镜和摄像机,载物台以及控制所述载物台作一定倾斜角度的旋转运动的运动控制设备,以及用来进行视觉处理的计算机系统,其特征在于,所述矫正方法包括步骤:(1)建立误差矫正模型,包括步骤:(1.1)采用标定板对所述摄像机进行标定,确定所述摄像机的内外参数,建立显微视觉系统的摄像机模型;(1.2)对标定板进行角点跟踪,得到N个角点的跟踪轨迹;(1.3)确定显微视觉系统载物台旋转轴;(1.4)建立理想的旋转运动轨迹;(1.5)构建误差矫正模型;其中,步骤(1.5)包括步骤:空间点P绕所述旋转轴旋转角度θ得到理想点坐标P’(x’,y’,z’)和实际点坐标<img file="FDA0000587271660000011.GIF" wi="202" he="75" />p’(u’,v’)和<img file="FDA0000587271660000012.GIF" wi="124" he="69" />分别为P’、<img file="FDA0000587271660000013.GIF" wi="46" he="66" />的二维图像点,那么满足下面两个公式:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>s</mi><mfenced open='[' close=']'><mtable><mtr><mtd><msup><mi>u</mi><mo>&prime;</mo></msup></mtd></mtr><mtr><mtd><msup><mi>v</mi><mo>&prime;</mo></msup></mtd></mtr><mtr><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mo>=</mo><mi>A</mi><mfenced open='[' close=']'><mtable><mtr><mtd><mi>R</mi></mtd><mtd><mi>T</mi></mtd></mtr></mtable></mfenced><mfenced open='[' close=']'><mtable><mtr><mtd><msup><mi>x</mi><mo>&prime;</mo></msup></mtd></mtr><mtr><mtd><msup><mi>y</mi><mo>&prime;</mo></msup></mtd></mtr><mtr><mtd><msup><mi>z</mi><mo>&prime;</mo></msup></mtd></mtr><mtr><mtd><mn>1</mn></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000587271660000014.GIF" wi="357" he="261" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>s</mi><mfenced open='[' close=']'><mtable><mtr><mtd><mover><mi>u</mi><mo>^</mo></mover></mtd></mtr><mtr><mtd><mover><mi>v</mi><mo>^</mo></mover></mtd></mtr><mtr><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mo>=</mo><mi>A</mi><mfenced open='[' close=']'><mtable><mtr><mtd><mi>R</mi></mtd><mtd><mi>T</mi></mtd></mtr></mtable></mfenced><mfenced open='[' close=']'><mtable><mtr><mtd><mover><mi>x</mi><mo>^</mo></mover></mtd></mtr><mtr><mtd><mover><mi>y</mi><mo>^</mo></mover></mtd></mtr><mtr><mtd><mover><mi>z</mi><mo>^</mo></mover></mtd></mtr><mtr><mtd><mn>1</mn></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000587271660000015.GIF" wi="344" he="267" /></maths>其中s是尺度因子,A为摄像机内参矩阵,[R T]为摄像机外参矩阵;将上述公式相减可得:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>s</mi><msubsup><mover><mi>d</mi><mo>~</mo></mover><mi>&theta;</mi><mi>T</mi></msubsup><mo>=</mo><msubsup><mi>ARD</mi><mi>&theta;</mi><mi>T</mi></msubsup><mo>,</mo></mrow>]]></math><img file="FDA0000587271660000016.GIF" wi="477" he="78" /></maths>其中<img file="FDA0000587271660000017.GIF" wi="388" he="67" />为世界坐标系上的位置偏移量,<img file="FDA0000587271660000018.GIF" wi="296" he="68" />为二维图像坐标系上的位置偏移量,而<img file="FDA0000587271660000019.GIF" wi="58" he="78" />为d<sub>θ</sub>的增0向量;其中<img file="FDA00005872716600000110.GIF" wi="56" he="76" />通过所述N个角点的跟踪轨迹与理想的旋转运动轨迹得到:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mover><mi>d</mi><mo>~</mo></mover><mi>&theta;</mi></msub><mo>=</mo><mrow><mo>(</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mrow><mo>(</mo><msub><mi>C</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow><mi>cos</mi><mi>&theta;</mi><mo>-</mo><msub><mover><mi>C</mi><mo>^</mo></mover><mi>n</mi></msub><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow><mi>cos</mi><mi>&theta;</mi><mo>)</mo></mrow><mo>,</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mrow><mo>(</mo><msub><mi>C</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow><mi>sin</mi><mi>&theta;</mi><mo>-</mo><msub><mover><mi>C</mi><mo>^</mo></mover><mi>n</mi></msub><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow><mi>sin</mi><mi>&theta;</mi><mo>)</mo></mrow><mo>,</mo><mn>0</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA00005872716600000111.GIF" wi="1326" he="153" /></maths>其中C<sub>n</sub>为第n个角点的理想运动轨迹的两次曲线,<img file="FDA00005872716600000112.GIF" wi="136" he="89" />为第n个角点的跟踪轨迹的两次曲线;经过变换得到误差矫正模型E<sub>θ</sub>为:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>E</mi><mi>&theta;</mi></msub><mo>=</mo><msubsup><mi>D</mi><mi>&theta;</mi><mi>T</mi></msubsup><mo>=</mo><msup><mi>sR</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><msup><mi>A</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><msubsup><mover><mi>d</mi><mo>~</mo></mover><mi>&theta;</mi><mi>T</mi></msubsup></mrow>]]></math><img file="FDA0000587271660000021.GIF" wi="471" he="79" /></maths>其中s是尺度因子,A为摄像机内参矩阵,[R T]为摄像机外参矩阵;(2)调整载物台位置以适应被观测物体,采用标定板对摄像机进行标定,计算出摄像机的内外参数;(3)根据所述摄像机录制的被观测物体的显微图像对应的旋转角度,将步骤(2)得到的摄像机的内外参数代入所述误差矫正模型,计算所述旋转角度对应的偏移矫正向量,对所述显微图像进行矫正。
地址 310014 浙江省杭州市下城区潮王路18号