激光与光电子学进展, 2018, 55 (1): 011003, 网络出版: 2018-09-10   

基于模糊边缘补足的自适应非局部均值图像去噪算法 下载: 877次

Self-Adaptive Non-Local Means Image Denoising Algorithm Based on Fuzzy Edge Complement
作者单位
中国矿业大学信息与控制工程学院, 江苏 徐州 221116
摘要
针对传统非局部均值图像去噪算法对纹理细节处理保持不足的问题,提出一种结合模糊边缘补足(FEC)的非局部均值图像去噪算法。利用FEC算法检测出图像的边缘纹理特征图像;基于此自适应选择相似性权重参数,有针对性地对边缘纹理区域和平坦区域进行不同程度的平滑,以防止边缘细节信息的过平滑;用边缘的结构相似性改进非局部图像块的相似性权重,加强相同区域像素的贡献而削弱不同区域像素的贡献,从而达到更好地保持纹理细节的作用。实验结果表明:该方法能够很好地达到去除噪声的效果,同时还能保持更多的纹理细节特征和几何结构特征。
Abstract
In view of the holding unsatisfactory effects of the traditional non-local means algorithm for texture details, an improved non-local means denoising algorithm combined with fuzzy edge complement (FEC) is proposed. The edge texture feature image is detected by the FEC algorithm. The similarity weight parameter is chosen adaptively according to the edge texture feature, and the edge texture region and the flat region are smoothed for different degrees pertinently, which prevents the edge texture information from being lost. The similarity weights of non-local image blocks are improved using the structural similarity of edges. The effects of pixels in the same area are increased, and those in different areas are reduced. Thus, the better texture hold effects can be achieved. Experimental result indicates that the image denosing can be effectively achieved by this method. Meanwhile, the more texture detail features and geometrical structural features are persevered.

1 引言

随着社会的进步和科技的发展,人们对图像质量的要求越来越高。但是,在采集和传输图像的过程中难免会受噪声的污染,因此去除图像噪声是图像处理领域内一个重要的课题。图像去噪的主要目的是要尽可能多地保留图像细节特征并尽可能多地去除噪声。传统的图像去噪方法主要分为空间域去噪法和变换域去噪法,其中中值滤波和双边滤波[1]是空间域局部去噪的经典方法,但他们都容易使图像模糊并丢失许多边缘纹理信息。基于小波变换去噪法是变换域去噪领域经典的去噪方法,但是变换域滤波的方法均依赖于固定的变换基,如曲波基、轮廓波或小波基等[2],并且仅能表示图像中的某些几何特征。基于字典学习[3]的图像去噪方法是变换域滤波领域的突破,其利用图像块字典作为变换基,使去噪效果得到提升。非局部均值(NLM)图像去噪[4]方法是空间域去噪方法的一大突破,其将结构特征引入到图像去噪中,由于其去噪效果良好、算法简单且易于拓展而得到广泛关注。

目前已经存在较多改进的NLM算法。Kervrann等[5]提出一种优化的空间自适应NLM方法,其基于局部方差度量图像片的相似性,并根据噪声水平自适应选取相似权重参数以及基于统计信息自适应选取搜索区域的大小,取得了很好的去噪效果。文献[ 6]提出了一种贝叶斯NLM滤波器,其采用局部自适应字典与噪声统计分布对图像之间的相似性进行衡量并做了相应的改进。Grewenig等[7]将主成分分析引入到图像块相似性度量中,选取一定维度的主成分用于计算相似性,因此会丢失掉一些结构信息。虽然这些算法都提高了去噪效果,但是普遍存在不能充分保持细节纹理信息的缺点。目前公认的最优图像去噪算法为三维块匹配(BM3D)算法[8],其将具有相似结构的二维图像块组成三维数组,然后用协同滤波的方法进行去噪,最后逆操作得到去噪图像。BM3D算法不仅能得到很高的峰值信噪比(PSNR),还具有很好的视觉效果,但是由于BM3D算法缺乏对图像整体信息的分析,因此,当噪声水平增高时,去噪性能快速下降,且该算法难度相对较大。

本文针对NLM纹理细节保持不充分的问题,提出了一种结合模糊边缘补足(FEC)修正图像块权重的NLM图像去噪方法。利用模糊边缘检测的方法得到边缘图像,结合边缘图像修正相似块权重和相似性权重参数,不仅能准确地检测相似图像块,而且还更加突出了边缘纹理信息。实验证明:PSNR值得到显著提高的同时保留了更多的边缘纹理信息。

