CN107194889B - 分块双边全变分正则化图像噪声消除方法 - Google Patents
分块双边全变分正则化图像噪声消除方法 Download PDFInfo
- Publication number
- CN107194889B CN107194889B CN201710330617.3A CN201710330617A CN107194889B CN 107194889 B CN107194889 B CN 107194889B CN 201710330617 A CN201710330617 A CN 201710330617A CN 107194889 B CN107194889 B CN 107194889B
- Authority
- CN
- China
- Prior art keywords
- image
- iteration
- total variation
- entering
- variation regularization
- 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 52
- 230000002146 bilateral effect Effects 0.000 title claims abstract description 47
- 230000008030 elimination Effects 0.000 title claims abstract description 13
- 238000003379 elimination reaction Methods 0.000 title claims abstract description 13
- 239000011159 matrix material Substances 0.000 claims abstract description 22
- 238000002945 steepest descent method Methods 0.000 claims abstract description 10
- 239000000470 constituent Substances 0.000 claims abstract description 7
- 238000011109 contamination Methods 0.000 claims abstract description 4
- 239000000126 substance Substances 0.000 claims description 9
- 230000000903 blocking effect Effects 0.000 claims description 3
- 238000009499 grossing Methods 0.000 claims description 3
- 238000001914 filtration Methods 0.000 description 10
- 230000001629 suppression Effects 0.000 description 6
- 230000000694 effects Effects 0.000 description 3
- 230000014759 maintenance of location Effects 0.000 description 3
- 230000015556 catabolic process Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
Images
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
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
Abstract
本发明涉及分块双边全变分正则化图像噪声消除方法,其包括以下步骤:(1)获取污染图像X0,并用污染图像X0初始化第1次迭代的去噪图像然后进入步骤(2);(2)计算第t次迭代去噪图像的分块双边结构相似距离矩阵DWt,然后进入步骤(3);(3)构建第t次迭代去噪图像的分块双边全变分正则项然后进入步骤(4);(4)构建由保真项和分块双边全变分正则项组成的能量泛函Et,进入步骤(5);(5)采用最陡下降法求解能量泛函Et的最小化问题,获得第t+1次迭代的去噪图像进入步骤(6);(6)判断迭代次数是否小于最大迭代次数N,如果迭代次数小于最大迭代次数N,令t=t+1并进入步骤(2);否则,输出第t+1次迭代的去噪图像结束本次操作。
Description
技术领域
本发明属于数字图像处理领域,具体涉及分块双边全变分正则化图像噪声消除方法,可用于数字图像的预处理过程。
背景技术
数字图像在获取和传输过程中,因受到电路自身和外界噪声源的干扰,不可避免地出现降质,从而严重影响到后续的特征提取与分析等处理。
图像去噪需要兼顾噪声的抑制和图像原始信息的保留两个方面。针对图像去噪问题,该领域的研究者进行了大量的探索和研究,提出了大量基于数字信号处理技术的噪声抑制方法。
现存的图像噪声抑制方法主要有空间域噪声抑制方法和变换域噪声抑制方法。空间域方法直接对像素输出进行处理,易于硬件实现,而获得广泛应用。传统的空域平均滤波法、中值滤波法、高斯滤波法、双边滤波法等,在滤除噪声的同时也会对图像边缘进行不同程度的平滑,Fasiu等人将边缘保持特性较好的全变分模型与双边滤波模型结合构建出双边全变分先验约束,取得了更好的去噪效果,输出图像具有更加锐利的边缘(Sina Farsiu,M.Dirk Robinson,Michael Elad et al.,Fast and Robust Multiframe SuperResolution,IEEE Transactions on Image Processing,Vol.13,No.10,2004,pp.1327-1344.)。但该方法中滤波权重仅涉及空间距离相似度,而忽略了亮度距离相似度,且滤波权重仍由单个像素亮度值决定,在对原图像局部结构特征的保持上,还有待提高。
综上前述,寻求一种在去噪同时保留更多的图像结构、细节信息的噪声消除方法显得尤为重要。
发明内容
发明目的:本发明针对上述现有技术存在的问题做出改进,即本发明的目的在于将基于块的双边滤波与全变分模型相结合,构建基于图像结构相似性的先验约束,在此基础上,提出一种新颖的分块双边全变分正则化噪声消除方法,以保留更多的图像结构信息,获得更高的信噪比。
技术方案:分块双边全变分正则化图像噪声消除方法,其用分块邻域相似度权取代传统双边滤波中的像素亮度距离权,构建分块双边全变分正则化项,并作为先验约束加入迭代回归估计过程中。
分块双边全变分正则化图像噪声消除方法,包括以下步骤:
||·||1表示L1范数。
α|m|+|l|为空间距离权;
dt为高斯平滑参数;
P为移位半径。
其中,λ为规整参数;
其中,Δt为迭代步长;
sign(·)代表符号函数。
有益效果:本发明公开的分块双边全变分正则化图像噪声消除方法与现有技术相比,具有以下优点:
本发明的方法在能量泛函中引入分块双边全变分正则项,与现有的双边全变分正则项图像去噪方法相比,充分利用了图像分块邻域的相似信息,可以在获得良好的去噪效果的同时,保留更多的图像细节信息。
附图说明
图1是本发明公开的分块双边全变分正则化图像噪声消除方法的流程图。
具体实施方式:
下面对本发明的具体实施方式详细说明。
参照图1,本发明的分块双边全变分正则化图像噪声消除方法,将分块双边全变分正则化引入数字图像消除中,具体包括如下步骤:
||·||1表示L1范数。
α|m|+|l|为空间距离权,取值为α=0.635;
dt为高斯平滑参数,取值为dt=1050;
P为移位半径,取值为P=2。
其中,λ为规整参数,取值为0.05<λ<0.2;
其中,Δt为迭代步长,取值为0.3<Δt<1.0;
sign(·)代表符号函数。
下面,分别采用峰值信噪比(PSNR)和结构相似度(SSIM)来量化对照评估本发明提出的分块双边全变分正则化图像噪声消除方法与现有全变分方法以及双边全变分方法的性能。
(1)、峰值信噪比
峰值信噪比PSNR的数学表达式为:
I(i,j)为未受噪声污染图像坐标(i,j)处的像素值;
I0(i,j)为校正后图像坐标(i,j)处的像素值;
M为图像的长度;
N为图像的宽度;
PSNR取值越大表明对噪声抑制越好。
(2)、结构相似度
结构相似度SSIM的数学表达式为:
SSIM(X,Y)=[I(X,Y)]α·[c(X,Y)]β·[s(X,Y)]γ,其中:
α为亮度比较项指数;
β为对比度比较项指数;
SSIM取值越大说明校正后图像与原始场景结构相似度越高,细节保持能力越好。
分别采用全变分方法、双边全变分方法和本发明方法对带有均值为0,方差为25的高斯白噪声的标准测试图像barbara进行100次迭代去噪。
对比测试的实验结果见表1。
表1对比测试结果的量化参数对比表
由表1可见:
(1)经本发明提出的方法去噪后的图像的峰值信噪比(PSNR)明显高于全变分法和双边全变分法的结果,说明经本发明方法去噪后的图像残留噪声更小。
(2)经本发明提出的方法去噪后的图像的结构相似性系数(SSIM)明显高于全变分法和双边全变分法的结果,说明经本发明方法去噪后的图像保留了原始场景更多的结构信息。
上述结果充分说明了,本发明的方法去噪效果更好,保留了原始场景更多的结构信息,具有更好的去噪性能和结构保持特性。
综上,本发明所提出的分块双边全变分正则化图像噪声去除方法可在有效去除噪声的同时更大限度地保留原始图像的边缘和细节等结构信息。
上面对本发明的实施方式做了详细说明。但是本发明并不限于上述实施方式,在所属技术领域普通技术人员所具备的知识范围内,还可以在不脱离本发明宗旨的前提下做出各种变化。
Claims (5)
1.分块双边全变分正则化图像噪声消除方法,其特征在于,包括以下步骤:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710330617.3A CN107194889B (zh) | 2017-05-11 | 2017-05-11 | 分块双边全变分正则化图像噪声消除方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710330617.3A CN107194889B (zh) | 2017-05-11 | 2017-05-11 | 分块双边全变分正则化图像噪声消除方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107194889A CN107194889A (zh) | 2017-09-22 |
CN107194889B true CN107194889B (zh) | 2021-02-02 |
Family
ID=59873443
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710330617.3A Active CN107194889B (zh) | 2017-05-11 | 2017-05-11 | 分块双边全变分正则化图像噪声消除方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107194889B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107945145A (zh) * | 2017-11-17 | 2018-04-20 | 西安电子科技大学 | 基于梯度置信变分模型的红外图像融合增强方法 |
CN108198139B (zh) * | 2017-12-01 | 2021-09-10 | 西安电子科技大学 | 基于传导全变分正则化的图像去噪方法 |
CN108550122B (zh) * | 2017-12-29 | 2021-10-29 | 西安电子科技大学 | 基于自主路径分块传导滤波的图像去噪方法 |
CN108510449B (zh) * | 2018-02-05 | 2020-06-30 | 西安电子科技大学 | 基于自适应核回归的全变分图像噪声消除方法 |
CN112561820B (zh) * | 2020-12-17 | 2023-10-27 | 三峡大学 | 一种适用于超声图像去噪的自适应加权混合总变分方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101807292A (zh) * | 2010-01-07 | 2010-08-18 | 清华大学 | 一种图像去噪方法 |
WO2012175010A1 (en) * | 2011-06-24 | 2012-12-27 | Technicolor (China) Technology Co., Ltd. | Method and device for processing of an image |
CN103868601A (zh) * | 2014-03-26 | 2014-06-18 | 西安电子科技大学 | Irfpa探测器非均匀响应的双边全变分正则化校正方法 |
CN104809461A (zh) * | 2015-05-08 | 2015-07-29 | 内蒙古科技大学 | 结合序列图像超分辨率重建的车牌识别方法及系统 |
CN105678715A (zh) * | 2016-02-15 | 2016-06-15 | 江苏大学 | 一种基于核回归全变分的去除图像噪声的方法 |
CN106254720A (zh) * | 2016-07-19 | 2016-12-21 | 四川大学 | 一种基于联合正则化的视频超分辨率重建方法 |
CN106251315A (zh) * | 2016-08-23 | 2016-12-21 | 南京邮电大学 | 一种基于全变分的图像去噪方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060291751A1 (en) * | 2004-12-16 | 2006-12-28 | Peyman Milanfar | Robust reconstruction of high resolution grayscale images from a sequence of low-resolution frames (robust gray super-resolution) |
-
2017
- 2017-05-11 CN CN201710330617.3A patent/CN107194889B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101807292A (zh) * | 2010-01-07 | 2010-08-18 | 清华大学 | 一种图像去噪方法 |
WO2012175010A1 (en) * | 2011-06-24 | 2012-12-27 | Technicolor (China) Technology Co., Ltd. | Method and device for processing of an image |
CN103868601A (zh) * | 2014-03-26 | 2014-06-18 | 西安电子科技大学 | Irfpa探测器非均匀响应的双边全变分正则化校正方法 |
CN104809461A (zh) * | 2015-05-08 | 2015-07-29 | 内蒙古科技大学 | 结合序列图像超分辨率重建的车牌识别方法及系统 |
CN105678715A (zh) * | 2016-02-15 | 2016-06-15 | 江苏大学 | 一种基于核回归全变分的去除图像噪声的方法 |
CN106254720A (zh) * | 2016-07-19 | 2016-12-21 | 四川大学 | 一种基于联合正则化的视频超分辨率重建方法 |
CN106251315A (zh) * | 2016-08-23 | 2016-12-21 | 南京邮电大学 | 一种基于全变分的图像去噪方法 |
Non-Patent Citations (5)
Title |
---|
An Iterative Super-Resolution Reconstruction of Image Sequences using Fast Affine Block-Based Registration with BTV Regularization;V. Patanavijit et al;《APCCAS 2006 - 2006 IEEE Asia Pacific Conference on Circuits and Systems》;20061207;全文 * |
Fast Gradient-Based Algorithms for Constrained Total Variation Image Denoising and Deblurring Problems;Amir Beck et al;《IEEE TRANSACTIONS ON IMAGE PROCESSING》;20091130;第18卷(第11期);全文 * |
Total variation regularized iterative back-projection method for single frame image super resolution;Rui Lai et al;《2012 IEEE 11th International Conference on Signal Processing》;20121025;全文 * |
一种改进全变差正则化的Shearlet自适应带钢图像去噪算法;韩英莉;《表面技术》;20141231;第43卷(第6期);全文 * |
基于Tikhonov和全变分正则化混合约束盲去模糊方法;李俭兵等;《南京邮电大学学报(自然科学版)》;20160630;第36卷(第3期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN107194889A (zh) | 2017-09-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107194889B (zh) | 分块双边全变分正则化图像噪声消除方法 | |
CN108921800B (zh) | 基于形状自适应搜索窗口的非局部均值去噪方法 | |
CN109671029B (zh) | 基于伽马范数最小化的图像去噪方法 | |
CN111369457B (zh) | 稀疏判别张量鲁棒性pca的遥感图像去噪方法 | |
Wen et al. | Joint adaptive sparsity and low-rankness on the fly: An online tensor reconstruction scheme for video denoising | |
CN108932699B (zh) | 基于变换域的三维匹配调和滤波图像去噪方法 | |
CN111681174A (zh) | 抑制椒盐/高斯混噪目标图像的联合滤波方法 | |
WO2016119745A1 (zh) | 一种视频去噪方法及装置 | |
CN113989168B (zh) | 一种针对椒盐噪声的自适应非局部均值滤波方法 | |
Zhang et al. | Kernel Wiener filtering model with low-rank approximation for image denoising | |
CN108460723B (zh) | 基于邻域相似性的双边全变分图像超分辨率重建方法 | |
CN108564544A (zh) | 基于边缘感知的图像盲去模糊组合稀疏优化方法 | |
CN105723416B (zh) | 图像降噪方法 | |
Li et al. | An in-loop filter based on low-complexity CNN using residuals in intra video coding | |
CN111915518B (zh) | 基于三重低秩模型的高光谱图像去噪方法 | |
CN111242854A (zh) | 一种图像去噪方法 | |
CN114820352A (zh) | 一种高光谱图像去噪方法、装置及存储介质 | |
CN110570404A (zh) | 一种基于Sobel算子的布料疵点检测方法 | |
Lv et al. | Convex blind image deconvolution with inverse filtering | |
US11915392B2 (en) | Image enhancement method and apparatus | |
CN110298803B (zh) | 一种基于压缩感知的图像条带噪声去除方法 | |
CN107644406B (zh) | 一种基于改进的正交匹配追踪的图像去噪方法 | |
CN112767261B (zh) | 一种基于广义非凸张量鲁棒主成分分析模型的针对彩色图像和视频的非局部去噪框架 | |
CN108198139B (zh) | 基于传导全变分正则化的图像去噪方法 | |
Zhang | Image motion deblurring based on salient structure selection and L0-2 norm kernel estimation |
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 |