CN104217405A - 融合局部与全局信息的图像椒盐噪声滤波方法 - Google Patents
融合局部与全局信息的图像椒盐噪声滤波方法 Download PDFInfo
- Publication number
- CN104217405A CN104217405A CN201410489377.8A CN201410489377A CN104217405A CN 104217405 A CN104217405 A CN 104217405A CN 201410489377 A CN201410489377 A CN 201410489377A CN 104217405 A CN104217405 A CN 104217405A
- Authority
- CN
- China
- Prior art keywords
- noise
- image
- pixel
- walkaway
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 44
- 239000006002 Pepper Substances 0.000 title claims abstract description 38
- 238000011084 recovery Methods 0.000 claims abstract description 14
- 238000001914 filtration Methods 0.000 claims description 25
- 239000011159 matrix material Substances 0.000 claims description 15
- 230000011218 segmentation Effects 0.000 claims description 10
- 230000001174 ascending effect Effects 0.000 claims description 4
- 241000196324 Embryophyta Species 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims description 3
- 238000009966 trimming Methods 0.000 claims description 3
- 230000004927 fusion Effects 0.000 claims 1
- 230000000694 effects Effects 0.000 abstract description 21
- 238000001514 detection method Methods 0.000 abstract description 6
- 241000209094 Oryza Species 0.000 description 13
- 235000007164 Oryza sativa Nutrition 0.000 description 13
- 235000009566 rice Nutrition 0.000 description 13
- 238000002474 experimental method Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 5
- 230000003044 adaptive effect Effects 0.000 description 4
- 238000012113 quantitative test Methods 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 230000009931 harmful effect Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 244000144992 flock Species 0.000 description 2
- 238000011158 quantitative evaluation Methods 0.000 description 2
- 150000003839 salts Chemical class 0.000 description 2
- 238000009333 weeding Methods 0.000 description 2
- 235000002566 Capsicum Nutrition 0.000 description 1
- 241001347978 Major minor Species 0.000 description 1
- 235000016761 Piper aduncum Nutrition 0.000 description 1
- 235000017804 Piper guineense Nutrition 0.000 description 1
- 244000203593 Piper nigrum Species 0.000 description 1
- 235000008184 Piper nigrum Nutrition 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 235000013339 cereals Nutrition 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000003909 pattern recognition Methods 0.000 description 1
- 230000001681 protective effect Effects 0.000 description 1
- 238000013138 pruning Methods 0.000 description 1
Landscapes
- Image Analysis (AREA)
- Image Processing (AREA)
Abstract
本发明涉及一种融合局部与全局信息的图像椒盐噪声滤波方法。该方法：首先，读入噪声图像，并进行噪声密度估计；其次，进行第一次噪声检测；再而，对所述步骤S2的噪声检测结果进行校正，并进行第一次噪声恢复；最后，进行第二次噪声检测，并进行第二次噪声恢复，得到去噪恢复的图像。本发明改善了图像去噪的效果,对自身含有灰度为0、255的非噪声极值点的图像改进幅度更明显。
Description
技术领域
本发明属于图像处理技术领域,主要应用于对被椒盐噪声污染的图像进行去噪恢复，特别是一种融合局部与全局信息的图像椒盐噪声滤波方法。
背景技术
椒盐噪声是一种常见的噪声，它的产生因素多种多样，比如，传感器硬件故障、通信信道受干扰。椒盐噪声表现为图像中与近邻像素点格格不入的黑色(灰度为0)或者白色 (灰度为255) 像素点。噪声点与近邻像素点间存在灰度突变，导致了较大的梯度值，使其与边缘点容易发生混淆，给图像分析尤其是边缘检测带来了很大的困难。
中值滤波是经典的椒盐去噪方法，其效果依赖于滤波窗口大小的选择。窗口小，去噪效果差，图像细节保留能力强；窗口大，去噪效果更佳，细节容易损失，图像模糊加重，如何平衡两者是个难题。为此，研究者提出了一系列改进型的中值滤波算法。比如，自适应中值滤波[1]根据局部窗口的内容自动选择滤波窗口的大小，改善了去噪效果，但其细节保护能力较弱。
开关滤波是图像椒盐噪声去除领域新近的研究热点。所谓开关是指先检测图像中的椒盐噪声点，然后仅对噪声点进行滤波恢复，有助于保护图像细节。SAWM[2]和MDWF[3]是去噪效果较好的两种开关滤波算法。SAWM和MDWF采用不同形式的方向灰度差异来检测图像中的椒盐噪声点,它们各自的噪声检测策略分别记为SAWM_D以及MDWF_D。SAWM_D进行噪声检测的具体过程如下：
(1) 任意给定一个像素点pi, j ，取以其为中心的Ld ×Ld 邻域灰度窗口Ω；
(2) 将窗口Ω内的像素点灰度按升序排列构成一个向量F，分别确定其第r个元素Fr 和第z-r+1个元素Fz -r+1，其中，z=Ld ×Ld ，1≤r≤(z-1)/2；
(3) 剔除掉窗口Ω内灰度小于等于Fr 或者大于等于Fz -r+1的像素点，构成新集合Ω1，其原理在于被剔除的像素点灰度位于窗口Ω内所有像素点所在灰度范围的两端，成为噪声点的概率比较大；
(4) 求窗口Ω1在方向k上的像素点灰度构成的集合
，这里只考虑水平、垂直、主副对角线4个方向；
(5) 根据下式计算中心像素点pi, j 在4个方向上与近邻像素点之间的加权绝对灰度差异之和：
(1)
(2)
式中，Φ代表空集合，ds ,t 代表了像素点pi, j 与其近邻像素点ps
,t 之间的绝对灰度差异；
(6) 计算4个方向上加权绝对灰度差异之和的最小值；
.
(3)
(7) 识别椒盐噪声点
(4)
其中P表示检测到的椒盐噪声点集合，T 1是一个参数；
MDWF_D进行噪声检测的具体过程如下：
(1) 任意给定一个像素点pi, j ，取以其为中心的m×m局部窗口Ω，fi , j 为pi , j 的灰度，求pi , j 与方向k上近邻像素点间绝对灰度差异之和：
(5)
(6)
(7)
其中1≤k≤4表示方向索引，Ω(k)为Ω在方向k上的内容，ws ,t 表示近邻像素点ps ,t 的权重；
（2）求方向绝对灰度差异之和的最小值：
(8)
（3）识别椒盐噪声点
(9)
其中P表示检测到的椒盐噪声点集合，T是一个参数。
SAWM_D和MDWF_D均利用了方向灰度差异来检测图像中的噪声点，其基本原理在于：噪声点没有方向性，而边缘像素点有一定的方向性，这种方向性的有无会造成方向加权灰度差异之和最小值的巨大差异，由此可以用来区分噪声点和边缘点。当然，这种方向加权灰度差异之和的最小值也可以用来区分噪声点和一般的非噪声点。因为噪声点与其周围邻域像素点往往存在较大的灰度差异，导致方向加权灰度差异之和的最小值也较大。而非噪声点与其周围邻域像素点间不存在灰度突变，灰度差异较小，方向加权灰度差异之和的最小值自然也较小。因此，方向加权灰度差异之和的最小值可以用来区分噪声点和非噪声点。
上述噪声检测的思路虽然基本正确，但在具体检测过程中仍然存在如下问题：
(1) 在SAWM_D的噪声检测过程中，通常设置参数r为1。此时，Fr 和Fz -r+1分别代表局部邻域窗口Ω的灰度极小值和极大值。这里的灰度极小值可能不为0，灰度极大值可能不为255。SAWM_D在剔除掉窗口Ω内灰度小于等于Fr 或者大于等于Fz -r+1的像素点过程中，可能将灰度不为0或255的非噪声点剔除掉，给方向加权灰度差异的计算带来不良影响。此外，即便窗口Ω内的灰度极小值和极大值刚好为0和255，直接剔除灰度小于等于0或者灰度大于等于255的像素点的做法也略显粗糙，容易将灰度刚好也为0或255的非噪声点一并删除掉，这也将给方向加权灰度差异的计算带来不良影响。方向加权灰度差异计算的不准确将直接影响噪声检测的准确性。实验结果表明SAWM_D往往存在较大的虚警率FA (False Alarm)，FA是指将非噪声点错误地识别为噪声点；
(2) MDWF_D在计算方向加权灰度差异的过程中未排除邻域内的任何像素点，如果邻域内存在噪声点，这些噪声点将对方向加权灰度差异的计算带来不良影响。方向加权灰度差异计算的不准确将直接影响噪声检测的准确性。实验结果表明MDWF_D往往存在较大的漏检率MD (Miss Detection)，MD是指将噪声点错误地识别为非噪声点；
(3) 现有的开关滤波算法包括SAWM和MDWF均只是利用像素点的局部邻域灰度信息来辨识噪声点，未考虑全局图像信息，对自身包含灰度为0或255的非噪声极值点块的图像，检测准确性急剧下降。而这类图像也经常出现，以著名的伯克利图像数据库BSDS300为例，300幅彩色图转换为灰度图后，有138幅图自身包含灰度为0或255的非噪声极值点块。
参考文献：
[1] Chan R H, Ho C W, Nikolova M. Salt-and-pepper noise removal by
median-type noise detectors and detail-preserving regularization [J]. IEEE
Transactions on Image Processing, 2005, 14(10): 1479~1485.
[2] Zhang X. Xiong Y. Impulse noise removal using directional difference
based noise detector and adaptive weighted mean filter. IEEE Signal Processing
Letters, 2009. 16(4), 295~298.
[3] Li Z, Liu G, Xu Y, Cheng Y. Modified directional weighted filter for
removal of salt & pepper noise, Pattern Recognition Letters, 2014, 40(1):
113~120.。
发明内容
本发明的目的在于提供一种能改善图像去噪的效果，特别是对自身含有灰度为0、255的非噪声极值点的图像改善幅度更明显的融合局部与全局信息的图像椒盐噪声滤波方法。
为实现上述目的，本发明的技术方案是：一种融合局部与全局信息的图像椒盐噪声滤波方法，包括如下步骤，
步骤S1：读入噪声图像，进行噪声密度估计；
步骤S2：进行第一次噪声检测；
步骤S3：对所述步骤S2的噪声检测结果进行校正，并进行第一次噪声恢复；
步骤S4：进行第二次噪声检测；
步骤S5：进行第二次噪声恢复，得到去噪图像。
在本发明实施例中，所述步骤S1中，所述噪声密度估计的具体过程如下：
步骤S11：将去噪图像的行、列分别分成w等份，计算每个等份对应的行、列数，若最后一个等份的行、列数不够，采用镜面反射的方式进行补足，最终将图像分成w×w个子块；
步骤S12：统计每个子块中灰度为0和255的像素点出现的概率，构造向量；
步骤S13：对向量中的元素按升序排序，排序后的向量记为；
步骤S14：估计整幅图像的噪声密度为：
相应地，整幅图像噪声点个数的估计值为：；
其中，mod代表取余，N表示图像中像素点的总数。
在本发明实施例中，所述步骤S2中，所述第一次噪声检测的具体过程如下：
步骤S21：采用算法SAWM噪声检测策略的改进版本SAWM_MD对噪声图像进行第一次噪声检测，即先计算出每个像素点pi,
j 对应的方向加权绝对灰度差异之和的最小值di ,j ，用所有像素点的di ,j 构建矩阵Md ；
步骤S22：对矩阵Md 中的元素按降序排列，构成新的向量，向量中的一个元素对应图像中的一个像素点；
步骤S23：取向量中前ne 个元素对应的像素点为噪声点，构成噪声点集合P，产生噪声检测结果标示矩阵Bm ：
其中，标示矩阵Bm 是一个二值化矩阵，1和0分别表示检测到的噪声点和非噪声点。
在本发明实施例中，所述步骤S3中，所述对噪声检测结果进行校正的具体过程如下：
步骤S31：根据椒盐噪声点灰度只可能为0或者255的事实，剔除掉噪声检测结果中灰度不为0或者255的噪声点，其形式化描述为：
步骤S32：若噪声密度估计结果De 最接近于9种噪声级别的前5种，则执行步骤S33~S36；否则，校正过程结束；
步骤S33：分别获取灰度0和255对应的阈值分割结果：
其中，1和0分别代表阈值分割结果中的目标像素和背景像素；
步骤S34：搜索8连通拓扑结构下，阈值分割结果seg 0和seg 255中目标像素构成的连通分量，将所有连通分量组成一个集合C=，其中，Ci 代表第i个连通分量对应像素点的集合，m代表连通分量的个数；
步骤S35：按如下条件移除掉可能由灰度为0或者255的噪声点构成的连通分量，得到新的连通分量集合：
其中，ni 代表连通分量Ci 包含的像素点个数；
参数：；其中，round代表四舍五入；
步骤S36：根据连通分量对Bm 进行如下校正：
。
在本发明实施例中，所述9种噪声级别分别为10%、20%、30%、40%、50%、60%、70%、80%、90%。
在本发明实施例中，所述步骤S4中，所述第二次噪声检测的具体过程如下：
步骤S41：采用算法MDWF的噪声检测策略MDWF_D对经过第一次噪声恢复的图像进行再次的噪声检测，判定椒盐噪声点的条件为：
其中，P表示检测到的椒盐噪声点集合，T为一参数；
步骤S42：对T值进行如下设置：
。
在本发明实施例中，所述第一次噪声恢复和第二次噪声恢复均采用算法SAWM的噪声点恢复策略的改进版本SAWM_MR进行噪声恢复，其具体过程如下：
首先，以噪声点为中心，初始化滤波窗口大小w 1=3，迭代地扩充滤波窗口大小为w 1=
w 1+2，直至满足滤波窗口内非噪声点的个数大于等于2或者w 1≥21；
然后，取噪声点为中心的w 1×w 1滤波窗口Ω，剔除掉其中的噪声点，将剩余非噪声点构成集合U；如果U为空集合，则将其重置为U=Ω；将集合U中像素点的加权灰度均值作为中心噪声点的滤波输出，其形式化描述为：
其中，fs ,t 代表像素点ps, t 的灰度，ω s, t 代表ps, t 的加权系数，加权函数的定义如下：
其中，Nu 代表集合U中像素点的个数，pm,n 和pk,l 代表集合U中的像素点，fm,n 和fk,l 分别代表像素点pm,n 和pk,l 的灰度。
相较于现有技术，本发明具有以下有益效果：
1、本发明算法明显改善了图像椒盐噪声去除的效果，特别是在含非噪声极值点块的图像上改善幅度更大；
2、本发明算法运行速度明显快于MDWF，后者的运行时间是前者的5~6倍。
附图说明
图1为本发明算法流程图。
图2为Lena图像和Rice图像的原图。
图3为30%噪声密度下Lena图像的噪声图及各算法去噪结果比较图。
图4为60%噪声密度下Lena图像的噪声图及各算法去噪结果比较图。
图5为90%噪声密度下Lena图像的噪声图及各算法去噪结果比较图。
图6为30%噪声密度下Rice图像的噪声图及各算法去噪结果比较图。
图7为60%噪声密度下Rice图像的噪声图及各算法去噪结果比较图。
图8为90%噪声密度下Rice图像的噪声图及各算法去噪结果比较图。
具体实施方式
下面结合附图，对本发明的技术方案进行具体说明。
如图1所示，本发明提出一种融合局部与全局信息的图像椒盐噪声滤波方法，包括如下步骤，
步骤S1：读入噪声图像，进行噪声密度估计；
步骤S2：进行第一次噪声检测；
步骤S3：对所述步骤S2的噪声检测结果进行校正，并进行第一次噪声恢复；
步骤S4：进行第二次噪声检测；
步骤S5：进行第二次噪声恢复，得到去噪图像。
在本发明一实施例中，所述步骤S1中，所述噪声密度估计的具体过程如下：
步骤S11：将待去噪图像的行、列分别分成w等份，计算每个等份对应的行、列数，若最后一个等份的行、列数不够，采用镜面反射的方式进行补足，最终将图像分成w×w个子块；
步骤S12：统计每个子块中灰度为0和255的像素点出现的概率，构造向量；
步骤S13：对向量中的元素按升序排序，排序后的向量记为；
步骤S14：估计整幅图像的噪声密度为：
相应地，整幅图像噪声点个数的估计值为：；
其中，mod代表取余，N表示图像中像素点的总数。
在本发明一实施例中，所述步骤S2中，所述第一次噪声检测的具体过程如下：
步骤S21：采用算法SAWM噪声检测策略的改进版本SAWM_MD对噪声图像进行第一次噪声检测，即先计算出每个像素点pi,
j 对应的方向加权绝对灰度差异之和的最小值di ,j ，用所有像素点的di ,j 构建矩阵Md ；
步骤S22：对矩阵Md 中的元素按降序排列，构成新的向量，向量中的一个元素对应图像中的一个像素点；
步骤S23：取向量中前ne 个元素对应的像素点为噪声点，构成噪声点集合P，产生噪声检测结果标示矩阵Bm ：
其中，标示矩阵Bm 是一个二值化矩阵，1和0分别表示检测到的噪声点和非噪声点。
在本发明一实施例中，所述步骤S3中，所述对噪声检测结果进行校正的具体过程如下：
步骤S31：根据椒盐噪声点灰度只可能为0或者255的事实，剔除掉噪声检测结果中灰度不为0或者255的噪声点，其形式化描述为：
步骤S32：若噪声密度估计结果De 最接近于9种噪声级别的前5种，则执行步骤S33~S36；否则，校正过程结束；所述9种噪声级别分别为10%、20%、30%、40%、50%、60%、70%、80%、90%；
步骤S33：分别获取灰度0和255对应的阈值分割结果：
其中，1和0分别代表阈值分割结果中的目标像素和背景像素；
步骤S34：搜索8连通拓扑结构下，阈值分割结果seg 0和seg 255中目标像素构成的连通分量，将所有连通分量组成一个集合C=，其中，Ci 代表第i个连通分量对应像素点的集合，m代表连通分量的个数；
步骤S35：按如下条件移除掉可能由灰度为0或者255的噪声点构成的连通分量，得到新的连通分量集合：
其中，ni 代表连通分量Ci 包含的像素点个数；
参数：；其中，round代表四舍五入；
步骤S36：根据连通分量对Bm 进行如下校正：
；
上述步骤S33~S36是把灰度为0和255的全局阈值分割结果作为全局图像信息对噪声检测结果进行校正，其主要作用对象是自身含有灰度为0和255的非噪声极值点块的图像。其基本原理在于：当噪声密度处于中下水平时，聚集在一起构成连通分量的灰度为0的噪声像素点个数应当不大，同理，聚集在一起构成连通分量的灰度为255的噪声像素点个数也应当不大。因此，我们把连通分量中像素点个数大于参数的视为由灰度为0或255构成的非噪声极值点块。这些非噪声极值点块是真实的图像块，而不是噪声，故将其从噪声检测结果中加以校正。
在本发明一实施例中，所述第一次噪声恢复和第二次噪声恢复均采用算法SAWM的噪声点恢复策略的改进版本SAWM_MR进行噪声恢复，其具体过程如下：
首先，以噪声点为中心，初始化滤波窗口大小w 1=3，迭代地扩充滤波窗口大小为w 1=
w 1+2，直至满足滤波窗口内非噪声点的个数大于等于2或者w 1≥21；附加的条件w 1≥21是为了防止高密度噪声下滤波窗口扩张太大导致局部性的丢失。
然后，取噪声点为中心的w 1×w 1滤波窗口Ω，剔除掉其中的噪声点，将剩余非噪声点构成集合U；如果U为空集合，则将其重置为U=Ω；将集合U中像素点的加权灰度均值作为中心噪声点的滤波输出，其形式化描述为：
其中，fs ,t 代表像素点ps, t 的灰度，ω s, t 代表ps, t 的加权系数，加权函数的定义如下：
其中，Nu 代表集合U中像素点的个数，pm,n 和pk,l 代表集合U中的像素点，fm ,n 和fk
,l 分别代表像素点pm,n 和pk,l 的灰度。
在本发明一实施例中，所述步骤S4中，所述第二次噪声检测的具体过程如下：
步骤S41：采用算法MDWF的噪声检测策略MDWF_D对经过第一次噪声恢复的图像进行再次的噪声检测，判定椒盐噪声点的条件为：
其中，P表示检测到的椒盐噪声点集合，T为一参数；
步骤S42：对T值进行如下设置：
其中，n代表迭代次数，其取值范围为5~10。
而在本发明设计的滤波算法中，MDWF的噪声检测策略只在第二次噪声检测时用到，无需迭代，故假定其最大迭代次数为10的基础上，将参数T的设置更新为：
。
上述流程中，第二次噪声恢复与第一次噪声恢复一样，采用算法SAWM噪声点恢复策略的改进版本SAWM_MR来实现，具体过程已经在第一次噪声恢复时阐述，这里不再赘述。
为了评价图像去噪算法的性能，我们采用客观的定量评价和主观的定性评价相结合的做法。峰值信噪比(PSNR)和平均结构相似性(MSSIM)被选定为定量评价指标，其中，参考图像为不含噪声的原始图像。PSNR的定义如下：
其中，N表示图像像素总数，xi,j 和ri,j 分别表示参考图像和去噪恢复后图像中像素点pi,j 的灰度。PSNR的值越大，说明去噪效果越好。MSSIM的定义如下：
其中，X和Y分别表示参考图像和去噪恢复后的图像，M代表图像分块窗口的数量，xj 和yj 代表X和Y中第j个分块窗口的内容，和代表各自的均值，和为各自的标准偏差，代表它们的协方差，C1 和C2 是参数。MSSIM的取值范围为0~1, 值越大说明去噪效果越好。
我们对自然图像进行了一系列的仿真实验，运用Matlab7.0编程，实验运行在1.7GHz英特尔酷睿i5-3317U CPU，4GB内存的联想笔记本上。本文算法与自适应中值滤波(AM)、基于边界判定噪声检测的开关中值滤波(BDND)、基于非对称裁剪策略的开关中值滤波(UTM)、改进的方向加权中值滤波(MDWM)、改进的方向加权均值滤波(MDWF)以及自适应加权均值开关滤波(SAWM)进行对比。
AM中表示最大窗口限制的参数wmax =39。MDWM和MDWF中迭代次数设置为10，为避免滤波窗口大小的不同给算法性能带来的影响，将3×3、5×5、7×7、9×9、11×11、13×13、15×15滤波窗口下最佳结果(PSNR最大)作为两种算法每次实验的最终结果。AM、BDND自适应地选择滤波窗口大小。UTM 采用固定的3×3的滤波窗口。SAWM中参数设置为Ld =7、r=1、T 1=2。本专利算法噪声密度估计阶段的参数w=9，第一次噪声检测阶段的参数Ld =7、r=1，第二次噪声检测阶段的参数m=7。
第一组实验由两幅常见的分辨率为256×256的图像构成。其中，Lena作为自身不包含灰度为0和255的非噪声极值点图像的代表，Rice作为自身包含灰度为0或255的非噪声极值点图像的代表。不含噪声的Lena图像中的最大、最小灰度分别为245和23。不含噪声的rice图像中的最大、最小灰度分别为254和0，其中灰度为0的像素点个数为3839。
两幅图分别加入密度范围为10%~90%，增量为10%的9种级别的噪声，表1和2分别给出了第一组实验下各椒盐噪声滤波算法去噪结果的PSNR和MSSIM度量结果。从表中数据可以观察到，对Lena图像，当噪声密度低于30%时，本发明算法的PSNR值最高。对于Rice图像，本发明算法在噪声密度高于30%时获得了最高的PSNR值。其他情况下，本发明算法的PSNR值略低于MDWF或者SAWM，但高于其他方法。MSSIM度量结果表明，本发明算法对噪声密度低于30%的Lena图像和所有噪声密度下的Rice图像获得了最高的MSSIM值。对于Lena图像，本发明算法的去噪效果与SAWM相当，两者的效果好于其它算法。与SAWM比较而言，本发明算法对Rice噪声图像获得的PSNR和MSSIM值有了明显的提高，特别是当噪声密度高于40%时，提高幅度更为显著。高的PSNR值意味着更佳的去噪效果。高的MSSIM值意味着经过去噪后的图像与参考的无噪声图像间保持了更大的结构相似性，去噪效果更好。
在定量比较的基础上，还提供了直观的定性比较。为了节省空间，图2为Lena图像和Rice图像的原图，图3至图8仅仅分别给出了低密度(30%)、中密度(60%)和高密度噪声(90%)下各种算法（从左到右，从上到下依次为噪声图、AM算法结果图、BDND算法结果图、UTM算法结果图、MDWM算法结果图、MDWF算法结果图、SAWM算法结果图及本发明算法结果图）的去噪结果。从图3至图8可以看到，AM和BDND的去噪效果最差。当噪声密度为60%和90%时，AM的去噪结果模糊现象严重，BDND的去噪结果有不少的白色噪声区域。UTM的去噪结果包含一定数量的白色和黑色的噪声区域。MDWM的去噪结果包含一些白色和黑色的噪声线条。本发明算法、MDWF和SAWM比其余算法拥有更佳的去噪效果。对Lena图像的去噪，本发明算法与SAWM效果相当，两者效果略好于MDWF。对噪声密度为30%和60%的Rice图像，本发明算法和MDWF的去噪效果明显好于SAWM，SAWM的去噪结果中含有不少的黑色噪声区域。对噪声密度90%的Rice图像，本发明算法的去噪效果明显好于SAWM和MDWF，因为我们的去噪结果包含了更完整和光滑的米粒边界。MDWF的去噪结果有少许的白色和黑色噪声区域。
第二组实验测试图像源至著名的伯克利图像数据库BSDS300中的300幅彩色图像。因为本发明算法仅对灰度图像进行去噪，故我们在保持每一幅图像分辨率不变的情况下运用Matlab的函数rgb2gray将它们统一转化为灰度图。每一幅灰度图像分别加入密度为10%到90%增量为10%的椒盐噪声。为了节省空间，表3和表4分别仅列出了各种算法在不同噪声密度下300幅图像以及138幅含有非噪声极值点块图像的平均PSNR和MSSIM值。
从表3可以看到，对300幅图像而言，每种噪声密度下，本发明算法均拥有最高的平均PSNR值，表明去噪效果最佳。定量测试结果的统计分析表明，与MDWF和SAWM相比较，本发明算法分别将9种噪声密度下2700幅噪声图像的平均PSNR值分别提高了0.5745556和0.354，最大提高幅度达到0.989和0.927。对138幅含有非噪声极值点块的图像而言，除开噪声密度90%时，本发明算法的平均PSNR值屈居第二，略低于MDWF，其余情况下均保持第一。定量测试结果的统计分析表明，与MDWF和SAWM相比较，本发明算法分别将9种噪声密度下1242幅噪声图像的平均PSNR值分别提高了0.4983333和0.7338889，最大提高幅度达到0.928和1.525。
从表4可以看到，本发明算法对300幅和138幅图像，所有噪声密度下均获得了最高的平均MSSIM值。对300幅图像定量测试结果的统计分析表明，与MDWF和SAWM相比较，本发明算法分别将9种噪声密度下2700幅噪声图像的平均MSSIM值分别提高了0.0312156和0.0016522，最大提高幅度达到0.05374和0.00228。对138幅含有非噪声极值点块的图像定量测试结果的统计分析表明，与MDWF和SAWM相比较，本发明算法分别将9种噪声密度下1242幅噪声图像的平均MSSIM值分别提高了0.0316756和0.0035944，最大提高幅度达到0.05487和0.00534。统计分析结果也表明，本发明算法在含非噪声极值点块的图像上获得了更大的去噪性能改进。
此外，实验也表明，本发明算法运行速度明显快于MDWF，后者的运行时间是前者的5~6倍。本发明算法运行速度慢于SAWM，运行时间大致是SAWM的2~3倍。MDWF速度慢的原因在于：1) 它是一个迭代算法，每次运行都需要进行10次的算法迭代。2) 为避免滤波窗口大小的不同给算法性能带来的影响，每次实验都对3×3、5×5、7×7、9×9、11×11、13×13、15×15总共9种滤波窗口下的MDWF进行了测试，取最佳结果(PSNR最大)作为此次实验的最终结果。
以上是本发明的较佳实施例，凡依本发明技术方案所作的改变，所产生的功能作用未超出本发明技术方案的范围时，均属于本发明的保护范围。
Claims (7)
1.一种融合局部与全局信息的图像椒盐噪声滤波方法，其特征在于：包括如下步骤，
步骤S1：读入噪声图像，进行噪声密度估计；
步骤S2：进行第一次噪声检测；
步骤S3：对所述步骤S2的噪声检测结果进行校正，并进行第一次噪声恢复；
步骤S4：进行第二次噪声检测；
步骤S5：进行第二次噪声恢复，得到去噪图像。
2.根据权利要求1所述的融合局部与全局信息的图像椒盐噪声滤波方法，其特征在于：所述步骤S1中，所述噪声密度估计的具体过程如下：
步骤S11：将噪声图像的行、列分别分成w等份，计算每个等份对应的行、列数，若最后一个等份的行、列数不够，采用镜面反射的方式进行补足，最终将图像分成w×w个子块；
步骤S12：统计每个子块中灰度为0和255的像素点出现的概率，构造向量
；
步骤S13：对向量中的元素按升序排序，排序后的向量记为；
步骤S14：估计整幅图像的噪声密度为：
相应地，整幅图像噪声点个数的估计值为：；
其中，mod代表取余，N表示图像中像素点的总数。
3.根据权利要求2所述的融合局部与全局信息的图像椒盐噪声滤波方法，其特征在于：所述步骤S2中，所述第一次噪声检测的具体过程如下：
步骤S21：采用算法SAWM噪声检测策略的改进版本SAWM_MD对噪声图像进行第一次噪声检测，即先计算出每个像素点pi,
j 对应的方向加权绝对灰度差异之和的最小值di ,j ，用所有像素点的di ,j 构建矩阵Md ；
步骤S22：对矩阵Md 中的元素按降序排列，构成新的向量，向量中的一个元素对应图像中的一个像素点；
步骤S23：取向量中前ne 个元素对应的像素点为噪声点，构成噪声点集合P，产生噪声检测结果标示矩阵Bm ：
其中，标示矩阵Bm 是一个二值化矩阵，1和0分别表示检测到的噪声点和非噪声点。
4.根据权利要求3所述的融合局部与全局信息的图像椒盐噪声滤波方法，其特征在于：所述步骤S3中，所述对噪声检测结果进行校正的具体过程如下：
步骤S31：根据椒盐噪声点灰度只可能为0或者255的事实，剔除掉噪声检测结果中灰度不为0或者255的噪声点，其形式化描述为：
步骤S32：若噪声密度估计结果De 最接近于9种噪声级别的前5种，则执行步骤S33~S36；否则，校正过程结束；
步骤S33：分别获取灰度0和255对应的阈值分割结果：
其中，1和0分别代表阈值分割结果中的目标像素和背景像素；
步骤S34：搜索8连通拓扑结构下，阈值分割结果seg 0和seg 255中目标像素构成的连通分量，将所有连通分量组成一个集合C=，其中，Ci 代表第i个连通分量对应像素点的集合，m代表连通分量的个数；
步骤S35：按如下条件移除掉可能由灰度为0或者255的噪声点构成的连通分量，得到新的连通分量集合：
其中，ni 代表连通分量Ci 包含的像素点个数；
参数：；其中，round代表四舍五入；
步骤S36：根据连通分量对Bm 进行如下校正：
。
5.根据权利要求4所述的融合局部与全局信息的图像椒盐噪声滤波方法，其特征在于：所述9种噪声级别分别为10%、20%、30%、40%、50%、60%、70%、80%、90%。
6.根据权利要求4所述的融合局部与全局信息的图像椒盐噪声滤波方法，其特征在于：所述步骤S4中，所述第二次噪声检测的具体过程如下：
步骤S41：采用算法MDWF的噪声检测策略MDWF_D对经过第一次噪声恢复的图像进行再次的噪声检测，判定椒盐噪声点的条件为：
其中，P表示检测到的椒盐噪声点集合，T为一参数；
步骤S42：对T值进行如下设置：
。
7.根据权利要求1至6任意一项所述的融合局部与全局信息的图像椒盐噪声滤波方法，其特征在于：所述第一次噪声恢复和第二次噪声恢复均采用算法SAWM的噪声点恢复策略是改进版本SAWM_MR进行噪声恢复，其具体过程如下：
首先，以噪声点为中心，初始化滤波窗口大小w 1=3，迭代地扩充滤波窗口大小为w 1= w 1+2，直至满足滤波窗口内非噪声点的个数大于等于2或者w 1≥21；
然后，取噪声点为中心的w 1×w 1滤波窗口Ω，剔除掉其中的噪声点，将剩余非噪声点构成集合U；如果U为空集合，则将其重置为U=Ω；将集合U中像素点的加权灰度均值作为中心噪声点的滤波输出，其形式化描述为：
其中，fs
,t 代表像素点ps,
t 的灰度，ω s, t 代表ps, t 的加权系数，加权函数的定义如下：
其中，Nu 代表集合U中像素点的个数，pm,n 和pk,l 代表集合U中的像素点，fm ,n 和fk ,l 分别代表像素点pm,n 和pk,l 的灰度。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410489377.8A CN104217405B (zh) | 2014-09-23 | 融合局部与全局信息的图像椒盐噪声滤波方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410489377.8A CN104217405B (zh) | 2014-09-23 | 融合局部与全局信息的图像椒盐噪声滤波方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104217405A true CN104217405A (zh) | 2014-12-17 |
| CN104217405B CN104217405B (zh) | 2017-01-04 |
Family
ID=
Cited By (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105243649A (zh) * | 2015-11-09 | 2016-01-13 | 天津大学 | 基于二次噪点检测的图像降噪方法 |
| CN105989578A (zh) * | 2015-02-26 | 2016-10-05 | 展讯通信（上海）有限公司 | 一种图像降噪方法和装置 |
| CN106920222A (zh) * | 2017-03-13 | 2017-07-04 | 苏州大学 | 一种图像平滑方法及装置 |
| CN108269243A (zh) * | 2018-01-18 | 2018-07-10 | 福州鑫图光电有限公司 | 一种图像信噪比的增强方法及终端 |
| CN108470344A (zh) * | 2017-02-23 | 2018-08-31 | 南宁市富久信息技术有限公司 | 自适应阈值的Sobel边缘检测方法 |
| CN109309862A (zh) * | 2018-07-26 | 2019-02-05 | 浠诲嘲 | 多媒体数据编辑系统 |
| CN109767389A (zh) * | 2019-01-15 | 2019-05-17 | 四川大学 | 基于局部和非局部联合先验的自适应加权双范数遥感图像盲超分辨重建方法 |
| CN110458783A (zh) * | 2019-08-22 | 2019-11-15 | 齐鲁工业大学 | 一种两步去除具有图像细节保护能力的非局部双边滤波图像椒盐噪声的方法 |
| CN111476537A (zh) * | 2020-03-31 | 2020-07-31 | 广州高新工程顾问有限公司 | 一种基于bim的工程造价动态控制系统及方法 |
| CN111476723A (zh) * | 2020-03-17 | 2020-07-31 | 哈尔滨师范大学 | 一种Landsat-7扫描线纠正器失效的遥感图像丢失像素恢复方法 |
| CN113591563A (zh) * | 2021-06-24 | 2021-11-02 | 金陵科技学院 | 图像固定值脉冲噪声去噪方法及其模型训练方法 |
| CN113989168A (zh) * | 2021-11-02 | 2022-01-28 | 华北电力大学（保定） | 一种针对椒盐噪声的自适应非局部均值滤波方法 |
| CN114051127A (zh) * | 2022-01-14 | 2022-02-15 | 深圳市艾科维达科技有限公司 | 网络机顶盒的图像传输降噪方法 |
| CN116630447A (zh) * | 2023-07-24 | 2023-08-22 | 成都海风锐智科技有限责任公司 | 一种基于图像处理的天气预测方法 |
| CN117437600A (zh) * | 2023-12-20 | 2024-01-23 | 山东海纳智能装备科技股份有限公司 | 基于图像识别技术的煤流监测系统 |
Non-Patent Citations (5)
| Title |
|---|
| WEIXIA LIU等: "Survey over Switching Filters for Salt & Pepper Noise Removal", 《2013 FIFTH INTERNATIONAL CONFERENCE ON INTELLIGENT HUMAN-MACHINE SYSTEMS AND CYBERNETICS》 * |
| ZUOYONG LI等: "Image Thresholding-based Switching Filter for Salt & Pepper Noise Removal", 《2013 SECOND INTERNATIONAL CONFERENCE ON ROBOT, VISION AND SIGNAL PROCESSING》 * |
| 张淑英等: "图像椒盐噪声的开关滤波算法综述", 《福建电脑》 * |
| 李佐勇等: "椒盐图像的方向加权均值滤波算法", 《中国图象图形学报》 * |
| 林亚明等: "分块策略实现图像椒盐噪声密度估计", 《中国图象图形学报》 * |
Cited By (24)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105989578A (zh) * | 2015-02-26 | 2016-10-05 | 展讯通信（上海）有限公司 | 一种图像降噪方法和装置 |
| CN105989578B (zh) * | 2015-02-26 | 2019-01-18 | 展讯通信（上海）有限公司 | 一种图像降噪方法和装置 |
| CN105243649A (zh) * | 2015-11-09 | 2016-01-13 | 天津大学 | 基于二次噪点检测的图像降噪方法 |
| CN108470344A (zh) * | 2017-02-23 | 2018-08-31 | 南宁市富久信息技术有限公司 | 自适应阈值的Sobel边缘检测方法 |
| CN106920222B (zh) * | 2017-03-13 | 2019-11-08 | 苏州大学 | 一种图像平滑方法及装置 |
| CN106920222A (zh) * | 2017-03-13 | 2017-07-04 | 苏州大学 | 一种图像平滑方法及装置 |
| CN108269243A (zh) * | 2018-01-18 | 2018-07-10 | 福州鑫图光电有限公司 | 一种图像信噪比的增强方法及终端 |
| CN109309862A (zh) * | 2018-07-26 | 2019-02-05 | 浠诲嘲 | 多媒体数据编辑系统 |
| CN109767389B (zh) * | 2019-01-15 | 2023-06-20 | 四川大学 | 基于局部和非局部联合先验的自适应加权双范数遥感图像盲超分辨重建方法 |
| CN109767389A (zh) * | 2019-01-15 | 2019-05-17 | 四川大学 | 基于局部和非局部联合先验的自适应加权双范数遥感图像盲超分辨重建方法 |
| CN110458783B (zh) * | 2019-08-22 | 2023-01-31 | 齐鲁工业大学 | 一种基于具有图像细节保护能力的非局部双边滤波的椒盐噪声去除方法 |
| CN110458783A (zh) * | 2019-08-22 | 2019-11-15 | 齐鲁工业大学 | 一种两步去除具有图像细节保护能力的非局部双边滤波图像椒盐噪声的方法 |
| CN111476723A (zh) * | 2020-03-17 | 2020-07-31 | 哈尔滨师范大学 | 一种Landsat-7扫描线纠正器失效的遥感图像丢失像素恢复方法 |
| CN111476537B (zh) * | 2020-03-31 | 2021-01-15 | 广州高新工程顾问有限公司 | 一种基于bim的工程造价动态控制系统及方法 |
| CN111476537A (zh) * | 2020-03-31 | 2020-07-31 | 广州高新工程顾问有限公司 | 一种基于bim的工程造价动态控制系统及方法 |
| CN113591563A (zh) * | 2021-06-24 | 2021-11-02 | 金陵科技学院 | 图像固定值脉冲噪声去噪方法及其模型训练方法 |
| CN113591563B (zh) * | 2021-06-24 | 2023-06-06 | 金陵科技学院 | 图像固定值脉冲噪声去噪方法及其模型训练方法 |
| CN113989168A (zh) * | 2021-11-02 | 2022-01-28 | 华北电力大学（保定） | 一种针对椒盐噪声的自适应非局部均值滤波方法 |
| CN113989168B (zh) * | 2021-11-02 | 2024-05-14 | 华北电力大学（保定） | 一种针对椒盐噪声的自适应非局部均值滤波方法 |
| CN114051127A (zh) * | 2022-01-14 | 2022-02-15 | 深圳市艾科维达科技有限公司 | 网络机顶盒的图像传输降噪方法 |
| CN116630447A (zh) * | 2023-07-24 | 2023-08-22 | 成都海风锐智科技有限责任公司 | 一种基于图像处理的天气预测方法 |
| CN116630447B (zh) * | 2023-07-24 | 2023-10-20 | 成都海风锐智科技有限责任公司 | 一种基于图像处理的天气预测方法 |
| CN117437600A (zh) * | 2023-12-20 | 2024-01-23 | 山东海纳智能装备科技股份有限公司 | 基于图像识别技术的煤流监测系统 |
| CN117437600B (zh) * | 2023-12-20 | 2024-03-26 | 山东海纳智能装备科技股份有限公司 | 基于图像识别技术的煤流监测系统 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Hendrycks et al. | Benchmarking neural network robustness to common corruptions and perturbations | |
| Zhou | Cognition and removal of impulse noise with uncertainty | |
| Jafar et al. | Efficient improvements on the BDND filtering algorithm for the removal of high-density impulse noise | |
| CN108765465B (zh) | 一种无监督sar图像变化检测方法 | |
| CN109741356B (zh) | 一种亚像素边缘检测方法及系统 | |
| CN103400357B (zh) | 一种去除图像椒盐噪声的方法 | |
| Hosseini et al. | Fast restoration of natural images corrupted by high-density impulse noise | |
| CN109978848B (zh) | 基于多光源颜色恒常模型检测眼底图像中硬性渗出的方法 | |
| CN109478329B (zh) | 图像处理方法和装置 | |
| WO2019026104A1 (ja) | 情報処理装置、情報処理プログラム及び情報処理方法 | |
| US6731806B1 (en) | Region growing based noise reduction method for digital images | |
| CN113822328A (zh) | 防御对抗样本攻击的图像分类方法、终端设备及存储介质 | |
| CN109102013A (zh) | 一种适于隧道环境特性的改进freak特征点匹配稳像方法 | |
| CN112329677B (zh) | 基于特征融合的遥感图像河道目标检测方法和装置 | |
| Chou et al. | A noise-ranking switching filter for images with general fixed-value impulse noises | |
| Flyckt et al. | Detecting ditches using supervised learning on high-resolution digital elevation models | |
| CN104217405A (zh) | 融合局部与全局信息的图像椒盐噪声滤波方法 | |
| CN104217405B (zh) | 融合局部与全局信息的图像椒盐噪声滤波方法 | |
| CN113963178A (zh) | 地空背景下红外弱小目标检测方法、装置、设备及介质 | |
| Gandikota et al. | Pixel Noise Localization Algorithm for Indian Satellite Data Quality Control: A Novel Approach | |
| JP4758272B2 (ja) | 移動物体追跡装置及びプログラム | |
| Huang et al. | Single image dehazing via color balancing and quad-decomposition atmospheric light estimation | |
| Kavitha et al. | Implementation of gray-level clustering algorithm for image segmentation | |
| CN117475157B (zh) | 基于无人机遥感的农业种植增强监测方法 | |
| Garg | An Improved Switching Median Filter Based on Local Outlier Factor |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20191217 Address after: 350000 Gulou Science and Technology Business Center Building, North Wushan Road (No. 451 Industrial Road), Gulou District, Fuzhou City, Fujian Province, Ninth Floor Patentee after: FUJIAN HUICHUAN INTERNET OF THINGS TECHNOLOGY SCIENCE AND TECHNOLOGY CO.,LTD. Address before: Minhou County of Fuzhou City, Fujian province 350108 Street Town Road No. 1 Patentee before: Minjiang University |
|
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20240119 Address after: 3/F, Building 2A, Fuzhou Internet of Things Industrial Park, Liandong U Valley, No. 12 Xingye West Road, Mawei District, Fuzhou City, Fujian Province, 350015 (within the Free Trade Zone) Patentee after: Huiyun data application (Fuzhou) Co.,Ltd. Address before: 350000 9th floor, Gulou science and technology business center building, north side of Wushan West Road, Gulou District, Fuzhou City, Fujian Province Patentee before: FUJIAN HUICHUAN INTERNET OF THINGS TECHNOLOGY SCIENCE AND TECHNOLOGY CO.,LTD. |