一种融合结构张量与非局域全变分的图像去噪方法技术

技术编号:15108908 阅读:53 留言:0更新日期:2017-04-09 00:19
本发明专利技术涉及一种融合结构张量与非局域全变分的图像去噪方法,属于图像处理技术领域。该方法包括以下步骤:1)对带噪图像f进行N*N的分块,N的取值可为N=3,5,7…;2)对图像去噪建模,建立融合结构张量与非局域全变分的目标函数;3)利用分裂Bregman算法求解目标函数得到去噪图像。本发明专利技术所述方法具有有效去噪同时保持图像的纹理特征与几何结构特征的优点,算法快速,同时能够减少阶梯效应的产生。

【技术实现步骤摘要】

本专利技术属于图像处理
,涉及一种融合结构张量与非局域全变分的图像去噪方法
技术介绍
图像去噪是图像处理领域的一个基础问题,旨在将图像中的噪声去除的同时尽可能的保留原始图像信息,例如边缘、纹理和细小图像结构等。图像去噪的实质是图像退化的逆过程,经典的去噪方法有频域处理的小波法,空域处理的高斯滤波、中值滤波等,这类方法通常在某几何领域内直接对像素点进行操作,取得了较好的平滑效果,但以丢失细节、纹理信息为代价。为了缓解这种矛盾,Yaroslavsky提出邻域滤波,以像素灰度相似性作为邻域间权重依据,双边滤波器进一步加入像素间的位置信息。但是仅考虑局域信息的滤波方法没有兼顾到图像的结构信息,因此容易模糊图像细节部分。Buades等提出的非局域滤波思想,将传统的空间邻域扩展到几何意义下的结构邻域,很好的利用了自然图像的自相似性,成为非局域方法在图像去噪领域的一个里程碑。Singer等认为非局域滤波实质是图像在块空间上的扩散,Milanfar等人提出了基于高阶核函数回归的图像去噪框架,建立起Bilateral滤波、非局部均值等方法的关系。Peter等用结构相似性代替L2距离作为度量因子,并以鲁棒统计权函数代替高斯函数,取得了较好的去噪效果。近年来,基于PDE(PartialDifferentialEquation)的图像处理方法引起研究者的广泛关注,这类方法通过偏微分方程建立起图像处理模型与数学方程之间的关系,已形成完整的理论体<br>系和数值方法。其中,变分PDE通过目标函数的最优化建立PDE方程,特别适用于不适定问题的求解。目标函数由保真项与正则项构成,正则项的定义通常基于先验知识,Rudin等人提出的经典图像去噪模型就是采用全变分(Totalvariation)正则项,该正则项基于图像分片光滑假设,允许存在不连续点,在抑制高频噪声的同时保持了边缘。之后,对正则项的改进不断推进:Shu等提出了增强梯度域的稀疏性和方向性的TV-L1模型,Chan等提出了高阶全变分模型,有效抑制了“阶梯效应”并很好的刻画了纹理。
技术实现思路
有鉴于此,本专利技术的目的在于提供一种融合结构张量与非局域全变分的图像去噪方法,具体包括以下技术方案:一种融合结构张量与非局域全变分的图像去噪方法,包括以下步骤:1)对带噪图像f进行N*N的分块,N的取值可为N=3,5,7...;2)对图像去噪建模,建立融合结构张量与非局域全变分的目标函数;3)利用分裂Bregman算法求解目标函数得到去噪图像。进一步,在步骤1)中,N的取值可为3-2M+1(M为自然数)的奇数。进一步,在步骤2)中,具体包括:21)建立基于全变分正则项的图像去噪模型E(u)=minu∈Rn∫Ω12||u-f||2+λ|▿u|]]>其中f是带噪图像,u是待恢复的去噪图像,图像去噪的目的在于从f中恢复出清晰的图像u;式中第一项为数据项,第二项为平滑项,λ是平衡数据项与平滑项的权重因子,Ω为整个图像空间,为n维实数空间,▽为梯度算子,|▽u|也称为全变分正则项;22)将传统的全变分正则项扩展到非局域空间,得到非局域全变分并作为正则项约束,进而得到基于非局域信息的全变分图像去噪模型其中非局域算子▽NLu由图像空间中任意两个像素点(x,y)∈Ω×Ω定义为:▿NLu(x,y)=(u(y)-u(x))wNL(x,y)]]>其中u(x),u(y)分别指u在像素点x,y处的值,wNL:Ω×Ω为像素点x,y之间的边缘权重,且满足对称性,即wNL(x,y)=wNL(y,x);用于衡量块Px和Py之间的相似性,其中Gσ是标准偏差为σ的高斯卷积核,h为尺度参数;23)融合图像的结构张量与非局域全变分:首先,定义结构张量矩阵为:Jp(▿Iσ)=Gρ*(▿Iσ⊗▿Iσ)=ΔJ11J12J21J22]]>其中,Iσ是经过方差为σ的高斯滤波器的图像I,▽是梯度算子,是张量积运算,*代表卷积运算,Gρ是方差为ρ的高斯滤波器,J11,J12,J21,J22分别为结构张量矩阵的四个元素;求解对称半正定的结构张量矩阵Jρ的两个正交特征向量:v1=2J12J22-J11+(J22-J11)2+4J122,v1=v1||v1||,v2=v1⊥]]>与对应的特征值:λ1,2=12(J11+J22±(J22-J11)2+4J122)]]>定义局部对比度为s(x)=|λ1-λ2|,则融合结构张量与非局域的权函数wSTNL(x,y)定义为:wSTNL(x,y)=exp(-Gσ*|Px-Py|2h2)·exp(-Gσ*|Psx-Psy|2h2)]]>其中和分别代表块Px和Py的局部对比度值;利用新的权函数wSTNL(x,y)修改正则项,得到融合结构变量与非局域全变分的目标函数其中|▽STNLu|称为融合结构变量与非局域信息的全变分正则项。进一步,在步骤3)中,引入变量目标函数改写为:minu,d∫Ω12||u-f||2+λ|d|s.t.d=▿STNLu]]>采用Bregman迭代法对上式进行求解得到:(uk+1,dk+1)=argminu,d∫Ω12||u-f||2+λ|d|+β2|d-▿STNLu-bk|2bk+1=bk+dk+1-▿STNLuk+1---(1)]]>其中uk+1,dk+1,bk+1分别为变量u,d,b在第k+1次迭代的值,bk为变量b在第k次迭代的值,β为权重因子;则uk+1的最小化方程由以下最优性条件给出:(uk+1-f)-βdivSTNL(dk-▽STNLuk+1-bk)=0其中div为散度算子,dk为变量d在第k次迭代的值;采用Gauss-Seidel迭代法得到uk+1在第i点的迭代方程:uik+1,n+1=11+βΣjwij(βΣjwijuik+1,n+fi-βΣjwij(dijk+1,n-djik+1,n-bijk+1,n+bjik+1,n))---(2)]]>uk+1,n=0=uk(3)其中下标i,j分别代表第i,j个像素点,上标k代表第k次外循环,上标n代表第n次内循环,n=0代表迭代的初始化状态;变量d的迭代方程由软阈值函数给出本文档来自技高网
...