2 基本原理

2.1 NLM图像去噪算法

NLM图像去噪算法[4]与均值去噪方法原理类似,后者将所有相似像素点相加求平均值以达到平滑去噪的效果。前者寻找相似像素点,并赋予相似像素点不同的权值以达到去噪的效果。设v为噪声图像,v(i)为噪声图像中的像素点,对像素点v(i)进行非均值平均估计:

NL[v(i)]=jΩ(i)w(i,j)v(j),(1)

式中Ω(i)为以像素点i为中心的搜索区域;w(i,j)为赋予像素点v(j)的权重,其值取决于像素点ij的相似性,其范围为[0,1],并且∑jw(i,j)=1。

分别以像素点ij为中心,取大小为w×w的图像块,分别记为NiNj。像素点ij的相似性依赖于图像块NiNj的高斯加权欧氏距离 GaNi-Nj22。则w(i,j)的计算式表示为

w(i,j)=exp-GaNi-Nj22/h2ZZ=jΩ(i)exp-GaNi-Nj22h2,(2)

式中h为相似权重参数,决定图像平滑的程度;Ga为高斯核,其作用是对图像进行预平滑处理,去除噪声对图像块相似性度量的影响,提高图像块相似性度量的准确性。

2.2 FEC算法

文献[ 9]提出FEC算法,对每个像素点进行处理:

θij=μijU,(3)

式中U为图像中最大的像素值;μij为图像中的像素点,ij分别为对应像素的横纵坐标。

以像素点(i,j)为中心点,取大小为w×w的窗口,记取出的窗口图像为P。基于FEC原理定义像素的锐化式为

θ̅ij=min1,τωxymin(θxy,1-θxy)×max(θxy)-min(θxy)U,x,y[1,w],(4)

式中τω均为缩放因子,(x,y)为窗口图像P内像素的坐标。

3 基于FEC的自适应NLM图像去噪算法

3.1 基于FEC改进的相似性权重

传统的NLM图像去噪算法中相似性权重依赖于相似图像块的高斯加权欧氏距离。(2)式中Ga的作用是减小噪声对相似图像块度量的影响,但Ga在减小噪声影响的同时也会平滑掉纹理信息。针对此问题,文献[ 10]提出了一种融合结构张量的非局部全变分的图像去噪(STNL)方法,用结构张量相似性权重优化非局部信息,能够很好地保持边缘轮廓信息,但是对于细节纹理信息保持的不够理想。如图1所示,对比图1 (a2)、(b2)和(a3)、(b3)可以看出,结构张量容易受到噪声的影响,在含有大量噪声的情况下不能很好地提取出细节纹理信息。基于结构张量的图像去噪方法在纹理信息的保持上,不具有很好的效果。

本文采用FEC方法修正相似图像块的权重,提高相似像素的权重,削弱不相似像素的权重,提高度量相似图像块的准确性,从而保持更多的纹理信息。图1 (a4)、(b4)为FEC方法处理得到的特征图像。可以看出,FEC方法能够很好地检测出边缘和纹理信息,同时也有一定的抗噪能力。

将2.2节中FEC得到的特征图像记作Ie,由 θ̅ij组成。在图像Ie上取以像素点i为中心的大小为w×w的图像块Nei,和以像素点j为中心的大小为w×w的图像块Nej。噪声图像v(i)和图像Ie上的像素点ij的位置保持一致。求取图像块Nei和图像块Nej的高斯加权欧氏距离 GaNei-Nej22,用此修正(2)式得到的计算表达为

wei,j=exp-GaNi-Nj22/h2·exp-GaNei-Nej22/h2ZeZe=jΩ(i)exp-GaNi-Nj22h2·exp-GaNei-Nej22h2(5)

图 1. 边缘纹理图。(a1)(b1)原始图像;(a2)(b2)原始图像结构张量的迹图;(a3)(b3)噪声图像结构张量的迹图;(a4)(b4)模糊边缘补足图

Fig. 1. Edge texture images. (a1)(b1) Original images; (a2)(b2) traces of structure tensor of original images; (a3)(b3) traces of structure tensor of noise images; (a4)(b4) images of fuzzy edge complement

下载图片 查看所有图片

3.2 基于FEC改进的相似权重参数

