CN103778611A - 利用边缘检测的开关加权矢量中值滤波方法 - Google Patents
利用边缘检测的开关加权矢量中值滤波方法 Download PDFInfo
- Publication number
- CN103778611A CN103778611A CN201410037489.XA CN201410037489A CN103778611A CN 103778611 A CN103778611 A CN 103778611A CN 201410037489 A CN201410037489 A CN 201410037489A CN 103778611 A CN103778611 A CN 103778611A
- Authority
- CN
- China
- Prior art keywords
- image
- pixel
- window
- weight
- noise
- 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
Images
Abstract
本发明涉及数字图像处理技术领域,为提高矢量中值滤波的降噪能力,同时保护图像细节。为此,本发明采用的技术方案是,利用边缘检测的开关加权矢量中值滤波方法,包括如下步骤:通过比较噪声图像中的像素点与对噪声图像的R、G、B分量分别进行标量中值滤波而得到的一幅参考图像中对应像素点的相似性来进行噪声监测,对疑似噪声点通过一种基于边缘检测的加权方法进行矢量中值滤波,从而提高矢量中值滤波的降噪能力,同时保护图像细节。本发明主要应用于数字图像处理。
Description
技术领域
本发明涉及数字图像处理技术领域,具体讲,涉及利用边缘检测的开关加权矢量中值滤波方法。
背景技术
标量中值滤波将给定窗口中各像素点灰度值的中值代替窗口中心像素,用于消除灰度图像矢量中值滤波(vector median filter,VMF)算法常用于消除彩色图像中产生的脉冲噪声。将每个像素的R(红色),G(绿色),B(蓝色)数据看作一个矢量,对于给定大小的窗口,计算每个像素与其他像素点的欧氏距离之和,然后对结果进行排序,最小值对应的像素矢量作为输出。矢量中值滤波能够有效的减小彩色图像中的噪声,并且不会引入图像中原本并不存在的色彩,因此广泛应用在各种数字图像处理系统中。滤波窗口越大,滤波器对脉冲噪声的滤除效果越强,但也更容易破坏图像细节。开关矢量中值滤波是通过某种方法进行噪声检测,只对噪声点进行滤波处理,对于非噪声点则保持不变。加权矢量中值滤波则是对滤波窗口内的各个像素点施以不同的权重,从而更精确的控制滤波过程,提高滤波效果,但目前未见相关成熟技术的报道。
发明内容
本发明旨在解决克服现有技术的不足,提高矢量中值滤波的降噪能力,同时保护图像细节。为此,本发明采用的技术方案是,利用边缘检测的开关加权矢量中值滤波方法,包括如下步骤:通过比较噪声图像中的像素点与对噪声图像的R、G、B分量分别进行标量中值滤波而得到的一幅参考图像中对应像素点的相似性来进行噪声监测,对疑似噪声点通过一种基于边缘检测的加权方法进行矢量中值滤波,从而提高矢量中值滤波的降噪能力,同时保护图像细节。
具体步骤为:
首先取一个n×n的窗口,令Xi=(Ri,Gi,Bi),i=1,2,…,N;N=n×n,N表示窗口内像素点的像素矢量;然后:
a.噪声检测
比较Xi与X'i的相似性,其中Q=R、G、B,k表示I’中各彩色分量图像的像素亮度范围,ρ表示矢量距离,如果
则认为Xi不是噪声点,直接输出,否则对该点进行加权矢量中值滤波处理,其中T为预先设定的阈值;
b.利用图像边缘检测的加权矢量中值滤波
加权矢量中值滤波的公式如下:
ρij(Xi,Xj)=||Xi-Xj||2=[(Ri-Rj)2+(Gi-Gj)2+(Bi-Bj)2]1/2 (3)
其中ρij表示两个像素的欧式距离,N表示窗口的大小,即滤波窗口中像素点的个数,ωj表示权重,用XWVMF代替窗口中心像素点即完成该点的滤波。
加权方法如下:
1)将待处理的彩色噪声图像转化为灰度图像;
2)对该灰度图像进行标量中值滤波;
3)对滤波后的灰度图像进行图像边缘检测,得到图像边缘的二值图像;
4)根据边缘检测结果以及当前像素点与图像边缘的相对位置信息,将滤波窗口内的像素点划分为不同的区域,并对不同区域的像素点施以不同的权重:
i.如果窗口中心像素在图像边缘上,则位于图像边缘上的像素点划分为同一组,将它们的权重置为ωe,窗口内其他像素点划分为另一组,权重为ω3,令ωe>ω3;
ii.如果窗口中心像素不在图像边缘上,但是滤波窗口内存在图像边缘,则将中心像素的权重设置为ωc,同时将窗口内的其他像素点分为三组:位于图像边缘上的像素点分为一组,权重为ω2;与中心像素位于图像边缘同一侧的分为另一组,权重为ω1;其余的像素点分为一组,权重为ω3,令ωc>ω1>ω2>ω3;
iii.如果滤波窗口内不存在图像边缘,则中心像素的权重置为ωc,其余点的权重为ω1,令ωc>ω1;
ρ通常是L1范数或者L2范数。
本发明的技术特点及效果:
本发明最主要的贡献是一方面设计了一种脉冲噪声检测方法,用于区分彩色图像中的噪声点和非噪声点,另一方面设计了一种基于图像边缘检测的加权方法,应用在矢量中值滤波中。通过本发明设计的方法,能够有效地抑制彩色图像中的脉冲噪声,同时能够很好地保护图像的细节,改善矢量中值滤波的处理效果。
附图说明
图1相邻滤波窗口。
图2加权原理示意图。
图3算法流程图。
图4图像边缘提取(a)噪声图像;(b)直接对噪声图像进行边缘提取后的结果;(c)先滤波,然后提取的图像边缘。
图5本发明方法的效果图(a)原始图像;(b)传统矢量中值滤波结果;(c)本发明方法滤波结果。
具体实施方式
本发明涉及一种基于图像边缘检测的开关加权矢量中值滤波器,用于消除彩色图像中的脉冲噪声同时保护图像细节。
本发明设计的目的是通过比较噪声图像中的像素点与一幅参考图像(对噪声图像的R、G、B分量分别进行标量中值滤波而得)中对应像素点的相似性来进行噪声监测,对疑似噪声点通过一种基于边缘检测的加权方法进行矢量中值滤波,从而提高矢量中值滤波的降噪能力,同时保护图像细节。
首先取一个n×n的窗口(n可以取3,5,7,9,11等),令Xi=(Ri,Gi,Bi),(i=1,2,…,N;N=n×n)表示窗口内像素点的像素矢量。图1表示一个5×5的滤波窗口,图中数字代表窗口中各像素点的编号。本发明采用的滤波过程如下:
c.噪声检测
比较Xi与X'i的相似性。其中Q=R、G、B。k表示I’中各彩色分量图像的像素亮度范围。ρ表示某种矢量距离(通常是L1范数或者L2范数)。如果
则认为Xi不是噪声点,直接输出,否则对该点进行加权矢量中值滤波处理。其中T为预先设定的阈值。
d.利用图像边缘检测的加权矢量中值滤波
加权矢量中值滤波的公式如下:
ρij(Xi,Xj)=||Xi-Xj||2=[(Ri-Rj)2+(Gi-Gj)2+(Bi-Bj)2]1/2 (3)
其中ρij表示两个像素的欧式距离(也称为L2范数),N表示窗口的大小,即滤波窗口中像素点的个数,ωj表示权重。用XWVMF代替窗口中心像素点即完成该点的滤波。本发明设计的加权方法如下:
5)将待处理的彩色噪声图像转化为灰度图像
6)对该灰度图像进行标量中值滤波
7)对滤波后的灰度图像进行图像边缘检测,得到图像边缘的二值图像。
8)根据边缘检测结果以及当前像素点与图像边缘的相对位置信息,如图2所示(图中深灰色部为图像边缘),将滤波窗口内的像素点划分为不同的区域,并对不同区域的像素点施以不同的权重:
iv.如果窗口中心像素在图像边缘上,如图2(a),则位于图像边缘上的像素点划分为同一组,将它们的权重置为ωe。窗口内其他像素点划分为另一组,权重为ω3。令ωe>ω3。
v.如果窗口中心像素不在图像边缘上,但是滤波窗口内存在图像边缘,如图2(b)(c),则将中心像素的权重设置为ωc,同时将窗口内的其他像素点分为三组:位于图像边缘上的像素点分为一组,权重为ω2;与中心像素位于图像边缘同一侧的分为另一组,权重为ω1;其余的像素点分为一组,权重为ω3。令ωc>ω1>ω2>ω3。
vi.如果滤波窗口内不存在图像边缘,如图2(d),则中心像素的权重置为ωc,其余点的权重为ω1。令ωc>ω1。
本发明的算法流程图如图3。
本发明的一个实现方式是窗口大小为5×5,如图1所示。在噪声检测过程中,计算相似性时的阈值T取值范围为0.75~0.95。在利用图像边缘检测过程中,使用Canny边缘检测算法。图4中,(a)为噪声图像,(b)为直接对噪声图像提取边缘的结果,(c)为先滤波然后使用Canny算子进行边缘检测的结果。进行加权时,各权重系数的关系为ωe/ω3=2.5~3.5,ωc/ω1=3~4,ωc/ω2=4~5,ωc/ω3=7~10。图5分别为原始图像、使用传统矢量中值滤波处理后的图像以及使用本发明方法处理后的图像。
Claims (3)
1.一种利用边缘检测的开关加权矢量中值滤波方法,其特征是,包括如下步骤:通过比较噪声图像中的像素点与对噪声图像的R、G、B分量分别进行标量中值滤波而得到的一幅参考图像中对应像素点的相似性来进行噪声监测,对疑似噪声点通过一种基于边缘检测的加权方法进行矢量中值滤波,从而提高矢量中值滤波的降噪能力,同时保护图像细节。
2.如权利要求1所述的利用边缘检测的开关加权矢量中值滤波方法,其特征是,具体步骤为:
首先取一个n×n的窗口,令Xi=(Ri,Gi,Bi),i=1,2,…,N;N=n×n,N表示窗口内像素点的像素矢量;然后:
a.噪声检测
比较Xi与X'i的相似性,其中Q=R、G、B,k表示I’中各彩色分量图像的像素亮度范围,ρ表示矢量距离,如果
则认为Xi不是噪声点,直接输出,否则对该点进行加权矢量中值滤波处理,其中T为预先设定的阈值;
b.利用图像边缘检测的加权矢量中值滤波
加权矢量中值滤波的公式如下:
ρij(Xi,Xj)=||Xi-Xj||2=[(Ri-Rj)2+(Gi-Gj)2+(Bi-Bj)2]1/2 (3)
其中ρij表示两个像素的欧式距离,N表示窗口的大小,即滤波窗口中像素点的个数,ωj表示权重,用XWVMF代替窗口中心像素点即完成该点的滤波。
3.如权利要求1所述的利用边缘检测的开关加权矢量中值滤波方法,其特征是,加权方法如下:
加权方法如下:
1)将待处理的彩色噪声图像转化为灰度图像;
2)对该灰度图像进行标量中值滤波;
3)对滤波后的灰度图像进行图像边缘检测,得到图像边缘的二值图像;
4)根据边缘检测结果以及当前像素点与图像边缘的相对位置信息,将滤波窗口内的像素点划分为不同的区域,并对不同区域的像素点施以不同的权重:
i.如果窗口中心像素在图像边缘上,则位于图像边缘上的像素点划分为同一组,将它们的权重置为ωe,窗口内其他像素点划分为另一组,权重为ω3,令ωe>ω3;
ii.如果窗口中心像素不在图像边缘上,但是滤波窗口内存在图像边缘,则将中心像素的权重设置为ωc,同时将窗口内的其他像素点分为三组:位于图像边缘上的像素点分为一组,权重为ω2;与中心像素位于图像边缘同一侧的分为另一组,权重为ω1;其余的像素点分为一组,权重为ω3,令ωc>ω1>ω2>ω3;
iii.如果滤波窗口内不存在图像边缘,则中心像素的权重置为ωc,其余点的权重为ω1,令ωc>ω1;
ρ是L1范数或者L2范数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410037489.XA CN103778611B (zh) | 2014-01-26 | 2014-01-26 | 利用边缘检测的开关加权矢量中值滤波方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410037489.XA CN103778611B (zh) | 2014-01-26 | 2014-01-26 | 利用边缘检测的开关加权矢量中值滤波方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103778611A true CN103778611A (zh) | 2014-05-07 |
CN103778611B CN103778611B (zh) | 2016-08-17 |
Family
ID=50570810
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410037489.XA Expired - Fee Related CN103778611B (zh) | 2014-01-26 | 2014-01-26 | 利用边缘检测的开关加权矢量中值滤波方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103778611B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104700372A (zh) * | 2015-03-19 | 2015-06-10 | 天津大学 | 基于边缘保护的加权各向异性扩散滤波方法 |
CN104766278A (zh) * | 2015-03-19 | 2015-07-08 | 天津大学 | 基于自适应平滑系数的各向异性滤波方法 |
CN107005623A (zh) * | 2014-10-16 | 2017-08-01 | 三星电子株式会社 | 用于图像处理的方法和装置 |
WO2018166083A1 (zh) * | 2017-03-13 | 2018-09-20 | 深圳市嘉和顺信息科技有限公司 | 一种对比度增强和视觉显著度优化方法在高尔夫球场图中的应用 |
CN109618076A (zh) * | 2018-08-07 | 2019-04-12 | 吴秋琴 | 歌手音乐自适应下载方法 |
CN110599422A (zh) * | 2019-09-12 | 2019-12-20 | 合肥中科离子医学技术装备有限公司 | 一种基于边缘保护的加权平均椒盐噪声降噪算法 |
CN113033574A (zh) * | 2021-02-26 | 2021-06-25 | 天津大学 | 一种基于fpga的图像数据降噪系统及方法 |
US11634175B2 (en) | 2021-07-09 | 2023-04-25 | Institute Of Geology And Geophysics, Chinese Academy Of Sciences | Dip angle-steering median filtering method based on a niche differential evolution algorithm |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1341246A (zh) * | 1999-02-05 | 2002-03-20 | 三星电子株式会社 | 彩色图像处理方法和装置 |
US20030219168A1 (en) * | 2002-05-27 | 2003-11-27 | Takuji Kawakubo | Edge emphasizing circuit |
CN101340594A (zh) * | 2007-07-06 | 2009-01-07 | 北京大学软件与微电子学院 | 基于矢量外积和中值滤波去彩色图像脉冲噪声方法及系统 |
CN102663708A (zh) * | 2012-04-27 | 2012-09-12 | 飞依诺科技(苏州)有限公司 | 基于方向加权中值滤波的超声图像处理方法 |
-
2014
- 2014-01-26 CN CN201410037489.XA patent/CN103778611B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1341246A (zh) * | 1999-02-05 | 2002-03-20 | 三星电子株式会社 | 彩色图像处理方法和装置 |
US20030219168A1 (en) * | 2002-05-27 | 2003-11-27 | Takuji Kawakubo | Edge emphasizing circuit |
CN101340594A (zh) * | 2007-07-06 | 2009-01-07 | 北京大学软件与微电子学院 | 基于矢量外积和中值滤波去彩色图像脉冲噪声方法及系统 |
CN102663708A (zh) * | 2012-04-27 | 2012-09-12 | 飞依诺科技(苏州)有限公司 | 基于方向加权中值滤波的超声图像处理方法 |
Non-Patent Citations (2)
Title |
---|
PEI-ENG NG, ET AL.: ""A Switching Median Filter With Boundary Discriminative Noise Detection for Extremely Corrupted Images"", 《IEEE TRANSACTIONS ON IMAGE PROCESS》 * |
张毛女等: ""一种基于边缘检测的改进的中值滤波去噪方法"", 《计算机与现代化》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107005623A (zh) * | 2014-10-16 | 2017-08-01 | 三星电子株式会社 | 用于图像处理的方法和装置 |
CN107005623B (zh) * | 2014-10-16 | 2020-04-28 | 三星电子株式会社 | 用于图像处理的方法和装置 |
CN104700372A (zh) * | 2015-03-19 | 2015-06-10 | 天津大学 | 基于边缘保护的加权各向异性扩散滤波方法 |
CN104766278A (zh) * | 2015-03-19 | 2015-07-08 | 天津大学 | 基于自适应平滑系数的各向异性滤波方法 |
CN104700372B (zh) * | 2015-03-19 | 2017-08-04 | 天津大学 | 基于边缘保护的加权各向异性扩散滤波方法 |
WO2018166083A1 (zh) * | 2017-03-13 | 2018-09-20 | 深圳市嘉和顺信息科技有限公司 | 一种对比度增强和视觉显著度优化方法在高尔夫球场图中的应用 |
CN109618076A (zh) * | 2018-08-07 | 2019-04-12 | 吴秋琴 | 歌手音乐自适应下载方法 |
CN110599422A (zh) * | 2019-09-12 | 2019-12-20 | 合肥中科离子医学技术装备有限公司 | 一种基于边缘保护的加权平均椒盐噪声降噪算法 |
CN110599422B (zh) * | 2019-09-12 | 2022-05-31 | 合肥中科离子医学技术装备有限公司 | 一种基于边缘保护的加权平均椒盐噪声降噪方法 |
CN113033574A (zh) * | 2021-02-26 | 2021-06-25 | 天津大学 | 一种基于fpga的图像数据降噪系统及方法 |
US11634175B2 (en) | 2021-07-09 | 2023-04-25 | Institute Of Geology And Geophysics, Chinese Academy Of Sciences | Dip angle-steering median filtering method based on a niche differential evolution algorithm |
Also Published As
Publication number | Publication date |
---|---|
CN103778611B (zh) | 2016-08-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103778611A (zh) | 利用边缘检测的开关加权矢量中值滤波方法 | |
CN107220988B (zh) | 基于改进canny算子的零部件图像边缘提取方法 | |
CN105096280A (zh) | 处理图像噪声的方法及装置 | |
CN104978715A (zh) | 一种基于滤波窗口及参数自适应的非局部均值图像去噪方法 | |
CN103337053B (zh) | 一种基于开关非局部全变分的椒盐噪声污染图像滤波方法 | |
CN105913396A (zh) | 一种噪声估计的图像边缘保持混合去噪方法 | |
CN102974551A (zh) | 一种基于机器视觉的多晶硅太阳能检测分选的方法 | |
CN109377450B (zh) | 一种边缘保护的去噪方法 | |
CN107665348B (zh) | 一种变电站数字仪表的数字识别方法和装置 | |
CN105528768A (zh) | 一种图像去噪方法 | |
CN109558908B (zh) | 一种给定区域最优边缘的确定方法 | |
CN105139391B (zh) | 一种雾霾天气交通图像边缘检测方法 | |
CN103455991A (zh) | 一种多聚焦图像融合方法 | |
CN104167005A (zh) | 一种基于相似度函数的自适应加权椒盐噪声滤波方法 | |
CN104580829A (zh) | 一种太赫兹图像增强方法及系统 | |
CN104517269A (zh) | 一种太赫兹图像条纹处理方法及系统 | |
CN106296763A (zh) | 一种金属材料工业ct图像质量快速校正方法 | |
CN111179186A (zh) | 一种保护图像细节的图像去噪系统 | |
CN104036461A (zh) | 一种基于联合滤波的红外复杂背景抑制方法 | |
CN102118547A (zh) | 图像加权滤波的方法 | |
CN106485663A (zh) | 一种车道线图像增强方法及系统 | |
CN104517270A (zh) | 一种太赫兹图像处理方法及系统 | |
CN105205829A (zh) | 基于改进的二维Otsu的变电站红外图像分割方法 | |
CN104899553A (zh) | 一种抗密集杂草干扰的田间作物行提取方法 | |
CN102509265A (zh) | 基于灰度值差和局部能量的数字图像去噪方法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160817 Termination date: 20210126 |
|
CF01 | Termination of patent right due to non-payment of annual fee |