CN101950414B - 自然图像非局部均值去噪方法 - Google Patents

Abstract

本发明公开了一种自然图像的非局部均值去噪方法，它属于图像处理技术领域，主要克服现有自然图像非局部均值去噪中相似性计算不准确的问题。其实现过程是：(1)对输入的含噪自然图像做小波变换；(2)对高频信息进行方差归一化处；(3)计算像素点x与搜寻区域内像素点y的相似性，得到搜寻区域内所有像素点的权值；(4)根据计算出的搜寻区域内所有像素点的权值，对搜寻区域内所有像素点进行加权平均，得到像素点修正后的灰度值；(5)用修正后像素点的灰度值取代输入的含噪自然图像中像素点的灰度值，得到去噪后的图像。本发明在总体性能上优于其它的去噪方法，能够更好地平滑噪声的同时保持自然图像的边缘和纹理等细节，可用于对自然图像的去噪处理。

Description

自然图像非局部均值去噪方法

技术领域

本发明属于图像处理技术领域，具体地说是一种去噪方法，可用于对自然图像的去噪处理。

背景技术

随着计算机和数码成像设备的日益普及，数字图像处理越来越受到人们的重视。然而由于成像设备及成像条件的限制，数字图像在采集，转换，以及运输过程中不可避免受到噪声的污染，因此图像去噪作为图像处理领域的基本技术之一，受到广泛的重视。许多实际的噪声可以近似的认为是高斯白噪声，去除图像中的高斯白噪声成为图像去噪领域中一个重要的方向。

传统的去噪方法大致可以分为两类，一类是基于空域的方法，一类是基于变换域的方法。空域去噪方法中比较经典的方法包括高斯滤波，中值滤波，双边滤波等。它们的共同特点就是利用局部窗口内像素灰度值的连续性来对当前像素进行灰度调整。这些方法大都在去除噪声的同时模糊了图像的细节信息，例如图像的边缘，纹理等。

由于自然图像，特别是纹理图像所含信息具有一定的冗余性，Buades等人提出了一种非局部均值的去噪方法。该方法以当前像素为中心取大小一定的窗口，在整幅图像内寻找与其具有相似结构的窗口，以窗口之间的相似度为权值对当前像素的灰度值进行调整。由于这种方法在去噪领域良好的性能，自提出以来迅速引起众多学者的广泛关注，但是它仍然存在以下问题：1：算法复杂性比较大；2：权值计算准确性欠佳；3：图像的边缘与细节仍存在一定程度的模糊。

基于变换域的去噪方法比较成熟的就是小波域的各种去噪方法，小波去噪的关键问题是对图像小波系数统计模型的研究，这类算法的基本思想是把统计模型作为小波系数的先验概率模型，然后利用这个先验信息，在贝叶斯框架下对原始图像的小波系数进行估计。然而，小波变换仍然具有以下不足：1：高维时系数的稀疏性较差；2：图像分解后方向信息有限，为了克服小波变换的不足，近几年来出现了Ridgelet，Curvelet，Contourlet，Brushlet，Bandelet等一系列图像分解的新工具。但是由于变换域的方法只是对高频图像的小波系数做萎缩阈值调整，对低频图像不做处理，所以其最终的去噪效果并不是很令人满意，而且常常会产生吉布斯现象。

发明内容

本发明的目的在于克服上述已有技术的不足，提出了自然图像非局部均值去噪方法，对图像块相似性度量的准确性进行改进，实现对自然图像去噪中边缘和平滑区域的兼顾，提高图像去噪效果。

为实现上述目的，本发明包括如下步骤：

(1)对输入的含噪自然图像做小波变换，将其分解为低频图像和高频图像；

(2)对高频图像进行方差归一化处理；

(3)利用如下公式对输入的含噪自然图像中待修正像素点x与搜寻区域内像素点y的相似性进行计算，得到搜寻区域内所有像素点的权值：

w(x，y)＝w₁(x，y)×w₂(x，y)

其中，w₁(x，y)表示低频图像中的相似性，

L₁(x)，L₂(y)分别表示低频图像中以x，y为中心的大小为n×n的图像块，h₁为平滑控制参数，h₁＝k₁nσ₀2^t-1，n＝7，σ₀是图像所含噪声的标准差，t指图像的分解层数，k₁是常数，实验中取k₁＝1.0；

w₂(x，y)表示高频图像中的相似性，

H₁(x)，H₂(y)分别表示所有高频图像中以x，y为中心的大小为n×n的图像块，h₂为平滑控制参数，h₂＝k₂nσ2^t-1，n＝7，σ是高频子带的噪声标准差，取值为1.0，k₂是常数，实验中取k₂＝0.6；

权值w(x，y)满足：0≤w(x，y)≤1，

(4)根据计算出的搜寻区域内所有像素点的权值，对搜寻区域内所有像素点进行加权平均，得到像素点修正后的灰度值：

$V^{\′}\left(x\right)=\underset{y\∈J}{\mathrm{\Σ}}V\left(y\right)w(x,y)$

