CN101877122A - 扩散程度可控的各向异性扩散图像去噪增强方法 - Google Patents
扩散程度可控的各向异性扩散图像去噪增强方法 Download PDFInfo
- Publication number
- CN101877122A CN101877122A CN2009102372525A CN200910237252A CN101877122A CN 101877122 A CN101877122 A CN 101877122A CN 2009102372525 A CN2009102372525 A CN 2009102372525A CN 200910237252 A CN200910237252 A CN 200910237252A CN 101877122 A CN101877122 A CN 101877122A
- Authority
- CN
- China
- Prior art keywords
- diffusion
- function
- tensor
- image
- pixel point
- 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
Landscapes
- Image Processing (AREA)
Abstract
本发明提供一种扩散程度可控的基于迹模型的各向异性扩散图像去噪增强方法,包括以下步骤:计算每个像元点的结构张量和Hessian矩阵;对每个像元点的结构张量进行特征值分解;构建每个像元点的扩散张量,使该扩散张量的特征向量是结构张量的特征向量,特征值即扩散率函数,是结构向量特征值的函数,而且该函数可通过参数调节使扩散程度可控;对迹模型迭代求解。采用该方法,可以有效对图像去噪增强。
Description
技术领域
本发明涉及图像处理技术,尤其是图像去噪技术,可用于处理有加性或乘性噪声的图像。
背景技术
图像去噪,尤其对于抑制斑点噪声的滤波算法,通常使用的有Lee滤波、Kuan滤波、Frost滤波等降噪方法,这类方法的共同之处是根据图像局部统计特征,选取合适窗口大小及调整滤波函数进行滤波,只是不同方法窗口选取准则和采取的滤波器不同。但共同的弱点是:这类方法在边缘方向上都是各向同性滤波,模糊了图像结构和细节信息。后来发展了一系列基于各向异性扩散的滤波方法,其中方法之一是Tschumperle的文章“Vector-Valued ImageRegularization with PDE’s:A Common Framework for Different Applications”(发表在国际会议IEEE International Conference on Computer Vision and Pattern Recognition(CVPR’03),Vol.1,pp.651--659,Madi son/USA,June 2003)所揭示的基于迹算子模型的各向异性滤波方法。该方法提出了基于局部几何结构的各向异性扩散模型-迹算子模型滤波方法,其中基于迹算子模型的偏微分方程是:
其中,trace(·)为矩阵的求迹运算,H为Hessian矩阵,T是扩散张量。扩散张量是驱动滤波过程的张量,它控制着局部的扩散过程,也是决定不同的各向异性扩散方法的本质变量。
已知的利用迭代方式求解迹模型的各向异性扩散方法包括以下步骤:
(a)读入二维(n=2)初始灰度图像,记为I(0)。符号I(m)表示图像序列,其中m表示迭代次数。
(b)计算图像I(m)上每个像元的一阶差分Ix,Iy;并计算对应每个像元点的2×2的对称矩阵,即结构张量S。
(c)计算I(m)每一像元点的二阶差分Ixx,Ixy,Iyy,得到每一像元点的2×2的Hessian矩阵H。
(e)使用公式(1)的离散方程I(m+1)=I(m)+c×trace(DH),从I(m)计算I(m+1)中的每一点。其中,c是预定义常数,0≤c≤1。
(f)重复步骤(b)至(e)共M次,获得结果图像。
上述步骤中Tschumperle基于平均曲率定义的扩散率函数g1(u1,u2)、g2(u1,u2)如下式:
这样的扩散率函数使得扩散滤波随着平均曲率的增加而减弱,即在局部结构强的地方,滤波减弱,而在局部结构弱的地方,滤波增强,且在梯度方向上扩散滤波减弱速度快于边缘方向。这是这类扩散率函数带来的扩散滤波的优势之处。尽管如此,由该扩散率函数生成的滤波器仍然具有不足之处:它缺少可控性和松弛性,即对于任何噪声图像,不管噪声水平如何,扩散率函数仅与平均曲率相关;而且不具备同时增强边缘等特征的能力。
发明内容
本发明针对现有技术中存在的缺陷和不足,在迹算子模型的偏微分方程(1)各向异性滤波的框架下提出一种新的用于扩散滤波保持边界的扩散率函数,该函数具有很好的可控性和松弛性,并具备能力在去噪同时增强图像。
本发明的技术方案如下:
一种基于图像局部几何结构的新的扩散率函数的各向异性扩散的图像降噪方法,其特征在于包含如下步骤:
(a)读入二维(n=2)初始灰度图像,记为I(0)。符号I(m)表示图像序列,其中m表示迭代次数。
(b)计算图像I(m)上每个像元的一阶差分Ix,Iy;并计算对应每个像元点的2×2的对称矩阵,即结构张量S。
(c)计算I(m)每一像元点的二阶差分Ixx,Ixy,Iyy,得到每一像元点的2×2的Hessian矩阵H。
(d)对每一像元点的结构张量进行矩阵的特征值分解,使得其中θ1、θ2分别是对应于特征值u1、u2的特征向量。θ1、θ2为图像灰度值变化的最大、最小方向;u1、u2则反映了特征方向上的灰度变化强弱程度。
(e)在每一像元点计算新的扩散率函数g1′(u1,u2)、g2′(u1,u2),其中u1,u2是步骤(d)中结构张量S矩阵分解得到的特征值。
(g)使用公式(1)的离散方程I(m+1)=I(m)+c×trace(DH),从I(m)计算I(m+1)中的每一点。其中,c是预定义常数,0≤c≤1。
(h)重复步骤(b)至(g)共M次,获得结果图像。
所述的矩阵分解是指图像上每个像元点对应的结构矩阵的特征值分解。
所述的扩散张量是一个2阶矩阵,它直接利用矩阵特征值分解的形式构建,构建的特征值与特征向量具有以下特点:
(1)扩散张量T的特征向量与局部结构张量的特征向量相同,即扩散张量的特征向量选取局部结构张量的特征向量θ1,θ2;
(2)扩散张量D的特征值分别为扩散率函数g1′(u1,u2)、g2′(u1,u2),分别表征在梯度方向和轮廓方向上的扩散滤波程度。
所述的新的扩散率函数具有形式:
其中α,β,a,b1,b2是五个控制参数,均是大于零的实数,可依据图像噪声特点选择。
所述的预定义常数c在区间[0,1],预先给出。
所述的迭代次数M预先给出。
本发明与现有技术相比有以下特点:本发明的扩散率函数有多个自由参数,这些参数从不同方面控制扩散滤波过程,使得可依据图像局部几何结构和噪声水平
附图说明
图1为迭代实现的基于偏微分方程图像扩散滤波的基本流程图
图2为基于迹算子模型的各向异性扩散滤波的基本流程图
图3分(a)、(b)、(c)图分别是控制参数a、b、β分别变化时扩散率函数变化的示意图
图4本发明去噪结果示例:(a)原噪声图像(含加性白噪声);(b)Lee滤波去噪结果;(c)本方法迭代20次去噪结果;(d)本方法迭代40次去噪结果
具体实施方式
下面结合附图通过实施例对本发明作进一步的详细说明。
图1是迭代实现的基于偏微分方程图像扩散滤波的基本流程图。迭代次数是M,是预先指定的。当扩散滤波次数小于指定次数M时,将滤波结果图像作为当前图像再次滤波,否则滤波结果图像作为最终滤波结果。其中滤波单元100是本发明基于偏微分方程扩散滤波的主要部分。
图2是本发明实现的基于迹算子模型的各向异性扩散滤波的流程图。图2中读入的图像经过下列6个单元处理后获得结果。
单元110对读入的图像计算对应每个像元的结构张量S。图像上每个像元点对应的结构张量是2阶对称矩阵,它是梯度算子的外积与某个低通核函数的卷积,梯度算子的外积具有形式其中,Ix、Iy是一阶微分或差分算子。卷积运算在降低结构张量噪声影响的同时揭示图像的局部几何结构。合适的低通核函数的例子是高斯核函数。2维的高斯核函数的形式是
单元130对每一像元点的结构张量S进行特征值分解,使得其中θ1、θ2分别是对应于特征值u1、u2的特征向量。θ1、θ2为图像灰度值变化的最大、最小方向;u1、u2则反映了特征方向上的灰度变化强弱程度。
单元140按照公式(3)计算新的扩散率函数g1′(u1,u2)、g2′(u1,u2)。新的扩散率函数中有五个调控参数α,β,a,b1,b2。这五个参数分别从不同方面控制扩散过程:α控制着平坦区域、边缘方向扩散滤波程度和梯度方向上最大的扩散滤波程度。在平坦区u1≈u2=0;g1≈g2=α这时梯度方向和边缘方向上的扩散滤波程度相同,相当于对平坦区域进行高斯滤波;对于完全结构区域u1<u2,g1≈0,g2≈α,这时在梯度方向不滤波,只进行边缘方向的滤波,相当于平均曲率运动扩散;对于不完全结构区域,a,b,β对扩散过程的控制可以由附图2的函数y(x)=1-xb/(xb+a+βb)的图示直观看出。附图3(a)是参数a,β固定b变化时的曲线的变化情况;附图2(b)是参数b,β固定a变化时的曲线的变化情况;附图2(c)是参数a,b固定β变化时的曲线的变化情况。由图2可见,参数b控制着扩散率函数从哪开始弯曲下降,作为一个软域值控制噪声信号区域的转换,控制转换的速度,b越大函数下降越快,随平均曲率增加,越快转换,滤波程度越快减弱;参数β控制从何处扩散率函数接近0;参数a控制着该扩散率函数最小值,即扩散率是否会过冲,给出小于0的值,扩散率小于零可增加边界对比度,起到图像增强的效果。实验中,令b1>b2使扩散率函数在梯度方向下降快于边缘方向,并选择参数控制不完全结构区域的滤波。
单元150对每一像元点重构新的扩散张量D,其中D的特征向量等于S的特征向量,D的特征值为公式(3)定义的新的扩散率函数g1′(u1,u2)、g2′(u1,u2)。
单元160使用公式(3)的离散方程I(m+1)=I(m)+c*trace(DH),从I(m)计算I(m+1)中的每一点。其中,c是预定义常数,0≤c≤1。
附图4是本发明应用的一个例子。图3(a)是添加了噪声的棋盘图;图3(b)是经典的Lee滤波的结果,显然去噪的同时模糊了边缘;图3(c)是本发明迭代20次处理的结果;图3(d)是本发明迭代40次处理的结果。图3(c)和图3(d)扩散滤波时选择的相关参数为:a=0.1,b1=3,b2=1,α=0.8,β=0.1,迭代步长c=0.05。图3结果表明对比Lee滤波结果,本发明方法在有效去噪的同时,边缘和细节信息保持得更充分。
Claims (5)
1.一种扩散程度可控的基于迹模型的各向异性扩散图像去噪增强方法,其特征在于包含以下步骤:计算每个像元点的结构张量和Hessian矩阵;对每个像元点的结构张量进行特征值分解;构建每个像元点的扩散张量,使该扩散张量的特征向量是结构张量的特征向量,特征值即扩散率函数,是结构向量特征值的函数,而且该函数可通过参数调节使扩散程度可控;最后对迹模型迭代求解,获得去噪增强的结果。
2.根据权利要求1所述的一种计算图像上每个像元点结构张量的方法,其特征在于:所述结构张量是该点梯度算子外积和一个合适的低通滤波核函数的卷积,合适的低通滤波核函数的一个例子是高斯核函数。
3.根据权利要求1所述的一种构建图像上每个像元点扩散张量的方法,其特征在于:使扩散张量的特征向量是结构张量的特征向量,特征值g′1(u1,u2)、g′2(u1,u2)是结构向量特征值u1、u2的函数,即扩散率函数,该扩散率函数可通过参数调节使扩散程度可控。
4.根据权利要求3所述的一种带有参数的扩散率函数,其特征在于:该函数具有形式
其中α,β,a,b1,b2是五个控制参数,均是大于零的实数,可以通过改变参数的大小改变函数的形状,从而控制扩散程度。
5.根据权利要求1所述的迹模型的求解方法,其特征在于:离散化迹模型方程后迭代求解,迭代次数预先给出。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 200910237252 CN101877122B (zh) | 2009-11-12 | 2009-11-12 | 扩散程度可控的各向异性扩散图像去噪增强方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 200910237252 CN101877122B (zh) | 2009-11-12 | 2009-11-12 | 扩散程度可控的各向异性扩散图像去噪增强方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101877122A true CN101877122A (zh) | 2010-11-03 |
CN101877122B CN101877122B (zh) | 2012-12-05 |
Family
ID=43019667
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 200910237252 Expired - Fee Related CN101877122B (zh) | 2009-11-12 | 2009-11-12 | 扩散程度可控的各向异性扩散图像去噪增强方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101877122B (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102095370A (zh) * | 2010-11-22 | 2011-06-15 | 北京航空航天大学 | 三x组合标记的检测识别方法 |
CN102572201A (zh) * | 2010-12-31 | 2012-07-11 | 北京大学 | 一种图像网纹去除方法及系统 |
CN103700089A (zh) * | 2013-12-01 | 2014-04-02 | 北京航空航天大学 | 一种三维医学图像多尺度异构特征的提取与分类方法 |
CN104657951A (zh) * | 2015-03-02 | 2015-05-27 | 桂林电子科技大学 | 图像乘性噪声移除方法 |
CN104700372A (zh) * | 2015-03-19 | 2015-06-10 | 天津大学 | 基于边缘保护的加权各向异性扩散滤波方法 |
CN105678770A (zh) * | 2016-01-07 | 2016-06-15 | 潘燕 | 一种轮廓识别滤波性能良好的墙体裂缝检测装置 |
CN107203976A (zh) * | 2017-04-19 | 2017-09-26 | 武汉科技大学 | 一种基于噪点检测的自适应非局部均值去噪方法及系统 |
CN108737686A (zh) * | 2018-06-11 | 2018-11-02 | 昆明理工大学 | 基于图像系数变更的边缘增强误差扩散数字半色调方法 |
CN108765332A (zh) * | 2018-05-23 | 2018-11-06 | 成都信息工程大学 | 一种椭圆搜索窗口和参数自适应的非局部均值去噪方法 |
CN110188614A (zh) * | 2019-04-30 | 2019-08-30 | 杭州电子科技大学 | 一种基于皮裂纹分割的nlm滤波指静脉去噪方法 |
US11593918B1 (en) | 2017-05-16 | 2023-02-28 | Apple Inc. | Gradient-based noise reduction |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1234091C (zh) * | 2002-07-24 | 2005-12-28 | 中国科学院自动化研究所 | 基于水平集和分水岭方法的医学图像分割方法 |
CN101493933B (zh) * | 2009-03-03 | 2010-12-08 | 北京科技大学 | 一种局部结构自适应的图像扩散去噪方法 |
CN101540042B (zh) * | 2009-04-24 | 2011-06-01 | 西安电子科技大学 | 基于第二代曲线波变换的sar图像相干斑抑制方法 |
-
2009
- 2009-11-12 CN CN 200910237252 patent/CN101877122B/zh not_active Expired - Fee Related
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102095370A (zh) * | 2010-11-22 | 2011-06-15 | 北京航空航天大学 | 三x组合标记的检测识别方法 |
CN102572201A (zh) * | 2010-12-31 | 2012-07-11 | 北京大学 | 一种图像网纹去除方法及系统 |
CN102572201B (zh) * | 2010-12-31 | 2015-01-28 | 北京大学 | 一种图像网纹去除方法及系统 |
CN103700089A (zh) * | 2013-12-01 | 2014-04-02 | 北京航空航天大学 | 一种三维医学图像多尺度异构特征的提取与分类方法 |
CN104657951A (zh) * | 2015-03-02 | 2015-05-27 | 桂林电子科技大学 | 图像乘性噪声移除方法 |
CN104700372A (zh) * | 2015-03-19 | 2015-06-10 | 天津大学 | 基于边缘保护的加权各向异性扩散滤波方法 |
CN104700372B (zh) * | 2015-03-19 | 2017-08-04 | 天津大学 | 基于边缘保护的加权各向异性扩散滤波方法 |
CN105678770A (zh) * | 2016-01-07 | 2016-06-15 | 潘燕 | 一种轮廓识别滤波性能良好的墙体裂缝检测装置 |
CN107203976A (zh) * | 2017-04-19 | 2017-09-26 | 武汉科技大学 | 一种基于噪点检测的自适应非局部均值去噪方法及系统 |
CN107203976B (zh) * | 2017-04-19 | 2019-07-23 | 武汉科技大学 | 一种基于噪点检测的自适应非局部均值去噪方法及系统 |
US11593918B1 (en) | 2017-05-16 | 2023-02-28 | Apple Inc. | Gradient-based noise reduction |
CN108765332A (zh) * | 2018-05-23 | 2018-11-06 | 成都信息工程大学 | 一种椭圆搜索窗口和参数自适应的非局部均值去噪方法 |
CN108765332B (zh) * | 2018-05-23 | 2022-05-06 | 成都信息工程大学 | 一种椭圆搜索窗口和参数自适应的非局部均值去噪方法 |
CN108737686A (zh) * | 2018-06-11 | 2018-11-02 | 昆明理工大学 | 基于图像系数变更的边缘增强误差扩散数字半色调方法 |
CN110188614A (zh) * | 2019-04-30 | 2019-08-30 | 杭州电子科技大学 | 一种基于皮裂纹分割的nlm滤波指静脉去噪方法 |
CN110188614B (zh) * | 2019-04-30 | 2021-03-30 | 杭州电子科技大学 | 一种基于皮裂纹分割的nlm滤波指静脉去噪方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101877122B (zh) | 2012-12-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101877122B (zh) | 扩散程度可控的各向异性扩散图像去噪增强方法 | |
CN110599409B (zh) | 基于多尺度卷积组与并行的卷积神经网络图像去噪方法 | |
Wang et al. | Dehazing for images with large sky region | |
Xu et al. | Deep edge-aware filters | |
CN103093433B (zh) | 基于区域划分和字典学习的自然图像去噪方法 | |
CN104715461B (zh) | 图像去噪方法 | |
CN103854262B (zh) | 基于结构聚类与稀疏字典学习的医学图像降噪方法 | |
CN102073999B (zh) | 基于双冗余字典学习的自然图像去噪方法 | |
CN102156975B (zh) | 基于支撑值变换和多尺度冗余字典学习的自然图像去噪方法 | |
CN105374007A (zh) | 融合骨架笔画和纹理特征的铅笔画生成方法和装置 | |
CN110136075B (zh) | 一种基于边缘锐化循环生成对抗网络的遥感影像去雾方法 | |
CN104285239A (zh) | 图像处理装置、图像处理方法、程序、印刷介质以及记录介质 | |
US8731318B2 (en) | Unified spatial image processing | |
CN111861886B (zh) | 一种基于多尺度反馈网络的图像超分辨率重建方法 | |
CN110782406B (zh) | 一种基于信息蒸馏网络的图像去噪方法及装置 | |
CN110503610A (zh) | 一种基于gan网络的图像雨雪痕迹去除方法 | |
CN105427262A (zh) | 基于双向增强扩散滤波的图像去噪方法 | |
CN103049891A (zh) | 基于自适应窗口选择的视频图像去模糊方法 | |
CN114723630A (zh) | 基于空洞双残差多尺度深度网络的图像去模糊方法及系统 | |
CN104616259B (zh) | 一种噪声强度自适应的非局部均值图像去噪方法 | |
CN109993701B (zh) | 一种基于金字塔结构的深度图超分辨率重建的方法 | |
CN104200438A (zh) | 红外图像多级细节增强处理方法及其处理装置 | |
Malik et al. | Image denoising by super neurons: Why go deep? | |
CN104299203A (zh) | 基于双树复小波变换的图像去噪方法 | |
CN114897884A (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 | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20121205 Termination date: 20131112 |