发明名称 基于折射图像偏差的三维信息视觉测量方法
摘要 本发明基于折射图像偏差的三维信息视觉测量方法属于计算机视觉测量技术领域,涉及一种基于折射图像偏差的三维信息视觉测量方法。测量方法采用单摄像机先拍摄一张图像A1,然后,将已知折射率的玻璃板以任意角度放入相机前,再利用摄像机拍摄一张图像A2,利用两张图像的偏差量,完成被测物体空间三维信息的测量。此方法只需单一摄像机搭配厚度已知的平行玻璃板,不需要辅助激光以及辅助投影仪,也不需要已知图像中的先验知识,即可对于图像内的特征进行三维信息的重建,降低了测量成本、增加了测量效率。实现了单目相机对于全视场内的快速测量。
申请公布号 CN105571518A 申请公布日期 2016.05.11
申请号 CN201610035785.5 申请日期 2016.01.19
申请人 大连理工大学 发明人 刘巍;马鑫;贾振元;鲁继文;李肖
分类号 G01B11/24(2006.01)I 主分类号 G01B11/24(2006.01)I
代理机构 大连理工大学专利中心 21200 代理人 关慧贞
主权项 一种基于折射图像偏差的三维信息视觉测量方法,采用单摄像机先拍摄一张图像A1,然后将已知折射率的玻璃板以任意角度放入相机前,再利用摄像机拍摄一张图像A2,利用两张图像的偏差量,完成被测物体空间三维信息的测量;测量方法的具体步骤如下:第一步进行单摄像机的标定测量方法采用单摄像机相对固定的方式,以精密加工的靶标板对摄像机进行标定的方法;以张正友等人提出的基于2维平面靶标的摄像机标定方法为基础,并针对测量过程中相对距离高精度的测量要求,实现单目系统在较大视场下的高精度标定,单摄像机模型如下所示:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>Z</mi><mi>c</mi></msub><mfenced open = "[" close = "]"><mtable><mtr><mtd><mi>u</mi></mtd></mtr><mtr><mtd><mi>v</mi></mtd></mtr><mtr><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>C</mi><mi>x</mi></msub></mtd><mtd><mn>0</mn></mtd><mtd><msub><mi>u</mi><mn>0</mn></msub></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><msub><mi>C</mi><mi>y</mi></msub></mtd><mtd><msub><mi>v</mi><mn>0</mn></msub></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mo>&lsqb;</mo><mtable><mtr><mtd><mi>R</mi></mtd><mtd><mi>T</mi></mtd></mtr></mtable><mo>&rsqb;</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>X</mi><mi>w</mi></msub></mtd></mtr><mtr><mtd><msub><mi>Y</mi><mi>w</mi></msub></mtd></mtr><mtr><mtd><msub><mi>Z</mi><mi>w</mi></msub></mtd></mtr><mtr><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mo>=</mo><mi>M</mi><mi>Q</mi><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>X</mi><mi>w</mi></msub></mtd></mtr><mtr><mtd><msub><mi>Y</mi><mi>w</mi></msub></mtd></mtr><mtr><mtd><msub><mi>Z</mi><mi>w</mi></msub></mtd></mtr><mtr><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000909374610000011.GIF" wi="1084" he="254" /></maths>其中,Z<sub>c</sub>为比例因子,(u<sub>0</sub>,v<sub>0</sub>)为图像的主点坐标,(C<sub>x</sub>,C<sub>y</sub>)为横纵方向的等效焦距,R,T分别为摄像机坐标系相对于世界坐标系的旋转和平移矩阵,(X<sub>w</sub>,Y<sub>w</sub>,Z<sub>w</sub>)为控制点在世界坐标系下的坐标,M为相机内参数矩阵,Q为相机外参数矩阵,(u,v)为控制点在摄像机平面上的像素坐标;第二步标记点特征的提取采用梯度重心法分别提取A1、A2图像上的特征,首先采用高斯一阶微分算子对图像进行卷积操作,以获得图像在各点的梯度,则对于标志图案所在子区域的梯度重心点坐标的求取方法如下:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>C</mi><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mo>-</mo><mi>h</mi></mrow><mi>h</mi></munderover><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mo>-</mo><mi>w</mi></mrow><mi>w</mi></munderover><mo>&lsqb;</mo><mo>|</mo><mi>G</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>|</mo><mo>&CenterDot;</mo><mi>P</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>&rsqb;</mo><mo>/</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mo>-</mo><mi>h</mi></mrow><mi>h</mi></munderover><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mo>-</mo><mi>w</mi></mrow><mi>w</mi></munderover><mo>|</mo><mi>G</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>|</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000909374610000012.GIF" wi="1461" he="167" /></maths>其中,C为应用梯度重心法所提取的标记点中心像素坐标值,(i,j)为标记点坐标,|G(i,j)|为(i,j)点的梯度幅值,w,h分别为目标图像的宽和高,P(i,j)为(i,j)点的图像坐标;第三步三维信息求解由于添加的玻璃板是平行板,所以通过以下方程求解平行玻璃板的法向量n;n(v0×b)=0  (3)其中,n为平行玻璃板的法向量,v0是图像上标记点坐标的空间方向矢量,b是同一特征点添加平行玻璃板前后的偏差矢量;在无平行板遮挡前,通过小孔成像模型得到公式(4),<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mfrac><mrow><msub><mi>v</mi><mn>0</mn></msub><msup><mi>P</mi><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></msup></mrow><mi>f</mi></mfrac><mo>=</mo><mi>P</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000909374610000021.GIF" wi="525" he="126" /></maths>在平行板遮挡后,通过多折射模型得到公式(5),<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>P</mi><mo>=</mo><mo>-</mo><mfrac><mrow><msup><mi>P</mi><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></msup><msub><mi>v</mi><mn>0</mn></msub></mrow><mrow><msup><msub><mi>v</mi><mn>0</mn></msub><mi>T</mi></msup><mi>n</mi></mrow></mfrac><mo>-</mo><mfrac><mrow><msub><mi>d</mi><mn>1</mn></msub><mrow><mo>(</mo><msub><mi>v</mi><mn>1</mn></msub><mo>-</mo><msub><mi>v</mi><mn>0</mn></msub><mo>)</mo></mrow></mrow><mrow><msup><msub><mi>v</mi><mn>1</mn></msub><mi>T</mi></msup><mi>n</mi></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000909374610000022.GIF" wi="686" he="127" /></maths>公式(4、5)中:P为空间一点在摄像机坐标下的三维坐标,P<sup>(3)</sup>为P的第三列,f为摄像机焦距,d<sub>1</sub>为玻璃板厚度,v<sub>0</sub>为靠近镜头端的光路的方向向量,v<sub>1</sub>光路在玻璃板中的方向向量,n为平行玻璃板的法向量;通过公式(4)和公式(5)的连立求解出P<sup>(3)</sup>,进而求解出P。
地址 116024 辽宁省大连市甘井子区凌工路2号