(2)式中相似权重参数h为常量,但是其对每个相似像素点所赋予的权重产生很大影响,起到平滑图像的作用。当赋予h较大值时,得到的w(i,j)值比较接近,因此对图像平滑作用较大。当赋予h较小值时,得到的w(i,j)值差距较大,因此对图像平滑作用较小。对于不同的噪声水平,应该采取不同程度的平滑。噪声水平大,则应赋予h较大的值,起到平滑噪声的作用;噪声水平小,则应赋予h较小的值[11]。此外,对于图像中平坦区域、纹理区域、边缘区域也应该赋予h不同大小的值。在平坦区域应该赋予h较大的值,而在边缘区域应该防止因过度平滑而使边缘模糊,所以应赋予h较小的值。

这里定义相似性权重参数h(i)为

h(i)=ασ,(6)

式中α值由图像内容决定,σ值由含噪图像噪声水平估计决定。

2.2节中求得的特征图像Ie突出了边缘纹理细节,因此基于特征图像Ie可以区分平坦区域、边缘区域、纹理区域。将特征图像Ie的像素值映射到[0,1]区间表示为

Iθij=θ̅ij-min(θ̅ij)i,jΩemax(θ̅ij)i,jΩe-min(θ̅ij)i,jΩe,(7

式中Ωe为FEC特征图像。取图像块中像素的均值作为区分图像内容的指标,表示为

eij=1SijsIθij,(8)

式中S为以像素Iθij为中心的搜索区域。

根据实验经验值,α值定义为[11]

α(Iθij)=α=2.8,0<eij13α=2.3,1/3<eij23α=1.8,2/3<eij1,(9)

式中α的值根据图像内容不同分别赋予不同的值,在内容丰富的区域赋予α较小的权值,其中平坦区域赋予α较大的权值,而纹理区域则赋予α中间值。

综上所述,采用(6)式对(5)式中的相似性权重系数h进行修正,得到(10)式,可以根据图像内容自适应匹配合适的相似性权重参数,从而实现对不同图像内容的合理平滑。采用文献[ 12]方法对含噪图像进行噪声水平估计:

we(i,j)=exp-GaNi-Nj22/hi2·exp-GaNei-Nej22/hi2ZeZe=jΩ(i)exp-GaNi-Nj22hi2·exp-GaNei-Nej22hi2,(10)

则基于FEC的自适应非局部均值图像去噪(FEC-NLM)算法定义为

NL[v(i)]=jΩ(i)we(i,j)v(j)(11)

3.3 算法整体流程

结合上述内容,其整体算法流程如下:

1) 初始化。读入噪声图像,初始化参数。

2) 评估噪声水平。

3) 求取边缘纹理特征图像。用FEC算法处理噪声图像得到边缘纹理特征图像。

4) 计算we(i,j)值。取一个像素点并以其为中心取图像块,在搜索区域内找出其相似块并取出边缘纹理特征图像中相同位置的图像块,用(6)式计算出相应的h值,然后用(10)式计算出we(i,j)值。

5) 计算去噪的像素值。用(11)式计算出去噪的像素值。

6) 循环步骤4)、5),直至计算出所有像素点。

4 实验结果

选取5张图片作为测试图片,实验中对图像分别加入方差为10和20的高斯白噪声,图像块大小为8 pixel×8 pixel,图像搜索区域大小设为15 pixel×15 pixel,缩放因子τ=9、ω=3。实验主要对FEC-NLM算法、NLM算法、STNL算法、全变分非局部均值去噪(NLTV)算法进行比较。为了评价去噪效果选用PSNR对其进行分析对比。实验操作系统为64位的Windows7,主要硬件为Intel Core i5处理器、内存为4G,软件环境为Matlab R2103a。表1给出了NLM、NLTV、STNL算法和本文算法对测试图像的PSNR。

表1数据可以看出,FEC-NLM算法的PSNR值均高于NLM算法、NLTV算法和STNL算法。当σ=10时,FEC-NLM算法的PSNR值比NLTV算法的PSNR值高出约为2.7 dB,比NLM算法高出约为3.5 dB,比STNL算法高出约为1.5 dB。当σ=20时,FEC-NLM算法的PSNR值比NLTV算法的PSNR值高出约为1.8 dB,比NLM算法高出约为2.5 dB,比STNL算法高出约为0.9 dB。纵向比对可以看出,本文算法对于Boat以及Lena这些边缘纹理信息丰富的图像改进效果比较大,而对于Peppers和House含有较少纹理信息的图像改进效果相对较小。

