基于模糊边缘补足的自适应非局部均值图像去噪算法 下载: 877次
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)进行非均值平均估计:
式中Ω(i)为以像素点i为中心的搜索区域;w(i,j)为赋予像素点v(j)的权重,其值取决于像素点i与j的相似性,其范围为[0,1],并且∑jw(i,j)=1。
分别以像素点i、j为中心,取大小为w×w的图像块,分别记为Ni和Nj。像素点i、j的相似性依赖于图像块Ni和Nj的高斯加权欧氏距离
式中h为相似权重参数,决定图像平滑的程度;Ga为高斯核,其作用是对图像进行预平滑处理,去除噪声对图像块相似性度量的影响,提高图像块相似性度量的准确性。
2.2 FEC算法
文献[ 9]提出FEC算法,对每个像素点进行处理:
式中U为图像中最大的像素值;μij为图像中的像素点,i、j分别为对应像素的横纵坐标。
以像素点(i,j)为中心点,取大小为w×w的窗口,记取出的窗口图像为P。基于FEC原理定义像素的锐化式为
式中τ和ω均为缩放因子,(x,y)为窗口图像P内像素的坐标。
3 基于FEC的自适应NLM图像去噪算法
3.1 基于FEC改进的相似性权重
传统的NLM图像去噪算法中相似性权重依赖于相似图像块的高斯加权欧氏距离。(2)式中Ga的作用是减小噪声对相似图像块度量的影响,但Ga在减小噪声影响的同时也会平滑掉纹理信息。针对此问题,文献[
10]提出了一种融合结构张量的非局部全变分的图像去噪(STNL)方法,用结构张量相似性权重优化非局部信息,能够很好地保持边缘轮廓信息,但是对于细节纹理信息保持的不够理想。如
本文采用FEC方法修正相似图像块的权重,提高相似像素的权重,削弱不相似像素的权重,提高度量相似图像块的准确性,从而保持更多的纹理信息。
将2.2节中FEC得到的特征图像记作Ie,由
图 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)为
式中α值由图像内容决定,σ值由含噪图像噪声水平估计决定。
2.2节中求得的特征图像Ie突出了边缘纹理细节,因此基于特征图像Ie可以区分平坦区域、边缘区域、纹理区域。将特征图像Ie的像素值映射到[0,1]区间表示为
式中Ωe为FEC特征图像。取图像块中像素的均值作为区分图像内容的指标,表示为
式中S为以像素Iθij为中心的搜索区域。
根据实验经验值,α值定义为[11]
式中α的值根据图像内容不同分别赋予不同的值,在内容丰富的区域赋予α较小的权值,其中平坦区域赋予α较大的权值,而纹理区域则赋予α中间值。
综上所述,采用(6)式对(5)式中的相似性权重系数h进行修正,得到(10)式,可以根据图像内容自适应匹配合适的相似性权重参数,从而实现对不同图像内容的合理平滑。采用文献[ 12]方法对含噪图像进行噪声水平估计:
则基于FEC的自适应非局部均值图像去噪(FEC-NLM)算法定义为
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. 图像去噪后的PSNR
Table 1. PSNR in denoised images dB
|
图 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)原图相似块;(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)式得到边缘检测特征图像,如
6 结论
提出了一种融合边缘检测和NLM的图像去噪算法。其在NLM的思想的基础上,利用边缘检测得到的信息,重新定义了NLM的相似性权重和相似性权重参数,充分利用了图像块间像素的相似性,同时考虑了低对比度纹理结构保持问题。对于图像内容不同区域赋予相似性权重参数,保证边缘纹理不会过度平滑。本文算法从相似图像块度量和相似性权重参数两个角度对NLM去噪算法的纹理细节保持不足问题进行了改进。实验表明:本文算法得到的去噪图像不仅具有很高的PSNR值,同时还具有很好的视觉效果。
[1] 杨永发, 李琦. 双边滤波算法的太赫兹共焦扫描图像去噪应用[J]. 激光与光电子学进展, 2015, 52(12): 121101.
[2] 李东明, 盖梦野, 李超然, 等. 基于小波域的Contourlet 变换法的自适应光学图像去噪算法研究[J]. 激光与光电子学进展, 2015, 52(11): 111001.
[3] 余临倩, 覃亚丽, 张晓帅. 基于字典学习的梯度重权非局部平均的强噪声图像去噪[J]. 激光与光电子学进展, 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.
[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.
[11] 许光宇. 非局部图像去噪方法及其应用研究[D]. 合肥: 合肥工业大学, 2013.
Xu GY. Research on nonlocal image denoising methods and their applications[D]. Hefei: Hefei University of Technology, 2013.
Article Outline
曹硕, 黄利萍, 侯倍倍, 陈刚. 基于模糊边缘补足的自适应非局部均值图像去噪算法[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.