留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

深度学习精确相位获取的离焦投影三维测量

赵洋 傅佳安 于浩天 韩静 郑东亮

赵洋, 傅佳安, 于浩天, 韩静, 郑东亮. 深度学习精确相位获取的离焦投影三维测量[J]. 红外与激光工程, 2020, 49(7): 20200012. doi: 10.3788/IRLA20200012
引用本文: 赵洋, 傅佳安, 于浩天, 韩静, 郑东亮. 深度学习精确相位获取的离焦投影三维测量[J]. 红外与激光工程, 2020, 49(7): 20200012. doi: 10.3788/IRLA20200012
Zhao Yang, Fu Jia'an, Yu Haotian, Han Jing, Zheng Dongliang. Defocus projection three-dimensional measurement based on deep learning accurate phase acquisition[J]. Infrared and Laser Engineering, 2020, 49(7): 20200012. doi: 10.3788/IRLA20200012
Citation: Zhao Yang, Fu Jia'an, Yu Haotian, Han Jing, Zheng Dongliang. Defocus projection three-dimensional measurement based on deep learning accurate phase acquisition[J]. Infrared and Laser Engineering, 2020, 49(7): 20200012. doi: 10.3788/IRLA20200012

深度学习精确相位获取的离焦投影三维测量

doi: 10.3788/IRLA20200012
基金项目: 国家自然科学基金(61727802,61971227);江苏省重点研发计划(BE2018126)
详细信息
    作者简介:

    赵洋(1996-),男,硕士生,主要从事三维测量等方面的研究工作。 Email: zhaoy_wxh@163.com

    郑东亮(1984-),男,副教授,硕士生导师,博士,主要从事三维测量和智能计算等方面的研究工作。 Email:dlzheng@njust.edu.cn

  • 中图分类号: TP391.4

