发明名称 一种基于质量图区域分割的快速相位展开方法
摘要 本发明属于计算机测量技术领域,具体涉及一种基于质量图区域分割的快速相位展开方法。本发明通过摄像机捕获处理得到包裹相位图;根据包裹相位图计算质量图;使用相位微分差分算法对质量图进行滤波;自动获得区域分割阈值;利用得到的阈值将包裹相位图分成优先级不同的区域;采用简单路径跟踪相位展开算法完成区域的相位展开;进行相位展开;按照优先级从高到低的顺序,将相位展开后的组合并。本发明基于区域分割的简单路径跟踪相位展开算法,能够在尽量不降低相位展开精度的前提下提高相位展开的速度。该方法满足结构光视觉三维测量的高精度,快速解缠,抗噪声干扰等要求。
申请公布号 CN103913128A 申请公布日期 2014.07.09
申请号 CN201410087942.8 申请日期 2014.03.12
申请人 哈尔滨工程大学 发明人 陆军;苏航;董东来;李积江;张鑫;宋景豪;彭仲涛;刘伟;华威;宁召可
分类号 G01B11/25(2006.01)I 主分类号 G01B11/25(2006.01)I
代理机构 代理人
主权项 一种基于质量图区域分割的快速相位展开方法,其特征在于:步骤1:使用投影仪将正弦光栅投影到目标物体上,通过摄像机捕获处理得到包裹相位图;步骤2:根据包裹相位图计算质量图,使用相位微分差分算法对质量图进行滤波,通过质量图来评价包裹相位精度的高低:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>q</mi><mi>min</mi></msub><mo>=</mo><mfrac><mrow><msqrt><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mi>m</mi><mo>-</mo><mfrac><mi>k</mi><mn>2</mn></mfrac></mrow><mrow><mi>m</mi><mo>+</mo><mfrac><mi>k</mi><mn>2</mn></mfrac></mrow></msubsup><msubsup><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mi>n</mi><mo>-</mo><mfrac><mi>k</mi><mn>2</mn></mfrac></mrow><mrow><mi>n</mi><mo>+</mo><mfrac><mi>k</mi><mn>2</mn></mfrac></mrow></msubsup><msup><mrow><mo>(</mo><msubsup><mi>&Delta;</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>x</mi></msubsup><mo>-</mo><msubsup><mover><mi>&Delta;</mi><mo>&OverBar;</mo></mover><mrow><mi>m</mi><mo>,</mo><mi>n</mi></mrow><mi>x</mi></msubsup><mo>)</mo></mrow><mn>2</mn></msup></msqrt><mo>+</mo><msqrt><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mi>m</mi><mo>-</mo><mfrac><mi>k</mi><mn>2</mn></mfrac></mrow><mrow><mi>m</mi><mo>+</mo><mfrac><mi>k</mi><mn>2</mn></mfrac></mrow></msubsup><msubsup><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mi>n</mi><mo>-</mo><mfrac><mi>k</mi><mn>2</mn></mfrac></mrow><mrow><mi>n</mi><mo>+</mo><mfrac><mi>k</mi><mn>2</mn></mfrac></mrow></msubsup><msup><mrow><mo>(</mo><msubsup><mi>&Delta;</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mi>y</mi></msubsup><mo>-</mo><msubsup><mover><mi>&Delta;</mi><mo>&OverBar;</mo></mover><mrow><mi>m</mi><mo>,</mo><mi>n</mi></mrow><mi>y</mi></msubsup><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mrow><mrow><mi>k</mi><mo>&times;</mo><mi>k</mi></mrow></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000475715120000011.GIF" wi="1178" he="185" /></maths>其中,<img file="FDA0000475715120000012.GIF" wi="73" he="80" />和<img file="FDA0000475715120000013.GIF" wi="67" he="92" />是包裹相位在x和y方向的差分,<img file="FDA0000475715120000014.GIF" wi="106" he="79" />和<img file="FDA0000475715120000015.GIF" wi="96" he="86" />是k×k窗口中x和y方向差分的平均值,步骤3:根据获得的质量图信息,按照质量值的由低到高将质量图分区域,根据质量值直方图的分布模式,自动获得区域分割阈值:对图像质量图I(x,y),高质量区域和低质量区域的分割阈值记作T,属于高质量区域的像素点数占整幅图像的比例记为ω<sub>0</sub>,其平均灰度μ<sub>0</sub>;低质量区域像素点数占整幅图像的比例为ω<sub>1</sub>,其平均灰度为μ<sub>1</sub>,图像的总平均灰度记为μ类间方差记为g,图像的大小为M*N,图像中像素的灰度值小于阈值T的像素个数记作N<sub>0</sub>,像素灰度大于阈值T的像素个数记作N<sub>1</sub>,<img file="FDA0000475715120000016.GIF" wi="482" he="103" />N<sub>0</sub>+N<sub>1</sub>=M*N,ω<sub>0</sub>+ω<sub>1</sub>=1,μ=ωx*μ<sub>0</sub>+ω<sub>1</sub>*μ<sub>1</sub>,g=ω<sub>0</sub>(μ<sub>0</sub>‑μ)<sup>2</sup>+ω<sub>1</sub>(μ<sub>1</sub>‑μ)<sup>2</sup>,g=ω<sub>0</sub>ω<sub>1</sub>(μ<sub>0</sub>‑μ<sub>1</sub>)<sup>2</sup>,采用遍历的方法得到使类间方差最大的阈值T,即为所求阈值;步骤4:利用得到的阈值将包裹相位图分成优先级不同的区域,优先级高的区域优先扩展,优先级低的区域随后扩展,按照展开的先后次序将得到的区域赋予组编号,在同一区域的组具有相同的优先级;步骤5:采用简单路径跟踪相位展开算法完成区域的相位展开:分别在不同优先级区域中选择一个起始点,将起始点按其左、上、右、下的次序扩展得到4个节点,如果这4个节点在本区域,将得到的4个节点依次放入open表,将起始点放入close表中,并令其绝对相位为零;步骤6:从open表中取出第一个节点,进行相位展开,<img file="FDA0000475715120000017.GIF" wi="1172" he="234" />式中φ<sub>A</sub>(r,c)为待展开处(r,c)的绝对相位,φ<sub>A</sub>(r′,c′)为(r′,c′)处已展开的绝对相位,(r′,c′)为(r,c)邻域中的一个点,且其在close表中,θ为阈值,此处取θ=π,Δφ=φ<sub>W</sub>(r,c)‑φ<sub>W</sub>(r′,c′),φ<sub>W</sub>(r,c)为(r,c)处的包裹相位,φ<sub>W</sub>(r′,c′)为(r′,c′)处的包裹相位;将展开的节点放入close表中并按其左、上、右、下的次序扩展得到4个节点,如果这4个节点在本区域中且不在close表中和open表中,将得到的这4个节点依次放入open表中;重复步骤6直到open表为空,完成区域的相位展开;步骤7:按照优先级从高到低的顺序,将相位展开后的组合并:将被低优先级组包围的孤立高优先级组同周围已经完成相位展开合并的低优先级组进行合并,得到最终绝对相位图。
地址 150001 黑龙江省哈尔滨市南岗区南通大街145号哈尔滨工程大学科技处知识产权办公室