图2(a)、(b)为原始图像Boat和Lena,图2(c)、(d)为对相应图像加入σ=20的高斯白噪声的噪声图像,图2(e)、(f)为NLM算法去噪图像,图2(g)、(h)为NLTV算法去噪图像,图2(i)、(j)为STNL算法去噪图像,图2(k)、(l)为FEC-NLM算法去噪图像。对比图2(e)、(g)、(i)和(k)可以看出,FEC-NLM算法对于船上的桅杆和天线保持得更好,而图2(e)很明显丢失了大量细小纹理特征,FEC-NLM算法相比前3种算法,对海面的处理没有过度平滑。对比图2(f)、(h)、(j)和(l)可以看出,采用NLM算法获得的图2(f)已基本区分不出帽尖和帽子上的纹理,采用NLTV算法获得的图2(h)帽尖和纹理比较模糊,采用STNL算法获得的图2(j)纹理比较模糊,而采用FEC-NLM算法获得的图2(l)能更好地恢复出了帽尖和帽子上的细小纹理。通过对比不同算法去噪后的PSNR值可知,FEC-NLM去噪算法能很好地去除噪声,同时由图2可知, FEC-NLM算法也能很好地保持图像的边缘和细节纹理特征。

表 1. 图像去噪后的PSNR

Table 1. PSNR in denoised images dB

σOriginalimageNoiseimageNLMalgorithmNLTValgorithmSTNLalgorithmFEC-NLMalgorithm
10Lena28.1330.8731.5632.3134.83
Peppers28.1429.4630.3831.0331.79
House28.1132.1432.9934.1135.61
Boat28.1428.0729.4331.2732.12
Barbara28.1429.6030.5431.9533.61
20Lena22.1130.6631.3732.0432.97
Peppers22.1328.7829.3730.3031.03
House22.1331.6532.2332.9533.82
Boat22.1127.7828.5129.6031.01
Barbara22.1229.1929.8830.8731.41

查看所有表

图 2. 去噪效果图。(a)(b)原始图像;(c)(d)噪声图像;(e)(f)传统NLM去噪;(g)(h) NLTV去噪;(i)(j) STNL去噪;(k)(l) FEC-NLM去噪

Fig. 2. Denoised images. (a)(b) Original images; (c)(d) noise images; (e)(f) traditional NLM denoising; (g)(h) NLTV denoising; (i)(j) STNL denoising; (k)(l) FEC-NLM denoising

下载图片 查看所有图片

5 分析与讨论

由于受噪声的影响,其细节纹理很难和背景区域区分出来,在NLM去噪的过程中,纹理区域和背景区域的高斯加权距离近似相等,所以背景区域的像素点也会被错认为是纹理区域像素点的相似像素点。图3(a)为没有噪声污染图像选取的图像相似块,红色方框为一组相似块,而蓝色方框为一组相似块。对于噪声图像图3(b),由于红色方框边缘明显,所以选取的相似块基本相同,但是蓝色方框标明相似块的纹理和背景区域区分度不大,导致其将背景区域也选为相似块,因此细节纹理被平滑掉。

图 3. 相似块。(a)原图相似块;(b)噪声图像相似块;(c)边缘纹理特征图像相似块

Fig. 3. Similar blocks. (a) Similar blocks on original image; (b) similar blocks on noise image; (c) similar blocks on edge texture futures

下载图片 查看所有图片

由此可以看出,准确选取相似性图像块是影响去噪效果的重要部分。根据(5)式得到边缘检测特征图像,如图3 (c)所示。图3(c)中边缘纹理区域和平坦区域有明显的差别,计算相应图像块中的高斯加权欧氏距离可知,纹理区域和纹理区域的高斯加权欧氏距离比较小,纹理区域与平坦区域的高斯加权欧氏距离比较大,而平坦区域之间的高斯加权欧氏距离比较小。因此,对应图像块的高斯加权欧氏距离修正NLM算法的相似度,可将图3(b)中错误分类的相似图像块的权值缩小,从而减小背景信息对纹理信息的影响。但是在高噪声的情况下,图3(b)中红色区域正确取出的相似图像块中其像素也被严重的污染,因此会减弱去噪效果。

6 结论

提出了一种融合边缘检测和NLM的图像去噪算法。其在NLM的思想的基础上,利用边缘检测得到的信息,重新定义了NLM的相似性权重和相似性权重参数,充分利用了图像块间像素的相似性,同时考虑了低对比度纹理结构保持问题。对于图像内容不同区域赋予相似性权重参数,保证边缘纹理不会过度平滑。本文算法从相似图像块度量和相似性权重参数两个角度对NLM去噪算法的纹理细节保持不足问题进行了改进。实验表明:本文算法得到的去噪图像不仅具有很高的PSNR值,同时还具有很好的视觉效果。