Defocus projection three-dimensional measurement based on deep learning accurate phase acquisition

  • 摘要: 数字光栅投影三维测量技术,通过离焦投影二值光栅条纹,生成三维测量所需正弦光栅条纹,能够实现超高投影速度,在快速三维测量领域具有极大潜力。但是,二值光栅条纹不可避免地包含高次谐波,导致计算所得相位包含相位误差,进而降低了快速三维测量精度。提出了一种基于深度学习精确相位获取的离焦投影三维测量方法,通过构建含噪声相位到精确相位的端到端深度卷积神经网络,降低高次谐波引入的相位误差,进而实现快速精确三维测量。首先,以理论分析证明所提方法的可行性,并以仿真和实验进一步验证了所提方法的有效性和精确性。与现有快速三维测量方法相比,所提方法在保证测量速度的同时保证测量精度。
  • 图  1  传统快速三维测量示意图

    Figure  1.  Schematic diagram of conventional high-speed 3D measurement

    图  2  精确相位获取的原理图。(a) 有噪声的相位图像;(b) 提出的PDNet;(c) 使用所提方法得到的结果

    Figure  2.  Schematic of accurate phase acquisition. (a) Phase with noise;(b) proposed PDNet;(c) result by the proposed method

    图  3  PDNet的结构

    Figure  3.  Structure of the PDNet

    图  4  仿真结果。(a) 有噪声的相位;(b) 使用所提方法得到的结果;(c) 图4(a)4(b)与理想相位对比第100行细节信息;(d) 图4(b)4(c)第100行误差曲线

    Figure  4.  Simulation results. (a) Phase map with noise;(b) result by the proposed method;(c) comparison of the phase in the 100th column from Fig. 4(a), 4(b) and the ground truth;(d) plotting error curves of the 100th column of Fig. 4(b) and 4(c)

    图  5  T=48,54,60,66,72,78时仿真得到包裹相位的相位误差

    Figure  5.  Phase error of simulated wrapped phase at T = 48,54,60,66,72,78

    图  6  实验结果。(a) 有噪声的相位;(b) 使用所提方法得到的结果;(c) 图6(a)6(b)第100行细节信息与理想相位对比;(d) 图6(b)6(c)第100行误差曲线

    Figure  6.  Experiment results. (a) Phase map with noise;(b) result by the proposed method;(c) comparison of the phase in the 100th column from Fig. 6(a), 6(b) and the ground truth;(d) plotting error curves of the 100th column of Fig. 6(b) and 6(c)

    图  7  T = 48,54,60,66,72,78时实际测量得到包裹相位的相位误差

    Figure  7.  Phase error of actual measurement of wrapped phase at T = 48, 54, 60, 66, 72, 78

    图  8  去噪前后包裹相位及其误差对比。(a) 有噪声的相位;(b) 使用所提方法得到的结果;(c) 去噪前相位误差;(d) 去噪后相位误差

    Figure  8.  Wrapping phase map before and after denoising and its error comparison. (a) Phase map with noise;(b) result by the proposed method;(c) phase error before denoising;(d) phase error after denoising

    图  9  重建效果对比。(a) 使用所提方法处理前的重建结果;(b) 理想情况下的重建结果;(c) 使用所提方法处理后的重建结果;(d) 重建误差曲线

    Figure  9.  Comparison of reconstruction results. (a) Reconstruction results before processing using the proposed method;(b) ideal reconstruction results;(c) reconstruction result by the proposed method;(d) reconstruction error curve

    表  1  网络主要模块及参数

    Table  1.   Main network modules and parameters

    TypeOut-FOut-Res
    ENCODER 1 Downsampler block 16 512×256
    2 Downsampler block 64 256×128
    3−7 5 × Non-bt-1D 64 256×128
    8 Downsampler block 128 128×64
    9 Non-bt-1D(dilated 2) 128 128×64
    10 Non-bt-1D(dilated 4) 128 128×64
    11 Non-bt-1D(dilated 8) 128 128×64
    12 Non-bt-1D(dilated 16) 128 128×64
    13 Non-bt-1D(dilated 2) 128 128×64
    14 Non-bt-1D(dilated 4) 128 128×64
    15 Non-bt-1D(dilated 8) 128 128×64
    16 Non-bt-1D(dilated 16) 128 128×64
    DECODER 17 Deconvolution (unsampling) 64 256×128
    18−19 2 × Non-bt-1D 64 256×128
    20 Deconvolution (unsampling) 16 512×256
    21−22 2 × Non-bt-1D 16 512×256
    23 Deconvolution (unsampling) C 1 024×512
    下载: 导出CSV
  • [1] Zhang S. Recent progresses on real-time 3D shape measurement using digital fringe projection techniques [J]. Optics and Lasers in Engineering, 2010, 48(2): 149−158. doi:  10.1016/j.optlaseng.2009.03.008
    [2] Hong Ziming, Ai Qingsong, Chen Kun. High precise 3D visual measurement based on fiber laser [J]. Infrared and Laser Engineering, 2018, 47(8): 0803011. (in Chinese) doi:  10.3788/IRLA201847.0803011
    [3] Zheng D, Da F, Kemao Q, et al. Phase error analysis and compensation for phase shifting profilometry with projector defocusing [J]. Applied Optics, 2016, 55(21): 5721. doi:  10.1364/AO.55.005721
    [4] Dai Meiling, Yang Fujun, He Xiaoyuan. Three-dimensional shape measurement of objects with discontinuities by dual-frequency color fringe projection [J]. Optics and Precision Engineering, 2013, 21(1): 11−16. (in Chinese)
    [5] Su X, Chen W. Fourier transform profilometry: a review [J]. Optics and Lasers in Engineering, 2001, 35(5): 263−284. doi:  10.1016/S0143-8166(01)00023-9
    [6] Zuo C, Feng S, Huang L, et al. Phase shifting algorithms for fringe projection profilometry: A review [J]. Optics and Lasers in Engineering, 2018, 109: 23−59. doi:  10.1016/j.optlaseng.2018.04.019
    [7] Song Z. Flexible 3D shape measurement using projector defocusing: extended measurement range [J]. Optics Letters, 2010, 35(7): 934−6. doi:  10.1364/OL.35.000934
    [8] Zuo C, Chen Q, Feng S, et al. Optimized pulse width modulation pattern strategy for three-dimensional profilometry with projector defocusing [J]. Applied Optics, 2012, 51(19): 4477−4490. doi:  10.1364/AO.51.004477
    [9] Zhao Liwei, Da Feipeng, Zheng Dongliang. Method for binary grating generation using defocused projection for three-dimensional measurement [J]. Acta Optica Sinica, 2016, 36(8): 0812005. (in Chinese) doi:  10.3788/AOS201636.0812005
    [10] Zheng D, Da F. Self-correction phase unwrapping method based on gray-code light [J]. Optics and Lasers in Engineering, 2012, 50(8): 1130−1139. doi:  10.1016/j.optlaseng.2012.01.019
    [11] Zhang Lei, Jiao Xiaoxue, Zhou Liqiu, et al. Three-dimensional shape acquisition method by integral imaging based on corresponding points [J]. Chinese Optics, 2015, 8(1): 45−50. (in Chinese) doi:  10.3788/co.20150801.0045
    [12] Sun Riming, Lin Tingting, Ji Lin, et al. A general linear imaging modeling method for space unstability targets and parameter optimization of linear laser radar [J]. Optics and Precision Engineering, 2018, 26(6): 1524−1532. (in Chinese) doi:  10.3788/OPE.20182606.1524
    [13] An Dong, Chen Li, Ding Yifei, et al. Optical system model and calibration of grating projection phase method [J]. Chinese Optics, 2015, 8(2): 248−254. (in Chinese) doi:  10.3788/co.20150802.0248
    [14] Romera E, Alvarez J M, Bergasa L M, et al. ERFNet: Efficient residual factorized ConvNet for real-time semantic segmentation [J]. IEEE Transactions on Intelligent Transportation Systems, 2018, 19(1): 263−272. doi:  10.1109/TITS.2017.2750080
  • [1] 朱坡, 张宗华, 高楠, 高峰, 王张颖.  彩色高反光物体表面三维形貌测量技术 . 红外与激光工程, 2023, 52(7): 20220761-1-20220761-7. doi: 10.3788/IRLA20220761
    [2] 何信欣, 刘斌, 王春柳, 吴冠豪.  用于高反表面测量的分块平滑自适应条纹投影方法 . 红外与激光工程, 2023, 52(5): 20220825-1-20220825-14. doi: 10.3788/IRLA20220825
    [3] 杨静雯, 张宗华, 付莉娜, 李雁玲, 高楠, 高峰.  利用抖动算法扩展深度范围的三维形貌测量术 . 红外与激光工程, 2023, 52(8): 20230059-1-20230059-10. doi: 10.3788/IRLA20230059
    [4] 刘今越, 马旭哲, 崔月盟, 薛路明, 李铁军.  基于自适应插值滤波的室内环境高速重建系统 . 红外与激光工程, 2022, 51(8): 20210625-1-20210625-8. doi: 10.3788/IRLA20210625
    [5] 朱硕, 郭恩来, 柏连发, 韩静.  高效学习的透过未知散射介质的相位恢复方法 . 红外与激光工程, 2022, 51(2): 20210889-1-20210889-9. doi: 10.3788/IRLA20210889
    [6] 杨程, 鄢秋荣, 祝志太, 王逸凡, 王明, 戴伟辉.  基于深度学习的压缩光子计数激光雷达 . 红外与激光工程, 2020, 49(S2): 20200380-20200380. doi: 10.3788/IRLA20200380
    [7] 左超, 张晓磊, 胡岩, 尹维, 沈德同, 钟锦鑫, 郑晶, 陈钱.  3D真的来了吗?— 三维结构光传感器漫谈 . 红外与激光工程, 2020, 49(3): 0303001-0303001-45. doi: 10.3788/IRLA202049.0303001
    [8] 姜宏志, 李宇曦, 赵慧洁.  单像素成像在三维测量中的应用 . 红外与激光工程, 2020, 49(3): 0303017-0303017-9. doi: 10.3788/IRLA202049.0303017
    [9] 宋晓凤, 李居朋, 陈后金, 李丰, 万成凯.  多场景下结构光三维测量激光中心线提取方法 . 红外与激光工程, 2020, 49(1): 0113004-0113004(8). doi: 10.3788/IRLA202049.0113004
    [10] 李红梅, 曹益平, 王亚品, 王璐, 李城梦, 安海华, 徐财.  基于灰度拓展复合光栅的单帧三维测量方法 . 红外与激光工程, 2020, 49(6): 20200034-1-20200034-8. doi: 10.3788/IRLA20200034
    [11] 伏燕军, 韩勇华, 陈元, 张鹏飞, 桂建楠, 钟可君, 黄采敏.  基于相位编码的三维测量技术研究进展 . 红外与激光工程, 2020, 49(3): 0303010-0303010-15. doi: 10.3788/IRLA202049.0303010
    [12] 冯世杰, 左超, 尹维, 陈钱.  深度学习技术在条纹投影三维成像中的应用 . 红外与激光工程, 2020, 49(3): 0303018-0303018-17. doi: 10.3788/IRLA202049.0303018
    [13] 钟锦鑫, 尹维, 冯世杰, 陈钱, 左超.  基于深度学习的散斑投影轮廓术 . 红外与激光工程, 2020, 49(6): 20200011-1-20200011-11. doi: 10.3788/IRLA20200011
    [14] 张钊, 韩博文, 于浩天, 张毅, 郑东亮, 韩静.  多阶段深度学习单帧条纹投影三维测量方法 . 红外与激光工程, 2020, 49(6): 20200023-1-20200023-8. doi: 10.3788/IRLA20200023
    [15] 周宏强, 黄玲玲, 王涌天.  深度学习算法及其在光学的应用 . 红外与激光工程, 2019, 48(12): 1226004-1226004(20). doi: 10.3788/IRLA201948.1226004
    [16] 张秀玲, 侯代标, 张逞逞, 周凯旋, 魏其珺.  深度学习的MPCANet火灾图像识别模型设计 . 红外与激光工程, 2018, 47(2): 203006-0203006(6). doi: 10.3788/IRLA201847.0203006
    [17] 唐聪, 凌永顺, 郑科栋, 杨星, 郑超, 杨华, 金伟.  基于深度学习的多视窗SSD目标检测方法 . 红外与激光工程, 2018, 47(1): 126003-0126003(9). doi: 10.3788/IRLA201847.0126003
    [18] 唐聪, 凌永顺, 杨华, 杨星, 郑超.  基于深度学习物体检测的视觉跟踪方法 . 红外与激光工程, 2018, 47(5): 526001-0526001(11). doi: 10.3788/IRLA201847.0526001
    [19] 刘巍, 高鹏, 张洋, 杨帆, 李晓东, 兰志广, 贾振元, 高航.  面向大型零部件的三维形面高精度测量方法 . 红外与激光工程, 2017, 46(3): 317003-0317003(9). doi: 10.3788/IRLA201746.0317003
    [20] 李彪, 吴海涛, 张建成, 伏燕军.  正弦脉冲宽度调制条纹结合相位编码条纹的三维测量方法 . 红外与激光工程, 2016, 45(6): 617006-0617006(6). doi: 10.3788/IRLA201645.0617006
  • 加载中
