CN107392864A - 一种去除高斯噪声和脉冲噪声的混合噪声滤波方法 - Google Patents

一种去除高斯噪声和脉冲噪声的混合噪声滤波方法 Download PDF

Info

Publication number
CN107392864A
CN107392864A CN201710528489.3A CN201710528489A CN107392864A CN 107392864 A CN107392864 A CN 107392864A CN 201710528489 A CN201710528489 A CN 201710528489A CN 107392864 A CN107392864 A CN 107392864A
Authority
CN
China
Prior art keywords
noise
mrow
pixel
point
value
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.)
Pending
Application number
CN201710528489.3A
Other languages
English (en)
Inventor
陈壮
胡雪冰
白良留
张健雨
张永
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University of Science and Technology
Original Assignee
Nanjing University of Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nanjing University of Science and Technology filed Critical Nanjing University of Science and Technology
Priority to CN201710528489.3A priority Critical patent/CN107392864A/zh
Publication of CN107392864A publication Critical patent/CN107392864A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/70Denoising; Smoothing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20024Filtering details

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Image Processing (AREA)

Abstract

本发明公开了一种去除高斯噪声和脉冲噪声的混合噪声滤波方法。该方法步骤如下:对于受电弓图像中高斯噪声和脉冲噪声的混合噪声,采用先进边界区分噪声检测算法即ABDND算法对脉冲噪声进行抑制,所述ABDND算法包含噪声像素检测、噪声像素核实、噪声像素滤除;针对对脉冲噪声进行抑制后的图像,采用滑动窗口对像素点进行处理,将窗口内所有点根据灰度值大小进行排序,去除极大值和极小值点,对剩余的像素点进行权值分配并进行加权计算,计算结果作为滤波窗口中心点的输出灰度值,代替该点原有的灰度值;重复上述步骤中两种算法进行递进滤波,最终滤除受电弓图像中的混合噪声。本发明原理简单、滤波效果好,对数字图像处理具有重要意义。

Description

