发明名称 基于条纹摄动的快速解相位方法
摘要 本发明公开了一种基于条纹摄动的快速解相位方法包括:首先获取加入条纹摄动的光栅的图像;然后从条纹图像求解出摄动度函数,利用摄动度函数描述的条纹掇动情况,解出相位主值;遍历条纹图,根据摄动度函数和条纹周期数的跳变以及相位主值的关系,得到摄动部分和其余部分的条纹周期数,解出完整的相位值。与现有技术相比,本发明的相移光栅中包含摄动条纹,通过对4幅相移图的处理可以解出完整的相位值,减少了解相位所需的投影光栅数量,提高了测量速度;通过摄动度函数描述摄动情况,进而求解出相位主值,使相位主值的求取有较高的精度;通过摄动度函数判断条纹的周期数,进而解出完整的相位值。
申请公布号 CN100487364C 申请公布日期 2009.05.13
申请号 CN200710192237.4 申请日期 2007.12.21
申请人 东南大学 发明人 达飞鹏;盖绍彦
分类号 G01B11/24(2006.01)I;G06T7/00(2006.01)I 主分类号 G01B11/24(2006.01)I
代理机构 南京经纬专利商标代理有限公司 代理人 陆志斌
主权项 1、一种基于条纹摄动的快速解相位方法,其特征在于:步骤1:获取光栅图像:由投影装置投影出4幅相移光栅条纹,摄像机拍摄条纹图像,获得的图像可以表示为:I<sub>1</sub>(x,y)=I′(x,y)+I″(x,y)cos[θ(x,y)],I<sub>2</sub>(x,y)=I′(x,y)+I″(x,y)cos[θ(x,y)+π/2],I<sub>3</sub>(x,y)=I′(x,y)+I″(x,y)cos[θ(x,y)+π],<maths num="0001"><![CDATA[<math><mrow><msub><mover><mi>I</mi><mo>~</mo></mover><mn>4</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mi>I</mi><mo>&prime;</mo><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>+</mo><mi>&beta;</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mi>I</mi><mo>&prime;</mo><mo>&prime;</mo><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mi>cos</mi><mrow><mo>[</mo><mi>&theta;</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>+</mo><mn>3</mn><mi>&pi;</mi><mo>/</mo><mn>2</mn><mo>]</mo></mrow><mo>,</mo></mrow></math>]]></maths>其中<maths num="0002"><![CDATA[<math><mrow><mi>&beta;</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mrow><mfenced open='{' close='' separators=','><mtable><mtr><mtd><mi>&lambda;</mi></mtd><mtd><mo>&ForAll;</mo><mi>&theta;</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>&Element;</mo><mrow><mo>[</mo><mrow><mo>(</mo><mn>2</mn><mi>Q</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mi>&pi;</mi><mo>,</mo><mrow><mo>(</mo><mn>2</mn><mi>Q</mi><mo>+</mo><mn>2</mn><mo>)</mo></mrow><mi>&pi;</mi><mo>]</mo></mrow></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mo>&ForAll;</mo><mi>&theta;</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>&NotElement;</mo><mrow><mo>[</mo><mrow><mo>(</mo><mn>2</mn><mi>Q</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mi>&pi;</mi><mo>,</mo><mrow><mo>(</mo><mn>2</mn><mi>Q</mi><mo>+</mo><mn>2</mn><mo>)</mo></mrow><mi>&pi;</mi><mo>]</mo></mrow></mtd></mtr></mtable></mfenced></mrow></mrow></math>]]></maths>其中,I<sub>1</sub>、I<sub>2</sub>、I<sub>3</sub>和<img file="C200710192237C00023.GIF" wi="34" he="58" />分别表示4幅光栅条纹图像的灰度值,θ(x,y)为待求的完整的相位值,I′(x,y)为条纹光强的背景值,I″(x,y)为调制强度,Q是位于中心的摄动光栅条纹的周期数;<img file="C200710192237C00024.GIF" wi="33" he="56" />是加入摄动的光栅条纹图像的灰度值,即在投影光栅时,调整投影装置:将[(2Q+1)π,(2Q+2)π]区间内的条纹光强的调制强度设置为原来的λ倍,λ是摄动光栅的常量参数,用于描述<img file="C200710192237C0002143847QIETU.GIF" wi="50" he="86" />摄动程度,取值范围在(0,1)之间;本发明中,定义上式中的β(x,y)为摄动度函数,用于描述光栅条纹的摄动情况;步骤2:求解摄动度函数:将条纹图像代入摄动度函数的求解公式,解出摄动度函数,该摄动度函数的求解公式为:<maths num="0003"><![CDATA[<math><mrow><mi>&beta;</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mrow><mfenced open='{' close=''><mtable><mtr><mtd><mi>&lambda;</mi></mtd><mtd><mo>&ForAll;</mo><mfrac><mrow><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mi>I</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mi>I</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>+</mo><msub><mover><mi>I</mi><mo>~</mo></mover><mn>4</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow><mrow><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mi>I</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>I</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mi>I</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow></mfrac><mo>&lt;</mo><mfrac><mi>&lambda;</mi><mn>2</mn></mfrac></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mo>&ForAll;</mo><mfrac><mrow><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mi>I</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mi>I</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>+</mo><msub><mover><mi>I</mi><mo>~</mo></mover><mn>4</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow><mrow><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mi>I</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>I</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mi>I</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow></mfrac><mo>&GreaterEqual;</mo><mfrac><mi>&lambda;</mi><mn>2</mn></mfrac></mtd></mtr></mtable></mfenced></mrow></mrow></math>]]></maths>其中,β(x,y)为摄动度函数,I<sub>1</sub>、I<sub>2</sub>、I<sub>3</sub>和<img file="C200710192237C00026.GIF" wi="33" he="61" />为光栅条纹图像的灰度值,λ为摄动光栅的常量参数;步骤3:求解相位主值:将条纹图像及上一步解出的β(x,y)代入相位主值的求解公式,解出相位主值,该相位主值的求解公式为:<img file="C200710192237C00031.GIF" wi="1200" he="349" />其中,<img file="C200710192237C00032.GIF" wi="133" he="44" />为相位主值,I<sub>1</sub>、I<sub>2</sub>、I<sub>3</sub>和<img file="C200710192237C00033.GIF" wi="32" he="58" />为光栅条纹图像的灰度值,β(x,y)为摄动度函数,λ为摄动光栅的常量参数;步骤4:解出完整的相位值:遍历图像,若点(x,y)的摄动度函数β(x,y)=λ,则点(x,y)位于摄动光栅条纹上,其条纹周期数k(x,y)=Q;然后根据k(x,y)的跳变和相位主值的关系,逐行遍历条纹图,解出其他点的条纹周期数k(x,y);将条纹周期数k(x,y)和相位主值代入求解公式,解出完整的相位值,该完整的相位值的求解公式为:<img file="C200710192237C00034.GIF" wi="555" he="43" />其中,θ(x,y)为完整的相位值,<img file="C200710192237C00035.GIF" wi="132" he="43" />为相位主值,k(x,y)为条纹周期数。
地址 210096江苏省南京市四牌楼2号