高级检索

ISSN1001-3806CN51-1125/TN 网站地图

留言板

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

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

基于时空域退偏的数字全息成像去噪研究

丁伟 孔勇 杜彤耀 钱心磊 王振伟

引用本文:
Citation:

基于时空域退偏的数字全息成像去噪研究

    作者简介: 丁伟(1996-),男,硕士研究生,从事数字全息及数字全息显微方面的研究.
    通讯作者: 孔勇, kkyy7757@aliyun.com
  • 基金项目:

    上海市自然科学基金资助项目 19ZR1421700

  • 中图分类号: O438.1

Research on digital holographic imaging denoising based on time-space domain depolarization

    Corresponding author: KONG Yong, kkyy7757@aliyun.com ;
  • CLC number: O438.1

  • 摘要: 为了抑制全息记录过程中的散斑噪声,采用一种通过降低相干光源相干性的方法,并进行了理论分析和实验验证。利用旋转二元波片进行退偏,在10s内记录70幅全息图进行平均化降噪; 同时利用毛玻璃对相干光进行空间退相干,降低了光源的相干性从而降低了散斑噪声。将单独采用时间退偏操作的重建图像与单独采用空间退偏操作的重建图像以及两者相结合的重建图像进行了对比。通过实验验证了所提方法的可行性,通过分析频谱比较证实了几种方法,利用信噪比直观地证实了几种方法的有效性。结果表明,时间、空间退偏结合对散斑噪声的抑制有着明显的优势,空间退偏的重建像的信噪比为0.4459,时间退偏的重建像信噪比为1.5155,而两者结合后的信噪比为1.7162,采用时间退偏以及空间退偏叠加的方式有效地抑制了相干噪声。该实验为研究数字全息成像的过程中由于光源相干度过高带来的散斑噪声干扰的问题提供了一个有效的抑制方法。
  • Figure 1.  Time domain depolarization device

    Figure 2.  Poincare ball depolarization effect simulation

    a—geometric representation of any point P b—a point P of fully polarized light on the surface of the bunga sphere c—the effect of arbitrary polarization states

    Figure 3.  Proposed method holography setup

    Figure 4.  Comparison of reconstruction maps of several methods

    a—reconstruction of original holography b—reconstruction of spatial depolarization hologram c—reconstruction of time domain depolarization hologram d—reconstruction of time-space depolarization hologram

    Figure 5.  The 3-D spectral distribution of Fig. 4a~Fig. 4d

    a—3-D spectrum distribution of the original hologram b—3-D spectrum distribution of the spatial depolarization c—3-D spectrum distribution of the temporal depolarization d—3-D spectrum distribution of the the proposed method

    Figure 6.  a~d sectional view of Fig. 5a~Fig. 5d, respectively

    a—spectral cross section of the original hologram b—spectral cross section of the spatial depolarizing c—spectral cross section of the time-domain de-polarizing d—spectral cross section of the proposed method

    Table 1.  Signal-to-noise ratio of different methods

    spatial-domain
    depolarization
    time-domain
    depolarization
    spatiotemporal
    depolarization
    RSNR/dB 0.4459 1.5155 1.7162
    下载: 导出CSV

    Table 2.  Parameters of noise suppression of different methods

    number of pixels between-1 and 0 image edges in x and y directions number of pixels between+1 and 0 image edges in x and y directions
    original hologram spectrum Δx=59, Δy=140 Δx=56, Δy=143
    spatial filtering spectrum Δx=180, Δy=172 Δx=182, Δy=176
    temporal filtering spectrum Δx=163, Δy=188 Δx=162, Δy=187
    spatiotemporal filtering spectrum Δx=85, Δy=164 Δx=112, Δy=168
    下载: 导出CSV
  • [1]

    MAKSIMOVA L A, RYABUKHO P V, MYSINA N Y, et al.Digital speckle photography of subpixel displacements of speckle structures bases on analysis of their spatial spectra [J]. Optics and Spectroscopy, 2018, 124(4): 549-559. doi: 10.1134/S0030400X18040124
    [2]

    JEONG W Y, SON K C, CHO J H, et al.Suppression algorithm of speckle noise for parallel phase-shift digital holography [J]. Optics and Laser Technology, 2019, 112:93-100. 
    [3]

    LIU Y, WANG Zh, HUANG J H. Recent progress on aberration compensation and coherent noise suppression in digital holography [J].Applied Sciences, 2018, 8(3) :11-15. 
    [4]

    BASANTA B, HOA P, MUSTAFA M, et al.Diffraction phase microscopy with white light[J].Optics Letters, 2012, 37(6):1094-1096. doi: 10.1364/OL.37.001094
    [5]

    YUTAKA M, TAKANORI N. Speckle reduction in hologram generation based on spherical waves synthesis using kow-coherence digital holograph[J]. Journal of Display Technology, 2015, 11(10):867-872. doi: 10.1109/JDT.2014.2342254
    [6]

    JORGE G S. Noise reduction in digital lensless holographic microscopy by engineering the light from a light-emitting diode [J]. Applied Optics, 2013, 52(1):A232-A239. 
    [7]

    FENG P, LI Z Y, WEN X.Coherent noise reduction in digital holographic microscopy by averaging multiple holograms recorded with a multimode laser [J]. Optics Express, 2017, 25(18):21815-21826. doi: 10.1364/OE.25.021815
    [8]

    CAO G Q, LIU Y T, WANG L L. Improved filtering of non-subsampled contour waves in speckle fringes [J].Laser Technology, 2019, 43(2):269-274(in Chinese).
    [9]

    TAKAHIKO F, YUTAKA M, TAKANORI N.Speckle reduction by spatial-domain mask in digital holography[J].Journal of Display Technology, 2015, 12(4):315-322. 
    [10]

    WANG J J, LIN H H, SUI Zh, et al.Study on a new liquid crystal depolarizer[J].Laser Technology, 2007, 31(3):311-313(in Chinese). 
    [11]

    SONG Sh X, SONG L K.Research on depolarizer of binary composite wave plate[J].Laser Technology, 2009, 33(6):654-656(in Chinese). 
    [12]

    SU H M, ZHENG X G, WANG X, et al.Computer simulation of the effect of polarization on laser holography[J].Acta Physica Sinica, 2002, 51(5):1044-1048(in Chinese). 
    [13]

    LIU Q, SONG L K, LIU C C.Research on the effect of polarimeter and λ/2 wave plate composite depolarization[J].Laser Journal, 2013, 34(1):13-14(in Chinese).
    [14]

    LI K T, WU F Q, LIU Q, et al.Muller matrix analysis of depolarization performance of ternary wave plate composite depolarizer[J].Laser Technology, 2015, 39(4):549-551(in Chinese). 
    [15]

    KONG W J, LI G H. Determining the ratio between the azimuth and the major and minor axes of elliptically polarized light using stokes vectors and Muller matrix[J].Journal of Applied Optics, 2002, 23(2):42-44(in Chinese).
    [16]

    HUANG X H, JIA Z H, ZHOU J L, et al. Speckle reduction of reconstructions of digital holograms using gamma-correction and filtering [J]. IEEE Access, 2018, 6: 5227-5235. 
    [17]

    SONG Sh X.Study on the effect of monochromatic optical depolarizer and wave plate depolarizer[D]. Qufu: Qufu Normal University, 2009: 31-33(in Chinese).
    [18]

    TU Q, YU Y J, ZHOU W J. Denoising of holographic system based on rotating frosted glass sheet[J].Optical Instruments, 2014, 36(4):337-341(in Chinese).
    [19]

    REN G J, YAO J Q, ZHAO J L. Research on the matrix of wave plate and rotator depolarization[J].Laser Technology, 2007, 31(3):314-316(in Chinese).
    [20]

    MENG F H, SONG L K.The poincare ball representation of the effect of polarized optical devices[J]. Journal of Qufu Normal University, 2005, 31(4):71-73(in Chinese).
  • [1] 秦怡王志鹏李根全 . 部分相干光数字全息消色差相移技术研究. 激光技术, 2012, 36(3): 424-427.
    [2] 王华英王广俊谢建军赵洁王大勇 . 数字全息显微中的准直光再现. 激光技术, 2008, 32(2): 131-133,136.
    [3] 姜雨甫周昕王运郭雪帆 . 采用二值化相位光瞳的光学扫描全息方法. 激光技术, 2013, 37(6): 723-727. doi: 10.7510/jgjs.issn.1001-3806.2013.06.004
    [4] 宋修法于梦杰王华英刘佐强高亚飞刘飞飞 . 物光与参考光强度比对数字全息再现像质的影响. 激光技术, 2014, 38(6): 859-862. doi: 10.7510/jgjs.issn.1001-3806.2014.06.029
    [5] 孙彩霞王素莲李若平路海唐道广黄明举 . 一种抗湿性光聚物在不同厚度下的全息特性. 激光技术, 2008, 32(5): 545-547,550.
    [6] 何飞黄明举 . 有机-无机纳米复合光聚物材料的研究进展. 激光技术, 2011, 35(4): 461-467. doi: 10.3969/j.issn.1001-3806.2011.04.007
    [7] 张巧鸽楼宇丽宋庆和桂进斌李重光 . 参考光为平面波和球面波时焦距检测对比研究. 激光技术, 2017, 41(3): 421-426. doi: 10.7510/jgjs.issn.1001-3806.2017.03.023
    [8] 钱国林李朝明陈新荣邹文龙吴建宏 . 全息拼接光栅的误差研究. 激光技术, 2013, 37(6): 747-751. doi: 10.7510/jgjs.issn.1001-3806.2013.06.009
    [9] 简献忠张晗范建鹏周志刚 . 基于FPGA技术的计算全息研究. 激光技术, 2012, 36(1): 26-28,32. doi: 10.3969/j.issn.1001-3806.2012.01.008
    [10] 周效东汤伟中周文 . 全光纤全息技术的研究与实现. 激光技术, 1995, 19(2): 115-118.
    [11] 简献忠周海乔静远李莹王佳 . 基于全变差重构算法的数字全息研究. 激光技术, 2014, 38(2): 236-239. doi: 10.7510/jgjs.issn.1001-3806.2014.02.019
    [12] 戴欣冉钱晓凡徐天杰 . 同轴菲涅耳全息中提取相位的算法. 激光技术, 2014, 38(2): 172-176. doi: 10.7510/jgjs.issn.1001-3806.2014.02.006
    [13] 王晓惠刘超陆英仕楼宇丽 . 基于彩色数字全息的透镜焦距检测. 激光技术, 2015, 39(4): 562-565. doi: 10.7510/jgjs.issn.1001-3806.2015.04.030
    [14] 魏亚运周昕王梦婷呼有军 . 基于光学扫描全息测量相位物体的相位分布. 激光技术, 2016, 40(6): 820-824. doi: 10.7510/jgjs.issn.1001-3806.2016.06.010
    [15] 郑小丹刘东陶世荃王大勇 . 正交读出方式体全息光栅通信波长衍射特性. 激光技术, 2009, 33(2): 113-116.
    [16] 陈翠茹王华英赵宝群王学朱巧芬王杰宇王文健雷家良 . 基于UMnet的数字全息相位解包裹. 激光技术, 2023, 47(1): 73-79. doi: 10.7510/jgjs.issn.1001-3806.2023.01.011
    [17] 程阳 . 1维全息光子晶体的偏振特性. 激光技术, 2010, 34(2): 279-281. doi: 10.3969/j.issn.1001-3806.2010.02.037
    [18] 王霄季进清张惠中刘会霞 . 基于全息技术的激光透射塑料焊接研究. 激光技术, 2008, 32(5): 510-512.
    [19] 向东桂进斌刘超郑立婷楼宇丽宋庆和 . 数字全息波前准确重建的实验研究. 激光技术, 2017, 41(3): 406-410. doi: 10.7510/jgjs.issn.1001-3806.2017.03.020
    [20] 张庆生吕晓旭余清婷刘赣勇 . 基于灰度相关的数字全息孔径合成方法. 激光技术, 2009, 33(3): 276-278,282.
  • 加载中
