发明名称 一种基于双频编码的快速三维测量方法
摘要 本发明公开了一种基于双频编码的快速三维测量方法,包括以下步骤:搭建三维测量系统;沿投影机纵向或横向编码生成5幅基于双频编码的结构光扫描图案;用照相机抓取所生成的结构光扫描图片,将照相机抓取回来的与<img file="DDA0000464237580000011.GIF" wi="81" he="95" />对应的图片进行处理;将高频相位φ<sub>h</sub>借助低频相位φ<sub>u</sub>进行相位展开,结合展开后的相位计算三维坐标。本发明采用双频策略的编码方式,可以产生同时含有高频与低频的编码图片,因此只需要使用5张图片,即可实现实时的三维重建,从而达到了更快的三维重建效果,并且还可以有效地抑制环境噪声的干扰。
申请公布号 CN103759650A 申请公布日期 2014.04.30
申请号 CN201410045452.1 申请日期 2014.02.07
申请人 成都天拓众成科技有限公司 发明人 刘凯;王帅军;胡子阳
分类号 G01B11/00(2006.01)I 主分类号 G01B11/00(2006.01)I
代理机构 北京天奇智新知识产权代理有限公司 11340 代理人 郭霞
主权项 1.一种基于双频编码的快速三维测量方法,其特征在于:包括以下步骤:(1)搭建三维测量系统:放置标靶于摄像头能够捕捉的范围,投影机射出的光照能够覆盖扫描范围;(2)根据以下公式(1)-(5),沿投影机纵向编码生成5幅基于双频编码的结构光扫描图案,也可沿投影机横向编码生成5幅双频编码策略的结构光扫描图案,则将公式(1)-(4)中的H和y<sup>p</sup>分别改为投影机分辨率的宽度W和对应于投影机中所在的列x<sup>p</sup>:<maths num="0001"><![CDATA[<math><mrow><msubsup><mi>I</mi><mn>0</mn><mi>p</mi></msubsup><mo>=</mo><msup><mi>A</mi><mi>p</mi></msup><mo>+</mo><msubsup><mi>B</mi><mn>1</mn><mi>p</mi></msubsup><mi>sin</mi><mrow><mo>(</mo><mfrac><mrow><mn>2</mn><mi>&pi;</mi><msub><mi>f</mi><mi>h</mi></msub><msup><mi>y</mi><mi>p</mi></msup></mrow><mi>H</mi></mfrac><mo>)</mo></mrow><mo>+</mo><msubsup><mi>B</mi><mn>2</mn><mi>p</mi></msubsup><mi>sin</mi><mrow><mo>(</mo><mfrac><mrow><mn>2</mn><mi>&pi;</mi><msub><mi>f</mi><mi>u</mi></msub><msup><mi>y</mi><mi>p</mi></msup></mrow><mi>H</mi></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><msubsup><mi>I</mi><mn>1</mn><mi>p</mi></msubsup><mo>=</mo><msup><mi>A</mi><mi>p</mi></msup><mo>+</mo><msubsup><mi>B</mi><mn>1</mn><mi>p</mi></msubsup><mi>cos</mi><mrow><mo>(</mo><mfrac><mrow><mn>2</mn><mi>&pi;</mi><msub><mi>f</mi><mi>h</mi></msub><msup><mi>y</mi><mi>p</mi></msup></mrow><mi>H</mi></mfrac><mo>)</mo></mrow><mo>+</mo><msubsup><mi>B</mi><mn>2</mn><mi>p</mi></msubsup><mi>cos</mi><mrow><mo>(</mo><mfrac><mrow><mn>2</mn><mi>&pi;</mi><msub><mi>f</mi><mi>u</mi></msub><msup><mi>y</mi><mi>p</mi></msup></mrow><mi>H</mi></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><msubsup><mi>I</mi><mn>2</mn><mi>p</mi></msubsup><mo>=</mo><msup><mi>A</mi><mi>p</mi></msup><mo>+</mo><msubsup><mi>B</mi><mn>1</mn><mi>p</mi></msubsup><mi>sin</mi><mrow><mo>(</mo><mfrac><mrow><mn>2</mn><mi>&pi;</mi><msub><mi>f</mi><mi>h</mi></msub><msup><mi>y</mi><mi>p</mi></msup></mrow><mi>H</mi></mfrac><mo>)</mo></mrow><mo>+</mo><msubsup><mi>B</mi><mn>2</mn><mi>p</mi></msubsup><mi>sin</mi><mrow><mo>(</mo><mfrac><mrow><mn>2</mn><mi>&pi;</mi><msub><mi>f</mi><mi>u</mi></msub><msup><mi>y</mi><mi>p</mi></msup></mrow><mi>H</mi></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><msubsup><mi>I</mi><mn>3</mn><mi>p</mi></msubsup><mo>=</mo><msup><mi>A</mi><mi>p</mi></msup><mo>+</mo><msubsup><mi>B</mi><mn>1</mn><mi>p</mi></msubsup><mi>cos</mi><mrow><mo>(</mo><mfrac><mrow><mn>2</mn><mi>&pi;</mi><msub><mi>f</mi><mi>h</mi></msub><msup><mi>y</mi><mi>p</mi></msup></mrow><mi>H</mi></mfrac><mo>)</mo></mrow><mo>+</mo><msubsup><mi>B</mi><mn>2</mn><mi>p</mi></msubsup><mi>cos</mi><mrow><mo>(</mo><mfrac><mrow><mn>2</mn><mi>&pi;</mi><msub><mi>f</mi><mi>u</mi></msub><msup><mi>y</mi><mi>p</mi></msup></mrow><mi>H</mi></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><msubsup><mi>I</mi><mn>4</mn><mi>p</mi></msubsup><mo>=</mo><msup><mi>A</mi><mi>p</mi></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths>公式(1)-(5)同时满足以下条件:<maths num="0006"><![CDATA[<math><mrow><msup><mi>A</mi><mi>p</mi></msup><mo>&GreaterEqual;</mo><mrow><mo>(</mo><msubsup><mi>B</mi><mn>1</mn><mi>p</mi></msubsup><mo>+</mo><msubsup><mi>B</mi><mn>2</mn><mi>p</mi></msubsup><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths><maths num="0007"><![CDATA[<math><mrow><msubsup><mi>B</mi><mn>1</mn><mi>p</mi></msubsup><mo>&GreaterEqual;</mo><msubsup><mi>B</mi><mn>2</mn><mi>p</mi></msubsup></mrow></math>]]></maths>且<maths num="0008"><![CDATA[<math><mrow><msup><mi>A</mi><mi>p</mi></msup><mo>+</mo><msubsup><mi>B</mi><mn>1</mn><mi>p</mi></msubsup><mo>+</mo><msubsup><mi>B</mi><mn>1</mn><mi>p</mi></msubsup><mo>=</mo><mn>255</mn></mrow></math>]]></maths>其中,f<sub>h</sub>和f<sub>u</sub>分别表示高频和低频,A<sup>p</sup>表示生成图片的直流分量,<img file="FDA0000464237550000019.GIF" wi="75" he="68" />和<img file="FDA00004642375500000110.GIF" wi="79" he="68" />分别表示对应于高频和低频的谐波调制参数,H表示投影机分辨率的高度,也是生成的结构光扫描图片的预设高度,y<sup>p</sup>对应于投影机中所在的行,<img file="FDA00004642375500000111.GIF" wi="72" he="75" />(i=0,1,2,3,4)表示结构光扫描图案;(3)用照相机抓取公式(1)-(5)所生成的图片,将照相机抓取回来的与<img file="FDA00004642375500000112.GIF" wi="75" he="73" />对应的图片按以下公式(6)-(10)进行处理:<maths num="0009"><![CDATA[<math><mrow><msubsup><mi>B</mi><mn>1</mn><mi>c</mi></msubsup><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msqrt><mfrac><mrow><msubsup><mi>I</mi><mn>0</mn><mi>c</mi></msubsup><mo>+</mo><msubsup><mi>I</mi><mn>2</mn><mi>c</mi></msubsup><mo>-</mo><mn>2</mn><msubsup><mi>I</mi><mn>4</mn><mi>c</mi></msubsup></mrow><mrow><msubsup><mi>I</mi><mn>1</mn><mi>c</mi></msubsup><mo>+</mo><msubsup><mi>I</mi><mn>3</mn><mi>c</mi></msubsup><mo>-</mo><mn>2</mn><msubsup><mi>I</mi><mn>4</mn><mi>c</mi></msubsup></mrow></mfrac></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0010"><![CDATA[<math><mrow><msubsup><mi>B</mi><mn>2</mn><mi>c</mi></msubsup><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msqrt><mfrac><mrow><msubsup><mi>I</mi><mn>0</mn><mi>c</mi></msubsup><mo>-</mo><msubsup><mi>I</mi><mn>2</mn><mi>c</mi></msubsup></mrow><mrow><msubsup><mi>I</mi><mn>1</mn><mi>c</mi></msubsup><mo>-</mo><msubsup><mi>I</mi><mn>3</mn><mi>c</mi></msubsup></mrow></mfrac></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0011"><![CDATA[<math><mrow><msup><mi>A</mi><mi>c</mi></msup><mo>=</mo><msubsup><mi>I</mi><mn>4</mn><mi>c</mi></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow></math>]]></maths><img file="FDA0000464237550000024.GIF" wi="1677" he="207" /><img file="FDA0000464237550000025.GIF" wi="1691" he="196" />其中,<img file="FDA0000464237550000027.GIF" wi="61" he="70" />(i=0,1,2,3,4)表示照相机抓取的对应于<img file="FDA0000464237550000028.GIF" wi="66" he="71" />编码的图片,A<sup>c</sup>对应于A<sup>p</sup>,表示抓取到的图片的直流分量,<img file="FDA00004642375500000211.GIF" wi="67" he="66" />与<img file="FDA00004642375500000212.GIF" wi="58" he="64" />分别对应于<img file="FDA0000464237550000029.GIF" wi="66" he="62" />和<img file="FDA00004642375500000210.GIF" wi="95" he="67" />分别表示抓取到的图片的高频和低频的调制参数,φ<sub>h</sub>和φ<sub>u</sub>分别表示高频和低频对应的相位信息;(4)根据公式(9)和(10),将高频相位φ<sub>h</sub>借助低频相位φ<sub>u</sub>进行相位展开,结合展开后的相位计算三维坐标。
地址 610041 四川省成都市高新区肖家河中街43号6幢1层