一种无透镜傅里叶变换数字全息的散斑降噪方法 下载: 1109次
1 引言
数字全息术可以快速记录物体的强度和相位信息,便于实现多种定量测量,再现时可通过数字图像处理技术,优化再现像的质量[1]。随着电荷耦合器件(CCD)等光电探测器件的快速发展,数字全息的研究也愈受关注[2]。基于其独特的优势,数字全息已经被广泛运用到多个领域,如:全息显微[3-4]、信息加密[5-6]、粒子测量[7-8]、3D显示[9]等。
当相干光束照射在物体的粗糙表面时,在不同表面点散射的光波会发生波动,这些散射波的叠加会在记录平面上形成一个固定的斑点图案[10],从而严重降低再现像的信噪比和空间分辨率。
目前,研究人员已提出了多种降低数字全息再现像散斑噪声的技术。例如:在数字信号处理方面,Uzan等[11-12]提出非局部均值(NLM)滤波法,通过多幅图像的加权平均强度来抑制散斑噪声,但该方法计算相似度的度量缺乏稳健性,且运算十分复杂;在低相干的方法上,Kemper等[13-14]提出用部分相干光源代替激光光源可以达到降低激光散斑的效果,但这要依赖于高亮度部分相干光源的采用;在复用全息图处理方面,Huang等[15-16]提出了合成孔径数字全息法,通过合成多个记录的全息图来增加全息图的尺寸,从而降低散斑的尺寸,但该方法也需要记录多幅全息图,不适合用于实时处理。
此外,Mori等[17-18]提出了一种通过空间域掩模(SDM)法抑制再现像的散斑强度来减小空间分辨率损失的方法。原始干涉图通过带有移动孔径的二值掩模进行过滤,得到多幅子干涉图。由于不同的子再现像是从干涉图的不同区域重建获得的,因此携带有不一样的散斑模型。将所有的子再现像叠加平均后可以有效降低散斑噪声,但该方法没有选择性地对所有的子再现像进行叠加,从而造成了图像信息的过度冗余,而且运行速度慢。
本文提出了一种基于无透镜傅里叶变换数字全息特性的方法,先将原始干涉图划分为多个区域,通过改变各区域的位置来获得多幅新干涉图,并对其进行再现叠加和平均,就可以得到散斑噪声更低的全息图。
2 基本原理
如
图 1. 无透镜傅里叶变换数字全息图记录光路图
Fig. 1. Lensless Fourier transform digital hologram record optical path diagram
根据菲涅耳衍射公式可知,干板处的光振幅由物光波和参考光波两部分组成[19],它们分别表示为
式中:F{}为快速傅里叶变换;
用发散球面波
从(4)式可以看出,无透镜傅里叶变换数字全息术可以通过直接对全息图进行傅里叶变换实现数字再现,从而使数字全息的再现得以简化。在再现光场中,除了零级光场外,还会得到清晰的原始像和共轭像。
3 数字全息降噪方法
3.1 空域掩模降噪法
如
式中:
式中:
3.2 干涉图重组降噪法
设
式中:rect(
由傅里叶变换定理可知,函数在空域中的平移会带来频域中的相移,即
所以,当
式中:
式中:
将干涉图划分为
式中:
4 实验分析
采用MATLAB软件对CCD记录的干涉图进行再现处理,得到的再现像如
再现像的降噪效果可以通过等效视数(ENL)来评估。等效视数是评价同质区域降噪质量的一个很好的指标,其值随噪声降低而增加。其表达式为
式中:
图 4. 无透镜傅里叶变换数字全息光路图
Fig. 4. Lensless Fourier transform digital holography optical path diagram
图 5. 无透镜傅里叶变换数字全息再现像。(a)硬币;(b)骰子
Fig. 5. Lensless Fourier transform digital holographic reconstructed images. (a) Coins; (b) dices
在本实验中,对于SDM方法,掩模窗口
经两种方法处理后,再现像的等效视数随着幅数的增加而改变,如
图 6. 等效视数与子再现像幅数的关系。(a)骰子;(b)硬币
Fig. 6. Equivalent number of looks changes with the increase of the number of sub-reconstructed images. (a) Dice; (b) coin
在
相比于SDM方法,所提方法能更快地提高再现像的等效视数,而且在子再现像幅数约为80时开始趋于稳定。此外,在SDM方法中,再现像等效视数的速度提高得较慢,耗时长,当子再现像幅数持续增长到一定值时,其等效视数也会比所提方法的等效视数高。但在这种情况下,SDM方法处理的再现像因将过多图像叠加而导致信息过度冗余,从而出现曝光现象。
在
两种方法处理再现像所需时间如
图 7. 所提方法与SDM方法得到的再现像。(a)所提方法,骰子;(b)所提方法,硬币;(c) SDM方法,骰子;(d) SDM方法,硬币
Fig. 7. Reconstructed images using proposed method and SDM method. (a) Proposed method, dices; (b) proposed method, coins; (c) SDM method, dices; (d) SDM method, coins
表 1. 两种方法的处理时间
Table 1. Processing time of two methodss
|
由
边缘保持系数(EPI)能表示图像经优化算法处理之后,垂直或水平方向上图像边缘的保持能力[20],其表达式为
式中:
表 2. 两种方法优化的再现像的边缘保持系数
Table 2. Edge preservation index of reconstructed image improved by two methods
|
由
5 结论
数字全息再现像中含有不可避免的噪声,这会严重影响再现像的质量。无透镜傅里叶变换数字全息只需进行一次傅里叶变换,就可以重建物光信息。本课题组基于无透镜傅里叶变换数字全息的特点,提出了一种新的降噪方法。该方法通过划分原始干涉图,随机摆布各个区域的位置,获得了多幅不同的干涉图,对这些重组的干涉图再现并进行累加平均,就可以得到噪声更低的再现像,其边缘保持系数明显提高。相比于SDM降噪方法,所提方法在MATLAB平台上处理的时间缩短了80%,再现像的边缘保持系数是SDM方法的2倍。此外,所提方法只需一幅干涉图,适用于实时处理。
[3] 曾雅楠, 卢钧胜, 刘源, 等. 基于数字全息显微的微粒三维位移跟踪测量技术[J]. 中国激光, 2017, 44(12): 1204001.
[6] 吴为民, 陈家祯. 基于双快速反应码双随机相位加密系统的多重水印方法[J]. 激光与光电子学进展, 2017, 54(12): 120901.
[8] 常新宇, 曾雅楠, 雷海, 等. 液态中交叠微球的三维纳米级位置测量方法[J]. 光学学报, 2017, 37(8): 0812002.
[16] 钟丽云, 张以谟, 吕晓旭. 合成孔径数字全息的分析模拟及多参考光合成孔径数字全息[J]. 光子学报, 2004, 33(11): 1343-1347.
[19] 陈家璧, 苏显渝. 光学信息技术原理及应用[M]. 北京: 高等教育出版社, 2009: 143- 144.
Chen JB, Su XY. Principles and applications of optical information technology[M]. Beijing: Higher Education Press, 2009: 143- 144.
[20] 张微, 孙蓉桦, 章孝灿. 基于改进的小波软阈值法的SAR图像去噪[J]. 遥感信息, 2004, 19(4): 4-6.
Article Outline
梁明大, 陈丽, 林伟涛, 陈永昊. 一种无透镜傅里叶变换数字全息的散斑降噪方法[J]. 激光与光电子学进展, 2018, 55(11): 110901. Mingda Liang, Li Chen, Weitao Lin, Yonghao Chen. A Speckle Noise Reduction Method for Lensless Fourier Transform Digital Holography[J]. Laser & Optoelectronics Progress, 2018, 55(11): 110901.