图(9) / 表(1)
计量
  • 文章访问数:  619
  • HTML全文浏览量:  144
  • PDF下载量:  78
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-04-13
  • 修回日期:  2020-05-22
  • 网络出版日期:  2020-07-23
  • 刊出日期:  2020-07-23

深度学习精确相位获取的离焦投影三维测量

doi: 10.3788/IRLA20200012
    作者简介:

    赵洋(1996-),男,硕士生,主要从事三维测量等方面的研究工作。 Email: zhaoy_wxh@163.com

    郑东亮(1984-),男,副教授,硕士生导师,博士,主要从事三维测量和智能计算等方面的研究工作。 Email:dlzheng@njust.edu.cn

基金项目:  国家自然科学基金(61727802,61971227);江苏省重点研发计划(BE2018126)
  • 中图分类号: TP391.4

摘要: 数字光栅投影三维测量技术,通过离焦投影二值光栅条纹,生成三维测量所需正弦光栅条纹,能够实现超高投影速度,在快速三维测量领域具有极大潜力。但是,二值光栅条纹不可避免地包含高次谐波,导致计算所得相位包含相位误差,进而降低了快速三维测量精度。提出了一种基于深度学习精确相位获取的离焦投影三维测量方法,通过构建含噪声相位到精确相位的端到端深度卷积神经网络,降低高次谐波引入的相位误差,进而实现快速精确三维测量。首先,以理论分析证明所提方法的可行性,并以仿真和实验进一步验证了所提方法的有效性和精确性。与现有快速三维测量方法相比,所提方法在保证测量速度的同时保证测量精度。

