一种基于噪声水平的梯度相关TV因子图像去噪去模糊方法技术

本发明专利技术公开了一种基于噪声水平的梯度相关TV因子图像去噪去模糊方法。日常拍摄所获取的图像都有一定程度的降质，即使是对静态目标的稳定成像，其获取图像也包含有像素级的模糊量，且噪声往往不可避免。对带噪模糊图像进行复原操作可以提升图像质量，然而，抑制噪声和去模糊两者间存在矛盾，图像复原时需兼顾两者的平衡。本发明专利技术方法提出一种基于总变分正则化的图像去噪去模糊方法，在对图像进行模糊复原的同时，结合Richardson‑Lucy算法的特点获取初步去模糊图像，通过若干正则化项进行约束，并利用噪声图像梯度的先验分布特性，对正则化权重因子进行梯度相关的分布式调整处理，能有效地抑制噪声同时保持良好的图像边缘特性，获得较高质量的图像。

An image denoising and deblurring method based on noise level and gradient correlation TV factor

The invention discloses an image denoising and deblurring method based on noise level and gradient correlation TV factor. The image acquired by daily shooting is degraded to a certain extent. Even for stable imaging of static targets, the acquired image also contains pixel-level blur, and noise is often unavoidable. Restoration of noisy blurred images can improve image quality. However, there is a contradiction between noise suppression and de-blurring. The balance between them should be taken into account in image restoration. The method of the invention proposes an image denoising and deblurring method based on total variation regularization. While restoring the blurred image, the preliminary deblurring image is obtained by combining the characteristics of Richardson_Lucy algorithm, and the regularization weight is obtained by restricting several regularization terms and utilizing the prior distribution characteristics of the noise image gradient. Multi-factor gradient-dependent distributed adjustment can effectively suppress noise while maintaining good image edge characteristics, and obtain high quality images.