图(6) / 表(2)
计量
  • 文章访问数:  5109
  • HTML全文浏览量:  3952
  • PDF下载量:  21
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-12-10
  • 录用日期:  2020-01-10
  • 刊出日期:  2021-01-25

基于时空域退偏的数字全息成像去噪研究

    通讯作者: 孔勇, kkyy7757@aliyun.com
    作者简介: 丁伟(1996-),男,硕士研究生,从事数字全息及数字全息显微方面的研究
  • 上海工程技术大学 电子电气工程学院 智能感知与控制系,上海 201620
基金项目:  上海市自然科学基金资助项目 19ZR1421700

摘要: 为了抑制全息记录过程中的散斑噪声,采用一种通过降低相干光源相干性的方法,并进行了理论分析和实验验证。利用旋转二元波片进行退偏,在10s内记录70幅全息图进行平均化降噪; 同时利用毛玻璃对相干光进行空间退相干,降低了光源的相干性从而降低了散斑噪声。将单独采用时间退偏操作的重建图像与单独采用空间退偏操作的重建图像以及两者相结合的重建图像进行了对比。通过实验验证了所提方法的可行性,通过分析频谱比较证实了几种方法,利用信噪比直观地证实了几种方法的有效性。结果表明,时间、空间退偏结合对散斑噪声的抑制有着明显的优势,空间退偏的重建像的信噪比为0.4459,时间退偏的重建像信噪比为1.5155,而两者结合后的信噪比为1.7162,采用时间退偏以及空间退偏叠加的方式有效地抑制了相干噪声。该实验为研究数字全息成像的过程中由于光源相干度过高带来的散斑噪声干扰的问题提供了一个有效的抑制方法。