English Abstract

    • 快速三维测量在先进制造、虚拟现实、安全防护等领域变得日益重要,数字光栅投影三维测量技术因其非接触、速度快、成本相对较低等优势,在快速三维测量领域表现了极大潜力[1-4]。传统光栅投影法投影符合正弦分布光栅条纹,通常由频域法或者相移法计算测量所需相位信息。频域法可依据单帧光栅条纹计算相位,但是较难保持待测物体边界、复杂面型等细节信息[5];相移法依据至少三帧以上相移光栅条纹,时域内计算得到相位信息,能较好地保持待测物体细节[6]

      数字投影设备利用数字微镜(Digital Micro-mirror Device,DMD)来投影光栅条纹,数字微镜中可以开和闭的微镜片表征对应像素,单个像素亮度强弱通过调整对应微镜片单位时间内开启时间长短来实现。投影“8”位正弦光栅条纹需要根据条纹灰度变化按照积分时间多次切换数字微镜开闭状态;而投影“1”位二值光栅条纹,数字微镜仅需保持开或闭两种状态,无需进行多次开闭状态切换,能够实现kHz(帧/秒)以上的投影速度。但投影二值光栅条纹不可避免会包含高次谐波,投影设备选择恰当离焦程度,离焦可由高斯低通滤波近似表征,即可滤除投影光栅包含的高次谐波,以生成测量所需正弦光栅条纹图像。基于二值光栅投影快速性,即可实现快速三维测量[7]

      然而,投影设备通常需要较大离焦程度来滤除二值光栅所包含高次谐波,以计算得到精确相位信息。但是,较大离焦程度会减少生成条纹对比度,降低三维测量动态范围;轻微离焦程度能够保证较高的动态测量范围,但是相应高次谐波会引入相位误差,降低三维测量精度。通常有两种思路解决上述问题:(1)优化设计二值光栅,降低所需离焦程度,主要包括:脉宽调制(PWM)和抖动(dithering)、方波等设计方法。其中,PWM [8]和dithering [9]各自适用于较窄和较宽周期二值光栅设计,投影高频方波二值光栅,操作便捷且在不同离焦程度保持良好的一致性。但是,实际测量过程中通常难以精确调整到恰当离焦程度,尤其当所选相移算法步数较低情况下,不可避免会包含高次谐波,降低三维测量精度。(2)通过附加多帧光栅或者进行频域分析进行相位误差补偿。附加多帧光栅条纹,需要投影多帧相移二值光栅,制约了三维测量速度[6];引入频域变换方法,相对难以保持待测物体边界[5]。基于上述,采用较低步数相移算法(例如:快速测量常用三步相移算法),投影设备处在不同离焦程度,均能实现精确相位获取,对快速三维测量具有重要意义。

      文中对传统快速三维测量方法做了改进,提出一种基于深度学习精确相位获取的离焦投影三维测量方法。具体而言,在快速测量过程中引入基于深度学习的PDNet(Phase Denoising Network)算法,使用含噪声相位到理想相位的端到端方式来训练PDNet模型,训练完成后使用模型就能够在测量过程中得到精确相位图像。文中将三帧相移二值光栅离焦投影,采用三步相移算法计算得到的包裹相位图像作为输入,对应的理想相位图像通过多步相移算法获得,所提方法能够实现不同离焦下的精确相位获取,进而实现快速精确三维测量。

    • 快速三维测量使用二值离焦技术来提高测量速度,二值方波凭借其获取简易以及一致性较好等一系列优点,在快速三维测量中得到广泛应用。三步相移算法因其所需投影光栅帧数较少,实际快速测量可通过离焦投影三帧相移二值方波光栅实现[7]。具体测量流程如图1所示。

      图  1  传统快速三维测量示意图

      Figure 1.  Schematic diagram of conventional high-speed 3D measurement

      图1(a)表示相机实际采集得到三帧相移光栅图像,具体可描述为[3]

      $$I_n^C(x,y) = {B_0}(x,y) + \sum\limits_{k = 1}^\infty {\left\{ {{B_k}(x,y)\cos [k{\phi _n}(x,y)]} \right\}} $$ (1)

      式中:n = 1,2,···,N表示相移步数;$(x,y)$为像素坐标;${B_0}$为背景光强;${B_k}$表示第k阶高次谐波的幅度信息;${\phi _n}$为耦合相位,可描述为:

      $${\phi _n}(x,y) = \phi (x,y) + {\delta _n}$$ (2)

      式中:$\phi $是理想相位信息;${\delta _n}$ 为第n帧条纹相移量,具体为:

      $${\delta _n} = 2\pi (n - 1)/N$$ (3)

      为了简化起见,下文省略像素坐标$(x,y)$

      图1(b)为采用相移算法计算所得真实相位${\phi ^C}$,通常采用最小二乘方法实现[6]

      $${\phi ^C} = \arctan \left[ {\frac{{ - \displaystyle\sum\limits_{n = 1}^\infty {(I_n^C\sin {\delta _n})} }}{{\displaystyle\sum\limits_{n = 1}^\infty {(I_n^C\cos {\delta _n})} }}} \right]$ .$ (4)

      由于采集条纹中含有高次谐波,真实相位信息中会引入相位误差形如[3]

      $$\Delta {\phi _N} = {\phi ^C} - \phi = \arctan \left\{ {\frac{{\displaystyle\sum\limits_{m = 1}^\infty {\left[ {B_{mN}^{ - 1}\sin (mN\phi )} \right]} }}{{{B_1} + \displaystyle\sum\limits_{m = 1}^\infty {\left[ {B_{mN}^{ + 1}\cos (mN\phi )} \right]} }}} \right\}$$ (5)

      式中:$B_{mN}^{ - 1}{\rm{ = }}B_{mN{\rm{ + }}1}^{} - B_{mN - 1}^{}$$B_{mN}^{{\rm{ + }}1}{\rm{ = }}B_{mN{\rm{ + }}1}^{}{\rm{ + }}B_{mN - 1}^{}$。该类型相位误差会给实际三维测量带来不可忽略的测量误差。

      上述过程求解得到的相位信息${\phi ^C}$是在(−π,π]范围内的一系列呈锯齿状不连续分布相位。为了获得物体的绝对相位$\varPhi$,需要对${\phi ^C}$进行相位展开,相位解包裹的方法有时域法和空域法,如图1(c)文中使用格雷编码法展开包裹相位[10]

      $$\varPhi {\rm{ }} = {\phi ^C}{\rm{ }} + K{\rm{ }} \times {\rm{ }}2{\text{π}} $$ (6)

      式中:$K$是对图像进行解码确定的相移条纹的级次信息,如图1(d)。最后,利用标定后的系统参数就可以从绝对相位中得到图1(f)的物体的三维信息[11-13]

    • 旨在轻微离焦程度下减少以致削弱高次谐波引入的测量误差,实现精确相位获取,进而实现快速精确三维测量,提出了基于深度学习精确相位获取的离焦投影三维测量方法,所提方法中使用深度学习实现精确相位获取的流程如图2所示。

      图  2  精确相位获取的原理图。(a) 有噪声的相位图像;(b) 提出的PDNet;(c) 使用所提方法得到的结果

      Figure 2.  Schematic of accurate phase acquisition. (a) Phase with noise;(b) proposed PDNet;(c) result by the proposed method

      与传统快速三维测量方法相比,所提方法在测量过程中引入基于深度学习的PDNet算法对包裹相位进行计算以得到精确相位信息。引入的PDNet算法构建含噪声相位${\phi ^C}$到近似无噪声理想相位$\phi $的端到端方式训练模型,训练过程中算法每次迭代计算出的精确相位$\phi _i^{out}$与理想相位图像$\phi $进行残差学习将产生的误差返回到网络中,使网络不断迭代,得到一个最优的网络参数,使其具有较好的去噪能力同时对不同程度的噪声具有较好的鲁棒性。文中使用的误差函数如下式:

      $$Loss = \frac{1}{m}\sum {\left\| {\phi _i^{{\rm{out}}} - \phi } \right\|^2}$$ (7)

      式中:$\phi _i^{{\rm{out}}}$为PDNet训练过程中每次迭代计算得到的相位图像,i为迭代次数,$\phi $为理想相位图像;m为相位图像中有效像素个数。最终,训练完成后利用得到的PDNet模型能够对不同场景下获取的相位图像${\phi ^C}$进行去噪,以得到测量所需精确相位图像${\phi ^{{\rm{out}}}}$

      投影设备轻微离焦程度下,实际测量过程中仅需30次以下高次谐波影响,根据方波误差分布规律,投影十五步高频二值方波条纹,即可得到理想相位$\phi $。对应的输入含噪声相位图像${\phi ^C}$即可选取三步相移算法求解得到。使用训练后得到的PDNet模型对输入相位图像${\phi ^C}$进行运算即可得到精确相位图像${\phi ^{{\rm{out}}}}$,精确相位图像${\phi ^{{\rm{out}}}}$与理想相位图像$\phi $之间存在的相位误差可用$\Delta \phi _{}^{{\rm{out}}}$表示,可简单用数学描述为:

      $$\Delta \phi _{}^{{\rm{out}}}{\rm{ }} = \phi - {\phi ^{{\rm{out}}}}$$ (8)

      据此式即可进行相位误差评估。依据去噪后的精确相位${\phi ^{{\rm{out}}}}$即可得到物体的准确三维信息。

      图2所示,PDNet算法训练完成后使用模型就能够在测量过程中将图2(a)含有噪声的相位图像${\phi ^C}$作为输入层计算,得到去噪后的精确相位图像${\phi ^{{\rm{out}}}}$,如图2(c)所示。图2(b)为引入的PDNet,具体结构见图3

      图  3  PDNet的结构

      Figure 3.  Structure of the PDNet

      图3所示,对于网络结构的设计,文中借鉴了高效残差分解网络(Efficient Residual Factorized Network,ERFNet)[14],在其基础上对网络结构进行了改进,增加了跳层连接并去掉了池化层,以更好地匹配文中的图像去噪任务。因为增加跳层可以更好地进行残差学习,利于噪声的去除,而池化层则可能导致图像丢失细节信息,内部的网络层均由特征图所构成。编码过程和解码过程都含有多尺寸的卷积层与相关的策略层,其主要模块及相关参数设置如表1所示。表1中,“Out-F”为输出图层通道数,“Out-Res”为假设输入尺寸为1 024×512时各层的输出分辨率。

      表 1  网络主要模块及参数

      Table 1.  Main network modules and parameters

      TypeOut-FOut-Res
      ENCODER 1 Downsampler block 16 512×256
      2 Downsampler block 64 256×128
      3−7 5 × Non-bt-1D 64 256×128
      8 Downsampler block 128 128×64
      9 Non-bt-1D(dilated 2) 128 128×64
      10 Non-bt-1D(dilated 4) 128 128×64
      11 Non-bt-1D(dilated 8) 128 128×64
      12 Non-bt-1D(dilated 16) 128 128×64
      13 Non-bt-1D(dilated 2) 128 128×64
      14 Non-bt-1D(dilated 4) 128 128×64
      15 Non-bt-1D(dilated 8) 128 128×64
      16 Non-bt-1D(dilated 16) 128 128×64
      DECODER 17 Deconvolution (unsampling) 64 256×128
      18−19 2 × Non-bt-1D 64 256×128
      20 Deconvolution (unsampling) 16 512×256
      21−22 2 × Non-bt-1D 16 512×256
      23 Deconvolution (unsampling) C 1 024×512

      该网络模型结构主要包含三个模块:Downsampler block,Non-bt-1D(dilated),Deconvolution (unsampling):

      (1) Downsampler block

      该模块包括:卷积核大小3×3的卷积层,并根据参数设置输入及输出的特征图通道数,用于获取噪声图的特征分布;BN层(Batch Normalization),可防止梯度消失或爆炸;ReLu激活层,也可克服梯度消失的问题并通过特征稀疏性加快网络速度。

      (2) Non-bt-1D(dilation rate)

      该模块包括:两层卷积核大小3×1的卷积层,两层卷积核大小1×3的卷积层,并可设置卷积核的间隔数量(dilation rate),形成空洞卷积以增大感受野,两层BN层,还包括dropout策略,即在前向传播的时候,让某个神经元的激活值以一定的概率停止工作,使得它不会太依赖某些局部的特征,这样可以使模型泛化性更强。

      (3) Deconvolution (unsampling)

      此模块主要包括一层卷积层、BN层和ReLu激活层,该模块主要目的是为了对特征进行上采样,最终恢复出与原图相同尺寸的目标输出图。

    • 该节在模拟环境下验证所提方法。选取的二值方波条纹尺寸为400×400像素,周期在[48,78]之间,使用高斯低通滤波来模拟离焦,标准差$\sigma $的范围在[2,5]之间,使用的高斯低通滤波器大小设定为$[6\sigma ] \times [6\sigma ]$。通过选取合适的方波条纹周期和高斯低通滤波器的标准差,以三步相移方波条纹仿真求解得到的包裹相位图像为输入,十五步相移方波条纹仿真求解得到的相位图像为近似理想相位,最终生成140组对应的包裹相位图像,其中100组作为训练集,验证集和测试集各20组。

      图4中,选择周期T = 60的二值方波条纹,即30个黑色像素加30个白色像素,对其进行仿真来验证所提方法。图4(a)为输入带有噪声的相位图像,图4 (b)是使用所提方法预测出的相位图像,从视觉对比可以看出高次谐波引起的噪声在相位图像中明显降低。为了更精确对比去噪效果,图4(c)展示了去噪前后以及理想相位图像第100行像素中部分细节信息对比,图4(d)展示了图4(a)图4(b)各自减去理想相位图像后第100行像素的部分误差曲线图。图4(c)展示了使用所提方法处理后的结果与理想相位图像更为接近,图4(d)展示了使用所提方法处理后的误差曲线更接近零,计算其MAE可知所提方法将相位噪声从0.108 6 rad降至0.020 2 rad,去噪效果显著。

      图  4  仿真结果。(a) 有噪声的相位;(b) 使用所提方法得到的结果;(c) 图4(a)4(b)与理想相位对比第100行细节信息;(d) 图4(b)4(c)第100行误差曲线

      Figure 4.  Simulation results. (a) Phase map with noise;(b) result by the proposed method;(c) comparison of the phase in the 100th column from Fig. 4(a), 4(b) and the ground truth;(d) plotting error curves of the 100th column of Fig. 4(b) and 4(c)

      为了验证所提方法在不同周期条纹下的表现,文中分别计算了模拟数据中周期T=48,54,60,66,72以及78的相位图像误差MAE数值并绘制误差曲线,如图5所示。由MAE误差曲线对比可得,文中所提出的方法对于不同周期相位图像具有良好的一致性效果。

      图  5  T=48,54,60,66,72,78时仿真得到包裹相位的相位误差

      Figure 5.  Phase error of simulated wrapped phase at T = 48,54,60,66,72,78

    • 该节在实际构建系统中验证所提方法。所构建测量系统主要包括DLP6500投影仪,带有8 mm焦距镜头的Basler acA800-510um CMOS相机,以及一台使用TITAN RTX(NVIDIA)进行运算加速的计算机。

      真实数据集共采集580组真实场景,选取500组真实场景作为训练集,验证集和测试集各40组。真实数据集中输入和近似理想相位图像分别为离焦投影三步和十五步二值方波条纹并使用相移算法求解得到的包裹相位图像,图像的大小为400 × 400像素,投影所用二值方波光栅的周期仍在[48,78]之间。

      首先选取了白色墙面作为目标物体,投影周期在[48,78]之间的二值方波条纹进行测量。文中选取投影周期T=60的二值方波条纹计算得到的包裹相位图像进行分析。图6(a)是输入带有噪声的相位图像,图6(b)是使用所提方法得到的精确相位图像。由视觉对比可以明显看出使用所提方法得到的精确相位图像,与含有噪声的相位图像相比更加平滑,计算其MAE数值可知,所提方法将相位噪声从0.156 5 rad降至0.018 8 rad。

      图  6  实验结果。(a) 有噪声的相位;(b) 使用所提方法得到的结果;(c) 图6(a)6(b)第100行细节信息与理想相位对比;(d) 图6(b)6(c)第100行误差曲线

      Figure 6.  Experiment results. (a) Phase map with noise;(b) result by the proposed method;(c) comparison of the phase in the 100th column from Fig. 6(a), 6(b) and the ground truth;(d) plotting error curves of the 100th column of Fig. 6(b) and 6(c)

      为了进一步对比去噪表现,图6(c)展示了图6(a)图6(b)与理想相位图像第100行部分细节对比,可以清楚看出使用所提方法得到结果与理想结果更为接近。图6(d)展示了图6(a)图6(b)各自减去理想相位后第100行部分误差曲线图。红色表示使用所提方法处理后的误差曲线,该方法能够将误差降低至随机噪声水平。

      为了分析所提方法在实际测量情况下投影不同周期条纹时的表现,同样绘制周期T=48,54,60,66,72以及78的相位图像MAE误差曲线,如图7所示。由MAE误差曲线对比可得,文中所提方法在实际测量中投影不同周期相位图像同样表现出良好的一致性效果。

      图  7  T = 48,54,60,66,72,78时实际测量得到包裹相位的相位误差

      Figure 7.  Phase error of actual measurement of wrapped phase at T = 48, 54, 60, 66, 72, 78

      其次选取玩具模型作为目标物体进行三维重建,使用所提方法对投影三步二值方波条纹计算得到的的包裹相位图像去噪前后效果如图8(a)~8(b)所示。为了更准确地评估文中所提方法的效果,对处理前后的包裹相位做了误差对比,如图8(c)~8(d)所示,由误差对比可以看出虽然在边缘以及细节处仍存在部分误差,但整体去噪效果显著。此外,还计算了误差去除前后的相位误差MAE数值,处理前为0.036 7 rad,处理后为0.024 6 rad。

      图  8  去噪前后包裹相位及其误差对比。(a) 有噪声的相位;(b) 使用所提方法得到的结果;(c) 去噪前相位误差;(d) 去噪后相位误差

      Figure 8.  Wrapping phase map before and after denoising and its error comparison. (a) Phase map with noise;(b) result by the proposed method;(c) phase error before denoising;(d) phase error after denoising

      图8中物体的相位误差对比结果可知,所提方法去噪效果显著。为了更进一步检验去噪效果,文中分别比较了使用所提方法处理前后相位图像的重建效果,如图9所示。图9(a)图9(b)分别为三步和十五步相移条纹求解相位的重建效果,图9(c)为使用所提方法对三步相移条纹处理后的重建结果。从视觉对比可以明显看出图9(a)的重建结果受高次谐波干扰较大,经过所提方法处理后的重建结果与近似理想重建结果较吻合。为了更准确说明所提方法的效果,将去噪前后物体三维坐标深度信息与对应近似理想结果相减得到误差曲线,如图9(d)所示。计算其MAE数值可知,使用所提方法将高次谐波引起的重建误差从0.161 2 mm降至0.089 6 mm。可以进一步得出结论,文中所提的快速测量方法能够有效较低高次谐波带来的干扰。综上,从模拟和真实数据均验证了所提快速测量方法的有效性和精确性。

      图  9  重建效果对比。(a) 使用所提方法处理前的重建结果;(b) 理想情况下的重建结果;(c) 使用所提方法处理后的重建结果;(d) 重建误差曲线

      Figure 9.  Comparison of reconstruction results. (a) Reconstruction results before processing using the proposed method;(b) ideal reconstruction results;(c) reconstruction result by the proposed method;(d) reconstruction error curve

    • 在传统快速三维测量基础上,提出一种降低高次谐波引起相位噪声的快速三维测量方法。所提方法的明显优势在于,能够在不增加额外投影光栅帧数前提下,实现轻微离焦程度下精确相位信息获取,结合已标定系统参数,即可进行精确快速三维重建。与现有快速三维测量方法相比,该方法能够在不牺牲测量速度的前提下保证测量精度。采用该方法对模拟数据和真实墙面进行验证分析。结果表明,该方法对于不同周期相移条纹计算得到的相位图像均能得到高质量相位图像,且去噪效果稳定,证明该方法的有效性。此外,文中还对真实玩具模型进行验证,对比处理前后的重建效果,证明该方法的精确性。另外,在实际测量中,仅需使用数据集训练后得到的模型,无需修改任何参数信息,即可实现不同离焦程度下的精确三维重建。

参考文献 (14)

目录

    /

    返回文章
    返回