发明名称 一种钢板表面三维重建快速实现方法
摘要 本发明涉及一种钢板表面三维重建的快速实现方法,采用红、绿、蓝三个单色光源沿不同角度同时照明钢板表面,其入射光与钢板轧制速度方向的角度在±50度范围,通过垂直于钢板表面的彩色CCD摄像机拍摄光源照射的钢板表面区域,分离其彩色图像的R、G、B通道,得到近似于红、绿、蓝光源单独照明下获得的三幅图像。本方法的复杂度与图像像素数量呈线性关系,而全局优化算法则呈平方关系,因此算法的耗时大大减少,且本算法的误差要略小于全局优化算法。
申请公布号 CN103886642A 申请公布日期 2014.06.25
申请号 CN201410137363.X 申请日期 2014.04.04
申请人 北京科技大学 发明人 徐科;王磊;徐金梧;周鹏;杨朝霖
分类号 G06T17/30(2006.01)I 主分类号 G06T17/30(2006.01)I
代理机构 北京金智普华知识产权代理有限公司 11401 代理人 皋吉甫
主权项 1.一种钢板表面三维重建的快速实现方法,其特征在于:采用红、绿、蓝三个单色光源沿不同角度同时照明钢板表面同一区域,其入射光与钢板轧制速度方向的角度在±50度范围,通过垂直于钢板表面的彩色CCD摄像机拍摄光源照射的钢板表面区域,分离其彩色图像的R、G、B通道,得到近似于红、绿、蓝光源单独照明下获得的三幅图像I<sub>R</sub>,I<sub>G</sub>,I<sub>B</sub>。I<sub>R</sub>,I<sub>G</sub>,I<sub>B</sub>的高h=2<sup>n</sup>,宽w=2<sup>n</sup>,图像坐标系(x,y)中某点的相对深度为Z(x,y),梯度矩阵P为深度矩阵Z沿x方向的差分,Q为深度矩阵Z沿y方向的差分,P、Q由I<sub>R</sub>,I<sub>G</sub>,I<sub>B</sub>和光源方向矩阵L得到,P(x,y),Q(x,y)的初始值用P<sup>1</sup>,Q<sup>1</sup>表示,深度矩阵Z当作小波分解的初始矩阵LL<sup>0</sup>,则Z可以通过式(1)所示的小波分解算法,式(2)所示的递推公式和式(3)所示的小波重构算法得到:<maths num="0001"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><msup><mi>LH</mi><mi>k</mi></msup><mo>=</mo><mi>D</mi><mrow><mo>(</mo><msup><mi>Q</mi><mi>k</mi></msup><mo>*</mo><msub><mi>L</mi><mi>D</mi></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msup><mi>HL</mi><mi>k</mi></msup><mo>=</mo><mi>D</mi><mrow><mo>(</mo><msup><mi>P</mi><mi>k</mi></msup><mo>*</mo><msup><msub><mi>L</mi><mi>D</mi></msub><mi>T</mi></msup><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msup><mi>HH</mi><mi>k</mi></msup><mo>=</mo><mi>D</mi><mrow><mo>(</mo><msup><mi>P</mi><mi>k</mi></msup><mo>*</mo><msup><msub><mi>H</mi><mi>D</mi></msub><mi>T</mi></msup><mo>)</mo></mrow><mo>=</mo><mi>D</mi><mrow><mo>(</mo><msup><mi>Q</mi><mi>k</mi></msup><mo>*</mo><msub><mi>H</mi><mi>D</mi></msub><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><msup><mi>P</mi><mi>k</mi></msup><mo>=</mo><mi>D</mi><mrow><mo>(</mo><msup><mi>P</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msup><mo>*</mo><msub><mi>L</mi><mi>D</mi></msub><mo>*</mo><msup><msub><mi>L</mi><mi>D</mi></msub><mi>T</mi></msup><mo>*</mo><msub><mi>L</mi><mi>D</mi></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msup><mi>Q</mi><mi>k</mi></msup><mo>=</mo><mi>D</mi><mrow><mo>(</mo><msup><mi>Q</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msup><mo>*</mo><msup><msub><mi>L</mi><mi>D</mi></msub><mi>T</mi></msup><mo>*</mo><msub><mi>L</mi><mi>D</mi></msub><mo>*</mo><msup><msub><mi>L</mi><mi>D</mi></msub><mi>T</mi></msup><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>LL<sup>k-1</sup>=U(LL<sup>k</sup>)*L<sub>R</sub><sup>T</sup>*L<sub>R</sub>+U(LH<sup>k</sup>)*H<sub>R</sub><sup>T</sup>*L<sub>R</sub>+U(HL<sup>k</sup>)*L<sub>R</sub><sup>T</sup>*H<sub>R</sub>+U(HH<sup>k</sup>)*H<sub>R</sub><sup>T</sup>*H<sub>R</sub>  (3)式(1)、式(2)、式(3)中:■“*”——二维卷积运算;■D(M)——二维下采样,抽取矩阵奇数行、列,D(M)<sub>(x,y)</sub>=M<sub>(2x-1,2y-1)</sub>;■U(M)——二维上采样,矩阵扩充为(2h+1)×(2w+1),其中偶数行列:U(M)<sub>(2x,2y)</sub>=M<sub>(x,y)</sub>,其余项填充0;■LL<sup>0</sup>——小波分解的初始矩阵,Z=LL<sup>0</sup>;■LL<sup>k</sup>、LH<sup>k</sup>、HL<sup>k</sup>、HH<sup>k</sup>——第k层分解的低频及高频子带(1≤k≤n),设LL<sup>n</sup>中元素的值为0;■P<sup>k</sup>,Q<sup>k</sup>——k层梯度矩阵(1≤k≤n),其中P<sup>1</sup>,Q<sup>1</sup>为P,Q,其余通过递推求得;■L<sub>D</sub>=(1,1),H<sub>D</sub>=(1,-1)——小波分解低通、高通滤波器;■L<sub>R</sub>=(0.5,0.5),H<sub>R</sub>=(-0.5,0.5)——小波重构低通、高通滤波器。
地址 100083 北京市海淀区学院路30号