一种块效应抑制的降质模糊图像非盲复原方法技术

本发明专利技术公布了一种块效应抑制的降质模糊图像非盲复原方法。经过低比特率块编码压缩重构后的降质图像边界会出现明显的块效应，常规的图像复原方法往往难以得到理想的复原结果。本发明专利技术提出了一种可以实现JPEG高压缩率下的模糊图像复原方法。在全变分(total variation TV)正则化约束的基础上，增加新的抑制块效应的正则化L2范数约束，保证了在复原图像的过程中逐步抑制块效应带来的负面影响。本发明专利技术能够得到兼顾消除块效应与保留图像细节的复原效果。

【技术实现步骤摘要】

本专利技术涉及计算机成像
，主要涉及模糊图像经过图像压缩后的图像复原问题。
技术介绍
在成像过程中，由于各种因素，比如光学器件缺陷，人为的抖动，大气非均匀性等会造成图像有某些程度的失真和不同程度的降质。相机的成像过程可以描述为理想清晰图像与其他外部影响而造成的模糊核的卷积。人们根据不同的图像退化成因研究了相应的图像复原方法。也相应取得了较为显著的图像质量提升结果。图像复原是二维卷积运算的逆过程，是一个病态问题。从复原方法上，分为已知模糊核的非盲复原算法和未知模糊核的盲复原算法。本专利技术主要针对已知模糊核的的非盲复原问题。同时，在很多情况下，受限于存储空间或者实时性等的要求，拍摄后的图像在保存时会为了节省空间而采用有损压缩技术。其中，被最为广泛使用的是JPEG保存格式，这种压缩格式经过低比特率块编码压缩重构后的图像8*8小块边界往往会出现明显的块状效应。这相当于在由于运动等造成的图像模糊退化基础上又进行了进一步的退化，由于块效应破坏了图像的先验信息，直接采用常规的图像复原方法往往难以得到理想的复原结果。对于这一类有块效应存在的模糊图像，无论是在复原前或复原后使用去块效应方法，获得的复原图像结果均有提升空间。针对这一问题尚未有人进行研究，本专利技术首次提出了将去块效应算法融合在图像复原的正则化优化框架中，得到了较为明显的复原效果提升。
技术实现思路
本专利技术公布了一种基于块效应抑制的压缩降质模糊图像非盲复原方法，其特征在于：在正则化约束的图像复原框架基础上，引入逐点形状自适应的DCT(离散余弦变换)去块效应滤波器。在复原迭代的过程中，增加新的抑制块效应的正则化L2范数约束，保证了在复原图像的过程中逐步抑制块效应带来的负面影响。具体采用以下步骤对具有压缩块的模糊图像进行图像复原：(1)输入一幅待复原的有压缩块效应的模糊图像b，并且输入该模糊图像对应的模糊核h。(2)获取压缩图像b的量化表Q。(3)对待求解的复原图f进行初次迭代求解。 f = argmin f μ 2 || h * f - b || 2 2 + || D x f || 1 + || D y f || 1 + ξ 2 || f - f 0 || 2 2 - - - ( 1 ) ]]>其中，Dx,Dy分别为水平方向与竖直方向的梯度算子，作用在f后求解出图像梯度。μ,ξ为正则化系数，f0表示迭代过程中的去块效应图像。分别表示图像的L1与L2范数。为了求解此问题，采用变量分离技术，将式(1)引入两个辅助变量，转化为可以方便求解的形式，如下式所示 argmin w 1 , w 2 , f μ 2 || h * f - b || 2 2 + Σ i = 1 2 || w i || 1 + β 2 || w 1 - D y f || 2 2 + β 2 || w 2 - D x f || 2 2 + ξ 2 || f - f 0 || 2 2 - - - ( 2 ) ]]>其中，w1,w2为引入的辅助变量。β为迭代过程中的正则化系数，根据本领域的公知常识可知，当β趋向于很大时，式(2)收敛于式(1)。因此，上述问题可以转化为轮转迭代求解，具体的求解过程为：(3.1)用步骤1输入的模糊图像b作为初始的清晰图像f输入。求解如下问题 argmin w 1 , w 2 Σ i = 1 本文档来自技高网...