一种去除高斯噪声和脉冲噪声的混合噪声滤波方法
技术领域
本发明属于数字图像滤波技术领域,特别是一种去除高斯噪声和脉冲噪声的混合噪声滤波方法。
背景技术
图像在产生和传输过程中常常会受到各种噪声的污染,致使图像的质量下降。在应用中,有两种噪声可以充分代表大部分图像中的噪音,即高斯噪声和脉冲噪声,这两种噪声和混合噪声在实际中也是难以避免的,对这种混合噪声的去除是图像处理中重要而又富有挑战性的课题。
随着图像处理技术的不断发展,混合噪声滤波问题一直是国内外学者关注的问题。方先成等针对高斯-脉冲混合噪声滤波问题提出了基于BM3D的两种不同去噪方法,当高斯噪声的方差和脉冲噪声浓度变化时,两种方法呈现出不同的优势。吴德刚先根据脉冲噪声的特点和局部能量信息提取出脉冲噪声,采用中值滤波加以去除,然后采用均值滤波对高斯噪声进行去噪。肖孟强等针对医学图像中的混合噪声滤除问题,综合考虑中值滤波算法和均值滤波算法的各自优点,提出了一种针对混合噪声的滤波算法,更完整的保存了图像细节信息。Liu等提出一个加权范数能量最小化模型以去除图像中的混合噪声,并通过实验结果验证该方法的可行性。Camarena等针对彩色图像引入一种基于模糊加权向量中值滤波器(FWMF)模型来检测和去除噪声,试验结果表明该方法能够有效地去除噪声并保留边缘信息。上述方法都能在一定程度上滤除图像中的混合噪声,但原理复杂,算法复杂度较高,图像处理速度较慢。
发明内容
本发明的目的在于提供一种简单方便、效果显著的去除高斯噪声和脉冲噪声的混合噪声滤波方法,针对高斯噪声和脉冲噪声,实现数字图像的混和噪声去噪。
实现本发明目的的技术解决方案是:一种去除高斯噪声和脉冲噪声的混合噪声滤波方法,步骤如下:
步骤1,对于受电弓图像中高斯噪声和脉冲噪声的混合噪声,采用先进边界区分噪声检测算法即ABDND算法对脉冲噪声进行抑制;
步骤2,针对步骤1中对脉冲噪声进行抑制后的图像,采用均值滤波算法对高斯噪声进行抑制;
步骤3,重复步骤1~2中两种算法进行递进滤波,最终滤除受电弓图像中的混合噪声。
进一步地,步骤1所述ABDND算法为基于自适应开关中值滤波算法ASMF的改进算法,包含噪声像素检测、噪声像素核实、噪声像素滤除三个阶段,具体步骤如下:
(1.1)噪声像素检测阶段:通过直方图处理得到噪声像素的范围,并以此为依据初步判断噪声像素和非噪声像素;
(1.2)噪声像素核实阶段:对临时噪声像素点附近的像素点进行局部信息统计,判断该像素点是否噪声像素点;
(1.3)噪声像素滤除阶段:采用ASMF算法滤除噪声像素,实现图像的脉冲噪声滤除。
进一步地,步骤(1.1)所述噪声像素检测阶段:通过直方图处理得到噪声像素的范围,并以此为依据初步判断噪声像素和非噪声像素,具体如下:
1)获取图像的全局灰度直方图,对相邻灰度值的直方图数值做差值处理,得到差值序列ΔH:
ΔHi=Hi+1-Hi,i=MIN,MIN+1,…,MAX-1 (1)
其中,Hi是灰度值i对应的直方图数值,ΔHi是灰度值i的前向差分值,MIN、MAX分别为灰度值的最小值、最大值;
2)根据差值序列ΔH,得到差值中的正最大值和负最大值,根据式(2)得到噪声像素的两个边界T1和T2
其中,Index2、Index1分别是正、负最大值在差值序列ΔH中的索引值;
3)采用大小为W×W的滤波窗口对图像像素点进行遍历,滤波窗口的中心点与待检测像素点重合,并统计该窗口内像素点的最大灰度值smax和最小灰度值smin,根据式(3)得到滤波图像的噪声标识矩阵r(i,j),初步判断该像素点是否是噪声像素:
其中,si,j是图像中点(i,j)的像素灰度值,r(i,j)=1则该点为临时噪声像素点,r(i,j)=0则该点不是噪声像素点。
进一步地,步骤(1.2)所述噪声像素核实阶段:对临时噪声像素点附近的像素点进行局部信息统计,判断该像素点是否噪声像素点,具体如下:
对于临时噪声像素做以下算法处理:
1)以临时噪声像素点为中心,取尺寸大小为w×w的窗口,其中w小于W,统计该窗口内非噪声像素点的数量nf,如果nf>w,计算窗口内非噪声像素点与中心临时噪声像素点的像素差值,所有像素差值之和值记为d,然后执行步骤3);如果nf≤w,执行步骤2);
2)令w=w+2,执行步骤1),滤波窗口的长度逐步递增2,且滤波窗口的长度的最大值小于等于W;
3)结合图像的最大像素灰度值MAX、最小像素灰度值MIN以及点(i,j)的灰度值si,j,确定阈值T3;若d小于T3,则该点不是噪声像素点,同时在噪声标识矩阵r(i,j)中将该点的值修改为0;反之,则该点是噪声像素点,矩阵r(i,j)该点的值保持不变;
T3=min(si,j-MIN,MAX-si,j) (4)
4)将w×w滤波窗口移至下一个临时噪声像素点,按照步骤1)~3)判断下一个临时噪声像素点是否是噪声像素;对所有临时噪声像素点逐一判断,排除非噪声像素点,得到最终的噪声像素点。
进一步地,步骤(1.3)所述噪声像素滤除阶段:采用ASMF算法滤除噪声像素,实现图像的脉冲噪声滤除,即利用滤波窗口中的非噪声像素点进行中值滤波,滤除噪声像素点,具体过程如下:
对图像任一噪声像素点(m,n),以该点为中心、采用大小为L×L的滤波窗口,将该窗口内所有非噪声像素点组成集合Em,n
对Em,n中灰度值进行升序排列,得到集合Fm,n,并提取集合Fm,n的排序中值Mm,n
在集合Fm,n中选取Q个最接近中值Mm,n的灰度值,和中值Mm,n组成新的集合,并计算该集合的灰度均值Qm,n,用Qm,n取代噪声点(m,n)的灰度值Fm,n;根据滤波窗口的大小对Q值进行选择,具体选择方法如式(6)所示:
进一步地,步骤2所述采用均值滤波算法对高斯噪声进行抑制,具体步骤如下:
(2.1)输入经过脉冲噪声抑制后的图像;
(2.2)以待滤波像素点为中心,采用滑动窗口W对像素点进行处理,将窗口内所有点根据灰度值大小进行排序,去除极大值和极小值点,将剩余点记为集合H;
(2.3)对集合H内的点进行权值分配并进行归一化处理,求出集合H内像素的均值Mean(∑H(i,j)),然后采用式(7)进行各像素权值的计算,并进行归一化处理;
其中,wk为集合H内各点的权值,k表示H内第k个点,N为H内点的总数,Dk为集合H内各点的灰度值与灰度均值Mean(∑H(i,j))的差的绝对值;
如式(8)所示,T为集合H内所有Dk的均值,作为各点权值选择的阈值,如式(9)所示,当集合H内一点的Dk值大于阈值T时,该点的权值由Dk决定,当其Dk小于阈值T时,点的权值由T决定;
Dk=|Hk-Mean(∑H(i,j))| (8)
(2.4)将集合H内的所有点与相对应的权值相乘并求和,和值作为滤波窗口中心点的输出灰度值,代替该点原有的灰度值,如式(10)所示:
(2.5)重复步骤(2.2)~(2.4)对图像的所有点进行遍历,实现图像的滤波处理。
本发明与现有技术相比,其显著优点在于:(1)滤波效果较好,能够有效滤除图像中高斯噪声和脉冲噪声;(2)能够对数字图像的混和噪声去噪,且算法简单,复杂度低。
附图说明
图1是本发明去除高斯噪声和脉冲噪声的混合噪声滤波方法的流程图。
图2是本发明中局部去极值加权均值滤波算法流程图。
图3是受电弓左半弓的示意图。
图4是不同算法对于高斯和脉冲混合噪声的滤除效果图,其中(a)为混合噪声的滤除效果图,(b)为均值滤波算法的滤除效果图,(c)为中值滤波算法滤除效果图,(d)为本发明混合滤波算法的滤除效果图。
具体实施方式
结合图1,本发明去除高斯噪声和脉冲噪声的混合噪声滤波方法,步骤如下:
步骤1,对于受电弓图像中高斯噪声和脉冲噪声的混合噪声,采用先进边界区分噪声检测算法(ABDND)对脉冲噪声进行抑制;所述ABDND算法为基于自适应开关中值滤波算法ASMF的改进算法,包含噪声像素检测、噪声像素核实、噪声像素滤除三个阶段,具体步骤如下:
(1.1)噪声像素检测阶段:通过直方图处理得到噪声像素的范围,并以此为依据初步判断噪声像素和非噪声像素,具体如下:
1)获取图像的全局灰度直方图,对相邻灰度值的直方图数值做差值处理,得到差值序列ΔH:
ΔHi=Hi+1-Hi,i=MIN,MIN+1,…,MAX-1 (1)
其中,Hi是灰度值i对应的直方图数值,ΔHi是灰度值i的前向差分值,MIN、MAX分别为灰度值的最小值、最大值;
2)根据差值序列ΔH,得到差值中的正最大值和负最大值,根据式(2)得到噪声像素的两个边界T1和T2
其中,Index2、Index1分别是正、负最大值在差值序列ΔH中的索引值;
3)采用大小为W×W的滤波窗口对图像像素点进行遍历,滤波窗口的中心点与待检测像素点重合,并统计该窗口内像素点的最大灰度值smax和最小灰度值smin,根据式(3)得到滤波图像的噪声标识矩阵r(i,j),初步判断该像素点是否是噪声像素:
其中,si,j是图像中点(i,j)的像素灰度值,r(i,j)=1则该点为临时噪声像素点,并通称为临时噪声像素点,r(i,j)=0则该点不是噪声像素点。
(1.2)噪声像素核实阶段:对临时噪声像素点附近的像素点进行局部信息统计,判断该像素点是否噪声像素点,具体如下:
对于临时噪声像素做以下算法处理:
1)以临时噪声像素点为中心,取尺寸大小为w×w的窗口,其中w远小于W,统计该窗口内非噪声像素点的数量nf,如果nf>w,计算窗口内非噪声像素点与中心临时噪声像素点的像素差值,所有像素差值之和值记为d,然后执行步骤3);如果nf≤w,执行步骤2);
2)令w=w+2,执行步骤1),滤波窗口的长度逐步递增2,且滤波窗口的长度的最大值小于等于W;
3)结合图像的最大像素灰度值MAX、最小像素灰度值MIN以及点(i,j)的灰度值si,j,确定阈值T3;若d小于T3,则该点不是噪声像素点,同时在噪声标识矩阵r(i,j)中将该点的值修改为0;反之,则该点是噪声像素点,矩阵r(i,j)该点的值保持不变;
T3=min(si,j-MIN,MAX-si,j) (4)
4)将w×w滤波窗口移至下一个临时噪声像素点,按照步骤1)~3)判断下一个临时噪声像素点是否是噪声像素;对所有临时噪声像素点逐一判断,排除非噪声像素点,得到最终的噪声像素点。
(1.3)噪声像素滤除阶段:采用ASMF算法滤除噪声像素,实现图像的脉冲噪声滤除采用ASMF算法滤除噪声像素,实现图像的脉冲噪声滤除,即利用滤波窗口中的非噪声像素点进行中值滤波,滤除噪声像素点;
经过阶段一和阶段二的处理后,整个图像的像素点分为噪声像素点和非噪声像素点(信号像素)两类。对于非噪声像素点,无需滤波,保持其灰度值不变即可;对于噪声像素点,需利用滤波窗口中的非噪声像素点进行中值滤波,具体过程如下:
对图像任一噪声像素点(m,n),以该点为中心、采用大小为L×L的滤波窗口,将该窗口内所有非噪声像素点组成集合Em,n
对Em,n中灰度值进行升序排列,得到集合Fm,n,并提取集合Fm,n的排序中值Mm,n;集合Em,n将噪声像素点排除在外,只包含非噪声像素点,这样就降低了噪声对滤波结果的影响。
为进一步验证阶段三的噪声滤除性能,在集合Fm,n中选取Q个最接近中值Mm,n的灰度值,和中值Mm,n组成新的集合,并计算该集合的灰度均值Qm,n,用Qm,n取代噪声点(m,n)的灰度值Fm,n;根据滤波窗口的大小对Q值进行选择,具体选择方法如式(6)所示:
步骤2,针对步骤1中对脉冲噪声进行抑制后的图像,采用均值滤波算法对高斯噪声进行抑制,具体步骤如下:
(2.1)输入经过脉冲噪声抑制后的图像;
(2.2)以待滤波像素点为中心,采用滑动窗口W对像素点进行处理,将窗口内所有点根据灰度值大小进行排序,去除极大值和极小值点,将剩余点记为集合H;
(2.3)对集合H内的点进行权值分配并进行归一化处理,求出集合H内像素的均值Mean(∑H(i,j)),然后采用式(7)进行各像素权值的计算,并进行归一化处理;
其中,wk为集合H内各点的权值,k表示H内第k个点,N为H内点的总数,Dk为集合H内各点的灰度值与灰度均值Mean(∑H(i,j))的差的绝对值;
如式(8)所示,T为集合H内所有Dk的均值,作为各点权值选择的阈值,如式(9)所示,当集合H内一点的Dk值大于阈值T时,该点的权值由Dk决定,当其Dk小于阈值T时,点的权值由T决定;
Dk=|Hk-Mean(∑H(i,j))| (8)
(2.4)将集合H内的所有点与相对应的权值相乘并求和,和值作为滤波窗口中心点的输出灰度值,代替该点原有的灰度值,如式(10)所示:
(2.5)重复步骤(2.2)~(2.4)对图像的所有点进行遍历,实现图像的滤波处理。
步骤3,重复步骤1~2中两种算法进行递进滤波,最终滤除受电弓图像中的混合噪声。
下面结合具体实施例对本发明作进一步详细说明。
实施例1
城轨列车通过受电弓滑板从接触线上接收电流作为列车行驶的动力,因接触线的高电压,受电弓周边存在一定的电磁干扰,所以在基于图像法的受电弓状态检测中,获取到的受电弓图像存在一定的混合噪声,主要是脉冲和高斯噪声。对于混合噪声的去噪问题,单一类型的滤波器去噪效果并不好,一般采用两种类型的滤波器递进地滤除噪声。本专利提出一种去除高斯噪声和脉冲噪声的混合噪声滤波算法,首先采用基于ABDND的中值滤波算法对脉冲噪声进行抑制,再采用一种改进的均值滤波算法对高斯噪声进行抑制,两种算法递进滤除受电弓图像中的混合噪声。评价主要采用均方差误差(MSE)和峰值信噪比(PSNR)来对图像滤波算法的性能进行评价,MSE越小,滤波效果越好,PSNR越大,滤波效果越好。以X和分别表示原始图像和滤波后的图像,图像大小为I×J,均方误差MSE和峰值信噪比PSNR的定义如式(11)和(22)。
X={x(i,j)|1≤i≤I,1≤j≤J} (11)
(1)均方误差
(2)峰值信噪比
其中,MAXI是表示图像点颜色的最大数值。
基于Matlab平台本文分别采用传统均值滤波算法、传统中值滤波算法和本文所提的混合滤波算法对不同密度的高斯和脉冲混合噪声图像进行滤波处理,结果如图4。
图4脉冲噪声与高斯噪声的方差分别为0.02、0.05的混合噪声图像算法处理结果其中图4(a)为混合噪声,图4(b)为均值滤波算法,图4(c)为中值滤波算法,图4(d)为本文混合滤波算法处理后图像。由图4可知,对于混合噪声,从视觉效果来看,本文混合滤波、传统中值、传统均值这三种滤波算法的去噪性能差别不大,为比较三者的优劣,对滤波后的图像进行MSE和PSNR的计算,计算结果如表1所示。
表1
表1是不同算法对于高斯和脉冲混合噪声图像处理后MES和PSNE指标,由表1可知,本文采用混合滤波算法的MSE值比传统中值、均值滤波的MSE值都小,PSNR值比传统中值、均值滤波的PSNR值都大。随着混合噪声密度的增加,三种算法的PSNR值都在减小,MSE值都在增大,表明在高密度噪声环境下,三种算法的滤波性能都在下降,但是本文采用算的法MSE值和PSNR值的变化趋势都比较平稳。综上所述,本文所采用的混合滤波算法的滤波性能优于传统中值、均值滤波算法。

