[发明专利]一种基于L0正则化的自然图像盲去运动模糊的方法

摘要

本发明提供一种基于L0正则化的自然图像盲去运动模糊的方法，利用自然图像的梯度信息以及运动模糊核的稀疏特性，在模糊核的求解模型中分别引入相应L0正则项。模型求解时，首先利用基于半二次分裂方法分别对中间复原图像以及模糊核进行求解，为增加鲁棒性，求解过程中采用金字塔模型逐层求解，并且模糊核的估计只用到图像的梯度信息，通过傅里叶变换，将求解过程转换到频率域上进行，避免了直接在空域上进行反卷积运算，从而达到快速求解目的。然后利用求得的模糊核，采用基于全变差分的非盲去卷积方法求得最终的复原图像。本发明提供的方法求解快速，鲁棒性高，且最终的复原图像具有很好的视觉效果。

主权项

1.一种基于L0正则化的自然图像盲去运动模糊的方法，其特征在于，包括以下步骤：S1：对于输入的原始模糊图像转换成单通道的灰度图像；S2：求解模糊核时，引入L0正则项，构造求解模型如公式(1)所示：据此求解模型再结合金字塔模型，采用基于半二次分裂方法对公式(1)中模型进行求解得到模糊核k；其中及γ为权重参数，y是模糊图像，x是清晰图像，k是模糊核，表示x的一阶偏导；S3：根据步骤S2求解得到的模糊核k，对原始模糊图像中每个通道的图像均采用基于全变差分的非盲去卷积方法进行复原，再将每个通道复原后的图像合并即求得最终复原图像；所述步骤S2的处理过程如下：S201：根据人工输入的模糊核大小k_size，下采样因子及规定的最小核尺寸k_min_size确定金字塔模型的层级数目；S202：根据公式(1)，采用基于半二次分裂方法求复原图像x，如公式(2)：引入辅助变量β及u，其中β初始化为γ，u初始化为0，将公式(2)写成如下公式(3)：当β接近于无穷时，公式(2)与公式(3)求解得到的x一致，而公式(3)能够通过交替的计算x和u来求解；S203：通过公式(4)求u：由于公式(4)是逐个像素的最小化问题，因此u能够通过公式(5)求解得：S204：根据S203求出u后，x能够通过解最小二乘问题(6)来求得：因此有：其中F(.)和F^‑1(.)分别代表傅里叶变换和反傅里叶变换，是F(.)的复共轭形式，和分别代表垂直以及水平差分，S205：调整β＝2β，若β未超过限定的最大值，则根据S204求得的x作为输入，执行S203，否则，执行S206；S206：根据公式(1)，采用基于半二次分裂方法求解k，在公式(8)：引入辅助变量θ及g，θ初始为g初始化为0，将公式(8)写成如下公式(9)：当θ接近于无穷时，公式(8)与公式(9)求解得到的k一致，而(9)式能够通过交替的计算k和g来求解；S207：通过公式(10)求g：由于公式(10)是逐个像素的最小化问题，因此g能够通过公式(11)求解得：S208：根据S207求出g后，k能够通过解最小二乘问题(12)来求得：因此有其中对k中小于规定阈值k_threshold的值赋0，并对求得的k归一化；S209：调整θ＝2θ，若θ未超过限定的最大值，则根据S208求得的k作为输入，执行S207，否则，执行S210；S210：若金字塔模型中所有层次求解完成，则保留k，执行S3，否则，调整γ：γ＝max{γ/1.1，1e‑4}，若金字塔模型中当前层次求解尚未完成，则执行S202，否则将k，x上采样到金字塔模型中下一层级的对应尺寸，执行S202。