【技术实现步骤摘要】
一种基于噪声水平的梯度相关TV因子图像去噪去模糊方法
本专利技术属于数字图像处理领域，涉及一种基于噪声水平的梯度相关TV因子图像去噪去模糊方法。
技术介绍
随着摄影设备和摄影技术的普及，摄影成为人们日常生活中重要的活动，相应地，人们对图像质量的要求也日趋提高。在自然图像形成的过程中，可能会受到如目标场景运动、成像设备抖动、成像传感器噪声等各方面的影响，最终获得一幅模糊降质的图像。一般来说，成像系统为线性移不变的前提下，带噪模糊图像的形成过程可以简单描述为：y(u,v)＝k(u,v)*x(u,v)+n(u,v)其中：k(u,v)为点扩散函数(图像模糊核，PSF)，*为卷积符号，x(u,v)表示清晰图像，n(u,v)代表加性噪声。图像模糊是摄像中一类常见的问题，即使是对于静止目标的稳定成像，也仍然存在小幅度的图像模糊，且噪声难以避免。图1所示为日常拍摄图像及其细节放大图。对带噪模糊图像进行复原处理是一种常见的手段。在进行图像复原的过程中，去模糊和去噪两者间存在矛盾，往往难以平衡。在一般的图像处理中，单独针对噪声图像的去噪方法不少，且能达到很好的效果，但以牺牲图像高频信息作为代价。本专利技术结合带噪图像的梯度分布先验特性，将其具体应用于图像复原中并取得显著效果。本专利技术提出的方法属于图像非盲复原范畴，在对图像进行模糊复原同时，结合Richardson-Lucy算法的特点获取初步去模糊图像，通过若干正则化项进行约束，并利用噪声图像梯度的先验分布特性，对正则化权重因子进行梯度相关的分布式调整处理，能有效地抑制噪声同时保持良好的图像边缘特性。
技术实现思路
图2所示为自然清晰图像和带噪清晰图像及其各自的梯度分布特性统计图对比，分析可得到一般带噪图像的梯度分布规律：梯度数据主要分布在低梯度范围，在较低梯度处存在极大值并向两边呈现不均等稀疏分布。本专利技术的目的是利用噪声梯度的不同响应，应用RL算法获取模糊图像得初步复原结果，并应用于随后的正则化变分方法图像复原操作中，通过添加图像约束条件，对正则化约束因子进行梯度相关的分布式调整处理，最终达到抑制噪声同时保持良好的图像边缘特性的复原效果。为达到以上目的，本专利技术采用以下技术方案：一种基于噪声水平的梯度相关TV因子图像去噪去模糊方法，该方法包括以下步骤：(1)利用Richardson-Lucy(RL)算法获取初步去模糊图像I0；(2)引入传统总变分正则化方法的能量方程模型：其中，I为潜在清晰图像，k为图像模糊核，λ为约束因子，为一阶求导操作，为卷积操作，E(I)为图像I的能量；(3)在步骤(2)能量方程模型的基础上，引入新变量w、v、θ、γ、β，加入梯度惩罚项和二次梯度平滑项，得到如下优化模型：其中Ω为图像积分域；将连续的优化模型转化为离散模型，该离散模型关于图像的解为：其中，F和F-1代表傅里叶变换与傅里叶逆变换，F*代表复共轭操作，x、y代表图像的两个积分方向，为梯度算子；同时，该离散模型分解为w子问题和v子问题：w子问题：v子问题：(4)梯度惩罚因子分布式处理，具体是：(4.1)引入梯度相关TV惩罚因子λs，其表达式如下：其中，变量s决定λs的整体约束力度，s∈[0,2]；表示梯度阈值，一般取值0.04；参数p1、b、d、p2、κ共同决定λs的整体范围，各自的范围是：p1∈[0,5]，b∈[50,300]，d∈[0,20]，p2∈[-10,50]，κ∈[0,0.5]；(4.2)利用λs替代约束因子λ，交替求解w子问题、v子问题和清晰图像I；(5)控制求解过程中的β、θ参数，反复迭代求解，直到得到清晰图像I。进一步地，将Richardson-Lucy算法与总变分正则化方法相结合，Richardson-Lucy算法获取初步去模糊图像作为总变分正则化方法的先验结果并应用于后续复原过程中。进一步地，引入梯度相关TV惩罚因子λs考虑的因素为：带噪图像的梯度数据主要分布在低梯度范围，在低梯度处存在极大值并向两边呈现不均等稀疏分布。进一步地，所述步骤(3)中，将连续的优化模型转化为离散模型：进一步地，所述步骤(4.2)中，w子问题的解为：本专利技术的有益效果：利用噪声梯度的不同响应，应用RL算法获取模糊图像得初步复原结果，并应用于随后的正则化变分方法图像复原操作中，通过添加图像约束条件，对正则化约束因子进行梯度相关的分布式调整处理，最终获取去噪去模糊图像。本专利技术适用于处理一般的带噪模糊图像，根据噪声梯度特性在去模糊过程中有效区别图像细节和噪声，能有效地抑制噪声同时保持良好的图像边缘特性，兼顾去模糊和去噪效果，同时由于RL算法的引入能一定程度上减少时间代价。附图说明图1为日常拍摄模糊图像及其细节放大示例，其中(a)为拍摄图，(b)为细节放大图。图2为带噪图像及其梯度统计图，其中(a)为清晰图像，(b)为添加均值0方差0.001的噪声方差后的图像结果，(c)为(a)的梯度分布统计特性图，(d)为(b)的梯度分布统计特性图。图3为本专利技术方法总体示意图。图4为本专利技术方法细节流程示意图。图5为参数设置为p1＝0.3、b＝90、d＝12、κ＝0.03、p2＝50、s＝0.5、时梯度相关的分布式因子λs的数值分布。图6复原效果对比图，其中，(a)为带噪声模糊图像，(b)为RL算法复原结果，(c)为FTVD算法复原结果，(d)为本专利技术方法复原结果。具体实施方式以下结合附图和具体实施例对本专利技术作进一步说明。本专利技术利用噪声梯度的不同响应，应用RL算法获取模糊图像得初步复原结果，并应用于随后的正则化变分方法图像复原操作中，通过添加图像约束条件，对正则化约束因子进行梯度相关的分布式调整处理，最终达到抑制噪声同时保持良好的图像边缘特性的复原效果。本专利技术的流程如附图3、图4所示，主要包括RL算法获取初步去模糊图像、图像正则化能量方程引入、新变量引入与优化方程求解、梯度惩罚因子分布式处理、迭代计算等几个步骤。步骤1.利用Richardson-Lucy(RL)算法获取初步去模糊图像。总的来说，RL算法假定图像噪声符合泊松分布，采用最大似然法估计清晰图像，是一种基于贝叶斯分析的迭代复原算法。对于泊松噪声模型，图像的似然概率可以表达为：为泊松过程。经过推导可得到RL算法的迭代式：其中，K*为K的伴随矩阵，t为迭代次数。虽然RL算法从泊松噪声模型的角度推导得到，但同样适用于其他类型的噪声模型。经过RL算法获取的初步较清晰图像将应用到随后的正则化变分模型中。步骤2.引入传统总变分正则化方法的能量方程模型，其为：上式中，I代表潜在图像，I0代表待处理图像，h代表图像模糊核，λx为约束因子，为一阶求导操作，为卷积操作。步骤3.新变量引入与优化方程求解，包括以下步骤：3-1在传统能量方程模型的基础上，引入新变量w、v、θ、γ，加入梯度惩罚项和二次梯度平滑项，有如下模型：其中Ω为图像积分域；当β趋近于0，θ趋近于0时，上述模型的解收敛为模型(3)的解。3-2模型(4)为连续的总变分模型，其可转化为离散模型如下：上式为I，w，v的最优化函数，其关于I有一闭合形式的最优解，作傅里叶变换后可等效求得：其中，F和F-1代表傅里叶变换与傅里叶逆变换，F*代表复共轭操作；3-3关于w和v的求解，由于w和v在各自的函数中独立，故其值可通过交替最小本文档来自技高网...