Claims (6)

1.一种去除高斯噪声和脉冲噪声的混合噪声滤波方法,其特征在于,步骤如下:
步骤1,对于受电弓图像中高斯噪声和脉冲噪声的混合噪声,采用先进边界区分噪声检测算法即ABDND算法对脉冲噪声进行抑制;
步骤2,针对步骤1中对脉冲噪声进行抑制后的图像,采用均值滤波算法对高斯噪声进行抑制;
步骤3,重复步骤1~2中两种算法进行递进滤波,最终滤除受电弓图像中的混合噪声。
2.根据权利要求1所述的去除高斯噪声和脉冲噪声的混合噪声滤波方法,其特征在于,步骤1所述ABDND算法为基于自适应开关中值滤波算法ASMF的改进算法,包含噪声像素检测、噪声像素核实、噪声像素滤除三个阶段,具体步骤如下:
(1.1)噪声像素检测阶段:通过直方图处理得到噪声像素的范围,并以此为依据初步判断噪声像素和非噪声像素;
(1.2)噪声像素核实阶段:对临时噪声像素点附近的像素点进行局部信息统计,判断该像素点是否噪声像素点;
(1.3)噪声像素滤除阶段:采用ASMF算法滤除噪声像素,实现图像的脉冲噪声滤除。
3.根据权利要求2所述的去除高斯噪声和脉冲噪声的混合噪声滤波方法,其特征在于,步骤(1.1)所述噪声像素检测阶段:通过直方图处理得到噪声像素的范围,并以此为依据初步判断噪声像素和非噪声像素,具体如下:
1)获取图像的全局灰度直方图,对相邻灰度值的直方图数值做差值处理,得到差值序列ΔH:
ΔHi=Hi+1-Hi,i=MIN,MIN+1,…,MAX-1 (1)
其中,Hi是灰度值i对应的直方图数值,ΔHi是灰度值i的前向差分值,MIN、MAX分别为灰度值的最小值、最大值;
2)根据差值序列ΔH,得到差值中的正最大值和负最大值,根据式(2)得到噪声像素的两个边界T1和T2
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>T</mi> <mn>1</mn> </msub> <mo>=</mo> <mi>I</mi> <mi>n</mi> <mi>d</mi> <mi>e</mi> <mi>x</mi> <mn>1</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>T</mi> <mn>2</mn> </msub> <mo>=</mo> <mi>M</mi> <mi>A</mi> <mi>X</mi> <mo>-</mo> <mi>I</mi> <mi>n</mi> <mi>d</mi> <mi>e</mi> <mi>x</mi> <mn>2</mn> <mo>+</mo> <mn>1</mn> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
其中,Index2、Index1分别是正、负最大值在差值序列ΔH中的索引值;
3)采用大小为W×W的滤波窗口对图像像素点进行遍历,滤波窗口的中心点与待检测像素点重合,并统计该窗口内像素点的最大灰度值smax和最小灰度值smin,根据式(3)得到滤波图像的噪声标识矩阵r(i,j),初步判断该像素点是否是噪声像素:
其中,si,j是图像中点(i,j)的像素灰度值,r(i,j)=1则该点为临时噪声像素点,r(i,j)=0则该点不是噪声像素点。
4.根据权利要求2所述的去除高斯噪声和脉冲噪声的混合噪声滤波方法,其特征在于,步骤(1.2)所述噪声像素核实阶段:对临时噪声像素点附近的像素点进行局部信息统计,判断该像素点是否噪声像素点,具体如下:
对于临时噪声像素做以下算法处理:
1)以临时噪声像素点为中心,取尺寸大小为w×w的窗口,其中w小于W,统计该窗口内非噪声像素点的数量nf,如果nf>w,计算窗口内非噪声像素点与中心临时噪声像素点的像素差值,所有像素差值之和值记为d,然后执行步骤3);如果nf≤w,执行步骤2);
2)令w=w+2,执行步骤1),滤波窗口的长度逐步递增2,且滤波窗口的长度的最大值小于等于W;
3)结合图像的最大像素灰度值MAX、最小像素灰度值MIN以及点(i,j)的灰度值si,j,确定阈值T3;若d小于T3,则该点不是噪声像素点,同时在噪声标识矩阵r(i,j)中将该点的值修改为0;反之,则该点是噪声像素点,矩阵r(i,j)该点的值保持不变;
T3=min(si,j-MIN,MAX-si,j) (4)
4)将w×w滤波窗口移至下一个临时噪声像素点,按照步骤1)~3)判断下一个临时噪声像素点是否是噪声像素;对所有临时噪声像素点逐一判断,排除非噪声像素点,得到最终的噪声像素点。
5.根据权利要求2所述的去除高斯噪声和脉冲噪声的混合噪声滤波方法,其特征在于,步骤(1.3)所述噪声像素滤除阶段:采用ASMF算法滤除噪声像素,实现图像的脉冲噪声滤除,即利用滤波窗口中的非噪声像素点进行中值滤波,滤除噪声像素点,具体过程如下:
对图像任一噪声像素点(m,n),以该点为中心、采用大小为L×L的滤波窗口,将该窗口内所有非噪声像素点组成集合Em,n
<mrow> <msub> <mi>E</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>n</mi> </mrow> </msub> <mo>=</mo> <mo>{</mo> <msub> <mi>s</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>n</mi> </mrow> </msub> <mo>|</mo> <mi>r</mi> <mrow> <mo>(</mo> <mi>m</mi> <mo>+</mo> <mi>p</mi> <mo>,</mo> <mi>n</mi> <mo>+</mo> <mi>q</mi> <mo>)</mo> </mrow> <mo>=</mo> <mn>0</mn> <mo>,</mo> <mi>p</mi> <mo>,</mo> <mi>q</mi> <mo>&amp;Element;</mo> <mo>&amp;lsqb;</mo> <mo>-</mo> <mfrac> <mi>L</mi> <mn>2</mn> </mfrac> <mo>,</mo> <mfrac> <mi>L</mi> <mn>2</mn> </mfrac> <mo>&amp;rsqb;</mo> <mo>}</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>
对Em,n中灰度值进行升序排列,得到集合Fm,n,并提取集合Fm,n的排序中值Mm,n
在集合Fm,n中选取Q个最接近中值Mm,n的灰度值,和中值Mm,n组成新的集合,并计算该集合的灰度均值Qm,n,用Qm,n取代噪声点(m,n)的灰度值Fm,n;根据滤波窗口的大小对Q值进行选择,具体选择方法如式(6)所示:
<mrow> <mi>Q</mi> <mo>=</mo> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mn>2</mn> </mtd> <mtd> <mtable> <mtr> <mtd> <mrow> <mi>L</mi> <mo>=</mo> <mn>3</mn> </mrow> </mtd> <mtd> <mrow> <mi>o</mi> <mi>r</mi> </mrow> </mtd> <mtd> <mrow> <mi>L</mi> <mo>=</mo> <mn>7</mn> </mrow> </mtd> </mtr> </mtable> </mtd> </mtr> <mtr> <mtd> <mn>4</mn> </mtd> <mtd> <mrow> <mi>L</mi> <mo>=</mo> <mn>5</mn> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> <mo>.</mo> </mrow>
6.根据权利要求1或2所述的去除高斯噪声和脉冲噪声的混合噪声滤波方法,其特征在于,步骤2所述采用均值滤波算法对高斯噪声进行抑制,具体步骤如下:
(2.1)输入经过脉冲噪声抑制后的图像;
(2.2)以待滤波像素点为中心,采用滑动窗口W对像素点进行处理,将窗口内所有点根据灰度值大小进行排序,去除极大值和极小值点,将剩余点记为集合H;
(2.3)对集合H内的点进行权值分配并进行归一化处理,求出集合H内像素的均值Mean(∑H(i,j)),然后采用式(7)进行各像素权值的计算,并进行归一化处理;
<mrow> <msub> <mi>w</mi> <mi>k</mi> </msub> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mn>1</mn> <mo>/</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mi>M</mi> <mi>a</mi> <mi>x</mi> <mo>(</mo> <mrow> <msub> <mi>D</mi> <mi>k</mi> </msub> <mo>,</mo> <mi>T</mi> </mrow> <mo>)</mo> <mo>)</mo> </mrow> </mrow> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <mn>1</mn> <mo>/</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mi>M</mi> <mi>a</mi> <mi>x</mi> <mo>(</mo> <mrow> <msub> <mi>D</mi> <mi>k</mi> </msub> <mo>,</mo> <mi>T</mi> </mrow> <mo>)</mo> <mo>)</mo> </mrow> <mi>T</mi> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>
其中,wk为集合H内各点的权值,k表示H内第k个点,N为H内点的总数,Dk为集合H内各点的灰度值与灰度均值Mean(∑H(i,j))的差的绝对值;
如式(8)所示,T为集合H内所有Dk的均值,作为各点权值选择的阈值,如式(9)所示,当集合H内一点的Dk值大于阈值T时,该点的权值由Dk决定,当其Dk小于阈值T时,点的权值由T决定;
Dk=|Hk-Mean(∑H(i,j))| (8)
<mrow> <mi>T</mi> <mo>=</mo> <mfrac> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msub> <mi>D</mi> <mi>k</mi> </msub> </mrow> <mi>N</mi> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>9</mn> <mo>)</mo> </mrow> </mrow>
(2.4)将集合H内的所有点与相对应的权值相乘并求和,和值作为滤波窗口中心点的输出灰度值,代替该点原有的灰度值,如式(10)所示:
<mrow> <mi>f</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msub> <mi>w</mi> <mi>k</mi> </msub> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>&amp;times;</mo> <mi>H</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow>
(2.5)重复步骤(2.2)~(2.4)对图像的所有点进行遍历,实现图像的滤波处理。
CN201710528489.3A 2017-07-01 2017-07-01 一种去除高斯噪声和脉冲噪声的混合噪声滤波方法 Pending CN107392864A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710528489.3A CN107392864A (zh) 2017-07-01 2017-07-01 一种去除高斯噪声和脉冲噪声的混合噪声滤波方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710528489.3A CN107392864A (zh) 2017-07-01 2017-07-01 一种去除高斯噪声和脉冲噪声的混合噪声滤波方法

