发明名称 一种应用于多频率三维测量的全局相位解相方法
摘要 本发明属于三维机器视觉领域,涉及一种应用于多频率三维测量的全局相位解相方法。该方法通过投射3种周期分别为P1、P2和P3的正弦或者余弦规律变化的光信息,每种光信息需要经过4-8步相移,通过每个周期的相移图案,分别解算每个周期的相位值、两个周期的合成相位值以及三个周期的最终合成相位值。最后通过合成相位分别计算每个周期的全局相位值,进而通过三维重建获得物体的全貌三维坐标。本发明所设计的全局相位解相方法,优于传统的格雷码相位解相方法,可以有效的解决三维测量中被测物体表面颜色不一致的测量难题,无需喷涂显影剂,更加绿色环保。
申请公布号 CN102589478A 申请公布日期 2012.07.18
申请号 CN201210048001.4 申请日期 2012.02.29
申请人 天津工业大学 发明人 宋丽梅;杨燕罡;董虓霄;张亮;陈昌曼
分类号 G01B11/25(2006.01)I 主分类号 G01B11/25(2006.01)I
代理机构 代理人
主权项 1.本发明所设计的全局相位解相方法,其特征是:对于图像中每一个点(x,y)的全局相位值的计算过程,工作步骤如下:步骤1:选取合适的P<sub>1</sub>、P<sub>2</sub>和P<sub>3</sub>的值,P<sub>1</sub>、P<sub>2</sub>和P<sub>3</sub>的取值均在0-L<sub>R</sub>之间,并且P<sub>1</sub>和P<sub>2</sub>的合成周期P<sub>12</sub>以及P<sub>2</sub>和P<sub>3</sub>的合成周期P<sub>23</sub>,以及P<sub>12</sub>和P<sub>23</sub>的合成周期P<sub>123</sub>,经过公式(1)的运算后,满足P<sub>123</sub>≥L<sub>R</sub>;<maths num="0001"><![CDATA[<math><mrow><msub><mi>P</mi><mn>12</mn></msub><mo>=</mo><mo>|</mo><mfrac><mrow><msub><mi>P</mi><mn>1</mn></msub><mo>&times;</mo><msub><mi>P</mi><mn>2</mn></msub></mrow><mrow><msub><mi>P</mi><mn>1</mn></msub><mo>-</mo><msub><mi>P</mi><mn>2</mn></msub></mrow></mfrac><mo>|</mo></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><msub><mi>P</mi><mn>23</mn></msub><mo>=</mo><mo>|</mo><mfrac><mrow><msub><mi>P</mi><mn>2</mn></msub><mo>&times;</mo><msub><mi>P</mi><mn>3</mn></msub></mrow><mrow><msub><mi>P</mi><mn>2</mn></msub><mo>-</mo><msub><mi>P</mi><mn>3</mn></msub></mrow></mfrac><mo>|</mo></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><msub><mi>P</mi><mn>123</mn></msub><mo>=</mo><mo>|</mo><mfrac><mrow><msub><mi>P</mi><mn>12</mn></msub><mo>&times;</mo><msub><mi>P</mi><mn>23</mn></msub></mrow><mrow><msub><mi>P</mi><mn>12</mn></msub><mo>-</mo><msub><mi>P</mi><mn>23</mn></msub></mrow></mfrac><mo>|</mo></mrow></math>]]></maths>公式(1)步骤2:利用光源投射装置,向物体投射周期为P<sub>1</sub>的相移图像,相移图像的步数应大于3步,最好为4-8步,采集每幅图像,并计算(x,y)在P<sub>1</sub>周期的相位<img file="FSA00000676370600014.GIF" wi="185" he="57" />步骤3:利用光源投射装置,向物体投射周期为P<sub>2</sub>的相移图像,相移图像的步数应大于3步,最好为4-8步,采集每幅图像,并计算(x,y)在P<sub>2</sub>周期的相位<img file="FSA00000676370600015.GIF" wi="188" he="57" />步骤4:利用光源投射装置,向物体投射周期为P<sub>3</sub>的相移图像,相移图像的步数应大于3步,最好为4-8步,采集每幅图像,并计算(x,y)在P<sub>3</sub>周期的相位<img file="FSA00000676370600016.GIF" wi="188" he="56" />步骤5:利用(x,y)点的<img file="FSA00000676370600017.GIF" wi="162" he="57" />和<img file="FSA00000676370600018.GIF" wi="189" he="58" />计算在P<sub>12</sub>周期的合成相位<img file="FSA00000676370600019.GIF" wi="196" he="57" />步骤6:利用(x,y)点的<img file="FSA000006763706000110.GIF" wi="167" he="58" />和<img file="FSA000006763706000111.GIF" wi="188" he="57" />计算在P<sub>23</sub>周期的合成相位<img file="FSA000006763706000112.GIF" wi="200" he="57" />步骤7:利用(x,y)点的<img file="FSA000006763706000113.GIF" wi="174" he="56" />和<img file="FSA000006763706000114.GIF" wi="201" he="57" />计算在P<sub>123</sub>周期的最终合成相位<img file="FSA000006763706000115.GIF" wi="208" he="57" />步骤8:对于在P<sub>1</sub>周期的每个点(x,y),利用其原始相位<img file="FSA000006763706000116.GIF" wi="215" he="57" />按公式(2)计算其全局相位<img file="FSA000006763706000117.GIF" wi="224" he="57" /><maths num="0004"><![CDATA[<math><mrow><msub><mi>&theta;</mi><mrow><mi>G</mi><mo>-</mo><msub><mi>P</mi><mn>1</mn></msub></mrow></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>&theta;</mi><mrow><mi>i</mi><mo>-</mo><msub><mi>P</mi><mn>1</mn></msub></mrow></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>+</mo><mn>2</mn><mi>&pi;</mi><mo>&times;</mo><mi>M</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow></math>]]></maths>公式(2)其中:<img file="FSA000006763706000119.GIF" wi="212" he="58" />为(x,y)点的原始相位;<img file="FSA000006763706000120.GIF" wi="225" he="58" />为(x,y)点的全局相位;M(x,y):为(x,y)点所在的相位周期在全局相位中的编码值,该值的计算如公式(3)所示;<maths num="0005"><![CDATA[<math><mrow><mi>M</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><mi>int</mi><mo>)</mo></mrow><mrow><mo>(</mo><mfrac><mrow><msub><mi>&theta;</mi><msub><mi>P</mi><mn>123</mn></msub></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow><mrow><mn>2</mn><mi>&pi;</mi></mrow></mfrac><mo>)</mo></mrow><mo>&times;</mo><msup><mrow><mo>(</mo><mo>|</mo><mrow><mo>(</mo><mfrac><mrow><mo>|</mo><mfrac><mrow><msub><mi>P</mi><mn>1</mn></msub><mo>&times;</mo><msub><mi>P</mi><mn>2</mn></msub></mrow><mrow><msub><mi>P</mi><mn>1</mn></msub><mo>-</mo><msub><mi>P</mi><mn>2</mn></msub></mrow></mfrac><mo>|</mo><mo>&times;</mo><mo>|</mo><mfrac><mrow><msub><mi>P</mi><mn>2</mn></msub><mo>&times;</mo><msub><mi>P</mi><mn>3</mn></msub></mrow><mrow><msub><mi>P</mi><mn>2</mn></msub><mo>-</mo><msub><mi>P</mi><mn>3</mn></msub></mrow></mfrac><mo>|</mo></mrow><mrow><mo>|</mo><mfrac><mrow><msub><mi>P</mi><mn>1</mn></msub><mo>&times;</mo><msub><mi>P</mi><mn>2</mn></msub></mrow><mrow><msub><mi>P</mi><mn>1</mn></msub><mo>-</mo><msub><mi>P</mi><mn>2</mn></msub></mrow></mfrac><mo>|</mo><mo>-</mo><mo>|</mo><mfrac><mrow><msub><mi>P</mi><mn>2</mn></msub><mo>&times;</mo><msub><mi>P</mi><mn>3</mn></msub></mrow><mrow><msub><mi>P</mi><mn>2</mn></msub><mo>-</mo><msub><mi>P</mi><mn>3</mn></msub></mrow></mfrac><mo>|</mo></mrow></mfrac><mo>)</mo></mrow><mo>&times;</mo><mfrac><mrow><msub><mi>P</mi><mn>1</mn></msub><mo>-</mo><msub><mi>P</mi><mn>2</mn></msub></mrow><mrow><msub><mi>P</mi><mn>1</mn></msub><mo>&times;</mo><msub><mi>P</mi><mn>2</mn></msub></mrow></mfrac><mo>|</mo><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mrow><mo>(</mo><mi>int</mi><mo>)</mo></mrow><mrow><mo>(</mo><mfrac><mrow><msub><mi>&theta;</mi><msub><mi>P</mi><mn>12</mn></msub></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow><mrow><mn>2</mn><mi>&pi;</mi></mrow></mfrac><mo>)</mo></mrow><mo>&times;</mo><mrow><mo>(</mo><mo>|</mo><mfrac><mrow><msub><mi>P</mi><mn>1</mn></msub><mo>&times;</mo><msub><mi>P</mi><mn>2</mn></msub></mrow><mrow><mrow><mo>(</mo><msub><mi>P</mi><mn>1</mn></msub><mo>-</mo><msub><mi>P</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>&times;</mo><msub><mi>P</mi><mn>1</mn></msub></mrow></mfrac><mo>|</mo><mo>)</mo></mrow></mrow></math>]]></maths>公式(3)(x,y)点的全局相位<img file="FSA000006763706000122.GIF" wi="201" he="57" />运算完毕。
地址 300160 天津市河东区成林道63号