【技术保护点】
一种块效应抑制的降质模糊图像非盲复原方法，其特征在于，具体可以采用以下步骤对存在较高压缩比的降质模糊图像进行复原：(1)输入一幅待复原的有压缩块效应的模糊图像b，并且输入该模糊图像对应的模糊核h。(2)获取压缩图像b的量化表Q。(3)对待求解的复原图f进行初次迭代求解。f=argfminμ2||h*f-b||22+||Dxf||1+||Dyf||1+ξ2||f-f0||22---(1)]]>其中，Dx,Dy分别为水平方向与竖直方向的梯度算子，作用在复原图f后表示图像两个方向的梯度。μ,ξ为正则化系数，f0表示迭代过程中的去块效应图像。分别表示图像的L1与L2范数。为了求解此问题，采用变量分离技术，将式(1)引入两个辅助变量，转化为可以方便求解的形式，如下式所示argminw1,w2,fμ2||h*f-b||22+Σi=12||wi||1+β2||w1-Dyf||22+β2||w2-Dxf||22+ξ2||f-f0||22---(2)]]>其中，w1,w2为引入的两个辅助变量。β为迭代过程中的正则化系数，求解过程如下：(3.1)用步骤1输入的模糊图像b作为初始的清晰图像f输入。求解如下问题argminw1,w2Σi=12||wi||1+β2||w1-Dyf||22+β2||w2-Dxf||22---(3)]]>该问题的解为w1=max(||Dxf||-1β)Dxf||Dxf||,w2=max(||Dyf||-1β)Dyf||Dyf||---(4)]]>(3.2)用步骤1输入的模糊图像b作为初始的去块效应图像f0输入，并用式(4)得到的wi(i＝1,2)，求解如下问题argminw1,w2,fμ2||h*f-b||22+β2||w1-Dyf||22+β2||w2-Dxf||22+ξ2||f-f0||22---(5)]]>该问题的解为其中，分别表示图像的傅里叶变换和傅里叶逆变换，“*”表示共轭运算，算符“ο”表示两个矩阵的点乘积。得到的f为下次迭代时(3.1)的输入。(3.3)利用(3.2)计算得到的f，结合步骤2得到的量化表Q，采用逐点形状自适应DCT滤波器，对其进行块效应消除，得到的图像为下次迭代时式(6)的f0。(4)利用得到的f,f0,w1,w2，并更新β值为其原值的2倍，重复步骤(3.1)～(3.3)，不断得到优化的清晰图像f，经过6～10次的轮转迭代后，得到最终的复原结果。...

【技术特征摘要】
1.一种块效应抑制的降质模糊图像非盲复原方法，其特征在于，具体可以采用以下步骤对存在较高压缩比的降质模糊图像进行复原：(1)输入一幅待复原的有压缩块效应的模糊图像b，并且输入该模糊图像对应的模糊核h。(2)获取压缩图像b的量化表Q。(3)对待求解的复原图f进行初次迭代求解。 f = arg f m i n μ 2 | | h * f - b | | 2 2 + | | D x f | | 1 + | | D y f | | 1 + ξ 2 | | f - f 0 | | 2 2 - - - ( 1 ) ]]>其中，Dx,Dy分别为水平方向与竖直方向的梯度算子，作用在复原图f后表示图像两个方向的梯度。μ,ξ为正则化系数，f0表示迭代过程中的去块效应图像。分别表示图像的L1与L2范数。为了求解此问题，采用变量分离技术，将式(1)引入两个辅助变量，转化为可以方便求解的形式，如下式所示 argmin w 1 , w 2 , f μ 2 | | h * f - b | | 2 2 + Σ i = 1 2 | | w i | | 1 + β 2 | | w 1 - D y f | | 2 2 + β 2 | | w 2 - D x f | | 2 2 + ξ 2 | | f - f 0 | | 2 2 - - - ( 2 ) ]]>其中，w1,w2为引入的两个辅助变量。β为迭代过程中的正则化系数，求解过程如下：(3.1)用步骤1输入的模糊图像b作为初始的清晰图像f输入。求解如下问题 argmin w 1 , w 2 Σ i = 1 2 | | w i | | 1 + β 2 | | w 1 - D y f | | 2 2 + β 2 | | w 2 - D x f | | 2 2 - - - ( 3 ) ]]>该问题的解为 w 1 = m a x ( | | D ...

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