Publications (1)

Publication Number Publication Date
CN107392864A true CN107392864A (zh) 2017-11-24

Family

ID=60334009

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710528489.3A Pending CN107392864A (zh) 2017-07-01 2017-07-01 一种去除高斯噪声和脉冲噪声的混合噪声滤波方法

Country Status (1)

Country Link
CN (1) CN107392864A (zh)

Cited By (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108198567A (zh) * 2018-02-22 2018-06-22 成都启英泰伦科技有限公司 一种新型语音除噪系统
CN108810320A (zh) * 2018-06-01 2018-11-13 深圳市商汤科技有限公司 图像质量提升方法和装置
CN108898139A (zh) * 2018-06-04 2018-11-27 上海大学 一种下雨环境下的激光雷达数据抗干扰处理方法及其实验装置
CN108937950A (zh) * 2018-05-30 2018-12-07 梧州井儿铺贸易有限公司 一种智能机器人
CN109003247A (zh) * 2018-07-26 2018-12-14 吉林大学 彩色图像混合噪声的去除方法
CN109188148A (zh) * 2018-09-26 2019-01-11 国网安徽省电力有限公司铜陵市义安区供电公司 应用于智能电网的输电线路可靠监测系统
CN109407664A (zh) * 2018-09-26 2019-03-01 东莞绿邦智能科技有限公司 一种变电站自动读表机器人
CN109657658A (zh) * 2019-02-19 2019-04-19 江苏邦融微电子有限公司 一种处理图像的硬件加速系统及方法
CN109859133A (zh) * 2019-01-30 2019-06-07 南京邮电大学 一种中值滤波图像去噪方法
CN110738621A (zh) * 2019-10-17 2020-01-31 内蒙古工业大学 线性结构滤波方法、装置、设备和存储介质
CN110988232A (zh) * 2019-11-28 2020-04-10 南通乐尔环保科技有限公司 一种色谱基线降噪方法及其装置
CN111175618A (zh) * 2019-12-31 2020-05-19 浙江黑卡电气有限公司 一种适用于同时处理红外和局放数据的局部相关性滤波方法
CN111951173A (zh) * 2020-06-16 2020-11-17 五邑大学 高自由度滤波算法的调整方法和存储介质
CN112053302A (zh) * 2020-09-03 2020-12-08 五邑大学 高光谱图像的去噪方法、装置及存储介质
CN112187318A (zh) * 2020-08-31 2021-01-05 南京理工大学 基于深度学习的脉冲噪声削减方法
CN112419190A (zh) * 2020-11-24 2021-02-26 海南大学 一种基于局部统计信息和几何距离的脉冲噪声滤波方法、系统
CN112785513A (zh) * 2020-08-25 2021-05-11 青岛经济技术开发区海尔热水器有限公司 一种用于过滤脉冲噪声的自适应中值滤波方法
CN112862708A (zh) * 2021-01-27 2021-05-28 牛津仪器科技(上海)有限公司 图像噪声的自适应识别方法、传感器芯片以及电子设备
CN113066389A (zh) * 2019-12-26 2021-07-02 杭州海康威视数字技术股份有限公司 发光二极管(led)灯板处理方法及装置
CN113077398A (zh) * 2021-04-09 2021-07-06 上海申瑞继保电气有限公司 断路器圆形分合指示灯图像噪声滤波方法
CN113570507A (zh) * 2020-04-29 2021-10-29 浙江宇视科技有限公司 一种图像降噪方法、装置、设备及存储介质
CN115310053A (zh) * 2022-10-10 2022-11-08 之江实验室 一种强鲁棒性的图像内容审查保护方法和系统
CN117745572A (zh) * 2023-12-14 2024-03-22 国网江苏省电力有限公司南通供电分公司 一种红外图像的去噪方法
CN117788570A (zh) * 2024-02-26 2024-03-29 山东济矿鲁能煤电股份有限公司阳城煤矿 一种基于机器视觉的斗轮机定位方法及系统
CN117970859A (zh) * 2024-04-02 2024-05-03 山西顺达胜业通信工程有限公司 一种低功耗广域物联网下的设备智能控制方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102897192A (zh) * 2012-10-18 2013-01-30 成都唐源电气有限责任公司 城市轨道交通接触轨检测系统及其检测方法
CN103634537A (zh) * 2013-09-27 2014-03-12 中国人民解放军空军工程大学 一种高强度脉冲噪声滤除方法
CN104050646A (zh) * 2014-07-04 2014-09-17 西安电子科技大学 基于迭代非局部均值抑制图像脉冲噪声的方法
CN104751426A (zh) * 2015-04-09 2015-07-01 西安电子科技大学 基于三维块匹配变换去除高密度脉冲噪声的方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102897192A (zh) * 2012-10-18 2013-01-30 成都唐源电气有限责任公司 城市轨道交通接触轨检测系统及其检测方法
CN103634537A (zh) * 2013-09-27 2014-03-12 中国人民解放军空军工程大学 一种高强度脉冲噪声滤除方法
CN104050646A (zh) * 2014-07-04 2014-09-17 西安电子科技大学 基于迭代非局部均值抑制图像脉冲噪声的方法
CN104751426A (zh) * 2015-04-09 2015-07-01 西安电子科技大学 基于三维块匹配变换去除高密度脉冲噪声的方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
A.K.TRIPATHI 等: "Switching media filter:advanced boundary discriminative noise detection algorithm", 《IET IMAGE PROCESSING》 *
孔凡震: "一种开关中值滤波器的改进", 《太原科技》 *
沈德海 等: "一种改进的加权均值滤波算法", 《现代电子技术》 *
肖孟强: "混合噪声图像滤波算法在医学图像中的应用研究", 《中国优秀硕士学位论文全文数据库信息科技辑》 *
贺东霞: "数字图像去噪算法的研究与应用", 《中国优秀硕士学位论文全文数据库信息科技辑》 *

Cited By (43)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108198567A (zh) * 2018-02-22 2018-06-22 成都启英泰伦科技有限公司 一种新型语音除噪系统
CN108937950A (zh) * 2018-05-30 2018-12-07 梧州井儿铺贸易有限公司 一种智能机器人
CN108810320A (zh) * 2018-06-01 2018-11-13 深圳市商汤科技有限公司 图像质量提升方法和装置
CN108810320B (zh) * 2018-06-01 2020-11-24 深圳市商汤科技有限公司 图像质量提升方法和装置
CN108898139A (zh) * 2018-06-04 2018-11-27 上海大学 一种下雨环境下的激光雷达数据抗干扰处理方法及其实验装置
CN108898139B (zh) * 2018-06-04 2022-06-10 上海大学 一种下雨环境下的激光雷达数据抗干扰处理方法及其实验装置
CN109003247A (zh) * 2018-07-26 2018-12-14 吉林大学 彩色图像混合噪声的去除方法
CN109003247B (zh) * 2018-07-26 2021-06-15 吉林大学 彩色图像混合噪声的去除方法
CN109407664A (zh) * 2018-09-26 2019-03-01 东莞绿邦智能科技有限公司 一种变电站自动读表机器人
CN109188148B (zh) * 2018-09-26 2019-10-25 国网安徽省电力有限公司铜陵市义安区供电公司 应用于智能电网的输电线路可靠监测系统
CN109407664B (zh) * 2018-09-26 2022-08-19 托马斯智通机器人技术(常州)有限公司 一种变电站自动读表机器人
CN109188148A (zh) * 2018-09-26 2019-01-11 国网安徽省电力有限公司铜陵市义安区供电公司 应用于智能电网的输电线路可靠监测系统
CN109859133A (zh) * 2019-01-30 2019-06-07 南京邮电大学 一种中值滤波图像去噪方法
CN109859133B (zh) * 2019-01-30 2022-08-02 南京邮电大学 一种中值滤波图像去噪方法
CN109657658A (zh) * 2019-02-19 2019-04-19 江苏邦融微电子有限公司 一种处理图像的硬件加速系统及方法
CN110738621A (zh) * 2019-10-17 2020-01-31 内蒙古工业大学 线性结构滤波方法、装置、设备和存储介质
CN110738621B (zh) * 2019-10-17 2022-05-17 内蒙古工业大学 线性结构滤波方法、装置、设备和存储介质
CN110988232A (zh) * 2019-11-28 2020-04-10 南通乐尔环保科技有限公司 一种色谱基线降噪方法及其装置
CN113066389A (zh) * 2019-12-26 2021-07-02 杭州海康威视数字技术股份有限公司 发光二极管(led)灯板处理方法及装置
CN111175618A (zh) * 2019-12-31 2020-05-19 浙江黑卡电气有限公司 一种适用于同时处理红外和局放数据的局部相关性滤波方法
CN111175618B (zh) * 2019-12-31 2022-06-24 浙江黑卡电气有限公司 一种适用于同时处理红外和局放数据的局部相关性滤波方法
CN113570507A (zh) * 2020-04-29 2021-10-29 浙江宇视科技有限公司 一种图像降噪方法、装置、设备及存储介质
CN113570507B (zh) * 2020-04-29 2024-03-26 浙江宇视科技有限公司 一种图像降噪方法、装置、设备及存储介质
CN111951173B (zh) * 2020-06-16 2024-01-12 五邑大学 高自由度滤波算法的调整方法和存储介质
CN111951173A (zh) * 2020-06-16 2020-11-17 五邑大学 高自由度滤波算法的调整方法和存储介质
CN112785513B (zh) * 2020-08-25 2023-04-18 青岛经济技术开发区海尔热水器有限公司 一种用于过滤脉冲噪声的自适应中值滤波方法
CN112785513A (zh) * 2020-08-25 2021-05-11 青岛经济技术开发区海尔热水器有限公司 一种用于过滤脉冲噪声的自适应中值滤波方法
CN112187318B (zh) * 2020-08-31 2022-02-18 南京理工大学 基于深度学习的脉冲噪声削减方法
CN112187318A (zh) * 2020-08-31 2021-01-05 南京理工大学 基于深度学习的脉冲噪声削减方法
CN112053302B (zh) * 2020-09-03 2023-10-17 五邑大学 高光谱图像的去噪方法、装置及存储介质
CN112053302A (zh) * 2020-09-03 2020-12-08 五邑大学 高光谱图像的去噪方法、装置及存储介质
CN112419190B (zh) * 2020-11-24 2023-06-20 海南大学 一种基于局部统计信息和几何距离的脉冲噪声滤波方法、系统
CN112419190A (zh) * 2020-11-24 2021-02-26 海南大学 一种基于局部统计信息和几何距离的脉冲噪声滤波方法、系统
CN112862708B (zh) * 2021-01-27 2024-02-23 牛津仪器科技(上海)有限公司 图像噪声的自适应识别方法、传感器芯片以及电子设备
CN112862708A (zh) * 2021-01-27 2021-05-28 牛津仪器科技(上海)有限公司 图像噪声的自适应识别方法、传感器芯片以及电子设备
CN113077398A (zh) * 2021-04-09 2021-07-06 上海申瑞继保电气有限公司 断路器圆形分合指示灯图像噪声滤波方法
CN115310053A (zh) * 2022-10-10 2022-11-08 之江实验室 一种强鲁棒性的图像内容审查保护方法和系统
CN115310053B (zh) * 2022-10-10 2022-12-13 之江实验室 一种强鲁棒性的图像内容审查保护方法和系统
CN117745572A (zh) * 2023-12-14 2024-03-22 国网江苏省电力有限公司南通供电分公司 一种红外图像的去噪方法
CN117745572B (zh) * 2023-12-14 2024-06-04 国网江苏省电力有限公司南通供电分公司 一种红外图像的去噪方法
CN117788570A (zh) * 2024-02-26 2024-03-29 山东济矿鲁能煤电股份有限公司阳城煤矿 一种基于机器视觉的斗轮机定位方法及系统
CN117788570B (zh) * 2024-02-26 2024-05-07 山东济矿鲁能煤电股份有限公司阳城煤矿 一种基于机器视觉的斗轮机定位方法及系统
CN117970859A (zh) * 2024-04-02 2024-05-03 山西顺达胜业通信工程有限公司 一种低功耗广域物联网下的设备智能控制方法

Similar Documents

Publication Publication Date Title
CN107392864A (zh) 一种去除高斯噪声和脉冲噪声的混合噪声滤波方法
Prabu et al. Design of cognitive image filters for suppression of noise level in medical images
US7187811B2 (en) Method for image resolution enhancement
CN102521813B (zh) 基于双平台直方图的红外图像自适应增强方法
CN106355561A (zh) 基于噪点先验约束的全变分图像去噪方法
CN108053374A (zh) 一种结合双边滤波与Retinex的水下图像增强方法
CN105046677B (zh) 一种用于交通视频图像的增强处理方法和装置
CN103942758A (zh) 基于多尺度融合的暗通道先验图像去雾方法
CN102256048B (zh) 一种自适应密度的图像椒盐噪声开关滤波方法
CN107871310A (zh) 一种基于模糊核精细化的单幅图像盲去运动模糊方法
CN110033417A (zh) 一种基于深度学习的图像增强方法
CN104899905B (zh) 人脸图像处理方法及装置
CN103440644B (zh) 一种基于最小描述长度的多尺度图像弱边缘检测方法
CN106251318A (zh) 一种序列图像的去噪装置及方法
Jin et al. An efficient color-impulse detector and its application to color images
CN105809630A (zh) 一种图像噪声过滤方法及系统
CN107085833A (zh) 基于梯度倒数自适应开关均中值融合的遥感图像滤波方法
CN106228528A (zh) 一种基于决策图与稀疏表示的多聚焦图像融合方法
CN105590301B (zh) 自适应正斜双十字窗均值滤波的脉冲噪声消除方法
CN102670251B (zh) 基于数值排序加权平均的超声图像空间滤波方法
CN103996177A (zh) 一种无参考检测消除雪花噪声算法
Bhardwaj et al. Adaptive Neuro-Fuzzy Inference System (ANFIS) Based Edge Detection Technique
CN110400274A (zh) 一种车载红外行人检测用红外图像增强方法
CN105719257B (zh) 图像超高密度椒盐噪声的降除方法
CN107103294A (zh) 一种基于移动智能手机的驾驶员疲劳检测系统

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20171124