发明名称 基于相位编码的混频结构光三维信息获取方法
摘要 本发明公开了一种基于相位编码的混频结构光三维信息获取方法,主要解决现有三维信息获取方法测量精度低、空间分辨率低、耗时长的问题。其实现步骤为:设计包含两种频率信息且强度变化的双颜色条纹模板;用投影仪将其投射到待测物体上,用摄像机记录变形条纹图像;求解变形条纹图像的颜色相位分布和强度相位分布;根据颜色相位分布和强度相位分布,计算变形条纹图像中像素点的截断相位展开值;根据截断相位展开值确定像素点在投影模板中的匹配点;用三角测距原理及匹配点坐标,求出变形条纹图像中每个像素点的三维坐标值。本发明具有空间分辨率高、测量精度高、测量速度快的优点,可用于对动态物体的三维信息获取。
申请公布号 CN103697815B 申请公布日期 2017.03.01
申请号 CN201410018554.4 申请日期 2014.01.15
申请人 西安电子科技大学 发明人 石光明;高山;李甫;石悦鑫;李芹;杨莉莉;李若岱
分类号 G01B11/00(2006.01)I 主分类号 G01B11/00(2006.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华;朱红星
主权项 一种基于相位编码的混频结构光三维信息获取方法,包括如下步骤:(1)设计包含两种频率信息且强度变化的双颜色条纹模板P:(1a)设定双颜色条纹模版P的颜色为红色和蓝色,设强度变化的周期为T<sub>i</sub>,颜色变化的周期为T<sub>c</sub>,且T<sub>i</sub>和T<sub>c</sub>为互质的整数;(1b)将双颜色条纹模板P中第x行、第y列的像素点(x,y)处的红色分量灰度值标记为P<sub>r</sub>(x,y),蓝色分量灰度值标记为P<sub>b</sub>(x,y),并按下式进行赋值,构成双颜色条纹模板P:<maths num="0001"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><msub><mi>P</mi><mi>r</mi></msub><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo><mo>=</mo><mfrac><mrow><mi>C</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>+</mo><mn>1</mn></mrow><mn>2</mn></mfrac><mi>P</mi><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mtd></mtr><mtr><mtd><msub><mi>P</mi><mi>b</mi></msub><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo><mo>=</mo><mo>-</mo><mfrac><mrow><mi>C</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>-</mo><mn>1</mn></mrow><mn>2</mn></mfrac><mi>P</mi><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0001045627710000011.GIF" wi="602" he="238" /></maths>其中,x、y分别表示像素点(x,y)的行、列坐标,P(x,y)=a+b cos(2πf<sub>i</sub>y)a、b为常量,分别表示模板的背景深度和调制强度,<img file="FDA0001045627710000012.GIF" wi="148" he="110" />C(x,y)=sign(l(x,y)‑l(x,y‑1)),l(x,y)=cos(2πf<sub>2</sub>y),<img file="FDA0001045627710000013.GIF" wi="148" he="110" />sign(·)为取正负操作,用1表示正,‑1表示负;(2)将投影仪T与摄像机V水平放置,并使两者光轴平行,再用投影仪T将双颜色条纹模板P投射到三维物体O上,用摄像机V拍摄经待测物体O调制过的变形条纹图像U;(3)将变形条纹图像U中灰度值大于10的像素点,标记为有效像素点(i,j),计算有效像素点(i,j)的颜色编码信息C(i,j)及强度编码信息I(i,j):<maths num="0002"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mi>C</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mrow><msub><mi>P</mi><mi>r</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>&GreaterEqual;</mo><msub><mi>P</mi><mi>b</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mo>-</mo><mn>1</mn></mrow></mtd><mtd><mrow><msub><mi>P</mi><mi>b</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>&gt;</mo><msub><mi>P</mi><mi>r</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced></mrow></mtd></mtr><mtr><mtd><mrow><mi>I</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mi>max</mi><mrow><mo>(</mo><msub><mi>P</mi><mi>r</mi></msub><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>P</mi><mi>b</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>)</mo></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0001045627710000014.GIF" wi="693" he="242" /></maths>其中,i、j分别表示有效像素点(i,j)的行、列坐标,P<sub>r</sub>(i,j),P<sub>b</sub>(i,j)分别为有效像素点(i,j)的红色分量和蓝色分量的灰度值,C(i,j)为1表示红色,C(i,j)为‑1表示蓝色;(4)计算变形条纹图像U中第i行、第j列的有效像素点(i,j)的颜色积分值:<maths num="0003"><math><![CDATA[<mrow><mi>T</mi><mi>r</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mi>&Sigma;</mi><mrow><msup><mi>j</mi><mo>&prime;</mo></msup><mo>=</mo><mi>s</mi><mi>t</mi><mi>a</mi><mi>r</mi><mi>t</mi></mrow><mrow><msup><mi>j</mi><mo>&prime;</mo></msup><mo>=</mo><mi>j</mi></mrow></msubsup><mi>C</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><msup><mi>j</mi><mo>&prime;</mo></msup><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001045627710000015.GIF" wi="492" he="92" /></maths>其中,start为第i行第一个有效像素点的列坐值;(5)将颜色积分值Tr(i,j)与波长为T<sub>c</sub>的一维Gabor滤波器进行卷积,求取像素点(i,j)的颜色相位分布<img file="FDA0001045627710000021.GIF" wi="169" he="63" />(6)将强度编码信息I(i,j)与波长为T<sub>i</sub>的一维Gabor滤波器进行卷积,求取像素点(i,j)的强度相位分布<img file="FDA0001045627710000022.GIF" wi="163" he="62" />(7)假定在深度为z<sub>R</sub>处存在虚拟参考平面R,根据空间几何关系,计算像素点(i,j)对应于虚拟参考平面R的参考颜色相位分布<img file="FDA0001045627710000023.GIF" wi="145" he="63" />及参考强度相位分布<img file="FDA0001045627710000024.GIF" wi="163" he="63" /><img file="FDA0001045627710000025.GIF" wi="726" he="303" />其中,<img file="FDA0001045627710000026.GIF" wi="1092" he="127" />f为摄像机的焦距,b为摄像机V光心和投影仪T光心的水平距离,θ<sub>V</sub>为摄像机的水平视场角,θ<sub>T</sub>为投影仪的水平视场角,Δd<sub>V</sub>为摄像机单个像素所代表的实际宽度,其值为<img file="FDA0001045627710000027.GIF" wi="417" he="135" />N<sub>V</sub>为图像U的总列数,Δd<sub>T</sub>为投影仪模版中单个像素所代表的实际宽度,其值为<img file="FDA0001045627710000028.GIF" wi="443" he="126" />N<sub>T</sub>为模版P的总列数,<img file="FDA0001045627710000029.GIF" wi="618" he="130" />T<sub>c</sub>为双颜色条纹模板P的颜色变化周期;(8)计算像素点(i,j)的强度截断相位差<img file="FDA00010456277100000210.GIF" wi="186" he="62" />和颜色截断相位差<img file="FDA00010456277100000211.GIF" wi="213" he="62" /><img file="FDA00010456277100000212.GIF" wi="628" he="156" />(9)根据空间几何关系及余数定理,利用像素点的强度截断相位差<img file="FDA00010456277100000213.GIF" wi="186" he="59" />和颜色截断相位差<img file="FDA00010456277100000214.GIF" wi="216" he="63" />求取像素点(i,j)的截断相位展开值Δφ(i,j);(10)利用像素点(i,j)的截断相位展开值Δφ(i,j),求取像素点(i,j)在投影模版P中的匹配点(x(i,j),y(i,j)),其中x(i,j)和y(i,j)分别表示匹配点在投影模板P中的行、列坐标值;(11)根据三角测距原理,利用像素点(i,j)与匹配点(x(i,j),y(i,j))的空间关系,计算待测物体的三维信息值。
地址 710071 陕西省西安市太白南路2号