CN106530254A - 一种基于小波阈值函数与改进中值滤波融合的抑制图像混合噪声算法 - Google Patents
一种基于小波阈值函数与改进中值滤波融合的抑制图像混合噪声算法 Download PDFInfo
- Publication number
- CN106530254A CN106530254A CN201611000703.XA CN201611000703A CN106530254A CN 106530254 A CN106530254 A CN 106530254A CN 201611000703 A CN201611000703 A CN 201611000703A CN 106530254 A CN106530254 A CN 106530254A
- Authority
- CN
- China
- Prior art keywords
- noise
- image
- algorithm
- threshold function
- 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
Links
- 238000001914 filtration Methods 0.000 title claims abstract description 27
- 230000004927 fusion Effects 0.000 title claims abstract description 6
- 230000002401 inhibitory effect Effects 0.000 title abstract 3
- 239000006002 Pepper Substances 0.000 claims description 25
- 230000001629 suppression Effects 0.000 claims description 15
- 238000000034 method Methods 0.000 claims description 13
- 230000000694 effects Effects 0.000 claims description 8
- 238000001514 detection method Methods 0.000 claims description 7
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 238000002372 labelling Methods 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims description 3
- 230000011664 signaling Effects 0.000 claims description 3
- 238000012790 confirmation Methods 0.000 claims description 2
- 235000002566 Capsicum Nutrition 0.000 description 3
- 241000722363 Piper Species 0.000 description 3
- 235000016761 Piper aduncum Nutrition 0.000 description 3
- 235000017804 Piper guineense Nutrition 0.000 description 3
- 235000008184 Piper nigrum Nutrition 0.000 description 3
- 150000003839 salts Chemical class 0.000 description 3
- 230000008676 import Effects 0.000 description 2
- 230000010355 oscillation Effects 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 241000196324 Embryophyta Species 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002146 bilateral effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
Classifications
-
- G06T5/70—
-
- 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
-
- 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/20048—Transform domain processing
- G06T2207/20064—Wavelet transform [DWT]
Landscapes
- Image Processing (AREA)
- Facsimile Image Signal Circuits (AREA)
Abstract
本发明公开了一种基于小波阈值函数与改进中值滤波融合的抑制图像混合噪声算法，属于图像处理领域。该算法的实现步骤如下：步骤一：导入一幅混有椒盐噪声和高斯噪声的灰度图像；步骤二：采用改进的中值滤波算法检测和抑制灰度图像中的椒盐噪声；步骤三：利用极限的思想选取合适的控制变量 m、β构建新的小波阈值函数，从而抑制灰度图像中的高斯噪声。该算法能够灵活有效的抑制混合噪声，能够良好的保留图像细节，从而有效提高了图像的质量。
Description
技术领域
本发明属于图像处理领域，尤其涉及一种基于小波阈值函数与改进中值滤波融合的抑制图像混合噪声算法。
背景技术
在计算机高度普及的今天，数字图像处理技术在各个领域都得到了广泛的应用，数字图像处理技术也日益系统化和完善化。但图像在拍摄、采样、量化和传输等过程中，常常受到外界各种噪声的干扰,给图像的分割、识别等后续处理过程带来影响，使图像不能真实反映景物,图像质量严重下降。因此，抑制图像噪声是数字图像处理的重要环节和步骤。
高斯噪声和椒盐噪声是图像中最常见的两种噪声，但多数情况下两种噪声同时存在。
目前抑制混合噪声的方法有：对传统的单一使用中值、均值算法进行改进，自适应中值-加权均值混合滤波器算法，改进中值滤波和改进加权均值结合的图像混合噪声滤波算法，双边滤波思想抑制图像中的混合噪声等。这些算法对图像混合噪声虽有一定的抑制效果，但不能很好地保留图像的细节和有效提高图像的质量。
为了克服现有的抑制图像混合噪声算法存在的缺陷，能够灵活的对噪声分类一一识别抑制，实现对图像混合噪声快速而高效的抑制，从而提高图像的质量，本发明提供了一种基于小波阈值函数与改进中值滤波融合的抑制图像混合噪声算法。
发明内容
本发明提供了一种基于小波阈值函数与改进中值滤波融合的抑制图像混合噪声算法，该算法的的实现步骤如下：
步骤一：导入一幅混有椒盐噪声和高斯噪声的灰度图像；
步骤二：采用改进的中值滤波算法检测和抑制灰度图像中的椒盐噪声；
步骤三：利用极限的思想选取合适的控制变量m、β构建新的小波阈值函数，从而抑制灰度图像中的高斯噪声。
所述步骤二中检测和抑制灰度图像中椒盐噪声的实现方法为：首先在小窗口内初选噪声点，设是3×3模版w3的中心像素灰度值，和分别表示窗口内所有像素的最大灰度值和最小灰度值，标记矩阵F的表达式为：
若，则认为是候选噪声点，标记为1；否则，标记为0；
为保证噪声尽可能完全被检测，采用7 × 7的窗口对可疑噪声进行二次检测确认，二次检测的模型如下：
上式中，F(i, j)=1表示位置(i, j)处的像素为噪声点，F(i, j)=1表示位置(i, j)处的像素为信号点。
所述步骤三中通过构建新的小波阈值函数，从而抑制灰度图像中的高斯噪声的实现方法为：本发明提出了一种新的阈值函数，其表达式如下：
上式中，sgn（x）是符号函数，其表达式为：
是阈值，其表达式为：
是噪声的标准方差，其表达式为：
W是原始小波系数，是阈值化小波分解系数，α、m、β均是大于零的控制变量，其中，0＜α＜10，10＜β，其作用是调节的大小，在一定阈值范围内尽可能多的缩小与W之间的偏差；当时，的取值为0，这种处理会误将图像中的有用信息滤除掉，影响图像细节；因此，当时，取，通过不断调节m的值，在保留图像细节的情况下，从而有效滤除噪声。
与现有技术相比，本发明的有益效果是：先检测图像中的椒盐噪声，并采用改进的中值滤波算法抑制椒盐噪声，比传统的中值滤波算法抑制噪声的效果更明显，噪声的抑制更彻底。针对图像中存在高斯噪声的问题，传统软、硬阈值函数在图像噪声抑制中存在过分平滑、边缘振荡和有恒定偏差的缺点，因此本发明提出了一种新的阈值函数，从而有效抑制灰度图像中的高斯噪声，其效果远远优于传统的软硬阈值函数处理高斯噪声。该发明算法对图像混合噪声的抑制快速而高效，能够有效提高图像的质量，且算法简单，易于实现。
附图说明
图1是本发明的流程示意图；
图2是抑制椒盐噪声的流程图；
图3(a)是原始图像，图3(b)是椒盐噪声图像，图3(c)是普通中值滤波抑制的椒盐噪声图像，图3 (d)是极值滤波算法抑制的椒盐噪声图像，图3(e)是本发明算法抑制的椒盐噪声图像；
图4(a)是传统硬阈值函数图像，图4(b)是传统软阈值函数图像，图4(c)是本发明阈值函数图像；
图5（a）为本发明函数m取不同值时的图像，图5(b)为本发明函数β取不同值时的图像；
图6（a）为原始图像，图6（b）为混有椒盐噪声和高斯噪声的图像，图6（c）为中值滤波噪声抑制图像，图6（d）为均值滤波噪声抑制图像，图6（e）为本发明算法噪声抑制图像。
具体实施方式
下面结合附图和具体实施方式对本发明作进一步的说明。
本发明提出了一种基于小波阈值函数与改进中值滤波融合的抑制图像混合噪声算法，该算法的实现步骤如图1所示，主要通过以下步骤实现：
步骤一：导入一幅混有椒盐噪声和高斯噪声的灰度图像。
本发明实例中导入的待处理图像是大小为的Lena图像，椒盐噪声为0.2，高斯噪声均值为0，方差为20，如图6(a)所示。
步骤二：采用改进的中值滤波算法检测和抑制灰度图像中的椒盐噪声。
首先在小窗口内初选噪声点，设是3×3模版w3的中心像素灰度值，和分别表示窗口内所有像素的最大灰度值和最小灰度值，标记矩阵F的表达式为：
若，则认为是候选噪声点，标记为1；否则，标记为0；
但是，局部极值点不一定是噪声点，若把局部极值点都作为噪声点进行滤除，必然会造成不必要的细节丢失，为保证噪声尽可能完全被检测，采用7 × 7的窗口对可疑噪声进行二次检测确认，二次检测的模型如下：
上式中，F(i, j)=1表示位置(i, j)处的像素为噪声点，F(i, j)=1表示位置(i, j)处的像素为信号点。
本发明采用改进的中值滤波处理图像中的椒盐噪声，针对椒盐噪声图像，滤波窗口的尺寸一般大于或等于3×3，为使图像中的椒盐噪声得到大程度的抑制，一次滤波很难得到好的滤波效果。因此本发明采取3×3的正方形模版进行反复迭代滤波。其迭代思想是用3×3的窗口在图像上滑过，若窗口内像素都为噪声，则先不进行处理，转而去处理后面那些窗口内不全为噪声的像素，反复如此操作，直至迭代结束，具体流程如图2所示。
步骤三：利用极限的思想选取合适的控制变量m、β构建新的小波阈值函数，从而抑制灰度图像中的高斯噪声。
针对图像中存在高斯噪声的问题，传统软、硬阈值函数在图像噪声抑制中存在过分平滑、边缘振荡和有恒定偏差的缺点，本发明提出了一种新的阈值函数，其表达式如下：
上式中，sgn（x）是符号函数，其表达式为：
是阈值，其表达式为：
是噪声的标准方差，其表达式为：
W是原始小波系数，是阈值化小波分解系数，α、m、β均是大于零的控制变量，其中，0＜α＜10，10＜β，其作用是调节的大小，在一定阈值范围内尽可能多的缩小与W之间的偏差；当时，的取值为0，这种处理会误将图像中的有用信息滤除掉，影响图像细节；因此，当时，取，通过不断调节m的值，在保留图像细节的情况下，从而有效滤除噪声。
传统软阈值函数如图3(a)所示，硬阈值函数如图3(b)所示，本发明阈值函数曲线如图3(c)所示，其中x表示W，y表示。由图可知，本发明算法克服了硬阈值函数不连续的缺点，并且改进了软阈值函数存在恒定偏差的缺点。
图5(a)为固定α、β，当m取不同值时本发明阈值函数的曲线图，其中x表示W，y表示。由图可知当m取值大于5时，函数曲线过分接近于x轴，当m取值小于1时，函数曲线过分接近于y轴。通过试验验证，本发明将m取值为1.8。
图5(b)为固定α、m，当β取不同值时本发明阈值函数的曲线图，其中x表示W，y表示。由图可知β值越大在一定阈值范围内W越无限接近于，本发明将β值定为100，当β＞100时，在一定阈值范围内曲线变化幅度极小。
由于变量的作用是保证本文阈值函数的连续性，改变α的值对图像的变化幅度影响不是太大，本发明确定其范围是0＜α＜10，取其值为2。
分别使用中值滤波、均值滤波以及本发明算法对噪声图像的抑制进行对比，图6(a)为原始图像，图6(b)为混有椒盐噪声和高斯噪声的图像，图6(c)为中值滤波噪声抑制图像，图6(d)为均值滤波噪声抑制图像，图6(e)为本发明算法噪声抑制图像。
本发明采用均方误差MSE与峰值信噪比PSNR作为评价参数，其中
上式中，f(x,y)为原始图像，为估计图像，M×N为图像尺寸，PSNR的单位为dB。将本发明噪声抑制方法与现有的噪声抑制方法的信噪比增益和均方差比较见下表：
| 参数 | 中值滤波 | 均值滤波 | 本发明算法 |
| PSNR | 23.7117 | 22.6925 | 25.8772 |
| MSE | 276.6344 | 349.8085 | 168.0198 |
通过上表的客观数据可知，本发明与其他算法相比提高了峰值信噪比，降低了最小均方误差。
本发明提出的一种基于小波阈值函数与改进中值滤波融合的抑制图像混合噪声算法，，在检测出椒盐噪声的前提下，能够有效抑制椒盐噪声，是处理后续图像噪声的前提。针对高斯噪声，本发明利用极限的思想采用新的小波阈值函数，并选取合适的控制变量m、β，有效的减少了小波系数和原小波系数之间的偏差，抑制了高斯噪声，本发明在有效抑制图像中的混合噪声的同时，在一定程度上也保留了图像的细节，有效提高了图像的质量。
Claims (3)
1.一种基于小波阈值函数与改进中值滤波融合的抑制图像混合噪声算法，其特征在于，该算法的实现步骤如下：
步骤一、导入一幅混有椒盐噪声和高斯噪声的灰度图像；
步骤二、根据图像所含的噪声种类，首先检测出图像中的椒盐噪声，并采用改进中值滤波算法抑制图像中的椒盐噪声；
步骤三：利用极限的思想选取合适的控制变量m、β构建新的小波阈值函数，从而抑制灰度图像中的高斯噪声。
2.根据权利要求1所述的一种基于小波阈值函数与改进中值滤波融合的抑制图像混合噪声算法，其特征在于，所述步骤二中检测灰度图像中的椒盐噪声的实现方法为：首先在小窗口内初选噪声点，设 是3×3模版w3的中心像素灰度值， 和 分别表示窗口内所有像素的最大灰度值和最小灰度值，标记矩阵F的表达式为：
若 ，则认为 是候选噪声点，标记为1；否则，标记为0；
为保证噪声尽可能完全被检测，采用7 × 7的窗口对可疑噪声进行二次检测确认，二次检测的模型如下：
上式中，F(i,j)=1表示位置(i,j)处的像素为噪声点，F(i,j)=1表示位置(i,j)处的像素为信号点。
3.根据权利要求1所述的一种基于小波阈值函数与改进中值滤波融合的抑制图像混合噪声算法，其特征在于，所述步骤三中通过构建新的小波阈值函数，从而抑制灰度图像中的高斯噪声的实现方法为：本发明提出了一种新的阈值函数，其表达式如下：
上式中，sgn（x）是符号函数，其表达式为：
W是原始小波系数， 是阈值化小波分解系数，α、m、β均是大于零的控制变量，其中，0＜α＜10，10＜β，其作用是调节 的大小，在一定阈值范围内尽可能多的缩小 与W之间的偏差；当 时， 的取值为0，这种处理会误将图像中的有用信息滤除掉，影响图像细节；因此，当 时，取 ，通过不断调节m的值，在保留图像细节的情况下，从而有效滤除噪声。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611000703.XA CN106530254A (zh) | 2016-11-14 | 2016-11-14 | 一种基于小波阈值函数与改进中值滤波融合的抑制图像混合噪声算法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611000703.XA CN106530254A (zh) | 2016-11-14 | 2016-11-14 | 一种基于小波阈值函数与改进中值滤波融合的抑制图像混合噪声算法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106530254A true CN106530254A (zh) | 2017-03-22 |
Family
ID=58351579
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611000703.XA Pending CN106530254A (zh) | 2016-11-14 | 2016-11-14 | 一种基于小波阈值函数与改进中值滤波融合的抑制图像混合噪声算法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106530254A (zh) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107066981A (zh) * | 2017-04-20 | 2017-08-18 | 上海博历机械科技有限公司 | 针对小型无人机自主降落的视觉分级地标定位识别系统 |
| CN109345475A (zh) * | 2018-09-19 | 2019-02-15 | 长安大学 | 一种无人机遥感山区公路图像融合滤波方法 |
| CN109596165A (zh) * | 2018-11-23 | 2019-04-09 | 湖南城市学院 | 一种智能地理信息动态预警布控系统及方法 |
| CN111681174A (zh) * | 2020-04-29 | 2020-09-18 | 西南电子技术研究所（中国电子科技集团公司第十研究所） | 抑制椒盐/高斯混噪目标图像的联合滤波方法 |
| WO2021258832A1 (zh) * | 2020-06-23 | 2021-12-30 | 青岛科技大学 | 基于自适应窗口滤波和小波阈值优化的水声信号去噪方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101860667A (zh) * | 2010-05-06 | 2010-10-13 | 中国科学院西安光学精密机械研究所 | 一种快速去除图像中混合噪声的方法 |
| CN102141403A (zh) * | 2010-12-17 | 2011-08-03 | 北京航空航天大学 | 基于小波阈值去噪、中值滤波和均值滤波的实时混合去噪方法 |
-
2016
- 2016-11-14 CN CN201611000703.XA patent/CN106530254A/zh active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101860667A (zh) * | 2010-05-06 | 2010-10-13 | 中国科学院西安光学精密机械研究所 | 一种快速去除图像中混合噪声的方法 |
| CN102141403A (zh) * | 2010-12-17 | 2011-08-03 | 北京航空航天大学 | 基于小波阈值去噪、中值滤波和均值滤波的实时混合去噪方法 |
Non-Patent Citations (5)
| Title |
|---|
| 李金伦 等: "基于改进中值滤波和提升小波变换的阈值去噪方法研究", 《应用光学》 * |
| 武英 等: "一种自适应图像去噪混合滤波方法", 《计算机工程与应用》 * |
| 田流芳: "基于中值滤波和小波变换的图像去噪算法研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
| 邱志祺: "基于中值滤波与小波变换的图像去噪研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
| 陈晓曦 等: "小波阈值去噪法的深入研究", 《激光与红外》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107066981A (zh) * | 2017-04-20 | 2017-08-18 | 上海博历机械科技有限公司 | 针对小型无人机自主降落的视觉分级地标定位识别系统 |
| CN109345475A (zh) * | 2018-09-19 | 2019-02-15 | 长安大学 | 一种无人机遥感山区公路图像融合滤波方法 |
| CN109345475B (zh) * | 2018-09-19 | 2021-07-23 | 长安大学 | 一种无人机遥感山区公路图像融合滤波方法 |
| CN109596165A (zh) * | 2018-11-23 | 2019-04-09 | 湖南城市学院 | 一种智能地理信息动态预警布控系统及方法 |
| CN111681174A (zh) * | 2020-04-29 | 2020-09-18 | 西南电子技术研究所（中国电子科技集团公司第十研究所） | 抑制椒盐/高斯混噪目标图像的联合滤波方法 |
| CN111681174B (zh) * | 2020-04-29 | 2023-02-03 | 西南电子技术研究所（中国电子科技集团公司第十研究所） | 抑制椒盐/高斯混噪目标图像的联合滤波方法 |
| WO2021258832A1 (zh) * | 2020-06-23 | 2021-12-30 | 青岛科技大学 | 基于自适应窗口滤波和小波阈值优化的水声信号去噪方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106530254A (zh) | 一种基于小波阈值函数与改进中值滤波融合的抑制图像混合噪声算法 | |
| CN106355561B (zh) | 基于噪点先验约束的全变分图像去噪方法 | |
| CN108416740B (zh) | 一种消除椒盐噪声的迭代自适应中值滤波方法 | |
| CN101094312A (zh) | 一种保持边缘的自适应图像滤波方法 | |
| CN109003233B (zh) | 一种基于自适应权重全变分模型的图像去噪方法 | |
| Zhang et al. | Decision-based non-local means filter for removing impulse noise from digital images | |
| CN104835127A (zh) | 一种自适应平滑滤波方法 | |
| CN113902638B (zh) | 结合卷积神经网络和轻型自适应中值滤波视频降噪方法 | |
| CN103971345A (zh) | 一种基于改进双边滤波的图像去噪方法 | |
| Patvardhan et al. | Denoising of document images using discrete curvelet transform for ocr applications | |
| Sunilkumar et al. | Removal of high density impulse noise through modified non-linear filter | |
| CN116862809A (zh) | 一种低曝光条件下的图像增强方法 | |
| Sankaran et al. | Non local image restoration using iterative method | |
| Mohammed | An improved median filter based on efficient noise detection for high quality image restoration | |
| RU2405200C2 (ru) | Способ и устройство быстрого фильтрования шума цифровых изображений | |
| CN113506259A (zh) | 一种基于融合网络的图像模糊判别方法及系统 | |
| CN108961290B (zh) | 一种基于Otsu的Ratio算子自适应图像边缘检测方法 | |
| Delon et al. | A patch-based approach for random-valued impulse noise removal | |
| Venkatachalam | A novel algorithm for image denoising using modified adaptive median filter | |
| Iqbal et al. | High density impulse noise reduction by denoising neighbor pixels | |
| Chen et al. | An impulse noise reduction method by adaptive pixel-correlation | |
| Vijaykumar et al. | Adaptive window based efficient algorithm for removing gaussian noise in gray scale and color images | |
| Metri et al. | Patch based wiener filter for image denoising | |
| Lee et al. | Neural-network-based impulse noise removal using group-based weighted couple sparse representation | |
| CN108765308B (zh) | 一种基于卷积模板的图像去噪方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WD01 | Invention patent application deemed withdrawn after publication | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20170322 |