基于旋滤波的散斑干涉图预处理方法.pdf

　第 2 1 卷第 2 期 强 激 光 与 粒 子 束 Vol . 2 1 , No . 2 　 　2009 年 2 月 H I GH PO W ER L A SER AND PA R TICL E B EA M S Feb . , 2009 　 文章编号 : 　100 14322 (2009) 基于旋滤波的散斑干涉图预处理方法 王 　峰, 　彭晓世 , 　刘慎业 , 　蒋小华 , 　丁永坤 ( 中国工程物理研究院 激光聚变研究中心, 四川 绵阳 62 1900) 　　摘 　要 : 　介绍了旋滤波算法中窗口大小变化对滤波效果的影响 ,得出旋滤波的窗口应该取在约 1/ 4 条纹 宽度的结论 。处理了大散斑噪声的干涉图 ,发现窗口尺寸必须超过散斑尺寸才能得到比较好的滤波效果 。确 定了旋滤波算法应用的范围是大于散斑尺寸 ,且在 1/ 4 条纹宽度处比较合适 。针对滤波后条纹对比度降低的 问题 ,使用灰度拉伸的方法进行对比度增强 。 　　关键词 : 　旋滤波 ; 　散斑噪声 ; 　中值滤波 ; 　干涉图像 ; 　预处理 　　中图分类号 : 　TN0 12 　　　　文献标志码 : 　A 　　现有的干涉条纹处理方法基本上可以分为两类 :一类是利用多幅干涉条纹进行处理 , 比较有代表性的是相 移法[ 12 ] ;另一类是利用单幅干涉条纹进行处理 , 比较有代表性的是傅里叶变换法[34 ] 。傅里叶变换法是干涉条 纹分析中一种比较常用的方法 ,它只需分析单一的一幅干涉图 ,适于处理动态和静态波面 ,且精度高 。在傅里 叶变换法中 ,噪声是通过选择适当的滤波窗来去除的。使用越窄的带通滤波器 ,就可以获得越平滑的相位 。然 而窄带滤波器也会使信号相位失真 ,虽然可以通过精确设计滤波器来减少这种失真 ,但这降低了条纹处理的自 动化程度 。这是该方法面临的最大的问题 ———频率泄露 。而且有些噪声的频谱比较宽 ,使得噪声频率和条纹 频率重叠在一起 ,因此 ,在频域中无法彻底去除这些噪声 。由于实际拍摄的干涉图经常有较严重的随机噪声 , ( ) 因此进行傅里叶变换分析之前必须对干涉图进行预处理 。一般抑制相干的随机噪声有两种方法 : 1 通过旋转 漫射器削弱光源的空间相干性[ 5 ] ,使用非相干光源[ 6 ] 或用其它设备移动光源[7 ] 等方法 ; (2) 干涉条纹预处理 。 本文将使用干涉条纹预处理中常用的旋滤波法对两种大小不同的散斑干涉图像进行处理 ,研究该方法的适用 性 ,从而确定旋滤波法窗口选择的范围。 1 　旋滤波法的基本原理 　　在条纹图的中间区域 ,条纹图的灰度值梯度分布有一定的规律 :在条纹的切线方向梯度值最小 ,而在条纹 ( ) ( ) 的法线方向梯度值最大 。在频域内分析此规律 , 图 1 a 和图 1 b 分别是条纹图中某一窗口内在条纹法线和切 ( ) ( ) ( ) ( ) 线方向的灰度分布 ; 图 1 c 和图 1 d 分别是图 1 a 和图 1 b 的频谱分布图。在条纹法线方向 ,灰度变化较 大 ,对应的频谱呈一条宽带 ,条纹信号和噪声的频谱是叠加在一起的 ,无法截然分开 。在条纹切线方向 ,灰度变 化很小 ,对应的频谱是在零频附近的一条窄带 ,而随机噪声仍然分布在高频带 , 因此可以用常规低通滤波器将 条纹和噪声分开 ,从而滤掉噪声 , 同时不损失条纹信号的任何信息 。旋滤波的基本思想是 :首先找到条纹灰度 值的等值线 ,然后在灰度值等值线上作低通滤波