非局部均值滤波

1、 图像去噪技术图像去噪技术 图像噪声图像噪声 n图像噪声：图像在摄取时或传输时所受到 的随机干扰信号。 椒盐噪声 高斯噪声 泊松噪声 n图像噪声 n空间域： 均值滤波、高斯滤波、形态学滤波、局 部滤波和非局部滤波等 n频域： 维纳滤波 和小波阈值收缩等 图像去噪的基本方法图像去噪的基本方法 均值滤波 原理 在图像上，对在图像上，对待处理的像素待处理的像素给定一个模板，给定一个模板， 该模板包括了其周围的邻近像素。将模板中该模板包括了其周围的邻近像素。将模板中 的全体像素的的全体像素的均值均值来替代原来的像素值的方来替代原来的像素值的方 法。法。 n非局部均值（NL-means）是近年来提出的

2、一项新型的去噪技术。该方法充分利用了图 像中的冗余信息，在去噪的同时能最大程度 地保持图像的细节特征。 n特点：算法简洁，性能优越，易于改进和扩 展 n基本思想是：当前像素的估计值由图像中与 它具有相似邻域结构的像素加权平均得到。 ( )( )( )v iu in i ( )v i ( )n i( )u i 2 2, 2 | ()()| ( , )exp() ij v Nv N w i j h 我们规定 为以i为中心的矩形邻域，那 么图像v中的像素i和像素j的相似性系数 为： ( , )w i j () i v N w(P,P2) = 1.4334e-11 w (P,P1) = 7.6567e

3、-04 Circuit 原始图像加入椒盐噪声图像 3X3均值滤波后的图像NLM滤波后的图像 Circuit 原始图像 加入高斯噪声图像 3X3均值滤波后的图像 NLM滤波后的图像 非局部均值去噪法存在的问题非局部均值去噪法存在的问题 相似性度量缺乏鲁棒性 高斯加权核各向同性性质影响 非相似像素块影响 运算量大 加权核系数选择 假设图像共N个像素点，搜索窗大小定义为D，邻域 窗口大小定义为d，计算每个矩形邻域间相似度的时 间为 ,对于每个像素点要计算他的搜索窗内 个像素点的相似度，故NL_Mean filter的复杂度为 2 ()O d 2 D 22 ()O ND d 基于积分图像的快速NLM

4、对图像整体处理，原图像与平移 后的 图像的欧氏距离为： 2 2 ( )() t sv xv xt ,tyxDs Ds 2 121122 12 (,):0,0 ( )( ),( ,)dt t zz zNzxzx Sxs z xx x 1212 1212 ( ,),1,1, ( )( )(1,)( ,1)ddd tttt xx xxx Sxs xSxxSx x 12 (1,1)dtSxx 如果我们先构造一个关于像素差值的积分图像： 上式在实际操作中可表达为： 那么对于不同区域的欧式距离可以写为： 2 1212 22, 1 ( )( )(,)(,)dd tsstss d v xv ySxd xdSx

5、d xd d 1212 (,)(,)dd tsstss SxdxdSxdxd 此时对于N个像素点的图像，搜索窗大小为D，计 算NL_Mean filter的复杂度为 2 ()O ND 积分图像算法的缺点积分图像算法的缺点： 它不直接允许的计算使用一个内核K加权范数,如 高斯。 当图像尺寸以及补丁距离很大, 积分图像的一些值 可能变得很大,即使使用双精度表示，最终结果的 准确性可能下降。 给出距离为t的两个相似区域的2范数的离散卷积 形式： 基于快速傅里叶变换的NLM加速算法 2 2 2 2, 2 : ( )()( )() s K z Zzd v xv xtK z v xzv xtz * t K

6、 s ( )()K zKz其中 表示卷积算子， * , 2 2 ( )() t sv xv xt 卷计算法可以用傅里叶变换求解。将上式做快速 二维傅里叶变换得到： 2 2, 1 ( )()( ) ( ) t K v xsv xtK FFF 128X128256X256512X512 原始算法计算时间130.54667.872743.24 FFT计算时间10.1144.77183.12 提高倍数12.9114.9214.98 表1. 原始算法与使用FFT加速算法计算时间比较 （单位：秒） pout 原始图像 h=10 滤波后图像 h=20 滤波后图像 h=2 滤波后图像 原始图像原始图像噪声图像

7、噪声图像NLM filterMean filter 0.81.42 10 pout25.25630.94636.194034.30034.532 Circuit 47.80744.22736.38435.330 32.944 20 pout30.65632.08032.38231.33629.358 Circuit 50.85741.04432.05830.92328.685 40 pout46.36127.96929.29629.13423.565 Circuit 61.54337.58428.13927.14223.434 n 表1. 去噪结果后PSNR比较 图图1. 不同噪声强度下不同平滑参数不同噪声强度下不同平滑参数h滤波效果滤波效果 表2. NLM 滤波优先参数选择参考 0,7131.5 7,9141.4 9,19151.3 20,28261.1 28,47371.0 47,70381.0 70,87381.0 s