一种融合结构张量与非局域全变分的图像去噪方法
【技术领域】
[0001] 本发明属于图像处理技术领域,涉及一种融合结构张量与非局域全变分的图像去 噪方法。
【背景技术】
[0002] 图像去噪是图像处理领域的一个基础问题,旨在将图像中的噪声去除的同时尽可 能的保留原始图像信息,例如边缘、纹理和细小图像结构等。图像去噪的实质是图像退化的 逆过程,经典的去噪方法有频域处理的小波法,空域处理的高斯滤波、中值滤波等,这类方 法通常在某几何领域内直接对像素点进行操作,取得了较好的平滑效果,但以丢失细节、纹 理信息为代价。为了缓解这种矛盾,Yaroslavsky提出邻域滤波,以像素灰度相似性作为邻 域间权重依据,双边滤波器进一步加入像素间的位置信息。但是仅考虑局域信息的滤波方 法没有兼顾到图像的结构信息,因此容易模糊图像细节部分。Buades等提出的非局域滤波 思想,将传统的空间邻域扩展到几何意义下的结构邻域,很好的利用了自然图像的自相似 性,成为非局域方法在图像去噪领域的一个里程碑。Singer等认为非局域滤波实质是图像 在块空间上的扩散,Milanfar等人提出了基于高阶核函数回归的图像去噪框架,建立起 Bilateral滤波、非局部均值等方法的关系。Peter等用结构相似性代替L2距离作为度量因 子,并以鲁棒统计权函数代替高斯函数,取得了较好的去噪效果。
[0003] 近年来,基于FOECPartial Differential Equation)的图像处理方法引起研究者 的广泛关注,这类方法通过偏微分方程建立起图像处理模型与数学方程之间的关系,已形 成完整的理论体系和数值方法。其中,变分toe通过目标函数的最优化建立toe方程,特别适 用于不适定问题的求解。目标函数由保真项与正则项构成,正则项的定义通常基于先验知 识,Rud i η等人提出的经典图像去噪模型就是采用全变分(To ta 1 var i at i on)正则项,该正 则项基于图像分片光滑假设,允许存在不连续点,在抑制高频噪声的同时保持了边缘。之 后,对正则项的改进不断推进:Shu等提出了增强梯度域的稀疏性和方向性的TV-L1模型, Chan等提出了高阶全变分模型,有效抑制了"阶梯效应"并很好的刻画了纹理。
【发明内容】
[0004] 有鉴于此,本发明的目的在于提供一种融合结构张量与非局域全变分的图像去噪 方法,具体包括以下技术方案:
[0005] -种融合结构张量与非局域全变分的图像去噪方法,包括以下步骤:
[0006] 1)对带噪图像f进行N*N的分块,N的取值可为N=3,5,7...;
[0007] 2)对图像去噪建模,建立融合结构张量与非局域全变分的目标函数;
[0008] 3)利用分裂Bregman算法求解目标函数得到去噪图像。
[0009] 进一步,在步骤1)中,N的取值可为3-2M+1 (M为自然数)的奇数。
[0010]进一步,在步骤2)中,具体包括:
[0011] 21)建立基于全变分正则项的图像去噪模型
[0012]
[0013] 其中f是带噪图像,u是待恢复的去噪图像,图像去噪的目的在于从f中恢复出清晰 的图像u;式中第一项为数据项,第二项为平滑项,λ是平衡数据项与平滑项的权重因子,Ω 为整个图像空间,;为η维实数空间,▽为梯度算子,| Vu |也称为全变分正则项;
[0014] 22)将传统的全变分正则项扩展到非局域空间,得到非局域全变分并作为正则项 约束,进而得到基于非局域信息的全变分图像去噪模型
[0015]
[0016] 其中非局域算子VnO!由图像空间中任意两个像素点(χ,5〇ΕΩ X Ω定义为:
[0017]
[0018] 其中u(x),u(y)分别指u在像素点x,y处的值,wnl: Ω X Ω为像素点x,y之间的边缘 权重,且满足对称性,即》1(1,7)=?^(7,1);
目于衡量块Px和 Py之间的相似性,其中G。是标准偏差为σ的高斯卷积核,h为尺度参数;
[0019] 23)融合图像的结构张量与非局域全变分:首先,定义结构张量矩阵为:
[0020]
[0021]其中,I。是经过方差为〇的高斯滤波器的图像I,▽是梯度算子,0是张量积运算,* 代表卷积运算,GP是方差为Ρ的高斯滤波器,如,如,加,加分别为结构张量矩阵的四个元 素;
[0022]求解对称半正定的结构张量矩阵JP的两个正交特征向量:
[0023]
[0024]与对应的特征值:
[0025]
[0026]定义局部对比度为s(x)= Ιλχ-λ〗I,则融合结构张量与非局域的权函数WSTNL(X,y) 定义为:
[0027]
[0028] 其中、和分别代表块Px和Py的局部对比度值;
[0029]利用新的权函数wSTNL(X,y)修改正则项,得到融合结构变量与非局域全变分的目 标函数
[0030]
[0031 ]其中I VstnOI I称为融合结构变量与非局域信息的全变分正则项。
[0032] 进一步,在步骤3)中,引入变量么-Cl,目标函数改写为:
[0033]
[0034] 米用Bregman迗代'/云对上式进仃豕觯得到:
[0035]
(1):
[0036] 其中uk+1,dk+1,bk+1分别为变量u,d,b在第k+Ι次迭代的值,b k为变量b在第k次迭代 的值,β为权重因子;则uk+1的最小化方程由以下最优性条件给出:
[0037] (uk+1 -f) -β(1 i vstnl (dk- V sTNLUk+1 -bk) = 0
[0038] 其中div为散度算子,dk为变量d在第k次迭代的值;
[0039] 采用Gauss-Seidel迭代法得到uk+1在第i点的迭代方程:
^ ρΛ >1 Λ ? (2)
[0041 ] uk+1>n=0 = uk (3)
[0042] 其中下标i,j分别代表第i,j个像素点,上标k代表第k次外循环,上标η代表第η次 内循环,η = 0代表迭代的初始化状态;
[0043] 变量d的迭代方程由软阈值函数给出:
[0044]
[0045] 变量d的离散化表达式:
[00461
(4)
[0047] 其中下标i,j分别代表第i,j个像素点,上标k代表第k次外循环;
[0048] 由式(1)得到变量b的迭代方稈为:
[0049]
(5):
[0050]联合式(2-5),得到模型求解的变量迭代方程组;迭代停止后,去噪图像为u = uR'N, 其中K为外循环最大迭代次数,N为内循环最大迭代次数。
[0051]本发明的有益效果在于:本发明在基于全变分的图像去噪模型中引入非局域描述 子,可以较好捕捉图像的细小纹理,同时考虑图像的结构信息,提出一种融合结构张量与非 局域全变分的图像去噪模型。该模型一方面考虑了图像的非局域特性,通过度量patch间的 灰度相似性,能够较好的保持图像的纹理特性;另一方面兼顾了几何结构特性,能够调节不 同对比度patch之间权重,尤其是提高了低对比度区域的权重值,同时较少阶梯效应的产 生。模型求解采用分裂Bregman算法,克服了传统梯度下降法中存在的非线性和不可微性带 来的最小化速度缓慢问题,本发明中采用的分裂Bregman算法快速有效。
【附图说明】
[0052]为了使本发明的目的、技术方案和有益效果更加清楚,本发明提供如下附图进行 说明:
[0053]图1为本发明所述方法的流程示意图。
【具体实施方式】
[0054]下面将结合附图,对本发明的优选实施例进行详细的描述。
[0055] 图1为本发明所述方法的流程示意图,如图所示,本方法包括以下步骤:1)对带噪 图像f进行N*N的分块,N的取值可为N=3,5,7...,在本实施例中为3-2M+UM为自然数)的奇 数;2)对图像去噪建模,建立融合结构张量与非局域全变分的目标函数;3)利用分裂 Bregman算法求解目标函数得到去噪图像。
[0056] 本实施例中,步骤2)中建立目标函数的详细步骤为:
[0057] 21)建立基于全变分正则项的图像去噪模型
[0058]
[0059] 其中f是带噪图像,u是待恢复的去噪图像,图像去噪的目的就是希望从f中恢复出 清晰的图像u。式中第一项为数据项,第二项为平滑项,λ是平衡数据项与平滑项的权重因 子,Ω为整个图像空间,X?为η维实数空间,▽为梯度算子,| Vu |也称为全变分正则项;
[0060] 22)将传统的全变分正则