基于L0梯度范数泛函优化和变换域阈值的图像去噪算法

浙江工贸职业技术学院学报JOURNAL OF ZHEJIANG INDUSTRY&TRADE VOCATIONAL COLLEGE第14卷第1期2014年3月V ol.14No.1Mar.2014基于L0梯度范数泛函优化和变换域阈值的图像去噪算法*赵秀芝（浙江工贸职业技术学院电子工程系，浙江温州325003）摘要：基于L0梯度泛函优化和变换域阈值法的图像去噪算法，首先利用L0梯度泛函的最优化算法把含有噪声的图像分解为显著边缘层和细节纹理细节层，然后对含有噪声的细节纹理层图像进行短时傅利叶转换，并在变换域中利用阈值法分离纹理细节层中的图像细节纹理和噪声，进而去除图像的噪声。

实验结果表明图像去噪算法获得了显著的去噪效果。

关键词：去噪；变换阈值；图像分解中图分类号：TP391文献标识码：A文章编号：1672-0105（2014）01-0036-04Image Denoising Algorithm Based on L0Gradient Minimization and Threshold ShrinkageZHAO Xiu-Zhi（Department of Electronic Engineering,Zhejiang Industry &Trade vocational College,Wenzhou,Zhejiang,325003,China ）Abstract:An image denoising algorithm based on L0gradient minimization and threshold shrinkage was proposed in this paper.In this algorithm,the L0gradient minimization algorithm was firstly used to decompose the noisy image into the large-scale edge lay-er and the fine-scale texture layer and then noises in the fine-scale texture layer was denoised by a threshold shrinkage algorithm on the short-time fourier transform space.The experiment results show that the proposed image denoising algorithm achieves significant improvement in its performance.Key words:denoising;threshold shrinkage;image decompositionDoi:10.3969/j.issn.1672-0105.2014.01.008收稿日期：2014-02-16基金项目：2013浙江省高等教育课堂教学改革研究项目（kg2013836）作者简介：赵秀芝(1978－)，女，浙江苍南人，讲师，硕士，主要研究方向：电子信息信号处理、图像彩色视觉模型。

0引言随着数字成像技术的飞速发展，成像设备（如相机）越来越先进，但数字图像却常常在成像的过程中受到设备噪声、周围环境噪声的干扰，以至于图像中含有噪声，为了图像的后续使用，必须对图像进行去噪处理。

目前，图像去噪算法大致分为两类：空间域滤波和变换域滤波。

早期典型的图像去噪算法均为空间域滤波算法，如维纳滤波、均值滤波器等典型的线性滤波器。

此类线性滤波算法简单却会模糊图像边缘和细节纹理。

近年，双边滤波器[1]、L0梯度范数泛函优化[2]等优秀算法被陆续提出。

这些去噪算法可以在去噪的同时保留图像主要边缘信息，但却无法在保留对比度较小的细节纹理信息和保留去除噪声之间做出平衡。

产生此现象的根本原因在于：在空间域中各图像特征并不能很好的表现出来，如小的细节纹理和噪声。

因此，对图像进行空间转换、分析图像特征，再通过阈值的方式去除噪声是另外一个去噪的主要思路。

例如，小波阈值算法[3,4]就是通过选择一个适合的阈值，对小波系数进行阈值处理，实现去噪的目的。

尽管变换域阈值算法在第14卷第1期去噪方面表现出了优越的性能，但却容易在显著边缘处产生振铃现象[5]。

L0梯度范数泛函优化虽然不能有效地分离噪声和图像的细节纹理，但却可将图像的边缘先分离出去，而变换阈值可以在无边缘下有效的去除噪声。

可先利用L0梯度范数泛函优化算法把图像的显著边缘先分离出来，余下的混有噪声的细节纹理图像再利用变换域阈值算法对其进行去噪处理，最后合并边缘和细节纹理重构图像。

此方法既去除了噪声又避免了振铃现象的产生。

1图像的去噪算法假定噪声图像的噪声为高斯白噪，其均值为0，方差为σ2，则噪声图像可表示为：y (i ,j )=x (i ,j )+η(i ,j )（1）式中，x (i ,j )表示没有噪声的原始图像，η(i ,j )表示噪声。

图像去噪的目标就是从噪声图像y (i ,j )中去除噪声η(i ,j )从而恢复原始图像x (i ,j )。

图1给出本文算法的框图，其中L0GF 表示空间域的L0梯度范数泛函优化，用以分解图像的边缘和细节纹理，因而C (i ,j )和c (i ,j )分别表示分解所得的边缘层图像（即高对比度图像）和细节纹理层图像（即低对比度图像），WTS 表示小波变化阈值，c͂(i ,j )表示c (i ,j )经过小波变化阈值算法去噪后的细节纹理层图像，x͂(i ,j )表示重构的、去噪后的图像。

在图1的图像分解过程中，边缘层图像是直接通过L0梯度范数泛函优化算法计算所得，而细节纹理层图像则是由噪声图像与边缘层图像的相减所得，因此噪声也被分解到细节纹理层图像中。