其中，J为搜索窗口像素点集合，V(y)为搜寻区域内像素点的灰度值；

(5)用修正后像素点的灰度值V′(x)取代输入的含噪自然图像中像素点的灰度值，得到去噪后的图像。

本发明由于使用了新提出的权值公式，与现有的技术相比具有以下优点：

1.本发明能够更准确地计算含噪自然图像中像素点之间的相似性，进而能更准确的计算出待修正像素点的灰度值；

2.本发明能够更准确计算出待修正像素点的灰度值，进而能够在较好平滑噪声的同时保持和恢复自然图像的边缘和纹理细节。

附图说明

图1是本发明的流程图；

图2是本发明使用的测试图像；

图3是本发明使用的含噪图像；

图4是用现有的BayesShrink小波阈值方法对图3进行去噪的结果图；

图5是用现有的双边滤波器方法对图3进行去噪的结果图；

图6是用现有的多分辨双边滤波器方法对图3进行去噪的结果图；

图7是用现有的非局部滤波器方法对图3进行去噪的结果图；

图8是用本发明方法对图3进行去噪的结果图。

具体实施方式

参照附图1，本发明给出自然图像非局部均值去噪方法，包括如下步骤：

步骤1，对输入的含噪自然图像做小波变换，将其分解为低频图像和高频图像。

数字图像由于成像设备以及成像条件的限制，不可避免的要受到噪声的污染，许多实际的噪声都可以认为是高斯加性白噪声，含噪的自然图像模型为：

v＝u+n

其中，v为含噪图像的灰度值，u是干净图像的灰度值，n是高斯加性白噪声，它主要利用自然图像中的冗余信息来达到去噪的目的，由于平稳小波分解出来的图像与原图是等大的，保持了图像信息的冗余性，本发明在实验中选择平稳小波对图像进行分解，使用的小波基是db1，但不局限于这种分解方法，例如也可以使用正交小波变换对图像进行分解，可将含噪自然图像分解为低频图像和高频图像两部分。

步骤2，对高频图像进行方差归一化处理。

2.1)利用Donoho提出的鲁棒中值法，估计各高频图像的噪声标准差：

${\mathrm{\σ}}_{\mathrm{ij}}=\frac{\mathrm{median}\left(\mathrm{abs}\left(V_{\mathrm{ij}}(:)\right)\right)}{0.6745}$

其中，V_ij是指第i层第j个方向的高频图像小波系数，V_ij(:)是指取出V_ij中含有的全部小波系数，abs是指取绝对值，median是指取中值；

2.2)利用2.1)求出的各高频图像的噪声标准差，对各高频图像的小波系数进行方差归一化处理：

$V_{\mathrm{ij}}^{\′}=\frac{V_{\mathrm{ij}}}{{\mathrm{\σ}}_{\mathrm{ij}}}.$

步骤3，利用如下公式计算像素点x与搜寻区域内像素点y的相似性，得到搜寻区域内所有像素点的权值：

w(x，y)＝w₁(x，y)×w₂(x，y)

其中，

表示低频图像中的相似性，L₁(x)，L₂(y)分别表示低频图像中以x，y为中心的大小为n×n的图像块，h₁为平滑控制参数，h₁＝k₁nσ₀2^t-1，n＝7，σ₀是图像所含噪声的标准差，t指图像的分解层数，k₁是常数，实验中取k₁＝1.0；

表示高频图像中的相似性，H₁(x)，H₂(y)分别表示所有高频图像中以x，y为中心的大小为n×n的图像块，h₂为平滑控制参数，h₂＝k₂nσ2^t-1，n＝7，σ是高频子带的噪声标准差，取值为1.0，t指图像的分解层数，k₂是常数，实验中取k₂＝0.6；

权值w(x，y)满足：0≤w(x，y)≤1，

步骤4，根据计算出的搜寻区域内所有像素点的权值，对搜寻区域内所有像素点进行加权平均，得到像素点修正后的灰度值：

$V^{\′}\left(x\right)=\underset{y\∈J}{\mathrm{\Σ}}V\left(y\right)w(x,y)$

其中，J为搜索窗口像素点集合，V(y)为搜寻区域内像素点的灰度值。

步骤5，用修正后像素点的灰度值V′(x)取代输入的含噪自然图像中像素点的灰度值，得到去噪后的图像。

本发明效果可以通过以下实验进一步证实：

一.实验条件和内容

实验条件：实验所使用的输入图像如图2和3所示，其中，图2(a)是测试图像lena，图2(b)是测试图像barbara，图2(c)是测试图像peppers，图2(d)是测试图像house，图3是对图2(a)加入噪声标准差为20的含噪的lena图像。

实验内容：在上述实验条件下，分别使用现有的BayesShrink小波阈值滤波方法，双边滤波器方法，多分辨双边滤波器方法以及非局部滤波器方法和本发明方法对图3进行实验。

二.实验结果

