白光衍射相位像光晕效应消除方法 下载: 1142次
1 引言
衍射相位显微术[1]是一种基于光栅衍射的共路离轴定量相位显微成像技术,可以精确地获取样品的高分辨率相位信息。按采用的照明光源的不同,相位成像系统可以分为激光衍射相位成像系统[2-3]和白光衍射相位成像系统[4-5]。后者采用白光作为照明光源,其相干长度非常短,可以很好地避免由光源强相干性导致的相干噪声,因此比前者具有更高的空间敏感性[6-7]。然而,白光衍射相位成像系统采用扩展的白光光源照明,照明光波不可能是完全的平面光波,它的空间相干面积一般远小于测量视场,从而使得相位像中出现与物体结构相依赖的光晕效应[8-9],阻碍对样品相位信息的准确测量。不仅是白光衍射相位成像系统,其他采用白光的相位成像技术,如Zernike相衬成像技术[10-11]、SLIM(Spatial Light Interference Microscopy)相位成像技术[12-13]、白光照明光强传输方程相位成像技术[14-15],均存在光晕效应影响。目前解决光晕效应影响的方法主要有两大类。一类是提高照明光源的空间相干性[8-9,16-18],一般通过减小照明光源的数值孔径,提高照明光源的空间相干面积来实现,但是这种方法会造成光通量降低,导致曝光时间增加,不利于实现动态测量。另一类方法是数值图像处理[19-21],其中最有效的方法是Nguyen等[19]提出的迭代数值重建方法,该方法在建立的部分相干相位成像物理模型基础上,结合硬件系统参数,采用最小二乘迭代算法,很好地消除了光晕效应,恢复了无光晕像,但是迭代运算费时,不利于实现实时消除光晕效应,且依赖于系统硬件参数。
本文基于白光衍射相位成像光晕效应的形成原理,提出了一种基于希尔伯特变换的新方法,可以快速有效地消除光晕效应。该方法不需要样品的任何先验知识,只需对图像进行一次希尔伯特变换,并混合正确的高频数据和低频数据,就可以有效地消除光晕效应,保留白光衍射相位成像系统的高分辨率信息。将该方法应用于实验获取的白光衍射相位像数据,包括标准微球数据和人体活体血红细胞数据。实验结果表明,提出的方法可以有效快速地消除白光衍射相位像的光晕效应。
2 白光衍射相位像光晕效应的形成原因
系统测量项为物光场和参考光场的互相干函数。由于物光场和参考光场经过相同的光学元器件,两光场的时间延迟为0,即τ=0,因此被测互相干函数可表示为
式中:Us表示物光波场函数;Ur表示参考光波场函数;<>t表示时间平均;r表示空间坐标;*表示共轭运算。假定光场是各态历经的稳态光场,其时间平均不随时间改变,即为时间统计平均。根据广义的维纳-欣钦定理,互相干函数Γs,r(r,r,τ)是交叉谱密度函数Ws,r(r,r,w)的傅里叶变换(w为时间角频率),则Γs,r(r,r,τ)又可表示为
式中:Us(r,w)、
根据光栅的衍射特性,0级光与+1级光都携带了样品的全部信息。滤波器允许+1级光全部通过矩形孔,则物光光场为
式中:Ui(r,w)表示被测样品的照明光场;T(r)为被测样品的透射率函数。
0级光经过针孔滤波后作为参考光场,参考光场函数的表达式为
式中:h0(r)为0级滤波孔径透过率函数的傅里叶变换;⊗表示二维卷积。将(3)式和(4)式代入(2)式得
式中:h(r)=
式中:φm(r)为测量相位;φ(r)为样品的真实相位;arg[(T⊗h)(r)]为平滑后样品透射率函数的相位,是一项低频项。正是该低频项导致样品测量相位值比真实相位值低,并使测量相位像的周边出现负值,这种现象称为光晕效应。
3 光晕效应实时消除方法
白光衍射相位成像系统中测量相位和真实相位的关系如(6)式所示。对(6)式进行求导,由于arg[(T⊗h)(r)]是一个低频函数,它的空间导数很小,因此可以近似忽略。以x维为例,求导可以得到
由(7)式可知,测量相位的导数近似等于真实相位的导数。将(7)式进行希尔伯特变换,可得
式中:P表示取x=x'时的柯西主值。由(8)式可知,对测量相位的导数进行希尔伯特变换得到的是真实相位与
数值处理时为了提高计算速度,可以通过傅里叶变换操作来进行希尔伯特变换,即
式中:kx为x方向的空间频率;ky为y方向的空间频率。
为了保留原白光衍射相位成像系统的高分辨率,本文将相位像中大于相干长度频率的空间频率完全保留,同时因为成像系统是圆对称的,(9)式又可以写成
式中:kr=
图 2. 基于希尔伯特变换的消除光晕效应算法示意图
Fig. 2. Schematic of the method of eliminating halo effect based on Hilbert transform
4 实验结果及分析
为了证实提出的基于希尔伯特变换的方法能准确消除白光衍射相位像的光晕效应,首先用该方法处理实验测量的标准聚苯乙烯微球像。微球直径为(2±0.05) μm,聚苯乙烯折射率为1.59。滴入奥林巴斯物镜浸油(折射率为1.518)浸没微球,并盖上盖玻片。将此样品放入载物台,调焦使之成像于CCD相机平面,记录的干涉图如
图 3. 标准聚苯乙烯微球实验结果。(a)测量干涉图;(b)测量厚度像;(c)利用所提方法重建的厚度像;(d)利用迭代方法重建的厚度像;(e)沿图3 (b)~(d)实线测得的微球厚度分布图
Fig. 3. Experimental results of standard polystyrene microspheres. (a) Measured interferogram; (b) measured thickness image; (c) reconstructed thickness image by the proposed method; (d) reconstructed thickness image by the iterative algorithm; (e) thickness profiles of the same standard polystyrene microsphere along the solid lines in Fig.3 (b)-(d)
对测量的血红细胞白光衍射相位像进行处理,结果如
图 4. 血红细胞实验结果。(a)测量得到的干涉图;(b)测量相位像;(c)利用所提方法重建的相位像;(d)沿图4 (b)和4(c)中实线测得的相位分布图
Fig. 4. Experimental results of red blood cells. (a) Measured interferogram; (b) measured phase image; (c) reconstructed phase image by the proposed method; (d) phase profiles of the red blood cell along the solid lines in Fig. 4 (b) and Fig. 4 (c)
5 结论
本文针对低相干光源的衍射相位成像系统存在的光晕效应问题,提出了一种基于希尔伯特变换的新方法。利用该方法处理了实验测量的标准聚苯乙烯微球和人体活体血红细胞的相位图,实验结果证明所提方法能有效地消除样品周围的光晕效应。虽然使用的初始相位图是基于白光衍射相位成像系统测得的,但所提方法也可以应用于其他受光晕影响的白光系统。该方法能有效快速地消除光晕效应,适用于实际无光晕白光相位成像。
[2] Popescu G, Ikeda T, Dasari R R, et al. Diffraction phase microscopy for quantifying cell structure and dynamics[J]. Optics Letters, 2006, 31(6): 775-777.
[3] Majeed H, Ma L H, Lee Y J, et al. Magnified Image Spatial Spectrum (MISS) microscopy for nanometer and millisecond scale label-free imaging[J]. Optics Express, 2018, 26(5): 5423-5440.
[4] Bhaduri B, Pham H, Mir M, et al. Diffraction phase microscopy with white light[J]. Optics Letters, 2012, 37(6): 1094-1096.
[5] 张佳恒, 马利红, 李勇, 等. 卤素灯照明光栅衍射共路数字全息显微定量相位成像[J]. 中国激光, 2018, 45(6): 0609003.
[7] Farrokhi H, Boonruangkan J, Chun B J, et al. Speckle reduction in quantitative phase imaging by generating spatially incoherent laser field at electroactive optical diffusers[J]. Optics Express, 2017, 25(10): 10791-10800.
[10] Zernike F. How I discovered phase contrast[J]. Science, 1955, 121(3141): 345-349.
[11] BornM, WolfE. Principles of optics: electromagnetic theory of propagation interference and diffraction of light[M]. Britain: Cambridge University Press, 1999.
[12] Majeed H, Nguyen T H, Kandel M E, et al. Label-free quantitative evaluation of breast tissue using Spatial Light Interference Microscopy (SLIM)[J]. Scientific Reports, 2018, 8: 6875.
[13] Min E, Kandel M E, Ko C J, et al. Label-free, multi-scale imaging of ex-vivo mouse brain using spatial light interference microscopy[J]. Scientific Reports, 2016, 6: 39667.
[14] Komuro K, Nomura T. Object plane detection and phase-amplitude imaging based on transport of intensity equation[J]. Optical Review, 2017, 24(5): 626-633.
[15] 张璐, 汤其剑, 邓定南, 等. 基于光强传输方程的双相机动态相位成像的视场校正[J]. 中国激光, 2019, 46(8): 0804005.
[16] Maurer C, Jesacher A, Bernet S, et al. Phase contrast microscopy with full numerical aperture illumination[J]. Optics Express, 2008, 16(24): 19821-19829.
[17] Zuo C, Sun J S, Li J J, et al. High-resolution transport-of-intensity quantitative phase microscopy with annular illumination[J]. Scientific Reports, 2017, 7: 7654.
[19] Nguyen T H, Kandel M, Shakir H M, et al. Halo-free phase contrast microscopy[J]. Scientific Reports, 2017, 7: 44034.
[20] Kandel M E, Fanous M, Best-Popescu C, et al. Real-time halo correction in phase contrast imaging[J]. Biomedical Optics Express, 2018, 9(2): 623-635.
Article Outline
丁海宁, 朱苗苗, 马利红, 李勇. 白光衍射相位像光晕效应消除方法[J]. 中国激光, 2020, 47(5): 0509001. Haining Ding, Miaomiao Zhu, Lihong Ma, Yong Li. Elimination Method of Halo Artifacts in White-Light Diffraction Phase Microscopy[J]. Chinese Journal of Lasers, 2020, 47(5): 0509001.