由上述可见，细节纹理层图像上基本都是小对比度的信息，且很多空间特性都具有相似性，而转换阈比较区细节纹理和噪声的特征。

因此，利用短时傅利叶转换（STFT ）[6]对细节纹理层图像进行转换得小波系数，再利用阈值方法区分噪声和细节纹理，最终从细节纹理中去除噪声，处理之后的两层图像重构得无噪声图像，公式如下所示：x ͂(i ,j )=c ͂(i ,j )+C ͂(i ,j )（2）式中，C ͂(i ,j )表示去噪后显著边缘层图像，在本文中等同于图像C (i ,j )。

1.1L0梯度范数泛函优化如上文假设所示y (i ,j )为噪声图像，C (i ,j )表示L0梯度范数泛函优化所得图像，(i ,j )表示图像中像素的位置，等同下文用p 表示。

像素p 与其相邻像素在i 和j 方向的图像梯度可以表示为∇C p =(∂i C p ,∂j C p )，并且在L0梯度范数泛函优化算法中以图像的L0梯度范数为泛函的正则项，如式（3）所示。

式（3）中符号#表示统计图像L0范数非零（即||∂i C p +||∂j C p 不为零）的个数。

最后，通过解式（4）得边缘层图像C (i ,j )（具体解法见参考文件[2]）。

G (C )=#{p |||∂i C p +||∂j C p ≠0}（3）min S ìíîüýþ∑p (C p -y p )2+λ⋅G (C )（4）式（4）中，λ表示正规化系数，平衡最小化变分泛函时正则项和保真函数∑p(C p -y p )2贡献，通常取λ≥0。

1.2小波变换阈值经过L0梯度范数泛函分解后，噪声只存在于细节纹理层中，因而利用变化域阈值进一步分离细节信号中的噪声。

如上所述，此阶段的处理中本文先利用短时傅利叶转换（short-time Fourier trans ⁃form,STFT ）[6]把细节纹理层图像由空间域转换到变换域中，再在软阈值算法[3,4]对STFT 系数进行阈值处理，以去除噪声。

STFT 是和傅立叶相关的一种数学转换关系，它是针对傅立叶变换的缺陷，给信号先加窗函数，再进行傅立叶变。

由于，固定窗函数的STFT 适合平稳信号的分析，而本文的细节纹理层中除了高斯图1本文去噪算法框架图赵秀芝：基于L0梯度范数泛函优化和变换域阈值的图像去噪算法37浙江工贸职业技术学院学报2014年3月白噪外还具有图像的细节纹理信息，其统计特性会随信号所在的位置有所不同。

因此，本根据参考文献[7]思想在STFT中采用一个随信号适应性变化得窗函数，即一个如同双边滤波器、随信号变化的窗函数，具体如下：S f,p =∑q∈N pe-i2πqf2r⋅[w p,q⋅c p,q]（5）式中，下表p表示当前像素位置，Np表示像素p的邻域范围，且范围的半径为r，下表q表示Np 邻域内像素位置。

cp,q表示细节纹理层信号，f表示STFT域中的频率，wp,q表示窗函数，具体如下：wp,q =e-||p-q22σs e-||||c p-c q22σr（6）式中，参数σs 和σr分别决定窗口函数的窗口大小和梯度大小，cp 和cq分别是细节纹理层图像p和q位置的像素值。

在重构图像之前，本文去噪算法的最后一步是利用软阈值算法对含有噪声的STFT系数阈值处理，通过去除噪声的STFT系数，达到分离噪声的目的，具体操作如下：S'f =sgn(Sf)⋅max(Sf-T,0)（7）在式（7）中，Sf 和S'f分别表示阈值前后的STFT系数。

参数T表示阈值，计算如下：T=σ2σ2x（8）σx =max(σ2y-σ2,0)（9）σ2y =1n2∑k,l=1ny2k,l（10）式（8）~（10）中，σ2y 、σ2x和σ2分别表示噪声图像、无噪声的原始图像和噪声的方差。

式（10）的参数n表示把图像分割多个区域，每个区域的大小为n×n。

2实验结果与分析实验选用包含显著边缘、小细节和噪声的一维信号，以及标准的灰度图像“Lena”、“Boats”和“Hill”进行仿真来验证本文的图像去噪算法的有效性。

本文实验中窗函数的参数σr 和σs分别去0.12和1%图像的大小，在计算噪声图像方差σ2y 中参数n选图像的宽度，而L0梯度滤波的参数λ则随噪声的水平有所不同，在噪声方差为10、20、30和40时λ分别为0.003、0.009、0.02和0.03。

图2给出只用L0梯度范数泛函优化算法和本文算法处理一维信号的效果图。

其结果发现：L0梯度范数泛函优化算法可以很好的保留边缘的锐度。

由于参数λ控制保留何种对比度于显著边缘层（即大对比度层）中，λ越小则越多小对比度保留在显著边缘层图像中，但也可能存在保留了一些噪声于显著边缘层中的问题（如图2(b)中圈的区域）。