CN110458783B - 一种基于具有图像细节保护能力的非局部双边滤波的椒盐噪声去除方法 - Google Patents
一种基于具有图像细节保护能力的非局部双边滤波的椒盐噪声去除方法 Download PDFInfo
- Publication number
- CN110458783B CN110458783B CN201910777064.5A CN201910777064A CN110458783B CN 110458783 B CN110458783 B CN 110458783B CN 201910777064 A CN201910777064 A CN 201910777064A CN 110458783 B CN110458783 B CN 110458783B
- Authority
- CN
- China
- Prior art keywords
- image
- noise
- formula
- value
- salt
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000001914 filtration Methods 0.000 title claims abstract description 64
- 239000006002 Pepper Substances 0.000 title claims abstract description 48
- 238000000034 method Methods 0.000 title claims abstract description 44
- 235000002566 Capsicum Nutrition 0.000 title claims abstract description 34
- 150000003839 salts Chemical class 0.000 title claims abstract description 34
- 235000016761 Piper aduncum Nutrition 0.000 title claims abstract description 31
- 235000017804 Piper guineense Nutrition 0.000 title claims abstract description 31
- 235000008184 Piper nigrum Nutrition 0.000 title claims abstract description 31
- 230000002146 bilateral effect Effects 0.000 title claims abstract description 29
- 244000203593 Piper nigrum Species 0.000 title 1
- 241000722363 Piper Species 0.000 claims abstract description 30
- 239000011159 matrix material Substances 0.000 claims abstract description 29
- 238000009499 grossing Methods 0.000 claims description 7
- 150000001875 compounds Chemical class 0.000 claims description 3
- 238000011524 similarity measure Methods 0.000 claims description 3
- 230000003044 adaptive effect Effects 0.000 description 11
- 230000000694 effects Effects 0.000 description 9
- 238000011156 evaluation Methods 0.000 description 7
- 238000013459 approach Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 5
- 238000012545 processing Methods 0.000 description 4
- 241000758706 Piperaceae Species 0.000 description 3
- 101001035951 Homo sapiens Hyaluronan-binding protein 2 Proteins 0.000 description 2
- 102100039238 Hyaluronan-binding protein 2 Human genes 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000035772 mutation Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005728 strengthening Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000003313 weakening effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20024—Filtering details
- G06T2207/20032—Median filtering
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
本发明涉及一种基于具有图像细节保护能力的非局部双边滤波的椒盐噪声去除方法，包括步骤如下：S1、遍历图像，利用mask矩阵记录噪声点位置；S2、利用NASMF算法和mask矩阵对噪声点的定位，对噪声点进行滤波将椒盐噪声转化为高斯噪声；S3、利用自适应非局部双边滤波方法对图像进行二次滤波，去除噪声。本发明利用NASMF将椒盐噪声转化为类高斯噪声，然后利用自适应非局部高斯滤波算法去除噪声，根据图像非局部细节特征，自适应的趋向于l1范数或者l2范数，从而达到保护边缘和纹理细节的结果。
Description
技术领域
本发明涉及基于一种两步去除具有图像细节保护能力的非局部双边滤波图像椒盐噪声的方法，属于计算机图像处理领域。
背景技术
图像椒盐噪声尤其是高频椒盐噪声给图像的后续处理工作带来极大的挑战。有效地去除椒盐噪声并尽可能多的回复图像的细节信息仍然是一个极具挑战性的工作。
现有的方法是关注于图像的局部信息，利用局部信息估计检测到的噪声点。但是，图像的局部极值点和非平稳点无法精确的恢复图像细节信息，甚至产生错误估计。为了解决这一问题，王等人提出了利用非局部信息去除图像椒盐噪声的方法(具体参见XiaotianWang,Shanshan Shen,Guangming Shi,Yuannan Xu,Peiyu Zhang.Iterative non-localmeans filter for salt and pepper noise removal.J.Vis.Commun.Image R.38(2016)440–450,2018)，该方法利用非局部区域的特征相似度寻找相似非局部区域块，通过更新权重参数提高图像恢复精度。该方法考虑了利用非局部均值迭代消除椒盐噪声，但是没有空间距离因素对于重建效果的影响。
双边滤波(Bilateral Filter)是一种非线性滤波，该方法考虑了图像的空间邻近度与像素值的相似度，被广泛的用于高斯噪声的去除。Veerakumar et al.提出了使用自适应双边滤波算法去除高斯噪声(具体参见T.Veerakumar,S.Esakkirajan,I.Vennila,Anapproach to minimize very high density salt and pepper noise through trimmedglobal mean.Int.J.Compt.Appl.39(12):29–33 2012)，该方法通过修改空间临近高斯函数和像素值相似度高斯函数实现。在像素值相似度函数计算中，像素值之间的差异性通过L2范数度量，该范数可以有效的保护图像纹理细节信息。该方法不能直接适用于椒盐噪声。
并且，现有的上述两种方法以L2范数度量区域之间的相似性，L2范数无法精确地逼近边缘区域，易造成边缘模糊过平滑的现象。
发明内容
针对现有技术的不足，本发明提供了一种两步去除具有图像细节保护能力的非局部双边滤波图像椒盐噪声的方法；
本发明提出了一种可以自适应保护图像边缘和纹理细节的去噪算法，该算法利用NASMF将椒盐噪声转化为类高斯噪声，然后利用自适应非局部高斯滤波算法去除噪声，根据图像非局部细节特征，自适应的趋向于l1范数或者l2范数，从而达到保护边缘和纹理细节的结果。
术语解释：
1、灰度化，在RGB模型中，如果R＝G＝B时，则彩色表示一种灰度颜色，其中R＝G＝B的值叫灰度值，因此，灰度图像每个像素只需一个字节存放灰度值(又称强度值、亮度值)，灰度范围为0-255。一般有分量法最大值法平均值法加权平均法四种方法对彩色图像进行灰度化。
2、图像细节保护能力，由于现有方法问题，在复原图像的过程中，采用临近点或者相似区域对噪声点进行估值，但是估计结果与图像整体风格有所偏差，导致图像失真。这些失真的原因并非所复原图像的结构存在问题，而是由于图像在局部区域内的边缘被过于弱化或者强化导致的。过于弱化会导致过平滑现象的产生，常见于纹理区域；过于强化会导致人工伪迹的产生，常见于强边缘区域。因此，本发明所提出的图像细节保护能力指的是对图像局部的边缘进行修复，避免所复原的图像具有人工伪迹、过平滑等现象的产生。
3、NASMF算法，NASMF算法的全称为非局部自适应中值滤波算法。参见论文：KennyKal Vin Toh,and Nor Ashidi Mat Isa.Noise Adaptive Fuzzy Switching MedianFilter for Salt-and-Pepper Noise Reduction.IEEE SIGNAL PROCESSING LETTERS,17(3):281-284MARCH 2010。
本发明的技术方案为：
一种两步去除具有图像细节保护能力的非局部双边滤波图像椒盐噪声的方法，包括步骤如下：
S1、遍历图像，利用mask矩阵记录噪声点位置；
S2、利用NASMF算法和mask矩阵对噪声点的定位，对噪声点进行滤波将椒盐噪声转化为高斯噪声；
S3、利用自适应非局部双边滤波方法对图像进行二次滤波，去除噪声。
根据本发明优选的，所述步骤S1，包括：
(1)图像灰度化，得到灰度图像；
(2)对于步骤(1)得到的灰度图像，Nsalt＝255,NPepper＝0；Nsalt表示盐噪声(白点)，Npepper表示椒噪声(黑点)，因此，可以得到一个二值蒙版矩阵M，该二值蒙版矩阵M的尺寸与该灰度图像相同，其元素用来标记灰度图像中相应位置的像素点是噪声点还是非噪声点，如式(Ⅰ)所示：
式(Ⅰ)中，I(i,j)表示灰度图像中像素点(i,j)的灰度值；M(i,j)表示二值蒙版矩阵中(i,j)位置的值，当M(i,j)＝1表示灰度图像中像素点(i,j)是非噪声点；M(i,j)＝0表示灰度图像中像素点(i,j)是噪声点。
根据本发明优选的，所述步骤S2，包括：
(3)计算滤波窗口的中值，并计算噪声点的估计值：
当M(i,j)＝0时，相应的图像像素点(i,j)要进行滤波处理，选取(2s+1)×(2s+1)的滤波窗口，s为预先设定的窗口参数；
滤波窗口形式化的表示为：W2s+1(i,j)＝{I(i+m,j+n)|m,n∈(-s,...,0,,,s)}，W2s+1(i,j)表示以(i,j)位置为中心的、边长为2s+1的窗口；m,n分别表示距离当前窗口中心点(i,j)的行、列距离；I(i+m,j+n)表示图像I中(i+m,j+n)处的像素值；
记：
T2s+1(i,j)表示蒙版矩阵M中以(i,j)位置为中心的、边长为2s+1的窗口内部非噪声点的数量，M(i+m,j+n)表示蒙版矩阵M中(i+m,j+n)处的值；
预先设定区域内非噪声点数量的阈值t＝3，若T2s+1(i,j)>t，说明以(i,j)为中心，(2s+1)×(2s+1)区域内的非噪声点具有足够的数量，则使用式(II)计算噪声点的所在区域的中值即滤波窗口的中值并将结果保存到Imed中，图像Imed的尺寸与含噪图像I的尺寸一致：
Imed(i,j)＝median{I(i+m,j+n)|m,n∈(-s,...,0...,s)} (II)
式(II)中，median{}指的是将图像内以(i,j)为中心，边长为2s+1的窗口内部所有的点从小到大排列，记录中间数的值，并记录到图像Imed(i,j)位置上；
若T2s+1(i,j)<t，这说明当前的滤波窗口内没有足够多的非噪声点,则扩展滤波窗口，即s→s+1，重新利用公式(II)计算滤波窗口的中值，并记录到Imed中；
当遍历完矩阵M中所有值为0的点后，利用式(Ⅲ)计算噪声点的估计值：
Ifinal＝(1-F(i,j))I(i,j)+F(i,j)Imed(i,j) (Ⅲ)
式(Ⅲ)中，Ifinal指的是去噪之后的图像；
D(i,j)为滤波窗口内局部灰度变化最大值：
D(i,j)＝max{|I(i+k,j+l)-I(i,j)|,k,l∈(-2,2),(i+k,j+l)≠(i,j)}；
max{}返回的是该集合中的最大值，并记录在D(i,j)中；T1和T2为预先设定的阈值；
(4)将若干副原始不含噪图像Ioriginal分别加入椒盐噪声，加入的椒盐噪声的强度为60％～90％，得到对应的若干副含噪图像I，通过步骤(3)得到对应的噪声点的估计值Ifinal，对含噪图像I和原始不含噪图像Ioriginal做逐点的差值运算，利用式(Ⅳ)得到误差图像Ierr：
Ierr＝Ioriginal-Ifinal (Ⅳ)
对误差图像Ierr进行直方图统计，统计结果见附图1，差值结果服从统计意义上均值为0的高斯分布，并对这些直方图进行拟合，该过程由Matlab软件自带的fit函数完成。得到一个均值为0的高斯分布，如式(Ⅴ)所示：
式(Ⅴ)中，
表示定义为，即将nGauss定义为Gauss(0,σ2)，Gauss(0,σ2)表示均值为0，方差为σ2的高斯分布，将式(Ⅴ)带入式(Ⅳ)中，得到式(Ⅵ)：
Ioriginal＝Ifinal+nGauss (Ⅵ)
由此，将原椒盐噪声去噪问题转化为高斯噪声去除问题，即通过图像Ifinal获取原始图像Ioriginal。
根据本发明优选的，s的初始值为1，即初始时，所选窗口大小为3×3。
T1＝10，T2＝30。
根据本发明优选的，所述步骤S3，包括：
(5)在本发明中ap,q,r,s的计算是根据图像的非局部信息自适应完成的。由于图像的局部特征各不相同，因此，在使用
来度量像素间差异性时，使用同一ap,q,r,s值，显然不符合实际。因此，度量像素间差异性，求取不同的ap,q,r,s值，如式(Ⅶ)所示：
式(Ⅶ)中，y(p,q)和y(r,s)分别表示以(i,j)和(k,l)为中心的图像块中的像素点的灰度值，ε是一个极小数；
下面分析函数
的性质。如果y(p,q)与y(r,s)灰度值相近，即(y(p,q)-y(r,s))值较小，则ap,q,r,s较大，此时，自适应范数效果趋近于L2范数，反之，图像区域突变区域，ap,q,r,s较小，使得自适应范数效果趋近于L1范数，
(6)计算空间临近度度量函数ωd(i,j,k,l)及块像素值相似度度量函数ωr(i,j,k,l)，分别如式(Ⅷ)、式(Ⅸ)所示：
式(Ⅷ)、式(Ⅸ)中,σd和σr为平滑系数，Θ(i,j)表示以(i,j)为中心、大小为3×3的图像块，y(p,q)表示Θ(i,j)中的点,Θ(k,l)表示以(k,l)为中心、大小为3×3的图像块，y(r,s)表示Θ(k,l)中的点，
且ap,q,r,s>0，exp()表示指数函数，即exp(·)＝e(·)；
(7)为了进一步对步骤S2的结果进行滤波,提出了如下非局部双边滤波定义，通过非局部双边滤波公式对图像进行二次滤波，去除噪声，得到输出图像Ioutput，用来替换噪声点的估计值Ifinal，非局部双边滤波公式如式(X)所示：
式(Ⅶ)中,Ioutput(i,j)表示在(i,j)位置上的输出估计,Ifinal(k,l)表示图像Ifinal中(k,l)的像素值,B(i,j)表示Ifinal图像中以(i,j)为中心点图像块，图像块的尺寸为奇数，该图像块区域内的点记为(k,l)，即(k,l)∈B(i,j)。ω(i,j,k,l)表示Ifinal(k,l)对于估计Ioutput(i,j)的贡献度，即权重参数，贡献度ω(i,j,k,l)的计算分为两部分，即ωd(i,j,k,l)和ωr(i,j,k,l)。
根据本发明优选的，ε为10-3，σd设置为80乘以噪声强度，σr设置为1，图像块的尺寸为7×7。
本发明的有益效果为：
1、本发明提出了一种两步去除椒盐噪声的算法，第一步利用自适应中值滤波算法去除椒盐噪声，第二步利用自适应非局部双边滤波算法进一步抑制噪声，从而达到复原图像的目的。本发明利用自适应非局部双边滤波算法来修复中值滤波算法带来的误差，该算法通过计算图像块之间的相似度和距离来自适应的求解贡献权重。
3、本发明提出了一种在自适应范数中计算尺度算子的方法，通过使用该发放，自适应范数可以根据图像非局部特征选择逼近真实值的范数，在平滑区域，逼近L2范数，在边缘区域，逼近L1范数，从而避免人工伪迹和图像过平滑现象的产生。
附图说明
图1为本发明椒盐噪声转化后的图像与原图像差值图的直方图；
图2为基于Lena图像的试验结果示意图；
图3为基于Peppers图像的试验结果示意图；
图4为Lena图像在加有90％强度的椒盐噪声条件下，不同算法的去噪重建结果示意图；
具体实施方式
下面结合说明书附图和实施例对本发明作进一步限定，但不限于此。
实施例1
一种两步去除具有图像细节保护能力的非局部双边滤波图像椒盐噪声的方法，包括步骤如下：
S1、遍历图像，利用mask矩阵记录噪声点位置；
S2、利用NASMF算法和mask矩阵对噪声点的定位，对噪声点进行滤波将椒盐噪声转化为高斯噪声；
S3、利用自适应非局部双边滤波方法对图像进行二次滤波，去除噪声。
实施例2
根据实施例1所述的一种两步去除具有图像细节保护能力的非局部双边滤波图像椒盐噪声的方法，其区别在于：
步骤S1，包括：
(1)图像灰度化，得到灰度图像；
(2)对于步骤(1)得到的灰度图像，Nsalt＝255,NPepper＝0；Nsalt表示盐噪声(白点)，Npepper表示椒噪声(黑点)，因此，我们可以得到一个二值蒙版矩阵M，该二值蒙版矩阵M的尺寸与该灰度图像相同，其元素用来标记灰度图像中相应位置的像素点是噪声点还是非噪声点，如式(Ⅰ)所示：
式(Ⅰ)中，I(i,j)表示灰度图像中像素点(i,j)的灰度值；M(i,j)表示二值蒙版矩阵中(i,j)位置的值，当M(i,j)＝1表示灰度图像中像素点(i,j)是非噪声点；M(i,j)＝0表示灰度图像中像素点(i,j)是噪声点。
步骤S2，包括：
(3)计算滤波窗口的中值，并计算噪声点的估计值：
当M(i,j)＝0时，相应的图像像素点(i,j)要进行滤波处理，选取(2s+1)×(2s+1)的滤波窗口，s为预先设定的窗口参数；
滤波窗口形式化的表示为：W2s+1(i,j)＝{I(i+m,j+n)|m,n∈(-s,...,0,,,s)}，W2s+1(i,j)表示以(i,j)位置为中心的、边长为2s+1的窗口；m,n分别表示距离当前窗口中心点(i,j)的行、列距离；I(i+m,j+n)表示图像I中(i+m,j+n)处的像素值；
记：
T2s+1(i,j)表示蒙版矩阵M中以(i,j)位置为中心的、边长为2s+1的窗口内部非噪声点的数量，M(i+m,j+n)表示蒙版矩阵M中(i+m,j+n)处的值；
预先设定区域内非噪声点数量的阈值t＝3，若T2s+1(i,j)>t，说明以(i,j)为中心，(2s+1)×(2s+1)区域内的非噪声点具有足够的数量，则使用式(II)计算噪声点的所在区域的中值即滤波窗口的中值并将结果保存到Imed中，图像Imed的尺寸与含噪图像I的尺寸一致：
Imed(i,j)＝median{I(i+m,j+n)|m,n∈(-s,...,0...,s)} (II)
式(II)中，median{}指的是将图像内以(i,j)为中心，边长为2s+1的窗口内部所有的点从小到大排列，记录中间数的值，并记录到图像Imed(i,j)位置上；
若T2s+1(i,j)<t，这说明当前的滤波窗口内没有足够多的非噪声点,则扩展滤波窗口，即s→s+1，重新利用公式(II)计算滤波窗口的中值，并记录到Imed中；
当遍历完矩阵M中所有值为0的点后，利用式(Ⅲ)计算噪声点的估计值：
Ifinal＝(1-F(i,j))I(i,j)+F(i,j)Imed(i,j) (Ⅲ)
式(Ⅲ)中，Ifinal指的是去噪之后的图像；
D(i,j)为滤波窗口内局部灰度变化最大值：
D(i,j)＝max{|I(i+k,j+l)-I(i,j)|,k,l∈(-2,2),(i+k,j+l)≠(i,j)}；
max{}返回的是该集合中的最大值，并记录在D(i,j)中；T1和T2为预先设定的阈值；
(4)将若干副原始不含噪图像Ioriginal分别加入椒盐噪声，加入的椒盐噪声的强度为60％～90％，得到对应的若干副含噪图像I，通过步骤(3)得到对应的噪声点的估计值Ifinal，对含噪图像I和原始不含噪图像Ioriginal做逐点的差值运算，利用式(Ⅳ)得到误差图像Ierr：
Ierr＝Ioriginal-Ifinal (Ⅳ)
对误差图像Ierr进行直方图统计，统计结果见附图1，Ierr图像统计结果，横坐标为误差统计范围，纵坐标为归一化后的统计值(百分比)，图例中“-”前部分表示图片名称，“-”后部分表示噪声强度。差值结果服从统计意义上均值为0的高斯分布，并对这些直方图进行拟合，该过程由Matlab软件自带的fit函数完成。得到一个均值为0的高斯分布，如式(Ⅴ)所示：
式(Ⅴ)中，
表示定义为，即将nGauss定义为Gauss(0,σ2)，Gauss(0,σ2)表示均值为0，方差为σ2的高斯分布，将式(Ⅴ)带入式(Ⅳ)中，得到式(Ⅵ)：
Ioriginal＝Ifinal+nGauss (Ⅵ)
由此，将原椒盐噪声去噪问题转化为高斯噪声去除问题，即通过图像Ifinal获取原始图像Ioriginal。
s的初始值为1，即初始时，所选窗口大小为3×3。
T1＝10，T2＝30。
步骤S3，包括：
(5)在本发明中ap,q,r,s的计算是根据图像的非局部信息自适应完成的。由于图像的局部特征各不相同，因此，在使用
来度量像素间差异性时，使用同一ap,q,r,s值，显然不符合实际。
因此，度量像素间差异性，求取不同的ap,q,r,s值，如式(Ⅶ)所示：
式(Ⅶ)中，y(p,q)和y(r,s)分别表示以(i,j)和(k,l)为中心的图像块中的像素点的灰度值，ε是一个极小数；
下面分析函数
的性质。如果y(p,q)与y(r,s)灰度值相近，即(y(p,q)-y(r,s))值较小，则ap,q,r,s较大，此时，自适应范数效果趋近于L2范数，反之，图像区域突变区域，ap,q,r,s较小，使得自适应范数效果趋近于L1范数，
(6)计算空间临近度度量函数ωd(i,j,k,l)及块像素值相似度度量函数ωr(i,j,k,l)，分别如式(Ⅷ)、式(Ⅸ)所示：
式(Ⅷ)、式(Ⅸ)中,σd和σr为平滑系数，Θ(i,j)表示以(i,j)为中心、大小为3×3的图像块，y(p,q)表示Θ(i,j)中的点,Θ(k,l)表示以(k,l)为中心、大小为3×3的图像块，y(r,s)表示Θ(k,l)中的点，
且ap,q,r,s>0，exp()表示指数函数，即exp(·)＝e(·)；
(7)为了进一步对步骤S2的结果进行滤波,提出了如下非局部双边滤波定义，通过非局部双边滤波公式对图像进行二次滤波，去除噪声，得到输出图像Ioutput，用来替换噪声点的估计值Ifinal，非局部双边滤波公式如式(X)所示：
式(Ⅶ)中,Ioutput(i,j)表示在(i,j)位置上的输出估计,Ifinal(k,l)表示图像Ifinal中(k,l)的像素值,B(i,j)表示Ifinal图像中以(i,j)为中心点图像块，图像块的尺寸为奇数，该图像块区域内的点记为(k,l)，即(k,l)∈B(i,j)。ω(i,j,k,l)表示Ifinal(k,l)对于估计Ioutput(i,j)的贡献度，即权重参数，贡献度ω(i,j,k,l)的计算分为两部分，即ωd(i,j,k,l)和ωr(i,j,k,l)。
ε为10-3，σd设置为80乘以噪声强度，σr设置为1，图像块的尺寸为7×7。
本实施例提出的算法在两幅常用图像(Lena和Peppers)上进行了测试，分别对两幅图像施加10％到90％强度的椒盐噪声，数值越大，噪声强度越强，然后利用本本实施例提出的算法对含噪图像进行去噪处理。并与现有的多种常用算法(如中值滤波及其改进算法(MF，NASMF，ACWMF，DAMF)，基于决策的算法(DBA)以及基于模糊集的去噪算法(FSAP)等)进行了比对，测试结果通过PSNR,SSIM和IEF等三种评价标准进行客观评价，表1是PSNR和SSIM两个指标的评价结果；
表1
表1为本实施例提出算法与传统算法的数值结果在lena和peppers图像上不同噪声强度下对比表。”/”前为PSNR值，后为SSIM值，这两个数值越大，说明重建效果越好。
图2和图3是IEF指标的评价结果，图2中，不同的线表示不同的对比算法，横坐标是不同的噪声强度，纵坐标不同强度噪声下，不同算法去噪结果的IEF评价值。图3中，不同的线表示不同的对比算法，横坐标是不同的噪声强度，纵坐标不同强度噪声下，不同算法去噪结果的IEF评价值。图4为Lena图像在加有90％强度的椒盐噪声条件下，不同算法的去噪重建结果示意图；图4中，第一行从左往右依次是Lena原图，MF方法,ACWMF方法,DBA方法和NAFSM方法；第二行从左往右依次是NASEPF方法,NASNLM方法,DAMF方法,FSAP方法以及本实施例提出的算法。其中，PSNR值、SSIM值和IEF值越高说明去噪效果越好。从评价结果(表1和图2，图3、图4)来看，本实施例提出的方法在客观评价指标上，均好于对比算法。
Claims (4)
1.一种两步去除具有图像细节保护能力的非局部双边滤波图像椒盐噪声的方法，其特征在于，包括步骤如下：
S1、遍历图像，利用mask矩阵记录噪声点位置；
S2、利用NASMF算法和mask矩阵对噪声点的定位，对噪声点进行中值滤波将椒盐噪声转化为高斯噪声；
所述步骤S2，包括：
(3)计算滤波窗口的中值，并计算噪声点的估计值：
当M(i,j)＝0时，相应的图像像素点(i,j)要进行滤波处理，选取(2s+1)×(2s+1)的滤波窗口，s为预先设定的窗口参数；
滤波窗口形式化的表示为：W2s+1(i,j)＝{I(i+m,j+n)|m,n∈(-s,...,0,,,s)}，W2s+1(i,j)表示以(i,j)位置为中心的、边长为2s+1的窗口；m,n分别表示距离当前窗口中心点(i,j)的行、列距离；I(i+m,j+n)表示图像I中(i+m,j+n)处的像素值；
记：
T2s+1(i,j)表示蒙版矩阵M中以(i,j)位置为中心的、边长为2s+1的窗口内部非噪声点的数量，M(i+m,j+n)表示蒙版矩阵M中(i+m,j+n)处的值；
预先设定区域内非噪声点数量的阈值t，若T2s+1(i,j)>t，说明以(i,j)为中心，(2s+1)×(2s+1)区域内的非噪声点具有足够的数量，则使用式(II)计算噪声点的所在区域的中值即滤波窗口的中值并将结果保存到Imed中，图像Imed的尺寸与含噪图像I的尺寸一致：
Imed(i,j)＝median{I(i+m,j+n)|m,n∈(-s,...,0...,s)} (II)
式(II)中，median{}指的是将图像内以(i,j)为中心，边长为2s+1的窗口内部所有的点从小到大排列，记录中间数的值，并记录到图像Imed(i,j)位置上；
若T2s+1(i,j)<t，则扩展滤波窗口，即s→s+1，重新利用公式(II)计算滤波窗口的中值，并记录到Imed中；
当遍历完矩阵M中所有值为0的点后，利用式(Ⅲ)计算噪声点的估计值：
Ifinal＝(1-F(i,j))I(i,j)+F(i,j)Imed(i,j) (Ⅲ)
式(Ⅲ)中，Ifinal指的是去噪之后的图像；
D(i,j)为滤波窗口内局部灰度变化最大值：
D(i,j)＝max{|I(i+k,j+l)-I(i,j)|,k,l∈(-2,2),(i+k,j+l)≠(i,j)}；
max{}返回的是集合中的最大值，并记录在D(i,j)中；T1和T2为预先设定的阈值；
(4)将若干副原始不含噪图像Ioriginal分别加入椒盐噪声，加入的椒盐噪声的强度为60％～90％，得到对应的若干副含噪图像I，通过步骤(3)得到对应的噪声点的估计值Ifinal，对含噪图像I和原始不含噪图像Ioriginal做逐点的差值运算，利用式(Ⅳ)得到误差图像Ierr：
Ierr＝Ioriginal-Ifinal (Ⅳ)
对误差图像Ierr进行直方图统计，差值结果服从统计意义上均值为0的高斯分布，并对这些直方图进行拟合，得到一个均值为0的高斯分布，如式(Ⅴ)所示：
式(Ⅴ)中，
表示定义为，即将nGauss定义为Gauss(0,σ2)，Gauss(0,σ2)表示均值为0，方差为σ2的高斯分布，将式(Ⅴ)带入式(Ⅳ)中，得到式(Ⅵ)：
Ioriginal＝Ifinal+nGauss (Ⅵ)
由此，将原椒盐噪声去噪问题转化为高斯噪声去除问题，即通过图像Ifinal获取原始图像Ioriginal；
S3、利用自适应非局部双边滤波方法对图像进行二次滤波，去除噪声；
所述步骤S3，包括：
(5)度量像素间差异性，求取不同的ap，q，r，s值，如式(Ⅶ)所示：
式(Ⅶ)中，y(p,q)和y(r,s)分别表示以(i,j)和(k,l)为中心的图像块中的像素点的灰度值，ε是一个极小数；
(6)计算空间临近度度量函数ωd(i,j,k,l)及块像素值相似度度量函数ωr(i,j,k,l)，分别如式(Ⅷ)、式(Ⅸ)所示：
式(Ⅷ)、式(Ⅸ)中,σd和σr为平滑系数，Θ(i,j)表示以(i,j)为中心、大小为3×3的图像块，y(p,q)表示Θ(i,j)中的点,Θ(k,l)表示以(k,l)为中心、大小为3×3的图像块，y(r,s)表示Θ(k,l)中的点，
且ap,q,r,s>0，exp()表示指数函数，即exp(·)＝e(·)；
(7)通过非局部双边滤波公式对图像进行二次滤波，去除噪声，得到输出图像Ioutput，用来替换噪声点的估计值Ifinal，非局部双边滤波公式如式(X)所示：
式(Ⅶ)中,Ioutput(i,j)表示在(i,j)位置上的输出估计,Ifinal(k,l)表示图像Ifinal中(k,l)的像素值,B(i,j)表示Ifinal图像中以(i,j)为中心点图像块，图像块的尺寸为奇数，该图像块区域内的点记为(k,l)，即(k,l)∈B(i,j)。
2.根据权利要求1所述的一种两步去除具有图像细节保护能力的非局部双边滤波图像椒盐噪声的方法，其特征在于，所述步骤S1，包括：
(1)图像灰度化，得到灰度图像；
(2)对于步骤(1)得到的灰度图像，Nsalt＝255,NPepper＝0；Nsalt表示盐噪声，Npepper表示椒噪声，得到一个二值蒙版矩阵M，该二值蒙版矩阵M的尺寸与该灰度图像相同，其元素用来标记灰度图像中相应位置的像素点是噪声点还是非噪声点，如式(Ⅰ)所示：
式(Ⅰ)中，I(i,j)表示灰度图像中像素点(i,j)的灰度值；M(i,j)表示二值蒙版矩阵中(i,j)位置的值，当M(i,j)＝1表示灰度图像中像素点(i,j)是非噪声点；M(i,j)＝0表示灰度图像中像素点(i,j)是噪声点。
3.根据权利要求1所述的一种两步去除具有图像细节保护能力的非局部双边滤波图像椒盐噪声的方法，其特征在于，s的初始值为1，即初始时，所选窗口大小为3×3；
T1＝10，T2＝30。
4.根据权利要求1所述的一种两步去除具有图像细节保护能力的非局部双边滤波图像椒盐噪声的方法，其特征在于，ε为10-3，σd设置为80乘以噪声强度，σr设置为1，图像块的尺寸为7×7。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910777064.5A CN110458783B (zh) | 2019-08-22 | 2019-08-22 | 一种基于具有图像细节保护能力的非局部双边滤波的椒盐噪声去除方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910777064.5A CN110458783B (zh) | 2019-08-22 | 2019-08-22 | 一种基于具有图像细节保护能力的非局部双边滤波的椒盐噪声去除方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110458783A CN110458783A (zh) | 2019-11-15 |
| CN110458783B true CN110458783B (zh) | 2023-01-31 |
Family
ID=68488517
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910777064.5A Active CN110458783B (zh) | 2019-08-22 | 2019-08-22 | 一种基于具有图像细节保护能力的非局部双边滤波的椒盐噪声去除方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110458783B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111242137B (zh) * | 2020-01-13 | 2023-05-26 | 江西理工大学 | 一种基于形态成分分析的椒盐噪声滤波方法、装置 |
| CN114155161B (zh) * | 2021-11-01 | 2023-05-09 | 富瀚微电子(成都)有限公司 | 图像去噪方法、装置、电子设备与存储介质 |
| CN113989168B (zh) * | 2021-11-02 | 2024-05-14 | 华北电力大学（保定） | 一种针对椒盐噪声的自适应非局部均值滤波方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104217405A (zh) * | 2014-09-23 | 2014-12-17 | 闽江学院 | 融合局部与全局信息的图像椒盐噪声滤波方法 |
| CN107038688A (zh) * | 2017-03-01 | 2017-08-11 | 武汉科技大学 | 基于Hessian矩阵的图像噪点检测及去噪方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7856150B2 (en) * | 2007-04-10 | 2010-12-21 | Arcsoft, Inc. | Denoise method on image pyramid |
| EA201200255A1 (ru) * | 2012-02-22 | 2013-08-30 | Закрытое Акционерное Общество "Импульс" | Способ подавления шума цифровых рентгенограмм |
-
2019
- 2019-08-22 CN CN201910777064.5A patent/CN110458783B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104217405A (zh) * | 2014-09-23 | 2014-12-17 | 闽江学院 | 融合局部与全局信息的图像椒盐噪声滤波方法 |
| CN107038688A (zh) * | 2017-03-01 | 2017-08-11 | 武汉科技大学 | 基于Hessian矩阵的图像噪点检测及去噪方法 |
Non-Patent Citations (1)
| Title |
|---|
| removal of salt-and pepper noise from bimodal images using a non-symmetry and anti-packing model;hu jiang.etc;《medical imaging and health information》;20190130;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110458783A (zh) | 2019-11-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110458783B (zh) | 一种基于具有图像细节保护能力的非局部双边滤波的椒盐噪声去除方法 | |
| Vijaykumar et al. | Fast and efficient algorithm to remove gaussian noise in digital images | |
| Harikiran et al. | Impulse noise removal in digital images | |
| Deepa et al. | Comparative analysis of noise removal techniques in MRI brain images | |
| CN103400357A (zh) | 一种去除图像椒盐噪声的方法 | |
| Erkan et al. | Improved adaptive weighted mean filter for salt-and-pepper noise removal | |
| Bhateja et al. | Performance improvement of decision median filter for suppression of salt and pepper noise | |
| Tiwari et al. | Image denoising using spatial gradient based bilateral filter and minimum mean square error filtering | |
| Kornilov et al. | Blind quality assessment for slice of microtomographic image | |
| Dosselmann et al. | A formal assessment of the structural similarity index | |
| Irum et al. | Salt and pepper noise removal filter for 8-bit images based on local and global occurrences of Grey levels as selection indicator | |
| CN111652810A (zh) | 一种基于小波域奇异值差分模型的图像去噪方法 | |
| Mohammed | An improved median filter based on efficient noise detection for high quality image restoration | |
| Biswas et al. | Speckle noise reduction from medical images using Gaussian fuzzy membership function | |
| Shrivastava et al. | Noise removal using first order neighborhood mean filter | |
| Bagchi et al. | Investigation of different spatial filters performance toward mammogram de-noising | |
| CN104966271A (zh) | 基于生物视觉感受野机制的图像去噪方法 | |
| Thanh et al. | Distorted image reconstruction method with trimmed median | |
| Zhu et al. | Two-stage non-local means filtering with adaptive smoothing parameter | |
| Sarangi et al. | A comparative study of filters for denoising mammograms | |
| Lyakhov et al. | Bilateral and Median Filter Combination for High-Quality Cleaning of Random Impulse Noise in Images | |
| Biswas | Impulse Noise Suppression in Color Images Using Median Filter and Deep Learning | |
| Kashyap et al. | Edge detection using sobel method with median filter | |
| Shrivastava et al. | Removal of impulse noise using first order neighborhood mean filter | |
| Ali et al. | A proposed de-noising algorithm |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |