CN104463811B - 基于能量泛函的图像平滑与锐化方法 - Google Patents
基于能量泛函的图像平滑与锐化方法 Download PDFInfo
- Publication number
- CN104463811B CN104463811B CN201410835528.0A CN201410835528A CN104463811B CN 104463811 B CN104463811 B CN 104463811B CN 201410835528 A CN201410835528 A CN 201410835528A CN 104463811 B CN104463811 B CN 104463811B
- Authority
- CN
- China
- Prior art keywords
- image
- energy functional
- diffusion
- formula
- algorithm
- 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.)
- Expired - Fee Related
Links
Landscapes
- Image Processing (AREA)
Abstract
本发明涉及一种基于能量泛函的图像平滑与锐化方法,包括以下步骤:(1)将各向异性扩散方程转化为最小化能量泛函;(2)将动态方程以内在坐标形式表示;(3)建立能量泛函(4)建立图像的梯度阈值函数k=e‑αt;(5)引入保真项
Description
技术领域
本发明涉及图像去噪算法领域,尤其是基于能量泛函的图像去噪方法。
背景技术
图像噪声的主要来源是图像在采集过程中的随机高斯噪声和图像传播过程中的椒盐噪声。传统的去噪方法有中值滤波、同态滤波、逆滤波等,这些方法在一定程度上可以达到去除噪声的目的,但它们有一个共同的弱点,在去噪的同时,也会使图像的边缘模糊化,甚至使图像的细节纹理信息丢失。相比于传统方法,基于偏微分方程的图像处理方法具有更强的局部自适应能力和更高的灵活性,在图像的去噪、分割、边缘检测、增强等方面都有重要应用。
基于偏微分方程的图像处理方法是在图像的连续数学模型基础上,令图像遵循某一指定的偏微分方程发生变化,基于变分模型最具代表性的算法是非线性偏微分方程去噪算法(TV算法),该算法将图像归类为变分有界函数空间,采用一次范数全变分作为其“平滑性”的度量,并沿着梯度垂直的方向进行约束,很好地保护图像的边缘,但是由于它并不完全符合图像处理的形态学原则,迭代多次后通常会产生阶梯效应。
发明内容
针对上述现有技术中的缺陷,本发明提供一种基于能量泛函的图像平滑与锐化方法,其应用于图像去噪时效性高、复杂度低、精确度高。
本发明采用以下技术方案:基于能量泛函的图像平滑与锐化方法,其特征在于,包括以下步骤:
步骤一、建立噪声图像模型I0=KI,为使去噪后的图像与原始图像充分接近,建立所述噪声图像模型的最小化模型其动态方程为引入Neumann边界,该动态方程等价于各向异性扩散方程即各向异性扩散方程转化为最小能量泛函,其中I是无噪图像,I0是降质图像,K是降质算子,Ω是图像区域;
步骤二、将步骤一中的动态方程以内在坐标形式表示为
步骤三、建立能量泛函式1)表示为
取代入式2)得
在图像边缘时,沿梯度方向逆扩散;在图像平坦区域时,沿边缘方向和梯度方向同时进行平滑处理,其中f(·)≥0,k是图像I的梯度阈值;
步骤四、建立图像的梯度阈值函数k=e-αt,其为随扩散时间和扩散次数变化的函数,其中t=n,t为扩散时间,n为扩散次数,max{W,H}为图像的最大灰度值;
步骤五、引入保真项建立基于能量泛函的图像平滑与锐化算法
其中λ是调整参数,在图像的边缘,在图像的平坦区域,则扩散程度达到最大;
步骤六、用中心差分数值算法对步骤五的结果进一步处理。
所述步骤六的中心差分数值算法具体为,令
式4)的离散化形式为
式中,n=0,1,2,······,表示分解尺度,则式4)的离散形式为
本发明达到的有益效果:在复杂度方面,利用内在坐标形式进行性能分析,得到扩散滤波算法,算法需要的信息量少,方法简单,实现图像去噪,使受污染图像更接近原始图像;在时效性方面,因为需要的信息量少,实施的复杂度低,从而降低了算法的处理时间;在去噪性能方面,图像的峰值信噪比大幅提高,受噪声污染的图像经本算法处理后更加接近原始图像。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
第一步:建立噪声图像模型I0=KI,为了使去噪后的图像与原始图像充分接近,现建立最小化模型求该式的拉格朗日方程,建立动态方程引入Neumann边界,从而动态方程等价于各向异性扩散方程这样从最小化问题出发得到了各向异性扩散方程,即各向异性扩散方程可以归结为最小能量泛函,其中I是无噪图像,I0是降质图像,K是降质算子,Ω是图像区域;
第二步:建立第一步中的内在坐标形式表示
由于第一步中的最小化模型有解,则故沿边缘方向的扩散系数在扩散过程中沿边缘方向正扩散,模糊图像;而梯度方向的扩散系数其正负号不确定,当时,沿梯度方向正扩散,模糊图像,当时,沿梯度方向逆扩散,锐化图像的边缘,能很好的保护图像的边缘。其中(η,ξ)为内在坐标系,η为图像的梯度方向,即垂直于图像特征(边缘)的方向;ξ为垂直于梯度的方向,即沿图像特征(边缘)的方向,
第三步:建立能量泛函由第一步解该式的最小能量泛函,建立数学模型由第二步建立数学模型的内在坐标形式现取代入数学模型中根据第二步,当在图像边缘时, 沿梯度方向逆扩散,锐化图像的边缘,改善了边缘模糊的现象,能很好的保护图像的边缘;当在图像平坦区域时,在平坦区域,沿着边缘方向和切线方向同时进行平滑,转化为各向同性的扩散方程,因此具有较强的去噪能力。其中Ω为图像区域,f(·)≥0,f'(·)>0,k是图像I的梯度阈值,下同;
第四步:建立图像的梯度阈值函数k。根据第三步本发明的整个扩散过程,实质上是通过图像的梯度与阈值大小进行比较实现的,故将梯度阈值k设计为一个随扩散时间和扩散次数变化的函数k=e-αt,其中扩散次数代表扩散时间,即t=n,max{W,H}为图像的最大灰度值。随着扩散次数增加,即扩散时间增加,k不断减小,即梯度阈值不断减小,这样就使得下一次的边缘判断更加准确,从而保留了更多的边缘信息;
第五步:建立保真项使图像不失真。其中λ是调整参数,一般取较小的值,有重要作用,在图像的边缘,此时 这样滤波结果会尽可能地接近原始图像,进一步加强了保边缘的效果;在图像的平坦区域,此时则扩散程度达到最大,尽可能的去除噪声。在此基础上,建立基于能量泛函的图像平滑与锐化算法
第六步:根据第五步得到的基于能量泛函的图像平滑与锐化算法,用中心差分数值算法进行数值计算。
为了实现本发明算法,我们采用有限差分法,令
式中,k为图像的梯度阈值。现采用中心差分法,则本发明算法的离散化形式为
式中,n=0,1,2,······,表示时间水平。其中
最终本发明算法的离散形式为
综上所述,本发明专利所述基于能量泛函的图像平滑与锐化算法复杂度低,时效性高去噪后的图像峰值信噪比与PM算法相比提高了15个dB左右,与TV算法相比提高了8个dB左右,能够增强与锐化图像中的边缘纹理等重要信息。
用受高斯随机噪声(σ=20)汚染的图进行仿真实验,把本发明算法与经典的PM算法与TV算法进行比较,由表1和表2,可以看出,本发明算法较PM算法的峰值信噪比提高了15个dB左右,较TV算法提高了8个dB左右,进一步证实了本发明算法的有效性与合理性。
表1各去噪算法的MSE和PSNR比较
表2各去噪算法不同噪声方差的MSE和PSNR比较
以上是本发明的较佳实施方式,但本发明的保护范围不限于此。任何熟悉本领域的技术人员在本发明所揭露的技术范围内,未经创造性劳动想到的变换或替换,都应涵盖在本发明的保护范围之内。因此本发明的保护范围应以权利要求所限定的保护范围为准。
Claims (2)
1.基于能量泛函的图像平滑与锐化方法,其特征在于,包括以下步骤:
步骤一、建立噪声图像模型I0=KI,为使去噪后的图像与原始图像充分接近,建立所述噪声图像模型的最小化模型其动态方程为引入Neumann边界,该动态方程等价于各向异性扩散方程即各向异性扩散方程转化为最小能量泛函,其中I是无噪图像,I0是降质图像,K是降质算子,Ω是图像区域;
步骤二、将步骤一中的动态方程以内在坐标形式表示为
步骤三、建立能量泛函式1)表示为
取代入式2)得
在图像边缘时,|▽I|>>k,沿梯度方向逆扩散;在图像平坦区域时,|▽I|<<k,沿边缘方向和梯度方向同时进行平滑处理,其中f(·)≥0,k是图像I的梯度阈值;
步骤四、建立图像的梯度阈值函数k=e-αt,其为随扩散时间和扩散次数变化的函数,其中t=n,t为扩散时间,n为扩散次数,max{W,H}为图像的最大灰度值;
步骤五、引入保真项建立基于能量泛函的图像平滑与锐化算法
其中λ是调整参数,在图像的边缘,|▽I|→∞,在图像的平坦区域,|▽I|→0,则扩散程度达到最大;
步骤六、用中心差分数值算法对步骤五的结果进一步处理。
2.根据权利要求1所述的基于能量泛函的图像平滑与锐化方法,其特征在于,所述步骤六的中心差分数值算法具体为,令
式4)的离散化形式为
式中,n=0,1,2,......,表示分解尺度,则式4)的离散形式为
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410835528.0A CN104463811B (zh) | 2014-12-29 | 2014-12-29 | 基于能量泛函的图像平滑与锐化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410835528.0A CN104463811B (zh) | 2014-12-29 | 2014-12-29 | 基于能量泛函的图像平滑与锐化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104463811A CN104463811A (zh) | 2015-03-25 |
CN104463811B true CN104463811B (zh) | 2017-08-04 |
Family
ID=52909799
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410835528.0A Expired - Fee Related CN104463811B (zh) | 2014-12-29 | 2014-12-29 | 基于能量泛函的图像平滑与锐化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104463811B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104851082B (zh) * | 2015-05-21 | 2017-10-27 | 南京信息工程大学 | 一种基于梯度与拉普拉斯算子的混合自适应图像去噪方法 |
CN105427262A (zh) * | 2015-12-15 | 2016-03-23 | 南京信息工程大学 | 基于双向增强扩散滤波的图像去噪方法 |
CN105913382B (zh) * | 2016-03-01 | 2018-11-23 | 南京信息工程大学 | 阈值寻优的高保真各向异性滤波方法 |
CN108428216B (zh) * | 2018-01-16 | 2021-09-24 | 辽宁师范大学 | 基于散布矩阵特征的二阶偏微分方程遥感图像去噪方法 |
CN110111261B (zh) * | 2019-03-28 | 2021-05-28 | 瑞芯微电子股份有限公司 | 图像的自适应平衡处理方法、电子设备和计算机可读存储介质 |
CN111435436B (zh) * | 2019-12-13 | 2021-01-08 | 珠海大横琴科技发展有限公司 | 一种基于目标位置的周界防入侵方法和装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101807292A (zh) * | 2010-01-07 | 2010-08-18 | 清华大学 | 一种图像去噪方法 |
CN101887576A (zh) * | 2010-06-04 | 2010-11-17 | 西安电子科技大学 | 基于偏微分方程滤波器的图像去噪方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7747045B2 (en) * | 2006-06-30 | 2010-06-29 | Fujifilm Corporation | Method and apparatus for diffusion based illumination normalization |
-
2014
- 2014-12-29 CN CN201410835528.0A patent/CN104463811B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101807292A (zh) * | 2010-01-07 | 2010-08-18 | 清华大学 | 一种图像去噪方法 |
CN101887576A (zh) * | 2010-06-04 | 2010-11-17 | 西安电子科技大学 | 基于偏微分方程滤波器的图像去噪方法 |
Non-Patent Citations (5)
Title |
---|
AN ADAPTIVE DIFFUSION SCHEME FOR IMAGE RESTORATION AND SELECTIVE SMOOTHING;V. R. SIJRYA PRASATH等;《International Journal of Image and Graphics》;20121231;第12卷(第1期);第1250003-1-125003-18页 * |
Modeling Textures with Total Variation Minimization and Oscillating Patterns in Image Processing;Luminita A. Vese等;《Journal of Scientific Computing》;20031231;第19卷(第1-3期);第553-572页 * |
一种改进的各向异性扩散去噪模型;赵海勇等;《计算机科学》;20131130;第40卷(第11A期);第147-149、169页 * |
基于偏微分方程的SAR图像去噪方法研究;王新楼;《中国优秀博硕士学位论文全文数据库(硕士)信息科技辑》;20050315;第18-52页 * |
基于自适应梯度阈值各向异性滤波抑制红外复杂背景;孙海江等;《光学 精密工程》;20140131;第22卷(第1期);第146-151页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104463811A (zh) | 2015-03-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104463811B (zh) | 基于能量泛函的图像平滑与锐化方法 | |
Jain et al. | A survey of edge-preserving image denoising methods | |
CN103116875B (zh) | 自适应双边滤波图像去噪方法 | |
Fang et al. | Single image dehazing and denoising: a fast variational approach | |
CN106127688B (zh) | 一种超分辨率图像重建方法及其系统 | |
CN104331869B (zh) | 梯度与曲率相结合的图像平滑方法 | |
CN104392418B (zh) | 基于各向异性扩散的图像平滑改进算法 | |
Xu et al. | A two-level domain decomposition method for image restoration | |
CN105427262A (zh) | 基于双向增强扩散滤波的图像去噪方法 | |
CN102750679B (zh) | 图像质量评估的盲去模糊方法 | |
CN107169932A (zh) | 一种适用于中子成像系统图像的基于高斯‑泊松混合噪声模型的图像复原方法 | |
CN105303538B (zh) | 一种基于nsct和pca的高斯噪声方差估计方法 | |
CN109816596A (zh) | 一种图像去噪方法及装置 | |
CN105243647A (zh) | 一种基于线性空间滤波的图像增强方法 | |
Yu et al. | Image denoising based on adaptive fractional order with improved PM model | |
CN104537613A (zh) | 一种用于改善图像视觉效果的分数阶I-divergence方法 | |
CN104463810B (zh) | 基于tv流的自适应扩散滤波图像去噪方法 | |
CN103679648B (zh) | 一种基于空间分割的矩匹配卫星影像条带噪声去除方法 | |
Chen et al. | Inexact alternating direction method based on proximity projection operator for image inpainting in wavelet domain | |
Youssif et al. | Adaptive algorithm for image denoising based on curvelet threshold | |
Das et al. | Histogram equalization techniques for contrast enhancement: a review | |
CN106127692A (zh) | 可用于灰度和彩色图像的双边回归滤波方法 | |
CN113139920B (zh) | 一种古籍图像修复方法、终端设备及存储介质 | |
Wang et al. | A total variation model based on edge adaptive guiding function for remote sensing image de-noising | |
Wang et al. | Multiplicative noise removal using deep CNN denoiser prior |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CP02 | Change in the address of a patent holder |
Address after: 210044 No. 219 Ningliu Road, Jiangbei New District, Nanjing City, Jiangsu Province Patentee after: Nanjing University of Information Science and Technology Address before: Zhongshan road Wuzhong District Mudu town of Suzhou city in Jiangsu province 215101 No. 70 Wuzhong Science Park Building 2 room 2310 Patentee before: Nanjing University of Information Science and Technology |
|
CP02 | Change in the address of a patent holder | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170804 Termination date: 20211229 |
|
CF01 | Termination of patent right due to non-payment of annual fee |