用BayesShrink小波阈值方法的去噪效果如图4所示，其中对图3做五层分解，噪声标准差采用Donoho的鲁棒中值法进行估计，从图4中可以看出，此方法的噪声抑制能力有限，而且边缘与细节存在模糊。

用双边滤波器方法的去噪结果如图5所示，其中σ_d＝1.8，σ_r＝2×σ，搜寻窗的大小为11×11，从图5可以看出，此方法的噪声抑制能力要优于BayesShrink小波阈值方法，但是它同样存在边缘和细节模糊的问题；

用多分辨双边滤波器方法的去噪结果如图6所示，其中对图3做两层小波分解，用BayesShrink方法对高频图像进行处理，用双边滤波器对低频图像进行处理，其中σ_d＝1.8，σ_r＝1.0×σ，搜寻窗的大小为11×11，从图6可以看出，此方法噪声抑制能力要优于上面两种方法，但是存在边缘和细节保持度欠佳的问题；

用非局部滤波器方法的去噪结果如图7所示，其中搜寻窗大小为21×21，相似窗大小为7×7，平滑参数h＝10σ，σ为图像所含噪声标准差，从图7中可以看出：此方法噪声抑制能力比较好，但是不能很好的保持图像的边缘和纹理信息。

用本发明方法的去噪结果如图8所示，从图8中可以看出：它的去噪效果要优于上面所提到的所有方法，同质区域也较平滑，图像的亮度保持效果较好，而且图像的边缘，细节也得到了很好的保持；

对图2中的测试图像分别加入噪声标准差为10，20，30的高斯加性白噪声，用PSNR作为去噪效果的评价指标，将上述四种现有的去噪方法和本发明的方法进行比较，各种方法的去噪效果PSNR值列在表1中。

表1各种去噪结果对比

表1结果均为10次平均后的结果，从表1中可以看出，本发明方法的去噪效果比BayesShrink小波阈值方法，双边滤波器方法，以及多分辨双边滤波器方法在PSNR值上都有很大的提高，与非局部均值滤波算法相比，本发明方法除了在Barbara图像的去噪效果上优势不是很明显，在其它几幅图像上都较非局部均值算法有很大的提高，从表1中我们还可以看出随着图像噪声标准差的增大，本发明方法的优势也越来越明显。

以上实验结果表明，本发明在总体性能上优于现有的其它去噪方法，能够更好地平滑噪声的同时保持自然图像的边缘和纹理等细节。

Claims (2)

1.一种自然图像的非局部均值去噪方法，包括如下步骤：

(1)对输入的含噪自然图像做小波变换，将其分解为低频图像和高频图像；

(2)对高频图像进行方差归一化处理；

(3)利用如下公式对输入的含噪自然图像中待修正像素点x与搜寻区域内像素点y的相似性进行计算，得到搜寻区域内所有像素点的权值：

w(x，y)＝w₁(x，y)×w₂(x，y)

其中，w₁(x，y)表示低频图像中的相似性，

L₁(x)，L₂(y)分别表示低频图像中以x，y为中心的大小为n×n的图像块，h₁为平滑控制参数，h₁＝k₁nσ₀2^t-1，1≤n≤5，σ₀是图像所含噪声的标准差，t指图像的分解层数，k₁是常数，取k₁＝1.0；

w₂(x，y)表示高频图像中的相似性，H₁(x)，

H₂(y)分别表示所有高频图像中以x，y为中心的大小为n×n的图像块，h₂为平滑控制参数，h₂＝k₂nσ2^t-1，σ是高频子带的噪声标准差，取值为1.0，k₂是常数，取k₂＝0.6；权值w(x，y)满足：0≤w(x，y)≤1，

(4)根据计算出的搜寻区域内所有像素点的权值，对搜寻区域内所有像素点进行加权平均，得到像素点修正后的灰度值：

$V^{\′}\left(x\right)=\underset{y\∈J}{\mathrm{\Σ}}V\left(y\right)w(x,y)$

其中，J为搜寻区域像素点集合，V(y)为搜寻区域内像素点的灰度值；

(5)用修正后像素点的灰度值V′(x)取代输入的含噪自然图像中像素点的灰度值，得到去噪后的图像。

2.根据权利要求1所述的自然图像非局部均值去噪方法，其特征在于步骤(2)中所述的对高频图像进行方差归一化处理，按如下步骤进行：

2a)对于输入的含噪自然图像进行小波变换，利用Donoho提出的鲁棒中值法，估计各高频图像的噪声标准差：

${\mathrm{\σ}}_{\mathrm{ij}}=\frac{\mathrm{median}\left(\mathrm{abs}\left(V_{\mathrm{ij}}(:)\right)\right)}{0.6745}$

其中，V_ij(：)是指取出V_ij中含有的全部小波系数，V_ij是指第i层第j个方向的高频图像小波系数，abs是指取绝对值，median是指取中值；

2b)利用如下公式对各高频图像的小波系数进行方差归一化：

$V_{\mathrm{ij}}^{\′}=\frac{V_{\mathrm{ij}}}{{\mathrm{\σ}}_{\mathrm{ij}}}.$

Images