CN107292836A - 基于外部图像块先验信息与稀疏表示的图像盲去模糊方法 - Google Patents
基于外部图像块先验信息与稀疏表示的图像盲去模糊方法 Download PDFInfo
- Publication number
- CN107292836A CN107292836A CN201710408056.4A CN201710408056A CN107292836A CN 107292836 A CN107292836 A CN 107292836A CN 201710408056 A CN201710408056 A CN 201710408056A CN 107292836 A CN107292836 A CN 107292836A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- msubsup
- munder
- mtd
- 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
- 238000000034 method Methods 0.000 title claims abstract description 31
- 239000011159 matrix material Substances 0.000 claims abstract description 16
- 230000003416 augmentation Effects 0.000 claims abstract description 7
- 238000012549 training Methods 0.000 claims abstract description 4
- 238000005457 optimization Methods 0.000 claims description 4
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 8
- 238000002474 experimental method Methods 0.000 description 12
- 238000013461 design Methods 0.000 description 2
- 206010019133 Hangover Diseases 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/73—Deblurring; Sharpening
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
- Cosmetics (AREA)
Abstract
本发明公开了一种基于外部图像块先验信息与稀疏表示的图像盲去模糊方法,具体涉及到对图像块的外部先验估计并将其应用于图像去模糊框架中。对清晰图像库中的图像训练得到高斯混合模型,用图像块期望对数似然估计来描述模糊图像块外部先验知识,并将其作为正则项加入到基于稀疏的去模糊框架中。基于稀疏去模糊的框架中,迭代恢复中间图像与求解模糊核两部分。首先,用每一类单高斯的协方差矩阵构建单个中间图像块的字典,利用字典得到稀疏系数,构建出中间清晰图像。其次,采用增广拉格朗日算法求解模糊核。最终,采用非盲去卷积中的超拉普拉斯算法求解恢复出最终的清晰图像。实验结果表明该方法在抑制振铃和减弱噪声方面效果更好。
Description
技术领域
本发明涉及一种图像去模糊的方法,具体涉及到基于外部图像块先验信息与稀疏表示的图像盲去模糊方法。
背景技术
由于相机和拍摄场景之间的相对运动,获得到的图像常常存在一定程度的运动模糊。图像退化的模型可以表示为如下的卷积过程:
其中,Y是获得的模糊图像,K是模糊核,X是清晰图像,N是噪声。
由于图像去模糊具有重要的应用价值,目前,受到了广泛关注。越来越多自然图像的先验知识应用到去模糊框架中,文献(R.Fergus,B.Singh,A.Hertzmann,S.T.Roweis,andW.T.Freeman,Removing camera shake from a single photograph,ACM Transactionson Graphics,vol.25,no.3,pp.787–794,2006)提出在变分贝叶斯理论的础上,通过采用零均值高斯混合模型来逼近自然图像的重拖尾分布特性,对模糊核和整体图像进行估计。文献(Q.Shan,J.Jia,and A.Agarwala.High-quality motion deblurring from a singleimage,ACM Trans.Graph.,27(3),2008)结合了图像和模糊核函数具有稀疏性的特点,通过使用交替优化方法同时更新模糊核和整体去模糊图像。文献(J.F.Cai,H.Ji,C.Liu,Z.Shen,Framelet based blind motion deblurring from a single image,IEEETrans.Image Process.,21(2),562–572,2012)采用多尺度几何方法对模糊核和整体图像进行稀疏表示。文献(D.Krishnan,T.Tay,R.Fergus,et al.,Blind deconvolution usinga normalized sparsity measure,CVPR 2011,2011,233-240)提出了基于规范稀疏测量的整体图像去模糊方法。上述的这些方法都是基于对整体的图像去模糊,并没有利用图像块间的特性。而本方法考虑到图像块间的特性,提出了基于外部图像块先验信息与稀疏表示的图像盲去模糊。
发明内容
本发明的目的是利用外部图像块与稀疏表示进行图像盲去模糊,从而提高图像的清晰度。
本发明的技术方案如下:
一种基于外部图像块先验信息与稀疏表示的图像盲去模糊方法,包括如下步骤:
(1)、通过对清晰图像库中的图像训练得到高斯混合模型表示外部图像块先验知识;
(2)、采用图像块期望对数似然估计,利用上述外部图像块先验知识对待处理图像块进行估计,并将其作为正则项加入到基于稀疏的去模糊框架模型中;
(3)、对去模糊模型框架进行分析,通过对每一类单高斯的协方差矩阵构建每一个待求中间图像块的字典,利用字典得到稀疏系数,构建出中间清晰图像;
(4)、采用增广拉格朗日算法求解出模糊核;采用非盲去卷积中的超拉普拉斯算法求解恢复出最终的清晰图像。
上述步骤(3)中构建中间清晰图像的步骤如下:
通过求解如下模型表达式式(1),迭代优化中间清晰图像和稀疏系数:
其中,Y是模糊图像,K是模糊核,X是待求中间清晰图像,ηi,γi和λi是权重系数,xi是第i个待求中间清晰图像块,Di表示图像块xi对应的字典,ai是图像块xi的稀疏系数,μi和∑i分别是对图像块xi进行高斯混合模型概率表示时对应概率最高的单高斯模型的均值和协方差矩阵;
对上述的模型表达式式(1)中构建系数稀疏进行一定修改,增加噪声控制因子σ,表达式如下:
将式(2)进行优化,得到如下的表达式:
其中,Λi是式(1)中协方差矩阵∑i对应的特征值矩阵;
通过图像块对应单高斯模型的协方差矩阵∑i构建图像块字典,其表达式为:
中间清晰图像的构建如下:
上述步骤(4)中最终的清晰图像的求解步骤如下:
(a)、估计模糊核,构建的模糊核的模型表达式如下:
其中β表示权重系数,f2(K)的表达式如下:
(b)、采用分裂Bregman算法,将式(6)转化为:
(c)、引入增广拉格朗日算法,将式(8)优化为:
(d)、得到模糊核之后,通过超拉普拉斯算法恢复出最终的清晰图像,其表达式为:
其中,||·||p表示为p范数,0<p<1,α表示为权重系数;
通过引入辅助变量W,将上式(10)转化为:
其中,α,α1是权重系数;
固定W,的求解表达式为:
最终采用FFT算法求解清晰图像的闭合式为:
本发明的有益效果:从本实验的效果图可以很明显的得出结论,对于去模糊图像本方法在抑制振铃和减弱噪声方面效果更好。
附图说明
图1是本发明图像去模糊整体框架;
图2是本方法采用的模糊图像;
图3是实施例1中对灰度图像的实验结果效果图像;
图4是实施例1中对彩色图像的实验结果效果图像。
具体实施方式
下面结合附图,对本发明的基于外部图像块稀疏表示的图像盲去模糊方法进一步阐述。
如图1、图2所示,一种基于外部图像块先验信息与稀疏表示的图像盲去模糊方法,包括如下步骤:
(1)、通过对清晰图像库中的图像训练得到高斯混合模型表示外部图像块先验知识;
(2)、采用图像块期望对数似然估计,利用上述外部图像块先验知识对待处理图像块进行估计,并将其作为正则项加入到基于稀疏的去模糊框架模型中;
(3)、对去模糊模型框架进行分析,通过对每一类单高斯的协方差矩阵构建每一个待求中间图像块的字典,利用字典得到稀疏系数,构建出中间清晰图像;
(4)、采用增广拉格朗日算法求解出模糊核;采用非盲去卷积中的超拉普拉斯算法求解恢复出最终的清晰图像。
上述步骤(3)中构建中间清晰图像的步骤如下:
通过求解如下模型表达式式(1),迭代优化中间清晰图像和稀疏系数:
其中,Y是模糊图像,K是模糊核,X是待求中间清晰图像,ηi,γi和λi是权重系数,xi是第i个待求中间清晰图像块,Di表示图像块xi对应的字典,ai是图像块xi的稀疏系数,μi和∑i分别是对图像块xi进行高斯混合模型概率表示时对应概率最高的单高斯模型的均值和协方差矩阵;
对上述的模型表达式式(1)中构建系数稀疏进行一定修改,增加噪声控制因子σ,表达式如下:
其中,σ的具体数值可通过一些已有图像噪声估计方法获得。
将式(2)进行优化,得到如下的表达式:
其中,Λi是式(1)中协方差矩阵∑i对应的特征值矩阵;
通过图像块对应单高斯模型的协方差矩阵∑i构建图像块字典,其表达式为:
中间清晰图像的构建如下:
上述步骤(4)中最终的清晰图像的求解步骤如下:
(a)、估计模糊核,构建的模糊核的模型表达式如下:
其中β表示权重系数,f2(K)的表达式如下:
(b)、采用分裂Bregman算法,将式(6)转化为:
(c)、引入增广拉格朗日算法,将式(8)优化为:
(d)、得到模糊核之后,通过超拉普拉斯算法恢复出最终的清晰图像,其表达式为:
其中,||·||p表示为p范数,0<p<1,α表示为权重系数;
通过引入辅助变量W,将上式(10)转化为:
其中,α,α1是权重系数;
固定W,的求解表达式为:
最终采用FFT算法求解清晰图像的闭合式为:
实施例1:
1)实验条件
本实验中具体参数设置为γi=2,β=1,α=3×10-4,p=2/3。其中模糊核的大小设置从15×15到35×35。为了便于实验比较,专利中所有算法都采用同一个大小的模糊核。
2)实验内容
根据具体的实验参数设置,对多种图像去模糊算法进行实验仿真,得出具体实验结果如下所示。图3共两组灰度图像,每组图像的左边的是(Z.Hu,B.J.Huang,H.M.Yang,"Single image deblurring with adaptive dictionary learning,"in IEEEInt.Conf.on Imaging Process.,pp.1169-1172,2010)的实验效果图。每组图像的右边是本方法的实验效果图。从塔顶部和女孩衣襟细节处看出,我们的实验恢复出来的纹理效果更好,减弱了噪声,抑制了振铃现象。本实验的结果在细节处和平滑处在都优于Hu的方法。
图4共四组彩色图像,每组分别是6种去模糊方法,从左到右依次是:Fergus,Pan,Krishnan,Cai,Filip,本方法。很明显的看出本方法在抑制振铃和减弱噪声方面比其他5种的算法效果更佳。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。
Claims (3)
1.一种基于外部图像块先验信息与稀疏表示的图像盲去模糊方法,其特征在于包括如下步骤:
(1)、通过对清晰图像库中的图像训练得到高斯混合模型表示外部图像块先验知识;
(2)、采用图像块期望对数似然估计,利用上述外部图像块先验知识对待处理图像块进行估计,并将其作为正则项加入到基于稀疏的去模糊框架模型中;
(3)、对去模糊模型框架进行分析,通过对每一类单高斯的协方差矩阵构建每一个中间图像块的字典,利用字典得到稀疏系数,构建出中间清晰图像;
(4)、采用增广拉格朗日算法求解出模糊核;采用非盲去卷积中的超拉普拉斯算法求解最终的清晰图像。
2.根据权利要求1所述的基于外部图像块先验信息与稀疏表示的图像盲去模糊方法,其特征在于所述步骤(3)中构建中间清晰图像的步骤如下:
通过求解如下模型表达式式(1),迭代优化中间清晰图像和稀疏系数:
<mrow>
<mo><</mo>
<mover>
<mi>X</mi>
<mo>~</mo>
</mover>
<mo>,</mo>
<mrow>
<mo>{</mo>
<msub>
<mover>
<mi>a</mi>
<mo>~</mo>
</mover>
<mi>i</mi>
</msub>
<mo>}</mo>
</mrow>
<mo>></mo>
<mo>=</mo>
<mi>arg</mi>
<munder>
<mi>min</mi>
<mrow>
<mi>X</mi>
<mo>,</mo>
<mrow>
<mo>{</mo>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<mo>}</mo>
</mrow>
</mrow>
</munder>
<mo>|</mo>
<mo>|</mo>
<mi>Y</mi>
<mo>-</mo>
<mi>K</mi>
<mo>&CircleTimes;</mo>
<mi>X</mi>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<munder>
<mi>&Sigma;</mi>
<mi>i</mi>
</munder>
<msub>
<mi>&eta;</mi>
<mi>i</mi>
</msub>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>D</mi>
<mi>i</mi>
</msub>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<munder>
<mi>&Sigma;</mi>
<mi>i</mi>
</munder>
<msub>
<mi>&lambda;</mi>
<mi>i</mi>
</msub>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<mo>|</mo>
<msub>
<mo>|</mo>
<mn>0</mn>
</msub>
<mo>+</mo>
<munder>
<mi>&Sigma;</mi>
<mi>i</mi>
</munder>
<msub>
<mi>&gamma;</mi>
<mi>i</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>D</mi>
<mi>i</mi>
</msub>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&mu;</mi>
<mi>i</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<msubsup>
<mi>&Sigma;</mi>
<mi>i</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>D</mi>
<mi>i</mi>
</msub>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&mu;</mi>
<mi>i</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,Y是模糊图像,K是模糊核,X是待求中间清晰图像,ηi,γi和λi是权重系数,xi是第i个待求中间清晰图像块,Di表示图像块xi对应的字典,ai是图像块xi的稀疏系数,μi和∑i分别是对图像块xi进行高斯混合模型概率表示时对应概率最高的单高斯模型的均值和协方差矩阵;
对上述的模型表达式式(1)中构建系数稀疏进行修改,增加噪声控制因子σ,表达式如下:
<mrow>
<msub>
<mover>
<mi>a</mi>
<mo>~</mo>
</mover>
<mi>i</mi>
</msub>
<mo>=</mo>
<mi>arg</mi>
<munder>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
</munder>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>D</mi>
<mi>i</mi>
</msub>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msub>
<mi>&gamma;</mi>
<mi>i</mi>
</msub>
<mi>&sigma;</mi>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>D</mi>
<mi>i</mi>
</msub>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&mu;</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<msubsup>
<mi>&Sigma;</mi>
<mi>i</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>D</mi>
<mi>i</mi>
</msub>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&mu;</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>&lambda;</mi>
<mi>i</mi>
</msub>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<mo>|</mo>
<msub>
<mo>|</mo>
<mn>0</mn>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
将式(2)进行优化,得到如下的表达式:
<mrow>
<msub>
<mover>
<mi>a</mi>
<mo>~</mo>
</mover>
<mi>i</mi>
</msub>
<mo>=</mo>
<mi>arg</mi>
<munder>
<mi>min</mi>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
</munder>
<mo>|</mo>
<mo>|</mo>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>&gamma;</mi>
<mi>i</mi>
<mo>&prime;</mo>
</msubsup>
<msubsup>
<mi>&Lambda;</mi>
<mi>i</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msubsup>
<msubsup>
<mi>D</mi>
<mi>i</mi>
<mi>T</mi>
</msubsup>
<msub>
<mi>&mu;</mi>
<mi>i</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<msub>
<mi>D</mi>
<mi>i</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>&gamma;</mi>
<mi>i</mi>
<mo>&prime;</mo>
</msubsup>
<msubsup>
<mi>&Lambda;</mi>
<mi>i</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msub>
<mi>&lambda;</mi>
<mi>i</mi>
</msub>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<mo>|</mo>
<msub>
<mo>|</mo>
<mn>0</mn>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,Λi是式(1)中协方差矩阵∑i对应的特征值矩阵;
通过图像块对应单高斯模型的协方差矩阵∑i构建图像块字典,其表达式为:
<mrow>
<msub>
<mi>&Sigma;</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<msub>
<mi>D</mi>
<mi>i</mi>
</msub>
<msub>
<mi>&Lambda;</mi>
<mi>i</mi>
</msub>
<msubsup>
<mi>D</mi>
<mi>i</mi>
<mi>T</mi>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
中间清晰图像的构建如下:
<mrow>
<mover>
<mi>X</mi>
<mo>~</mo>
</mover>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<munder>
<mo>&Sigma;</mo>
<mi>i</mi>
</munder>
<msubsup>
<mi>R</mi>
<mi>i</mi>
<mi>T</mi>
</msubsup>
<msub>
<mi>R</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<munder>
<mo>&Sigma;</mo>
<mi>i</mi>
</munder>
<msubsup>
<mi>R</mi>
<mi>i</mi>
<mi>T</mi>
</msubsup>
<msub>
<mi>D</mi>
<mi>i</mi>
</msub>
<msub>
<mover>
<mi>a</mi>
<mo>~</mo>
</mover>
<mi>i</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
3.根据权利要求1所述的基于外部图像块先验信息与稀疏表示的图像盲去模糊方法,其特征在于所述步骤(4)中最终的清晰图像的求解步骤如下:
(a)、估计模糊核,构建的模糊核的模型表达式如下:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mover>
<mi>K</mi>
<mo>~</mo>
</mover>
<mo>=</mo>
<mi>arg</mi>
<munder>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
<mi>K</mi>
</munder>
<mo>|</mo>
<mo>|</mo>
<mi>Y</mi>
<mo>-</mo>
<mi>K</mi>
<mo>&CircleTimes;</mo>
<mi>X</mi>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msub>
<mi>&beta;f</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>K</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>f</mi>
<mn>2</mn>
</msub>
<mo>(</mo>
<mi>K</mi>
<mo>)</mo>
<mo>=</mo>
<munder>
<mi>&Sigma;</mi>
<mi>j</mi>
</munder>
<mi>&psi;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>k</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
其中β表示权重系数,f2(K)的表达式如下:
<mrow>
<msub>
<mi>f</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>K</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munder>
<mi>&Sigma;</mi>
<mi>j</mi>
</munder>
<mi>&psi;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>k</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mi>&psi;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>k</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<msub>
<mi>k</mi>
<mi>j</mi>
</msub>
</mtd>
<mtd>
<mrow>
<msub>
<mi>k</mi>
<mi>j</mi>
</msub>
<mo>&GreaterEqual;</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mi>&infin;</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>k</mi>
<mi>j</mi>
</msub>
<mo><</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
(b)、采用分裂Bregman算法,将式(6)转化为:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mover>
<mi>K</mi>
<mo>~</mo>
</mover>
<mo>=</mo>
<mi>arg</mi>
<munder>
<mi>min</mi>
<mi>K</mi>
</munder>
<mo>|</mo>
<mo>|</mo>
<mi>Y</mi>
<mo>-</mo>
<mi>K</mi>
<mo>&CircleTimes;</mo>
<mi>X</mi>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msub>
<mi>&beta;f</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>H</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>H</mi>
<mo>=</mo>
<mi>K</mi>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
(c)、引入增广拉格朗日算法,将式(8)优化为:
<mrow>
<mo><</mo>
<mover>
<mi>K</mi>
<mo>~</mo>
</mover>
<mo>,</mo>
<mover>
<mi>H</mi>
<mo>~</mo>
</mover>
<mo>></mo>
<mo>=</mo>
<mi>arg</mi>
<munder>
<mi>min</mi>
<mrow>
<mi>K</mi>
<mo>,</mo>
<mi>H</mi>
</mrow>
</munder>
<mo>|</mo>
<mo>|</mo>
<mi>Y</mi>
<mo>-</mo>
<mi>K</mi>
<mo>&CircleTimes;</mo>
<mi>X</mi>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msub>
<mi>&beta;f</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>H</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>&beta;</mi>
<mn>1</mn>
</msub>
<mo>|</mo>
<mo>|</mo>
<mi>K</mi>
<mo>-</mo>
<mi>H</mi>
<mo>-</mo>
<mi>B</mi>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
(d)、得到模糊核之后,通过超拉普拉斯算法恢复出最终的清晰图像,其表达式为:
<mrow>
<mover>
<mi>X</mi>
<mo>~</mo>
</mover>
<mo>=</mo>
<munder>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
<mi>X</mi>
</munder>
<mo>|</mo>
<mo>|</mo>
<mi>Y</mi>
<mo>-</mo>
<mi>K</mi>
<mo>&CircleTimes;</mo>
<mi>X</mi>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<mi>&alpha;</mi>
<mo>|</mo>
<mo>|</mo>
<mo>&dtri;</mo>
<mi>X</mi>
<mo>|</mo>
<msub>
<mo>|</mo>
<mi>p</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,||.||p表示为p范数,0<p<1,α表示为权重系数;
通过引入辅助变量W,将上式(10)转化为:
<mrow>
<mo><</mo>
<mover>
<mi>X</mi>
<mo>~</mo>
</mover>
<mo>,</mo>
<mover>
<mi>W</mi>
<mo>~</mo>
</mover>
<mo>></mo>
<mo>=</mo>
<mi>arg</mi>
<munder>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
<mrow>
<mi>X</mi>
<mo>,</mo>
<mi>W</mi>
</mrow>
</munder>
<mo>|</mo>
<mo>|</mo>
<mi>Y</mi>
<mo>-</mo>
<mi>K</mi>
<mo>&CircleTimes;</mo>
<mi>X</mi>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<mi>&alpha;</mi>
<mo>|</mo>
<mo>|</mo>
<mi>W</mi>
<mo>|</mo>
<msub>
<mo>|</mo>
<mi>p</mi>
</msub>
<mo>+</mo>
<msub>
<mi>&alpha;</mi>
<mn>1</mn>
</msub>
<mo>|</mo>
<mo>|</mo>
<mo>&dtri;</mo>
<mi>X</mi>
<mo>-</mo>
<mi>W</mi>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,α,α1是权重系数;
固定W,的求解表达式为:
<mrow>
<mover>
<mi>X</mi>
<mo>~</mo>
</mover>
<mo>=</mo>
<mi>arg</mi>
<munder>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
<mi>X</mi>
</munder>
<mo>|</mo>
<mo>|</mo>
<mi>Y</mi>
<mo>-</mo>
<mi>K</mi>
<mo>&CircleTimes;</mo>
<mi>X</mi>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msub>
<mi>&alpha;</mi>
<mn>1</mn>
</msub>
<mo>|</mo>
<mo>|</mo>
<mo>&dtri;</mo>
<mi>X</mi>
<mo>-</mo>
<mi>W</mi>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
最终采用FFT算法求解清晰图像的闭合式为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710408056.4A CN107292836B (zh) | 2017-06-02 | 2017-06-02 | 基于外部图像块先验信息与稀疏表示的图像盲去模糊方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710408056.4A CN107292836B (zh) | 2017-06-02 | 2017-06-02 | 基于外部图像块先验信息与稀疏表示的图像盲去模糊方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107292836A true CN107292836A (zh) | 2017-10-24 |
CN107292836B CN107292836B (zh) | 2020-06-26 |
Family
ID=60094985
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710408056.4A Expired - Fee Related CN107292836B (zh) | 2017-06-02 | 2017-06-02 | 基于外部图像块先验信息与稀疏表示的图像盲去模糊方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107292836B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110648286A (zh) * | 2019-08-02 | 2020-01-03 | 广东工业大学 | 一种毫米波图像非盲目去模糊的方法 |
CN111640077A (zh) * | 2020-05-29 | 2020-09-08 | 扆亮海 | 一种简单高效的模糊文本图片清晰化处理方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102354395A (zh) * | 2011-09-22 | 2012-02-15 | 西北工业大学 | 基于稀疏表示的模糊图像盲复原方法 |
CN105046664A (zh) * | 2015-07-13 | 2015-11-11 | 广东工业大学 | 一种基于自适应epll算法的图像去噪方法 |
CN105184742A (zh) * | 2015-08-07 | 2015-12-23 | 河海大学常州校区 | 一种基于拉普拉斯图特征向量的稀疏编码的图像去噪方法 |
CN105184744A (zh) * | 2015-08-24 | 2015-12-23 | 西安电子科技大学 | 基于标准化稀疏度量图像块先验的模糊核估计方法 |
CN105279740A (zh) * | 2015-11-17 | 2016-01-27 | 华东交通大学 | 一种基于稀疏正则化的图像去噪方法 |
CN105957024A (zh) * | 2016-04-20 | 2016-09-21 | 西安电子科技大学 | 基于图像块先验与稀疏范数的盲去模糊方法 |
CN105957025A (zh) * | 2016-04-21 | 2016-09-21 | 天津大学 | 基于稀疏表示的非一致图像盲复原方法 |
US20160335747A1 (en) * | 2015-04-22 | 2016-11-17 | Adobe Systems Incorporated | Scale adaptive blind deblurring |
CN106339996A (zh) * | 2016-09-09 | 2017-01-18 | 江南大学 | 一种基于超拉普拉斯先验的图像盲去模糊方法 |
-
2017
- 2017-06-02 CN CN201710408056.4A patent/CN107292836B/zh not_active Expired - Fee Related
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102354395A (zh) * | 2011-09-22 | 2012-02-15 | 西北工业大学 | 基于稀疏表示的模糊图像盲复原方法 |
US20160335747A1 (en) * | 2015-04-22 | 2016-11-17 | Adobe Systems Incorporated | Scale adaptive blind deblurring |
CN105046664A (zh) * | 2015-07-13 | 2015-11-11 | 广东工业大学 | 一种基于自适应epll算法的图像去噪方法 |
CN105184742A (zh) * | 2015-08-07 | 2015-12-23 | 河海大学常州校区 | 一种基于拉普拉斯图特征向量的稀疏编码的图像去噪方法 |
CN105184744A (zh) * | 2015-08-24 | 2015-12-23 | 西安电子科技大学 | 基于标准化稀疏度量图像块先验的模糊核估计方法 |
CN105279740A (zh) * | 2015-11-17 | 2016-01-27 | 华东交通大学 | 一种基于稀疏正则化的图像去噪方法 |
CN105957024A (zh) * | 2016-04-20 | 2016-09-21 | 西安电子科技大学 | 基于图像块先验与稀疏范数的盲去模糊方法 |
CN105957025A (zh) * | 2016-04-21 | 2016-09-21 | 天津大学 | 基于稀疏表示的非一致图像盲复原方法 |
CN106339996A (zh) * | 2016-09-09 | 2017-01-18 | 江南大学 | 一种基于超拉普拉斯先验的图像盲去模糊方法 |
Non-Patent Citations (6)
Title |
---|
JEREMIAS SULAM 等: "Expected Patch Log Likelihood with a Sparse Prior", 《ENERGY MINIMIZATION METHODS IN COMPUTER VISION AND PATTERN RECOGNITION》 * |
JINSHAN PAN 等: "Blind Image Deblurring Using Dark Channel Prior", 《THE IEEE CONFERENCE ON COMPUTER VISION AND PATTERN RECOGNITION (CVPR)》 * |
WEISHENG DONG 等: "Image Restoration via Simultaneous Sparse Coding: Where Structured Sparsity Meets Gaussian Scale Mixture", 《INTERNATIONAL JOURNAL OF COMPUTER VISION 》 * |
廖斌 等: "基于空间约束高斯混合模型的EPLL自然图像复原", 《量子电子学报》 * |
蔡德生: "基于字典稀疏表示和梯度稀疏的图像盲去模糊", 《燕山大学学报》 * |
詹曙 等: "基于Gabor特征和字典学习的高斯混合稀疏表示图像识别", 《电子学报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110648286A (zh) * | 2019-08-02 | 2020-01-03 | 广东工业大学 | 一种毫米波图像非盲目去模糊的方法 |
CN111640077A (zh) * | 2020-05-29 | 2020-09-08 | 扆亮海 | 一种简单高效的模糊文本图片清晰化处理方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107292836B (zh) | 2020-06-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zhang et al. | Learning deep CNN denoiser prior for image restoration | |
Remez et al. | Class-aware fully convolutional Gaussian and Poisson denoising | |
CN106683067B (zh) | 一种基于残差子图像的深度学习超分辨率重建方法 | |
Bigdeli et al. | Image restoration using autoencoding priors | |
CN104599242B (zh) | 使用多尺度非局部正则的模糊核估计方法 | |
Matakos et al. | Accelerated edge-preserving image restoration without boundary artifacts | |
CN106991440B (zh) | 一种基于空间金字塔的卷积神经网络的图像分类方法 | |
Sun et al. | Lightweight image super-resolution via weighted multi-scale residual network | |
CN106709877B (zh) | 一种基于多参数正则优化模型的图像去模糊方法 | |
CN107133923B (zh) | 一种基于自适应梯度稀疏模型的模糊图像非盲去模糊方法 | |
CN107464217B (zh) | 一种图像处理方法及装置 | |
CN112541877B (zh) | 基于条件生成对抗网络的去模糊方法、系统、设备及介质 | |
CN107590779B (zh) | 一种基于图像块聚类字典训练的图像去噪去模糊方法 | |
CN110533579B (zh) | 基于自编码结构与梯度保序的视频风格转换方法 | |
CN113689517B (zh) | 一种多尺度通道注意力网络的图像纹理合成方法及系统 | |
Sahu et al. | Blind deblurring using deep learning: A survey | |
CN113781308A (zh) | 图像超分辨率重建方法、装置、存储介质及电子设备 | |
CN107067373A (zh) | 一种基于0范数的梯度最小化二值图像复原方法 | |
CN107292836A (zh) | 基于外部图像块先验信息与稀疏表示的图像盲去模糊方法 | |
CN112801899B (zh) | 基于互补结构感知的内外循环驱动图像盲去模糊方法和装置 | |
Xiao et al. | A dive into sam prior in image restoration | |
CN114626984A (zh) | 中文文本图像的超分辨率重建方法 | |
Umer et al. | Deep iterative residual convolutional network for single image super-resolution | |
CN112734763B (zh) | 基于卷积与k-svd字典联合稀疏编码的图像分解方法 | |
CN111986121B (zh) | 基于Framelet l0范数约束的模糊图像非盲复原方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200626 |