CN111242854B - 一种图像去噪方法 - Google Patents
一种图像去噪方法 Download PDFInfo
- Publication number
- CN111242854B CN111242854B CN202010005060.8A CN202010005060A CN111242854B CN 111242854 B CN111242854 B CN 111242854B CN 202010005060 A CN202010005060 A CN 202010005060A CN 111242854 B CN111242854 B CN 111242854B
- Authority
- CN
- China
- Prior art keywords
- image
- representing
- noise
- wavelet coefficients
- wavelet
- 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
- 238000000034 method Methods 0.000 title claims abstract description 44
- 238000001914 filtration Methods 0.000 claims abstract description 27
- 230000009466 transformation Effects 0.000 claims abstract description 15
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 8
- 230000000694 effects Effects 0.000 claims description 14
- 238000010606 normalization Methods 0.000 claims description 3
- 239000004973 liquid crystal related substance Substances 0.000 abstract 2
- 230000002708 enhancing effect Effects 0.000 abstract 1
- 238000012545 processing Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000002146 bilateral effect Effects 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000013441 quality evaluation Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000004800 variational method Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000004377 microelectronic Methods 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/10—Internal combustion engine [ICE] based vehicles
- Y02T10/40—Engine management systems
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
Abstract
本发明提供一种图像去噪方法。所述图像去噪方法包括以下步骤：S1：根据图像噪声服从高斯分布：Yi,j＝Xi,j+εi,j；式中，i,j＝1,2…n，Xi,j代表原始图像，εi,j是噪声，且独立分布于N(0,σ2)，且与原始图像相互独立，图像经小波变换后：其中，代表带噪图像小波系数，代表原始图像小波系数，代表噪声小波系数；S2：对各子带各方向小波系数进行计算：S21：对各部分噪声方差进行估计：S＝1,2,…n代表水平方向HLJ、垂直方向LHJ和对角线方向系数HHJ，J＝1,2…L，代表分解层数，L代表最大分解层数；S22：对各部分图像小波系数进行估计。本发明提供的图像去噪方法具有能有效地区分图像信号与噪声，增强了局部滤波的自适应性，同时能够更大程度上保留图像边缘细节信息部分的优点。
Description
技术领域
本发明涉及图像去噪技术领域，尤其涉及一种图像去噪方法。
背景技术
数字图像去噪算法研究涉及了光学、微电子技术、计算机科学、数学分析等学科领域，是一门综合性很强的学科，如今的理论体系已经十分完善，经过几十年的发展研究，已经实践应用在军事、医疗、农业、等多种方向，形成了以下几大类主要的传统图像去噪算法：
1)空间域滤波
空域滤波是直接在原始数字图像直接对像素的灰度值进行滤波处理，主要可分为线性滤波和非线性滤波方法两类。常用的线性滤波方法有：均值滤波、高斯滤波、空域维纳滤波等。线性滤波最明显的缺点是会丢失图像的边缘细节信息，使得去噪后图像边缘变得模糊。而非线性滤波方法，有效地保护了图像细节和边缘信息，广泛应用到了工业、医学等控制领域，典型的有中值滤波和双边滤波、自适应中值滤波。
2)变换域滤波
变换域滤波方法是对图像或信号从空间域转换到变换域空间中，再进行降噪处理的方法。经典方法如傅立叶变换(Fourier transform)、余弦变换(cosine transform)、K-L变换(Karhunen-Loeve Transform)、小波变换(wavelet transform,WT)等等。
这其中最具代表性的是傅里叶变换和小波变换。傅立叶变换是将时域信号转化为频域信号，进行分析，对于信号x(t)，其傅立叶变换形式为：
3)偏微分方程
偏微分方程是针对低层图像处理的一种去噪方法，取得了一定的去噪效果。该方法对于系数的处理较为简单，在噪声密度较低情况下，降噪性能较好，但是对于高噪声密度的情况下，降噪性能较低，具有良好的平滑图像和边缘尖锐化的能力。
4)变分法
变分法[2]的核心步骤就是找到匹配的能量函数方程，增强方法的鲁棒性能，去噪的效果才能达到最好，具有代表性的方法有全变分TV模型。
5)形态学噪声滤除
形态学(morphology)通常表示生物学分支，我们使用同一词语表示数学形态学的内容。形态学中的开操作是为了平滑物体的轮廓、断开较窄的狭颈并消除细的突出物。与开操作相反的，闭操作虽然也是平滑轮廓的一部分，但通常会弥合较窄的简短和细长的沟壑，消除小的孔洞，填补轮廓线中的断裂，开、闭操作是形态学中最重要的部分。
然而传统的图像去噪声算法区分图像信号与噪声较差，且局部滤波的自适应性较差，保留图像边缘细节信息较差。
因此，有必要提供一种新的图像去噪方法解决上述技术问题。
发明内容
本发明解决的技术问题是提供一种能有效地区分图像信号与噪声，增强了局部滤波的自适应性，同时能够更大程度上保留图像边缘细节信息部分的图像去噪方法。
为解决上述技术问题，本发明提供的图像去噪方法包括：以下步骤：
S1：根据图像噪声服从高斯分布：
Yi,j＝Xi,j+εi,j；
式中，i,j＝1,2…n，Xi,j代表原始图像，εi,j是噪声，且独立分布于N(0,σ2)，且与原始图像相互独立，图像经小波变换后：
其中，代表带噪图像小波系数，/>代表原始图像小波系数，/>代表噪声小波系数；
S2：对各子带各方向小波系数进行计算：
S21：对各部分噪声方差进行估计：
S＝1,2,…n代表水平方向HLJ、垂直方向LHJ和对角线方向系数HHJ，J＝1,2…L，代表分解层数，L代表最大分解层数；
S22：对各部分图像小波系数进行估计：
由于服从高斯分布，所以：
由得到：
这样就得到不同方向不同层级的阈值：
S23：针对阈值函数的选取：
采用软硬阈值折衷函数[i]：
S3：图像经过小波变换后，结合非局部均值算法处理小波变换后的低频系数：
非局部均值算法(NLM)输出图像定义：
其中，I为搜索区域，ω(i,j)代表权重，该值由匹配块的相似度决定；
块的相似度定义如下：
该值表示点i和j邻域差值平方卷积高斯核，Z(i)表示权重归一化系数，h为衰减因子。
优选的，所述S23中，为了获得更好的滤波效果：避免硬阈值函数门限不连续，使处理后的图像产生伪Gibbs失真效应；软阈值函数人为的减小边缘的小波系数，使图像变得模糊，使图像失去细节信息，采用软硬阈值折衷函数。
优选的，所述S3中，非局部均值算法处理小波变换后的低频系数用于解决图像经过小波变换后造成边缘模糊的问题，用于保护图像的边缘和细节信息。
优选的，所述S3中，利用积分图像技术对非局部均值算法进行加速，用以解决非局部均值算法的复杂度高的问题。
与相关技术相比较，本发明提供的图像去噪方法具有如下有益效果：
本发明提供一种图像去噪方法，所得到结果的峰值信噪比(PSNR)和结构相似性(SSIM)指数均高于其他算法上，各项图像质量评价指标都得到了有效地提升，图像去噪效果也得到了验证，在恢复图像轮廓信息方面有更好的效果，能有效地区分图像信号与噪声，增强了局部滤波的自适应性，在滤除图像高频噪声部分的同时，在更大程度上保留图像边缘细节信息部分，因此，量子去噪算法相比较其他几种方法的去噪效果和图像评价指标水平都得到了很大提高。
附图说明
图1为传统图像处理方法与本发明提供的图像去噪方法对图像处理后的对比图；
图2为本发明提供的图像去噪方法与其他传统方法得到的PSNR值；
图3为本发明提供的图像去噪方法与其他传统方法得到的SSIM值。
具体实施方式
下面结合附图和实施方式对本发明作进一步说明。
图像去噪方法包括以下步骤：
S1：根据图像噪声服从高斯分布：
Yi,j＝Xi,j+εi,j；
式中，i,j＝1,2…n，Xi,j代表原始图像，εi,j是噪声，且独立分布于N(0,σ2)，且与原始图像相互独立，图像经小波变换后：
其中，代表带噪图像小波系数，/>代表原始图像小波系数，/>代表噪声小波系数；
S2：对各子带各方向小波系数进行计算：
S21：对各部分噪声方差进行估计：
S＝1,2,…n代表水平方向HLJ、垂直方向LHJ和对角线方向系数HHJ，J＝1,2…L，代表分解层数，L代表最大分解层数；
S22：对各部分图像小波系数进行估计：
由于服从高斯分布，所以：
由可以得到：
这样就可以得到不同方向不同层级的阈值：
S23：针对阈值函数的选取：
采用软硬阈值折衷函数[ii]：
S3：图像经过小波变换后，结合非局部均值算法处理小波变换后的低频系数：
非局部均值算法(NLM)输出图像定义：
其中，I为搜索区域，ω(i,j)代表权重，该值由匹配块的相似度决定；
块的相似度定义如下：
该值表示点i和j邻域差值平方卷积高斯核，Z(i)表示权重归一化系数，h为衰减因子。
所述S23中，为了获得更好的滤波效果：避免硬阈值函数门限不连续，使处理后的图像产生伪Gibbs失真效应；软阈值函数人为的减小边缘的小波系数，使图像变得模糊，使图像失去细节信息，采用软硬阈值折衷函数。
所述S3中，非局部均值算法处理小波变换后的低频系数用于解决图像经过小波变换后造成边缘模糊的问题，用于保护图像的边缘和细节信息。
所述S3中，利用积分图像技术对非局部均值算法进行加速，用以解决非局部均值算法的复杂度高的问题。
结合非局部均值的小波去噪图像数据去噪算法在峰值信噪比(PSNR)和结构相似性(SSIM)两个衡量指标全部高于VishShrink算法、Bayes算法、中值滤波算法和结合双边滤波的小波去噪滤波算法。
请参阅图1-图3在添加均值为0，方差为0.5高斯噪声图像中，结合非局部均值的小波去噪算法在边缘保持指数和结构相似性(SSIM)指数高于其他算法。
与相关技术相比较，本发明提供的图像去噪方法具有如下有益效果：
本发明提供一种图像去噪方法，所得到结果的峰值信噪比(PSNR)和结构相似性(SSIM)指数均高于其他算法上，各项图像质量评价指标都得到了有效地提升，图像去噪效果也得到了验证，在恢复图像轮廓信息方面有更好的效果，能有效地区分图像信号与噪声，增强了局部滤波的自适应性，在滤除图像高频噪声部分的同时，在更大程度上保留图像边缘细节信息部分，因此，量子去噪算法相比较其他几种方法的去噪效果和图像评价指标水平都得到了很大提高。
以上所述仅为本发明的实施例，并非因此限制本发明的专利范围，凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换，或直接或间接运用在其它相关的技术领域，均同理包括在本发明的专利保护范围内。
Claims (4)
1.一种图像去噪方法，其特征在于，包括以下步骤：
S1：根据图像噪声服从高斯分布：
Yi,j＝Xi,j+εi,j；
式中，i,j＝1,2…n，Xi,j代表原始图像，εi,j是噪声，且独立分布于N(0,σ2)，且与原始图像相互独立，图像经小波变换后：
其中，代表带噪图像小波系数，/>代表原始图像小波系数，/>代表噪声小波系数；
S2：对各子带各方向小波系数进行计算：
S21：对各部分噪声方差进行估计：
S＝1,2,…n代表水平方向HLJ、垂直方向LHJ和对角线方向系数HHJ，J＝1,2…L，代表分解层数，L代表最大分解层数；
S22：对各部分图像小波系数进行估计：
由于服从高斯分布，所以：
由得到：
这样就得到不同方向不同层级的阈值：
S23：针对阈值函数的选取：
采用软硬阈值折衷函数[i]：
S3：图像经过小波变换后，结合非局部均值算法处理小波变换后的低频系数：
非局部均值算法(NLM)输出图像定义：
其中，I为搜索区域，ω(i,j)代表权重，该值由匹配块的相似度决定；
块的相似度定义如下：
该值表示点i和j邻域差值平方卷积高斯核，Z(i)表示权重归一化系数，h为衰减因子。
2.根据权利要求1所述的图像去噪方法，其特征在于，所述S23中，为了获得更好的滤波效果：避免硬阈值函数门限不连续，使处理后的图像产生伪Gibbs失真效应；软阈值函数人为的减小边缘的小波系数，使图像变得模糊，使图像失去细节信息，采用软硬阈值折衷函数。
3.根据权利要求1所述的图像去噪方法，其特征在于，所述S3中，非局部均值算法处理小波变换后的低频系数用于解决图像经过小波变换后造成边缘模糊的问题，用于保护图像的边缘和细节信息。
4.根据权利要求1所述的图像去噪方法，其特征在于，所述S3中，利用积分图像技术对非局部均值算法进行加速，用以解决非局部均值算法的复杂度高的问题。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010005060.8A CN111242854B (zh) | 2020-01-03 | 2020-01-03 | 一种图像去噪方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010005060.8A CN111242854B (zh) | 2020-01-03 | 2020-01-03 | 一种图像去噪方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111242854A CN111242854A (zh) | 2020-06-05 |
| CN111242854B true CN111242854B (zh) | 2023-09-01 |
Family
ID=70877641
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010005060.8A Active CN111242854B (zh) | 2020-01-03 | 2020-01-03 | 一种图像去噪方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111242854B (zh) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112183259B (zh) * | 2020-09-17 | 2022-11-04 | 昆明理工大学 | 一种基于ceemd与峭度加权平均阈值去噪的滚动轴承故障诊断方法 |
| CN112348031A (zh) * | 2020-11-17 | 2021-02-09 | 安徽理工大学 | 一种去除指纹图像混合噪声的改进小波阈值去噪方法 |
| CN112634164B (zh) * | 2020-12-29 | 2022-09-30 | 哈尔滨工业大学芜湖机器人产业技术研究院 | 一种基于多尺度哈尔小波变换的图像去噪方法 |
| CN112907482A (zh) * | 2021-03-17 | 2021-06-04 | 浙江理工大学 | 一种美图图像去噪的方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104318527A (zh) * | 2014-10-21 | 2015-01-28 | 浙江工业大学 | 基于小波变换和引导滤波器的医学超声图像去噪方法 |
| CN105550997A (zh) * | 2015-12-08 | 2016-05-04 | 天津津航计算技术研究所 | 一种基于多个变换域的三维匹配图像去噪方法 |
| CN105913393A (zh) * | 2016-04-08 | 2016-08-31 | 暨南大学 | 一种自适应小波阈值图像去噪算法及装置 |
| CN106023103A (zh) * | 2016-05-16 | 2016-10-12 | 济南大学 | 一种基于精确局部方差先验建模的自适应正交小波图像去噪方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| AU2014202322A1 (en) * | 2014-04-29 | 2015-11-12 | Canon Kabushiki Kaisha | Wavelet denoising of fringe image |
-
2020
- 2020-01-03 CN CN202010005060.8A patent/CN111242854B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104318527A (zh) * | 2014-10-21 | 2015-01-28 | 浙江工业大学 | 基于小波变换和引导滤波器的医学超声图像去噪方法 |
| CN105550997A (zh) * | 2015-12-08 | 2016-05-04 | 天津津航计算技术研究所 | 一种基于多个变换域的三维匹配图像去噪方法 |
| CN105913393A (zh) * | 2016-04-08 | 2016-08-31 | 暨南大学 | 一种自适应小波阈值图像去噪算法及装置 |
| CN106023103A (zh) * | 2016-05-16 | 2016-10-12 | 济南大学 | 一种基于精确局部方差先验建模的自适应正交小波图像去噪方法 |
Non-Patent Citations (1)
| Title |
|---|
| 基于轮廓波变换的图像混合噪声去除算法研究;许友雷等;《现代计算机》;20191125(第33期);第47-53页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111242854A (zh) | 2020-06-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111242854B (zh) | 一种图像去噪方法 | |
| Zhou et al. | Retinex-based laplacian pyramid method for image defogging | |
| Padmavathi et al. | Comparison of filters used for underwater image pre-processing | |
| Wang et al. | Enhanced ridge structure for improving fingerprint image quality based on a wavelet domain | |
| Nikpour et al. | Using diffusion equations for improving performance of wavelet-based image denoising techniques | |
| CN111696061A (zh) | 空间域与变换域相融合的噪声图像重构方法 | |
| Raj et al. | Medical image denoising using multi-resolution transforms | |
| Vijayakumar | Selective image enhancement and restoration for skin cancer identification | |
| Qian | Image denoising algorithm based on improved wavelet threshold function and median filter | |
| CN111192204A (zh) | 一种图像增强方法、系统和计算机可读存储介质 | |
| CN112750090A (zh) | 一种改进小波阈值的水下图像降噪方法及系统 | |
| Zhang et al. | Diffusion scheme using mean filter and wavelet coefficient magnitude for image denoising | |
| Sun et al. | Color image denoising based on guided filter and adaptive wavelet threshold | |
| CN111461999A (zh) | 一种基于超像素相似性测量的sar图像相干斑抑制方法 | |
| Bhonsle et al. | De-noising of CT images using combined bivariate shrinkage and enhanced total variation technique | |
| Iqbal et al. | Reduction of speckle noise in medical images using stationary wavelet transform and fuzzy logic | |
| Bhargava et al. | An Effective Method for Image Denoising Using Non-local Means and Statistics based Guided Filter in Nonsubsampled Contourlet Domain. | |
| Dhiman et al. | An improved threshold estimation technique for image denoising using wavelet thresholding techniques | |
| Singh | Wavelet transform in image processing: Denoising, segmentation and compression of digital, images | |
| Yan et al. | A Power Thresholding Function-based Wavelet Image Denoising Method. | |
| Ehsaeyan | A new shearlet hybrid method for image denoising | |
| Xu | Color enhancement algorithm for visual communication posters based on homomorphic filtering | |
| Jung et al. | Adaptive image denoising in scale-space using the wavelet transform | |
| Shi et al. | Image sharpening via image denoising in the complex wavelet domain | |
| Parihar et al. | External Filtering and Wavelet Domain Thresholding-based Denoising Method for AWGN corrupted images |
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 |