发明名称 一种降低摄像机多路径效应影响的相位融合方法
摘要 本发明公开了一种降低摄像机多路径效应影响的相位融合方法,包括以下步骤:搭建三维测量系统,将标靶置于摄像头能够捕捉的范围和投影机能够投影的范围;对标靶分别进行沿X方向和沿Y方向的扫描并计算得到两组相位值φ<sub>x</sub>和φ<sub>y</sub>;将X方向的相位值转换到Y方向,转换后的相位为φ<sub>xy</sub>;通过分别投射一副全黑和全白的图片,然后计算出标靶表面的反射率,并依据该反射率计算出标靶表面每一像素点反射率的变化方向,记为θ;利用相位融合算法求出新的相位值φ;对该新生成的相位值φ进行三维重建。本发明仅仅通过增加不同算法对扫描得到的相位进行修复,再通过修复后的相位值对标靶进行三维重建,得到更加精确的三维信息,尤其是对于高对比度区域,能显著地减小测量误差。
申请公布号 CN103837094B 申请公布日期 2016.06.15
申请号 CN201410086658.9 申请日期 2014.03.11
申请人 四川大学 发明人 刘凯;许晓畅
分类号 G01B11/24(2006.01)I 主分类号 G01B11/24(2006.01)I
代理机构 北京天奇智新知识产权代理有限公司 11340 代理人 杨春
主权项 一种降低摄像机多路径效应影响的相位融合方法,其特征在于:包括以下步骤:(1)搭建三维测量系统,将标靶置于摄像头能够捕捉的范围和投影机能够投影的范围;如果三维测量系统已校准并已有摄像机的参数矩阵M<sub>wc</sub>和投影仪的参数矩阵M<sub>wp</sub>,则进入步骤(2),如果三维测量系统未校准,则先调整投影仪的焦距,再调整摄像头的焦距、光圈大小以及曝光时间的参数,然后对三维测量系统进行校准,求得摄像机的参数矩阵M<sub>wc</sub>和投影仪的参数矩阵M<sub>wp</sub>,之后再进入步骤(2);(2)对标靶分别进行沿X方向和沿Y方向的扫描并计算得到两组相位值,其中沿X方向的相位值为φ<sub>x</sub>,沿Y方向的相位值为φ<sub>y</sub>;(3)利用如下公式(1)将X方向的相位值转换到Y方向,转换后的相位为φ<sub>xy</sub>:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>&phi;</mi><mrow><mi>x</mi><mi>y</mi></mrow></msub><mo>=</mo><mfrac><mrow><mo>&lsqb;</mo><mtable><mtr><mtd><msubsup><mi>m</mi><mn>21</mn><mrow><mi>w</mi><mi>p</mi></mrow></msubsup></mtd><mtd><msubsup><mi>m</mi><mn>22</mn><mrow><mi>w</mi><mi>p</mi></mrow></msubsup></mtd><mtd><msubsup><mi>m</mi><mn>23</mn><mrow><mi>w</mi><mi>p</mi></mrow></msubsup></mtd></mtr></mtable><mo>&rsqb;</mo><msup><msub><mi>C</mi><mi>x</mi></msub><mrow><mo>-</mo><mn>1</mn></mrow></msup><msub><mi>D</mi><mi>x</mi></msub><mo>+</mo><msubsup><mi>m</mi><mn>24</mn><mrow><mi>w</mi><mi>p</mi></mrow></msubsup></mrow><mrow><mo>&lsqb;</mo><mtable><mtr><mtd><msubsup><mi>m</mi><mn>11</mn><mrow><mi>w</mi><mi>p</mi></mrow></msubsup></mtd><mtd><msubsup><mi>m</mi><mn>12</mn><mrow><mi>w</mi><mi>p</mi></mrow></msubsup></mtd><mtd><msubsup><mi>m</mi><mn>13</mn><mrow><mi>w</mi><mi>p</mi></mrow></msubsup></mtd></mtr></mtable><mo>&rsqb;</mo><msup><msub><mi>C</mi><mi>x</mi></msub><mrow><mo>-</mo><mn>1</mn></mrow></msup><msub><mi>D</mi><mi>x</mi></msub><mo>+</mo><msubsup><mi>m</mi><mn>14</mn><mrow><mi>w</mi><mi>p</mi></mrow></msubsup></mrow></mfrac><msub><mi>&phi;</mi><mi>x</mi></msub><mo>+</mo><mi>t</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000917769400000011.GIF" wi="1838" he="206" /></maths>其中,(i,j)是对应相位的坐标,公式(1)中的各变量分别由下述公式(2)~(6)所得:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>t</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><msub><mi>C</mi><mn>1</mn></msub></mfrac><msub><mi>&Sigma;&phi;</mi><mi>y</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>&phi;</mi><mi>x</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000917769400000012.GIF" wi="1844" he="174" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>M</mi><mrow><mi>w</mi><mi>c</mi></mrow></msub><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msubsup><mi>m</mi><mn>11</mn><mrow><mi>w</mi><mi>c</mi></mrow></msubsup></mtd><mtd><msubsup><mi>m</mi><mn>12</mn><mrow><mi>w</mi><mi>c</mi></mrow></msubsup></mtd><mtd><msubsup><mi>m</mi><mn>13</mn><mrow><mi>w</mi><mi>c</mi></mrow></msubsup></mtd><mtd><msubsup><mi>m</mi><mn>14</mn><mrow><mi>w</mi><mi>c</mi></mrow></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>m</mi><mn>21</mn><mrow><mi>w</mi><mi>c</mi></mrow></msubsup></mtd><mtd><msubsup><mi>m</mi><mn>22</mn><mrow><mi>w</mi><mi>c</mi></mrow></msubsup></mtd><mtd><msubsup><mi>m</mi><mn>23</mn><mrow><mi>w</mi><mi>c</mi></mrow></msubsup></mtd><mtd><msubsup><mi>m</mi><mn>24</mn><mrow><mi>w</mi><mi>c</mi></mrow></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>m</mi><mn>31</mn><mrow><mi>w</mi><mi>c</mi></mrow></msubsup></mtd><mtd><msubsup><mi>m</mi><mn>32</mn><mrow><mi>w</mi><mi>c</mi></mrow></msubsup></mtd><mtd><msubsup><mi>m</mi><mn>33</mn><mrow><mi>w</mi><mi>c</mi></mrow></msubsup></mtd><mtd><msubsup><mi>m</mi><mn>34</mn><mrow><mi>w</mi><mi>c</mi></mrow></msubsup></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000917769400000013.GIF" wi="1854" he="391" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>M</mi><mrow><mi>w</mi><mi>p</mi></mrow></msub><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msubsup><mi>m</mi><mn>11</mn><mrow><mi>w</mi><mi>p</mi></mrow></msubsup></mtd><mtd><msubsup><mi>m</mi><mn>12</mn><mrow><mi>w</mi><mi>p</mi></mrow></msubsup></mtd><mtd><msubsup><mi>m</mi><mn>13</mn><mrow><mi>w</mi><mi>p</mi></mrow></msubsup></mtd><mtd><msubsup><mi>m</mi><mn>14</mn><mrow><mi>w</mi><mi>p</mi></mrow></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>m</mi><mn>21</mn><mrow><mi>w</mi><mi>p</mi></mrow></msubsup></mtd><mtd><msubsup><mi>m</mi><mn>22</mn><mrow><mi>w</mi><mi>p</mi></mrow></msubsup></mtd><mtd><msubsup><mi>m</mi><mn>23</mn><mrow><mi>w</mi><mi>p</mi></mrow></msubsup></mtd><mtd><msubsup><mi>m</mi><mn>24</mn><mrow><mi>w</mi><mi>p</mi></mrow></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>m</mi><mn>31</mn><mrow><mi>w</mi><mi>p</mi></mrow></msubsup></mtd><mtd><msubsup><mi>m</mi><mn>32</mn><mrow><mi>w</mi><mi>p</mi></mrow></msubsup></mtd><mtd><msubsup><mi>m</mi><mn>33</mn><mrow><mi>w</mi><mi>p</mi></mrow></msubsup></mtd><mtd><msubsup><mi>m</mi><mn>34</mn><mrow><mi>w</mi><mi>p</mi></mrow></msubsup></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000917769400000021.GIF" wi="1869" he="350" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>C</mi><mi>x</mi></msub><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><msubsup><mi>m</mi><mn>11</mn><mrow><mi>w</mi><mi>c</mi></mrow></msubsup><mo>-</mo><msubsup><mi>m</mi><mn>31</mn><mrow><mi>w</mi><mi>c</mi></mrow></msubsup><msup><mi>x</mi><mi>c</mi></msup></mrow></mtd><mtd><mrow><msubsup><mi>m</mi><mn>12</mn><mrow><mi>w</mi><mi>c</mi></mrow></msubsup><mo>-</mo><msubsup><mi>m</mi><mn>32</mn><mrow><mi>w</mi><mi>c</mi></mrow></msubsup><msup><mi>x</mi><mi>c</mi></msup></mrow></mtd><mtd><mrow><msubsup><mi>m</mi><mn>13</mn><mrow><mi>w</mi><mi>c</mi></mrow></msubsup><mo>-</mo><msubsup><mi>m</mi><mn>33</mn><mrow><mi>w</mi><mi>c</mi></mrow></msubsup><msup><mi>x</mi><mi>c</mi></msup></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mi>m</mi><mn>21</mn><mrow><mi>w</mi><mi>c</mi></mrow></msubsup><mo>-</mo><msubsup><mi>m</mi><mn>31</mn><mrow><mi>w</mi><mi>c</mi><mi>c</mi></mrow></msubsup><mi>y</mi></mrow></mtd><mtd><mrow><msubsup><mi>m</mi><mn>22</mn><mrow><mi>w</mi><mi>c</mi></mrow></msubsup><mo>-</mo><msubsup><mi>m</mi><mn>32</mn><mrow><mi>w</mi><mi>c</mi></mrow></msubsup><msup><mi>y</mi><mi>c</mi></msup></mrow></mtd><mtd><mrow><msubsup><mi>m</mi><mn>23</mn><mrow><mi>w</mi><mi>c</mi></mrow></msubsup><mo>-</mo><msubsup><mi>m</mi><mn>33</mn><mrow><mi>w</mi><mi>c</mi></mrow></msubsup><msup><mi>y</mi><mi>c</mi></msup></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mi>m</mi><mn>11</mn><mrow><mi>w</mi><mi>p</mi></mrow></msubsup><mo>-</mo><msubsup><mi>m</mi><mn>31</mn><mrow><mi>w</mi><mi>p</mi></mrow></msubsup><msub><mi>&phi;</mi><mi>x</mi></msub></mrow></mtd><mtd><mrow><msubsup><mi>m</mi><mn>12</mn><mrow><mi>w</mi><mi>p</mi></mrow></msubsup><mo>-</mo><msubsup><mi>m</mi><mn>32</mn><mrow><mi>w</mi><mi>p</mi></mrow></msubsup><msub><mi>&phi;</mi><mi>x</mi></msub></mrow></mtd><mtd><mrow><msubsup><mi>m</mi><mn>13</mn><mrow><mi>w</mi><mi>p</mi></mrow></msubsup><mo>-</mo><msubsup><mi>m</mi><mn>33</mn><mrow><mi>w</mi><mi>p</mi></mrow></msubsup><msub><mi>&phi;</mi><mi>x</mi></msub></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000917769400000022.GIF" wi="1846" he="332" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>D</mi><mi>x</mi></msub><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><msubsup><mi>m</mi><mn>34</mn><mrow><mi>w</mi><mi>c</mi></mrow></msubsup><msup><mi>x</mi><mi>c</mi></msup><mo>-</mo><msubsup><mi>m</mi><mn>14</mn><mrow><mi>w</mi><mi>c</mi></mrow></msubsup></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mi>m</mi><mn>34</mn><mrow><mi>w</mi><mi>c</mi></mrow></msubsup><msup><mi>y</mi><mi>c</mi></msup><mo>-</mo><msubsup><mi>m</mi><mn>24</mn><mrow><mi>w</mi><mi>c</mi></mrow></msubsup></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mi>m</mi><mn>34</mn><mrow><mi>w</mi><mi>c</mi></mrow></msubsup><msub><mi>&phi;</mi><mi>x</mi></msub><mo>-</mo><msubsup><mi>m</mi><mn>14</mn><mrow><mi>w</mi><mi>p</mi></mrow></msubsup></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000917769400000023.GIF" wi="1861" he="367" /></maths>其中,t(i,j)即为坐标(i,j)的t,是补偿系数,(x<sup>c</sup>,y<sup>c</sup>)代表摄像机平面的坐标,<img file="FDA0000917769400000024.GIF" wi="1812" he="87" />和<img file="FDA0000917769400000025.GIF" wi="105" he="79" />分别是摄像机的参数矩阵M<sub>wc</sub>的参数,<img file="FDA0000917769400000026.GIF" wi="888" he="79" /><img file="FDA0000917769400000027.GIF" wi="796" he="79" />和<img file="FDA0000917769400000028.GIF" wi="105" he="78" />分别是投影机的参数矩阵M<sub>wp</sub>的参数;(4)通过分别投射一副全黑和全白的图片,然后计算出标靶表面的反射率,并依据该反射率计算出标靶表面每一像素点反射率的变化方向,记为θ;(5)根据θ利用下述公式(7)求出新的相位值φ(x<sup>c</sup>,y<sup>c</sup>):<img file="FDA0000917769400000029.GIF" wi="1918" he="512" />(6)对该新生成的相位值φ(x<sup>c</sup>,y<sup>c</sup>)进行三维重建。
地址 610000 四川省成都市一环路南一段24号