English Abstract

    • 数字全息技术(digital holography, DH)由于其高效性和稳定性应用广泛,可用于3维动态成像,测量生物细胞及金属瑕疵等。通过物光与参考光的相互干涉产生的干涉条纹组成全息图,并由CMOS/CCD记录。而在数字全息中,记录全息图时会引入多种噪声,包括散斑噪声、拍摄噪声,以及额外的环境噪声[1-5]。相干噪声或散斑噪声是具有相干性的激光照射在粗糙物体表面引起的。拍摄噪声是由于光学器件的来回反射引起的多重噪声。而环境噪声则是由不完全的暗室环境引起的,如日光、照明反射灯光等。全息技术中散斑去噪的研究一直是一个重点。降低散斑噪声的方法一般分为3种[3]。第1种是对光源进行空间时间上的去相干操作。BASANTA等人利用白光光源提高了成像的质量[4]。JORGE等人引入LED作为光源,降低了数字全息显微(digital holographic microscopy, DHM)成像的噪声[6]。第2种是对全息图进行时域叠加复用。YUTAKA等人提出利用低相干光源根据不同位置物体的清晰细节不同拍摄多幅图并进行组合的方法,降低了散斑噪声[5]; FENG等人提出利用多模光纤记录多幅光程差不同的全息图,通过平均记录的全息图来降低散斑噪声[7]。第3种是利用数字处理的方法即利用各种图像处理的算法矩阵对全息图进行处理,通过软件处理图像获得降噪的效果。传统的掩模法将全息图的频谱分割成多个小块进行重建,最后进行叠加去噪。CAO等人针对电子散斑干涉中的条纹图的散斑噪声,提出一种改进非下采样轮廓波滤波算法,更好地恢复了干涉条纹的细节信息[8]。TAKAHIKO等人提出利用空域掩模将全息图分割之后进行重建,得到不同散斑噪声的重建图再进行叠加平均,有效抑制了1082pixel×1082pixel大小图像的散斑噪声[9]。JEONG等人提出将插值算法与平行相移相结合抑制散斑噪声,提高了信噪比[2]。然而上述的方法都有明显的缺陷。使用非相干光源较复杂且对记录的物体的大小及像素有着严格的要求; 使用多幅全息图去噪的方法需要记录多张的全息图,增加了噪声的不确定性,而算法处理进行图像去噪则比较耗时。

      近年来,学者们越来越多地关注退偏器对于光学成像质量的影响。光的不同偏振态对于非均匀物质的照射有不同的结果,这就使得退偏器可以应用于不同的领域,如医学上用于测量细胞等。有人曾提出利用双液晶片组合成新型的单色退偏器消除了激光驱动器系统的偏振化影响[10]。有人研究出旋转二元波片形成理想的夹角得到合适的退偏器[11]。有人在曾用计算机模拟偏振对激光全息的影响, 找出了平面内3束光形成最佳干涉图案的偏振状态[12]。还有人将旋光器与半波片的组合,并得到了很好的退偏效果[13]。进一步的,三元波片复合退偏器的理想效果也被验证[14]。上述方法的公式皆使用Muller矩阵和斯托克斯参量来进行分析[13-15]

      本文中提出一种通过降低相干光源的相干性来抑制散斑噪声的方法。利用旋转二元波片对光源实现时间退偏,利用毛玻璃对相干光进行空间退相干的操作降低了光源的相干性从而降低了散斑噪声。将单独采用时间退偏操作的重建图像与单独采用空间退偏操作的重建图像以及两者相结合的重建图像进行了对比。实验结果表明, 时间、空间退偏结合对散斑噪声的抑制有着明显的优势。

    • 数字全息的过程就是利用相干光源对物体进行照射形成物光场EO,同时引入同样的一束无物体的参考光场与物光进行干涉,形成干涉条纹,通过CMOS/CCD进行记录并利用计算机进行重建的过程。记录的全息图为:

      $I = {E_R}^2 + {E_O}^2 + {E_O}^*{E_R} + {E_R}^*$

      (1)

      式中, *代表共轭。而在数字全息的实验过程中,由于激光的高度相干特性以及照射物体表面的不规则,粗糙的表面会引起物体相位的骚乱,导致最后CMOS上接收的光信号携带散斑噪声。事实上,散斑噪声可以看成是许多随机相位点的集合,任一散斑点的强度可以表示为[16]

      $I_{\mathrm{SP}}=-2 \sigma^{2} \ln \left[2 \sigma^{2} P_{I}(I)\right]$

      (2)

      式中,PI(I)代表强度的概率密度,表达式为:

      $P_{I}(I)=\frac{1}{2 \sigma^{2}} \exp \left(-\frac{I}{2 \sigma^{2}}\right)$

      (3)

      式中, σ2为方差,I≥0。

      本文中重点研究对散斑噪声的去除,通过时间退偏以及空间退偏叠加的方式有效地抑制了相干噪声。接下来讨论具体的操作方法。

    • 图 1所示为时间退偏器的结构示意图。使用λ/4和λ/2波片调制物光,QWP(quarter wave-plate)和HWP(half wave-plate)分别为λ/4, λ/2波片。以MULLER矩阵的形式表示任意波片对光的调制[17]

      Figure 1.  Time domain depolarization device

      $\boldsymbol{M}=\left[\begin{array}{cccc}1 & 0 & 0 & 0 \\ 0 & \cos ^{2}(2 \alpha)+\sin ^{2}(2 \alpha) \cos \delta & (1-\cos \delta) \cos (2 \alpha) \sin (2 \alpha) & -\sin (2 \alpha) \sin \delta \\ 0 & (1-\cos \delta) \cos (2 \alpha) \sin (2 \alpha) & \sin ^{2}(2 a)+\cos ^{2}(2 \alpha) \cos \delta & \cos (2 \alpha) \sin \delta \\ 0 & \sin (2 \alpha) \sin \delta & -\cos (2 \alpha) \sin \delta & \cos \delta\end{array}\right]$

      (4)

      式中, α为波片的旋转角度,δ为相位调制量。λ/4波片的相位调制量为π/2;而λ/2波片的相位调制量为π。实际波片的相位调制量为[17]

      $\delta(\lambda)=\frac{2 \pi d\left(n_{\mathrm{e}}-n_{\mathrm{o}}\right)}{\lambda}$

      (5)

      式中, d为波片厚度,neno分别为波片的折射率及空气折射率。研究表明[3],当二元复合波片λ/4波片的角速度ω1λ/2波片的角速度ω2之比为1:2时可达到理论上的理想退偏, 即ω1:ω2=1:2。此时,二元波片的复合Muller矩阵[3]为:

      $\begin{aligned} \boldsymbol{M}_{t}=& \boldsymbol{M}_{\lambda / 4} \times \boldsymbol{M}_{\lambda / 2}=\boldsymbol{M}_{\delta=\pi / 2} \times \boldsymbol{M}_{\delta=\pi / 4}=\\ &\left[\begin{array}{llll}1 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0\end{array}\right] \end{aligned}$

      (6)
    • 研究表明,毛玻璃可以减少光学系统的衍射噪声[18]。假设通过毛玻璃后的光场为G(x, y),毛玻璃引入的相位为φg,入射光的相位为φ0,则G(x, y)的表达式为:

      $G(x, y)=G_{0} \exp \left[\mathrm{i} \varphi_{0}(x, y)\right] \exp \left[\mathrm{i} \varphi_{g}(x, y)\right]$

      (7)

      式中, G0为常数。

    • 本文中将时间退偏装置与空间退偏装置相结合,有效地降低了光源的相干性,并进而降低了全息图的散斑噪声。通过在物光侧将二元复合波片的装置与毛玻璃空间退偏装置相结合,对物光的相位进行了多重调制并抑制了散斑噪声。最终获得的相位结果为:

      $\varphi=\delta(\lambda) \varphi_{\mathrm{g}}$

      (8)
    • 本文中研究时空域退偏对于相干光的效果,采用邦加球进行仿真。邦加球由于其独特的偏振敏感特性, 一直以来用作分析偏振态的测试仪器。如图 2a所示,单色平面偏振光的斯托克斯矢量的邦加球表示为[19-20]

      Figure 2.  Poincare ball depolarization effect simulation

      $\boldsymbol{S}=\left[\begin{array}{cccc}S_{0} & S_{1} & S_{2} & S_{3}\end{array}\right]^{\mathrm{T}}= \\ \left[\begin{array}{cc}\mathrm{i} & \operatorname{icos}(2 \theta) \cos (2 \Phi) & \operatorname{icos}(2 \theta) \sin (2 \Phi) & \operatorname{isin}(2 \theta)]^{\mathrm{T}}\end{array}\right.$

      (9)

      式中, S1, S2, S3为球面上某点P的坐标,θΦ分别为对应球坐标的两个角度。邦加球的不同区域模拟不同的偏振态,坐标轴无量纲。当满足S02=S12+S22+S32时, 为全偏振光,代表球体表面的点; S02>S12+S22+S32, 为部分偏振光,代表在球体内部的点。如图 2b所示,邦加球上各点偏振态皆不同,邦加球水平最大圆周上的点表示线偏振光,从上极点到下极点右旋椭圆偏振光向左旋椭圆偏振光过度,球面为全偏振光,球心代表未经偏振调制的自然光,球体内为部分偏振光。采用二元波片旋转退偏与毛玻璃空间退偏时会产生任意偏振态的光,如图 2c所示,其中部分偏正光的坐标点未标出。

    • 图 3为所提方法的全息退偏装置。其中激光源为5mW,532nm绿光激光器,BS(beam spliter)为分光器,OBJ(objective)为20×倍率显微镜用于扩束,A为小孔,用于小孔滤波,QWP为λ/4波片,HWP为λ/2波片,两波片以ω1:ω2=1:2的角速度旋转, L为透镜,GG为毛玻璃,D为骰子,近似一个边长1cm的圆角立方体。CMOS用于记录全息图,其分辨率为1280pixel×1022pixel且像素大小为5.3μm×5.3μm,记录距离为45cm。将波片去除时为空间退偏装置,将毛玻璃去掉时为时间退偏装置。实验结果如图 4所示。

      Figure 3.  Proposed method holography setup

      Figure 4.  Comparison of reconstruction maps of several methods

      图 4中将上述提到的几种方法的实验结果进行了对比。图 4a为初始全息图重建,图 4b为空间退偏重建图,图 4c为时域退偏重建图,图 4d为时域空域退偏重建图。显然,图 3b空域退偏的效果最差,毛玻璃将散斑打乱,于是导致了整体质量的下降,同时将部分高频信息淹没,导致物体的边缘信息难以感知,但0级像的噪声显著降低,验证了毛玻璃降低相干度的作用。图 4c中采用时域退偏,拍摄的时间为10s,期间一共拍摄了70幅全息图,得出一个叠加平均的全息图。可以看到,对比初始全息图,时域退偏的重建图无论是±1级像还是0级像的噪声都得到了大幅度的降低。图 4c比起图 4b,其全局噪声也得到抑制,同时散斑噪声也远小于图 4a中所示。图 4d图 4c的基础上添加毛玻璃进行时域空域叠加退偏,起到了进一步抑制噪声的作用,且从感知质量上来看,图 4d的去噪效果是最好的。用图像的信噪比来评价去噪的效果。信噪比定义为:

      $R_{\mathrm{SNR}}=10 \lg \left\{\frac{\sum\limits_{x=1}^{N} \sum\limits_{y=1}^{N}\left[f_{0}(x, y)-\hat{f}_{O}(x, y)\right]^{2}}{\sum\limits_{x=1}^{N} \sum\limits_{y=1}^{N}\left[f_{\mathrm{p}}(x, y)-\hat{f}_{\mathrm{p}}(x, y)\right]^{2}}\right\}$

      (10)

      式中, fO为原始图像的像素值, $\hat{f}_{o}$为原始图像的像素平均值,fp为处理后的像素值, $\hat{f}_{\mathrm{p}}$为处理后的像素平均值,N为像素数。计算所得信噪比记录在表 1中。

      Table 1.  Signal-to-noise ratio of different methods

      spatial-domain
      depolarization
      time-domain
      depolarization
      spatiotemporal
      depolarization
      RSNR/dB 0.4459 1.5155 1.7162

      为进一步验证上述方法的去噪效果,取各实验结果的频谱图进行了分析。图 5a~图 5d分别为图 4a~图 4d的频谱图,其中x, y轴为像素坐标,z轴为灰度值。

      Figure 5.  The 3-D spectral distribution of Fig. 4a~Fig. 4d

      显然,对比图 5b~图 5d的频谱和图 5a中的频谱可以看出,上述方法都在一定程度上抑制了全息散斑噪声。以频谱的灰度值作为分析指标,从0级像可以直观地得出,空间退偏以及时间退偏对于0级像的散斑抑制效果比较明显,两种方法都将±1级像与0级像之间的噪声抑制了一定的程度; 同时,也将±1级像中的噪声也抑制了一部分,将时空退偏效果叠加之后可以看到,时空域退偏将噪声抑制了一部分也保留了许多+1级像中的物体信息。此外,实验过程中的噪声抑制与物体信息的丢失是同步的,因此无法精确地比较噪声抑制的同时物体信息的损失,对此,通过感知质量比较最后的去噪效果,并在之后采用信噪比来直观地比较时空退偏的效果。

      为了直观地比较各种方法,进一步对频谱进行2维曲线的分析。图 6a~图 6d分别为图 5a~图 5dx=0时的截面图,其中x轴为像素坐标,y轴为灰度值。

      Figure 6.  a~d sectional view of Fig. 5a~Fig. 5d, respectively

      图 6中可以直观地看出不同方法降噪的效果。频谱图降噪的评价指标为±1级像与0级像边缘之间的空隙质量。显然,从图 6a可以看出,初始全息图在记录过程中会受到多重噪声的影响,±1级像频谱与0级像频谱之间包含许多噪声,导致了如图 4a所示的模糊的重建像; 图 6b为空域退偏的效果,可以看到,±1级像与0级像大大分离,这是利用毛玻璃降低了激光的相干度, 从而抑制了散斑噪声,使得噪声下降,但是缺陷在于由于毛玻璃的不均匀,导致噪声去除的不理想,如图中品红色虚线框所示; 图 6c为时域退偏效果,经过多幅全息图的记录平均降低了噪声,但同时也会将少量高频细节叠加,会产生图中品红色虚线框所示的毛刺,观察±1级图像可以看出,此方法在去噪的同时也去除了许多物体信息; 图 6d为时空域退偏的效果,可以看出,比起初始全息图,此方法显著地去除了噪声,对比图 6b图 6c, 此方法保留了许多物体信息。

      表 2所示,将图 6中的4种情况的参量全部都记录下来, 便于直观地比较不同方法的降噪效果。

      Table 2.  Parameters of noise suppression of different methods

      number of pixels between-1 and 0 image edges in x and y directions number of pixels between+1 and 0 image edges in x and y directions
      original hologram spectrum Δx=59, Δy=140 Δx=56, Δy=143
      spatial filtering spectrum Δx=180, Δy=172 Δx=182, Δy=176
      temporal filtering spectrum Δx=163, Δy=188 Δx=162, Δy=187
      spatiotemporal filtering spectrum Δx=85, Δy=164 Δx=112, Δy=168
    • 将时间域的退偏域空间域的退偏相结合,对数字全息的散斑噪声进行了抑制并获得了显著的好效果。通过毛玻璃对物光进行空间上的退偏,同时利用二元复合波片系统进行时域退偏,电机控制旋转波片进行时间退偏,记录了70幅全息图平均噪声,结合空间退偏有效地抑制了散斑噪声,相比于之前单独采用时间退偏或空间退偏的装置,信噪比得到提高。时空退偏的重建图信噪比达到了1.7162,相比空间退偏的重建像的信噪比为0.4459,时间退偏的重建像信噪比为1.5155, 有了相对的提高。此研究对于今后的数字全息研究具有促进的意义。

      下一步研究的重点是:将时空退偏的方法应用到数字全息显微术中,以及利用去噪算法获得更好的降噪效果。

参考文献 (20)

目录

    /

    返回文章
    返回