发明名称 基于范数比值正则化的快速图像盲去模糊方法
摘要 本发明公开了一种基于范数比值正则化的快速图像盲去模糊方法,包括利用范数比值l<sub>1/2</sub>/l<sub>2</sub>的正则化项作为先验知识,采用多尺度的方法估计出模糊核;利用得到的模糊核k矩阵,通过封闭式的阈值公式的非盲去卷积方法,快速和高质量地恢复出原始清晰图像u。本发明求模糊核的过程由粗尺度到细尺度逐渐进行,多尺度的算法保证了模糊核函数计算的准确性和鲁棒性;采用范数比值的先验模型更逼近自然图像的梯度分布,使得复原结果更准确,同时计算效率高,性能要优于传统算法;采用范数比值正则化先验作为光滑项,保证求解时能量是下降的;在估计出模糊核后,利用封闭式的阈值公式的非盲去卷积的方法,能够快速和高质量地得到清晰图像。
申请公布号 CN104112261A 申请公布日期 2014.10.22
申请号 CN201410342554.X 申请日期 2014.07.17
申请人 五邑大学 发明人 余义斌;彭念;甘俊英
分类号 G06T5/00(2006.01)I 主分类号 G06T5/00(2006.01)I
代理机构 广州嘉权专利商标事务所有限公司 44205 代理人 谭志强
主权项 基于范数比值正则化的快速图像盲去模糊方法,其特征在于,包括以下步骤:(1)输入模糊图像;(2)初始化模糊核k矩阵大小为3×3;(3)利用范数比值l<sub>1/2</sub>/l<sub>2</sub>的正则化项作为先验知识,采用多尺度的方法,估计出模糊核:(3a)构建求解模型:给出模糊函数f,利用离散滤波器▽<sub>x</sub>=[1,‑1],▽<sub>y</sub>=[1;‑1],产生高频的图像信息y=[▽<sub>x</sub>f,▽<sub>y</sub>f],构建求解模型:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><munder><mi>min</mi><mi>x</mi></munder><mi>&lambda;</mi><msub><mrow><mo>|</mo><mo>|</mo><mi>x</mi><mo>&CircleTimes;</mo><mi>k</mi><mo>-</mo><mi>y</mi><mo>|</mo><mo>|</mo></mrow><mn>1</mn></msub><mo>+</mo><mfrac><msub><mrow><mo>|</mo><mo>|</mo><mi>x</mi><mo>|</mo><mo>|</mo></mrow><mi>q</mi></msub><msub><mrow><mo>|</mo><mo>|</mo><mi>x</mi><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msub></mfrac></mtd></mtr><mtr><mtd><mrow><munder><mi>min</mi><mi>k</mi></munder><msub><mi>&lambda;</mi><mn>1</mn></msub></mrow><msubsup><mrow><mo>|</mo><mo>|</mo><mi>x</mi><mo>&CircleTimes;</mo><mi>k</mi><mo>-</mo><mi>y</mi><mo>|</mo><mo>|</mo></mrow><mn>2</mn><mn>2</mn></msubsup><mo>+</mo><mi>&psi;</mi><msub><mrow><mo>|</mo><mo>|</mo><mi>k</mi><mo>|</mo><mo>|</mo></mrow><mn>1</mn></msub></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000539471870000011.GIF" wi="739" he="290" /></maths>其中,k≥0,Σ<sub>i</sub>k<sub>i</sub>=1,k<sub>i</sub>表示模糊核k中的元素,x表示高频领域的未知清晰图像,q取1/2,<img file="FDA0000539471870000012.GIF" wi="50" he="49" />表示二维的卷积操作;(3b)估计模糊核:采用多尺度估计,每一尺度的运算,都涉及到清晰图像和模糊核的插值与更新;插值过程是将粗尺度更新出的清晰图像x和模糊核k进行上采样来作为细尺度的初始值;更新过程包括清晰图像x的更新和模糊核k的更新的两个过程;最终迭代得到模糊核k矩阵;(4)利用得到的模糊核k矩阵,恢复出原始清晰图像u:(4a)采用超拉普拉斯模型作为先验知识,对自然图像的梯度分布进行近似模拟,建立图像非盲求解模型:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><munder><mi>min</mi><mi>u</mi></munder><mfrac><msub><mi>&lambda;</mi><mn>2</mn></msub><mn>2</mn></mfrac><msubsup><mrow><mo>|</mo><mo>|</mo><mi>u</mi><mo>&CircleTimes;</mo><mi>k</mi><mo>-</mo><mi>f</mi><mo>|</mo><mo>|</mo></mrow><mn>2</mn><mn>2</mn></msubsup><mo>+</mo><msubsup><mrow><mo>|</mo><mo>|</mo><mo>&dtri;</mo><mi>u</mi><mo>|</mo><mo>|</mo></mrow><mi>q</mi><mi>q</mi></msubsup><mo>,</mo></mrow>]]></math><img file="FDA0000539471870000021.GIF" wi="735" he="145" /></maths>(4b)引入辅助项d=▽u,将变量u与梯度变换分离,进行分步求解,定义<img file="FDA0000539471870000027.GIF" wi="430" he="67" />其中i=1,2,且f<sub>1</sub>=[1,‑1],f<sub>2</sub>=[1;‑1],并加入权衡参数β,整理模型为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><munder><mi>min</mi><mrow><mi>u</mi><mo>,</mo><mi>d</mi></mrow></munder><mfrac><msub><mi>&lambda;</mi><mn>2</mn></msub><mn>2</mn></mfrac><msubsup><mrow><mo>|</mo><mo>|</mo><mi>u</mi><mo>&CircleTimes;</mo><mi>k</mi><mo>-</mo><mi>f</mi><mo>|</mo><mo>|</mo></mrow><mn>2</mn><mn>2</mn></msubsup><mo>+</mo><mfrac><mi>&beta;</mi><mn>2</mn></mfrac><msubsup><mrow><msubsup><mrow><mo>|</mo><mo>|</mo><msub><mi>F</mi><mi>i</mi></msub><mi>u</mi><mo>-</mo><mi>d</mi><mo>|</mo><mo>|</mo></mrow><mn>2</mn><mn>2</mn></msubsup><mo>+</mo><mrow><mo>|</mo><mo>|</mo><mi>d</mi><mo>|</mo><mo>|</mo></mrow></mrow><mi>q</mi><mi>q</mi></msubsup><mo>,</mo></mrow>]]></math><img file="FDA0000539471870000022.GIF" wi="1059" he="145" /></maths>将上式对u求微分,并用二维快速傅里叶法求得u的最优解:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>u</mi><mo>=</mo><msup><mi>F</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><mfrac><mrow><mi>F</mi><mrow><mo>(</mo><mfrac><mi>&beta;</mi><msub><mi>&lambda;</mi><mn>2</mn></msub></mfrac><msup><msub><mi>F</mi><mi>i</mi></msub><mi>T</mi></msup><mi>d</mi><mo>)</mo></mrow><mo>+</mo><mi>F</mi><msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>*</mo></msup><mi>eF</mi><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow></mrow><mrow><mfrac><mi>&beta;</mi><msub><mi>&lambda;</mi><mn>2</mn></msub></mfrac><mrow><mo>(</mo><mi>F</mi><mrow><mo>(</mo><msup><msub><mi>F</mi><mi>i</mi></msub><mi>T</mi></msup><msub><mi>F</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow><mo>+</mo><mi>F</mi><msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>*</mo></msup><mi>eF</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow></mfrac><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000539471870000023.GIF" wi="1044" he="360" /></maths>将上式对d求微分,利用封闭式的阈值公式去卷积的方法,选取q=1/2时的阈值公式:<img file="FDA0000539471870000024.GIF" wi="1318" he="488" />其中,<img file="FDA0000539471870000028.GIF" wi="740" he="233" /><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>p</mi><mrow><mo>(</mo><msub><mi>&lambda;</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mroot><mn>54</mn><mn>3</mn></mroot><mn>4</mn></mfrac><msup><mrow><mo>(</mo><msub><mi>&lambda;</mi><mn>2</mn></msub><mo>)</mo></mrow><mfrac><mn>2</mn><mn>3</mn></mfrac></msup><mo>,</mo><mi>a</mi><mo>=</mo><mo>&dtri;</mo><mi>u</mi><mo>;</mo></mrow>]]></math><img file="FDA0000539471870000029.GIF" wi="710" he="231" /></maths>变量u和v的交互迭代,迭代条件为:β的初始值设为1,β的最大值设为256,β在迭代过程中以<img file="FDA0000539471870000026.GIF" wi="91" he="86" />的倍数递增,直到β>256结束,最终恢复出清晰图像。
地址 529000 广东省江门市蓬江区东成村22号