发明名称 实时光学三维测量中动态场景无效点的自动识别与排除方法
摘要 本发明提出一种实时光学三维测量中动态场景无效点的自动识别与排除方法。本发明利用N步相移法和双频时间相位法求出的相移图、光栅条纹、相位调制度图以及高斯滤波器来侦察无效点。本发明方法高效地排除了由于物体运动和不稳定的测量环境所产生了无效点,显著提高了测量结果的精度。
申请公布号 CN103697834A 申请公布日期 2014.04.02
申请号 CN201310732419.1 申请日期 2013.12.26
申请人 南京理工大学 发明人 陈钱;冯世杰;顾国华;左超;孙佳嵩;喻士领;申国辰;李如斌
分类号 G01B11/25(2006.01)I 主分类号 G01B11/25(2006.01)I
代理机构 南京理工大学专利中心 32203 代理人 朱显国
主权项 1.实时光学三维测量中动态场景无效点的自动识别与排除方法,其特征在于,包括以下步骤:步骤一:利用投影仪对待测物表面各投影一组高频正弦光栅条纹和一组低频正弦光栅条纹,两组正弦光栅条纹均包含N幅相移为2πn/N的正弦光栅条纹,使用摄像机采集包含了待测物表面信息的高频正弦光栅条纹和低频正弦光栅条纹;分别在高频正弦光栅条纹和低频正弦光栅条纹里使用N步相移算法计算获得待测物表面每个被测点对应的像素点被包裹的高频相位值φ(x,y)和低频相位值φ′(x,y),然后使用双频时间相位去包裹算法计算获得去包裹以后的相位值Φ<sub>1</sub>(x,y);同时计算待测物表面同一被测点对应像素点的高频正弦光栅条纹所包含的相位调制度值B<sub>High</sub>(x,y)和低频正弦光栅条纹所包含的相位调制度值B<sub>Low</sub>(x,y);其中(x,y)为摄像机图像像素点坐标;步骤二:将每个像素点去包裹以后的相位值Φ<sub>1</sub>(x,y)代入公式(1),判断公式(1)是否成立,若公式(1)不成立,则对应像素点为无效点,予以排除;若公式(1)成立,则对应像素点予以保留,<![CDATA[<math><mrow><mi>E</mi><mo>></mo><msqrt><mfrac><mrow><msubsup><mi>&Sigma;</mi><mi>n</mi><mi>N</mi></msubsup><msup><mrow><mo>(</mo><msubsup><mi>K</mi><mi>n</mi><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msubsup><mo>-</mo><msubsup><mi>K</mi><mi>n</mi><mo>&prime;</mo></msubsup><mo>)</mo></mrow><mn>2</mn></msup></mrow><mi>N</mi></mfrac></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>式(1)中,n取值依次为1、2、3…N,N为相移算法的总相移步数,E为公式(6)的判断阈值,取值范围在0至1之间;式(1)中,中间量K<sub>n</sub>′=cos[Φ<sub>1</sub>(x,y)+2πn/N];式(1)中,中间量<img file="FDA0000447598660000012.GIF" wi="385" he="141" />其中,(x,y)为摄像机图像像素坐标,A(x,y)为摄像机采集图像的平均光强,B(x,y)为摄像机采集图像的相位调制度值,且<![CDATA[<math><mrow><mi>B</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>2</mn><mi>N</mi></mfrac><msqrt><msup><mrow><mo>[</mo><msubsup><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></msubsup><msubsup><mi>I</mi><mi>n</mi><mi>c</mi></msubsup><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mi>sin</mi><mrow><mo>(</mo><mn>2</mn><mi>&pi;n</mi><mo>/</mo><mi>N</mi><mo>)</mo></mrow><mo>]</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>[</mo><msubsup><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></msubsup><msubsup><mi>I</mi><mi>n</mi><mi>c</mi></msubsup><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mi>cos</mi><mrow><mo>(</mo><mn>2</mn><mi>&pi;n</mi><mo>/</mo><mi>N</mi><mo>)</mo></mrow><mo>]</mo></mrow><mn>2</mn></msup></msqrt><mo>;</mo></mrow></math>]]></maths>步骤三:将每个像素点的高频正弦光栅条纹所包含的相位调制度值B<sub>High</sub>(x,y)和低频正弦光栅条纹所包含的相位调制度值B<sub>Low</sub>(x,y)代入公式(2),判断公式(2)是否成立,若公式(2)不成立,则对应像素点为无效点,予以排除;若公式(2)成立,则对应像素点予以保留,<![CDATA[<math><mrow><mfrac><mrow><mi>abs</mi><mo>[</mo><msub><mi>B</mi><mi>High</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>B</mi><mi>Low</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>]</mo></mrow><mrow><mn>0.5</mn><mo>&times;</mo><mo>[</mo><msub><mi>B</mi><mi>High</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>B</mi><mi>Low</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>]</mo></mrow></mfrac><mo>&lt;</mo><mi>&sigma;</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>式(2)中,abs为求取绝对值运算符,σ为公式(2)的判断阈值,取值范围在0至1之间;步骤四:将每一个像素点(x,y)以及其右侧像素点(x+1,y)的相位值Φ<sub>1</sub>(x,y)和Φ<sub>1</sub>(x+1,y)代入公式(3),判断公式(3)是否成立,若公式(3)不成立,则对应像素点为无效点,予以排除;若公式(3)成立,则对应像素点予以保留,Δ<sub>1</sub><Φ<sub>1</sub>(x+1,y)-Φ<sub>1</sub>(x,y)<Δ<sub>2</sub>    (3)式(3)中,Δ<sub>1</sub>和Δ<sub>2</sub>均为公式(3)的判断阈值,Δ<sub>1</sub>与Δ<sub>2</sub>的取值范围在-π至π之间;步骤五:用高斯滤波器对每一个像素点的相位值Φ<sub>1</sub>(x,y)做高斯滤波获得滤波后的相位值<img file="FDA0000447598660000021.GIF" wi="186" he="80" />将每一个像素点的相位值Φ<sub>1</sub>(x,y)和滤波后相位值<img file="FDA0000447598660000022.GIF" wi="148" he="79" />代入公式(4),判断公式(4)是否成立,若公式(4)不成立,则对应像素点为无效点,予以排除;若公式(4)成立,则对应像素点予以保留,<![CDATA[<math><mrow><mi>abs</mi><mrow><mo>(</mo><msub><mi>&Phi;</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>-</mo><mover><mi>&Phi;</mi><mo>~</mo></mover><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>&lt;</mo><mi>V</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths>式(4)中,V公式(4)的判断阈值,取值范围在0至1之间。
地址 210094 江苏省南京市孝陵卫200号