发明名称 整数线性规划搜索法的双频正弦光栅绝对相位解包裹方法
摘要 本发明提出一种整数线性规划搜索法的双频正弦光栅绝对相位解包裹方法,根据传统的双频解包裹原理,提出了关于相对相位、条纹级数的整数线性规划法,建立对应的目标函数和约束条件;以变形的正弦光栅图像为计算对象,采用该方法对其进行遍历搜索,可以快速准确地找到条纹级数的整数对,从而达到解包裹的目的。本发明不要求高频与低频之间成一定的比例关系,扩展了方法的使用范围;不分先后顺序,可以实现对高频与低频的同时解包裹,增强了方法的可操作性;不引入任何中间过度量和划分子区间,简化了计算过程,大大减少了计算量,提高了运算效率。
申请公布号 CN104034285B 申请公布日期 2017.01.04
申请号 CN201410292664.X 申请日期 2014.06.25
申请人 西北工业大学 发明人 莫蓉;卢津;孙惠斌;聂寇准;蒋超峰
分类号 G01B11/25(2006.01)I 主分类号 G01B11/25(2006.01)I
代理机构 西北工业大学专利中心 61204 代理人 陈星
主权项 一种整数线性规划搜索法的双频正弦光栅绝对相位解包裹方法,其特征在于:采用以下步骤:步骤1:设计编码,经过N步相移,生成高频f<sub>high</sub>、低频f<sub>low</sub>互为质数的标准双频正弦光栅条纹图像2N幅,高频正弦光栅条纹图像<img file="FDA0000935924410000011.GIF" wi="226" he="71" />低频正弦光栅条纹图像<img file="FDA0000935924410000012.GIF" wi="190" he="69" />各N幅;(x,y)为图像中水平方向和垂直方向上的像素序列;步骤2:利用投影仪将步骤1中生成的正弦光栅条纹图像<img file="FDA0000935924410000013.GIF" wi="194" he="69" />和<img file="FDA0000935924410000014.GIF" wi="195" he="70" />依次投影到被测物体表面,并用摄像机采集被物体表面调制变形的正弦光栅图像<img file="FDA0000935924410000015.GIF" wi="234" he="71" />和<maths num="0001"><math><![CDATA[<mrow><msubsup><mi>I</mi><mi>n</mi><mrow><mi>C</mi><mo>_</mo><mi>high</mi></mrow></msubsup><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000935924410000016.GIF" wi="267" he="70" /></maths>步骤3:对步骤2中变形正弦光栅图像<img file="FDA0000935924410000017.GIF" wi="235" he="64" />和<img file="FDA0000935924410000018.GIF" wi="252" he="71" />分别进行高斯平滑滤波,并通过建立调制灰度函数,剔除变形的正弦光栅图像中的背景无效像素信息,保留变形正弦光栅图像的有效像素信息;步骤4:将步骤3处理得到的变形正弦光栅图像代入N步相移公式:<img file="FDA0000935924410000019.GIF" wi="902" he="270" />分别计算得到低频f<sub>low</sub>和高频f<sub>high</sub>对应的相对相位值<img file="FDA00009359244100000110.GIF" wi="193" he="62" />和<img file="FDA00009359244100000111.GIF" wi="222" he="70" />N步相移公式中<img file="FDA00009359244100000112.GIF" wi="165" he="69" />根据计算的低频或高频对应取步骤3处理得到的变形正弦光栅图像;步骤5:建立关于频率f<sub>low</sub>、f<sub>high</sub>,相对相位值<img file="FDA00009359244100000113.GIF" wi="429" he="70" />和条纹级数n<sub>low</sub>(x,y)、n<sub>high</sub>(x,y)的目标函数以及约束条件的整数线性规划法则:<maths num="0002"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mi>min</mi></mtd><mtd><mrow><mo>|</mo><msub><mi>n</mi><mrow><mi>h</mi><mi>i</mi><mi>g</mi><mi>h</mi></mrow></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>-</mo><mfrac><msub><mi>f</mi><mrow><mi>h</mi><mi>i</mi><mi>g</mi><mi>h</mi></mrow></msub><msub><mi>f</mi><mrow><mi>l</mi><mi>o</mi><mi>w</mi></mrow></msub></mfrac><msub><mi>n</mi><mrow><mi>l</mi><mi>o</mi><mi>w</mi></mrow></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>-</mo><mfrac><mrow><mi>I</mi><mi>N</mi><mi>t</mi></mrow><msub><mi>f</mi><mrow><mi>l</mi><mi>o</mi><mi>w</mi></mrow></msub></mfrac><mo>|</mo></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA00009359244100000114.GIF" wi="874" he="148" /></maths><img file="FDA00009359244100000115.GIF" wi="1005" he="253" />其中floor()为函数符号,floor(f<sub>low</sub>/2)表示不大于f<sub>low</sub>/2的最大整数,INt为图像矩阵的整数集合;INt采用以下步骤得到:步骤5.1:建立变形的正弦光栅图像整数矩阵:<img file="FDA0000935924410000021.GIF" wi="908" he="142" />其中round表示四舍五入函数;步骤5.2:对整数矩阵INT进行统计:INt=tabulate(INT),其中tabulate表示统计函数,得到整数矩阵INT中出现的元素值;步骤6:采用步骤5中整数线性规划法则对步骤4中的相对相位值图像进行遍历搜索,确定出条纹级数n<sub>low</sub>(x,y)、n<sub>high</sub>(x,y),并计算出绝对相位值Φ<sub>low</sub>(x,y)和Φ<sub>high</sub>(x,y),完成绝对相位展开。
地址 710072 陕西省西安市友谊西路127号