【技术保护点】
一种融合结构张量与非局域全变分的图像去噪方法,其特征在于:包括以下步骤:1)对带噪图像f进行N*N的分块,N的取值可为N=3,5,7...;2)对图像去噪建模,建立融合结构张量与非局域全变分的目标函数;3)利用分裂Bregman算法求解目标函数得到去噪图像。

【技术特征摘要】
1.一种融合结构张量与非局域全变分的图像去噪方法,其特征在于:包括以下步骤:
1)对带噪图像f进行N*N的分块,N的取值可为N=3,5,7...;
2)对图像去噪建模,建立融合结构张量与非局域全变分的目标函数;
3)利用分裂Bregman算法求解目标函数得到去噪图像。
2.根据权利要求1所述的一种融合结构张量与非局域全变分的图像去噪方法,其特征在
于:在步骤1)中,N的取值可为3-2M+1(M为自然数)的奇数。
3.根据权利要求1所述的一种融合结构张量与非局域全变分的图像去噪方法,其特征在
于:在步骤2)中,具体包括:
21)建立基于全变分正则项的图像去噪模型
其中f是带噪图像,u是待恢复的去噪图像,图像去噪的目的在于从f中恢复出清晰的图像u;
式中第一项为数据项,第二项为平滑项,λ是平衡数据项与平滑项的权重因子,Ω为整个图像
空间,为n维实数空间,为梯度算子,也称为全变分正则项;
22)将传统的全变分正则项扩展到非局域空间,得到非局域全变分并作为正则项约束,
进而得到基于非局域信息的全变分图像去噪模型
其中非局域算子由图像空间中任意两个像素点(x,y)∈Ω×Ω定义为:
▿NLu(x,y)=(u(y)-u(x))wNL(x,y)]]>其中u(x),u(y)分别指u在像素点x,y处的值,wNL:Ω×Ω为像素点x,y之间的边缘权重,且
满足对称性,即wNL(x,y)=wNL(y,x);用于衡量块Px和Py之间
的相似性,其中Gσ是标准偏差为σ的高斯卷积核,h为尺度参数;
23)融合图像的结构张量与非局域全变分:首先,定义结构张量矩阵为:
Jρ(▿Iσ)=Gρ*(▿Iσ⊗▿Iσ)=ΔJ11J12J21J22]]>其中,Iσ是经过方差为σ的高斯滤波器的图像I,是梯度算子,是张量积运算,*代表卷积
运算,Gρ是方差为ρ的高斯滤波器,J11,J12,J21,J22分别为结构张量矩阵的四个元素;
求解对称半正定的结构张量矩阵Jρ的两个正交特征向量:
v1=2J12J22-J11+(J22-J11)2+4J122,v1=v1||v1||,v2=v1⊥]]>与对应的特征值:
λ1,2=12(J11+J22±(J22-J11)2+4J122)]]>定义局部对比度为s(x)=|λ1-λ2|,则融合结构张量与非局域的权函数wSTNL(x,y)定义
为:
wSTNL(x,y)=exp(-Gσ*|Px-Py|2h2)·exp(-Gσ*|Psx-Psy|2h2)]...

【专利技术属性】
技术研发人员:王诗言
申请(专利权)人:重庆邮电大学
类型:发明
国别省市:重庆;50

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1