【技术保护点】
1.一种基于噪声水平的梯度相关TV因子图像去噪去模糊方法，其特征在于，该方法包括以下步骤：(1)利用Richardson‑Lucy(RL)算法获取初步去模糊图像I0；(2)引入传统总变分正则化方法的能量方程模型：

【技术特征摘要】
1.一种基于噪声水平的梯度相关TV因子图像去噪去模糊方法，其特征在于，该方法包括以下步骤：(1)利用Richardson-Lucy(RL)算法获取初步去模糊图像I0；(2)引入传统总变分正则化方法的能量方程模型：其中，I为潜在清晰图像，k为图像模糊核，λ为约束因子，为一阶求导操作，为卷积操作，E(I)为图像I的能量；(3)在步骤(2)能量方程模型的基础上，引入新变量w、v、θ、γ、β，加入梯度惩罚项和二次梯度平滑项，得到如下优化模型：其中Ω为图像积分域；将连续的优化模型转化为离散模型，该离散模型关于图像的解为：其中，F和F-1代表傅里叶变换与傅里叶逆变换，F*代表复共轭操作，x、y代表图像的两个积分方向，为梯度算子；同时，该离散模型分解为w子问题和v子问题：w子问题：v子问题：(4)梯度惩罚因子分布式处理，具体是：(4.1)引入梯度相关TV惩罚因子λs，其表达式如下：其中，变量s决定λs的整体约束力度，s∈[0,2]；表示梯度阈值，参数p1、b、d、p2、κ共同决定λs的整体范围，各自的范围是：p1∈[0,5]，b∈[50,300]，d...

【专利技术属性】
技术研发人员：冯华君，黄加紫，徐之海，李奇，陈跃庭，
申请(专利权)人：浙江大学，
类型：发明
国别省市：浙江,33