参考文献

[1] 杨永发, 李琦. 双边滤波算法的太赫兹共焦扫描图像去噪应用[J]. 激光与光电子学进展, 2015, 52(12): 121101.

    Yang Y F, Li Q. Application of bilateral filtering algorithm on terahertz confocal scanning image denoising[J]. Laser & Optoelectronics Progress, 2015, 52(12): 121101.

[2] 李东明, 盖梦野, 李超然, 等. 基于小波域的Contourlet 变换法的自适应光学图像去噪算法研究[J]. 激光与光电子学进展, 2015, 52(11): 111001.

    Li D M, Gai M Y, Li C R, et al. Research on adaptive optics image denoising algorithm based on the wavelet-based Contourlet transform[J]. Laser & Optoelectronics Progress, 2015, 52(11): 111001.

[3] 余临倩, 覃亚丽, 张晓帅. 基于字典学习的梯度重权非局部平均的强噪声图像去噪[J]. 激光与光电子学进展, 2016, 53(11): 111002.

    Yu L Q, Qin Y L, Zhang X S. Denoising of strong noisy image via gradient reweighted non-local averaging over learned dictionaries[J]. Laser & Optoelectronics Progress, 2016, 53(11): 111002.

[4] BuadesA, CollB, Morel JM. A non-local algorithm for image denoising[C]// 2005 IEEE Computer Society Conference on Computer Vision and Pattern Recognition, 2005, 2: 60- 65.

[5] Kervrann C, Boulanger J. Optimal spatial adaptation for patch-based image denoising[J]. IEEE Transactions on Image Processing, 2006, 15(10): 2866-2878.

[6] KervrannC, BoulangerJ, CoupéP. Bayesian non-local means filter, image redundancy and adaptive dictionaries for noise removal[C]// International Conference on Scale Space and Variational Mehtods in Computer Vision, 2007, 4485: 520- 532.

[7] Grewenig S, Zimmer S, Weickert J. Rotationally invariant similarity measures for nonlocal image denoising[J]. Journal of Visual Communication & Image Representation, 2011, 22(2): 117-130.

[8] Dabov K, Foi A, Katkovnik V, et al. Image denoising by sparse 3-D transform-domain collaborative filtering[J]. IEEE Transactions on Image Processing, 2007, 16(8): 2080-2095.

[9] Chang C C, Lee J S. Le T H N. Hybrid wet paper coding mechanism for steganography employing n-indicator and fuzzy edge detector[J]. Digital Signal Processing, 2010, 20(4): 1286-1307.

[10] Wang S Y. Image denoising model based on structure tensor and non-local total variation[J/OL].Application Research of Computers, 2018( 2017- 03- 15). . http://kns.cnki.net/kcms/detail/51.1196.TP.20170315.1537.124.html

    王诗言. 一种融合结构张量与非局域全变分的图像去噪方法[J/OL]. 计算机应用研究, 2018( 2017- 03- 15). . http://kns.cnki.net/kcms/detail/51.1196.TP.20170315.1537.124.html

[11] 许光宇. 非局部图像去噪方法及其应用研究[D]. 合肥: 合肥工业大学, 2013.

    Xu GY. Research on nonlocal image denoising methods and their applications[D]. Hefei: Hefei University of Technology, 2013.

[12] Amer A, Dubois E. Fast and reliable structure-oriented video noise estimation[J]. IEEE Transactions on Circuits & Systems for Video Technology, 2005, 15(1): 113-118.

曹硕, 黄利萍, 侯倍倍, 陈刚. 基于模糊边缘补足的自适应非局部均值图像去噪算法[J]. 激光与光电子学进展, 2018, 55(1): 011003. Cao Shuo, Huang Liping, Hou Beibei, Chen Gang. Self-Adaptive Non-Local Means Image Denoising Algorithm Based on Fuzzy Edge Complement[J]. Laser & Optoelectronics Progress, 2018, 55(1): 011003.

本文已被 1 篇论文引用
被引统计数据来源于中国光学期刊网
引用该论文: TXT   |   EndNote

相关论文

加载中...

关于本站 Cookie 的使用提示

中国光学期刊网使用基于 cookie 的技术来更好地为您提供各项服务,点击此处了解我们的隐私策略。 如您需继续使用本网站,请您授权我们使用本地 cookie 来保存部分信息。
全站搜索
您最值得信赖的光电行业旗舰网络服务平台!