CN103413277A - 基于l0稀疏先验的盲相机抖动去模糊方法 - Google Patents
基于l0稀疏先验的盲相机抖动去模糊方法 Download PDFInfo
- Publication number
- CN103413277A CN103413277A CN2013103629330A CN201310362933A CN103413277A CN 103413277 A CN103413277 A CN 103413277A CN 2013103629330 A CN2013103629330 A CN 2013103629330A CN 201310362933 A CN201310362933 A CN 201310362933A CN 103413277 A CN103413277 A CN 103413277A
- Authority
- CN
- China
- Prior art keywords
- old
- fuzzy core
- new
- camera shake
- lambda
- 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
Images
Landscapes
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
基于L0稀疏先验的盲相机抖动去模糊方法属于数字图像处理技术领域,是一种对相机抖动造成的模糊图像进行去模糊的方法,能够估计各种空间不变的相机抖动模糊核(点扩散函数);该方法不仅避免了当前变分贝叶斯估计方法计算复杂度高的问题,而且解决了当前最大后验估计方法缺乏严格最优化理论支撑的问题;首先,引入基于L0范数的显著边缘稀疏先验,利用迭代硬阈值收缩法实现显著边缘特征的隐性自动预测;其次,引入基于L0范数的相机抖动模糊核稀疏先验,利用迭代重新加权最小二乘法实现模糊核的快速估计;最后,利用基于超拉普拉斯先验的图像非盲去模糊方法,获得高质量去模糊图像;该方法流程图如附图1所示。
Description
技术领域
本发明属于数字图像处理技术领域,具体涉及对相机抖动造成的模糊图像进行去模糊的领域。
背景技术
盲相机抖动去模糊是指去除或减轻相机拍摄数字图像过程中发生的图像运动模糊现象,是近年来数字图像处理中非常重要又极具挑战性的研究内容,其核心是估计对应各种相机抖动的模糊核(点扩散函数)。当前,研究人员已经提出了多种不同方法。
最经典的盲相机抖动去模糊方法是的最大似然估计法,已集成到数值计算软件MATLAB中的图像处理工具箱,但是该方法在实际应用中具有明显的局限性,不仅要求模糊图像具有较高的信噪比,而且要求模糊核的尺寸较小。该方法的另一个缺点是,去模糊图像中的边缘等细节信息不能很好地恢复,常伴有明显的振铃效应。
目前,国际上提出了多种盲相机抖动去模糊的新方法。根据模糊核的估计准则,这些方法主要分为两大类:变分贝叶斯估计方法和最大后验估计方法。变分贝叶斯估计方法的计算复杂度较高,其中的两种代表性方法是:Fergus等人提出了基于高斯混合模型的变分贝叶斯盲相机抖动去模糊方法,参见文献《Removing camera shake from a single photograph》,ACM Trans.Graph.,2006,vol.25,no.3,pp.787–794;之后,Levin等人基于类似的图像建模思想降低了上述方法的计算复杂度,参见文献《Efficient marginallikelihood optimization in blind deconvolution》,IEEE Conf.Computer Vision and Pattern Recognition,2011,pp.2657–2664。相比变分贝叶斯估计方法,盲相机抖动去模糊最大后验估计方法的计算复杂度要低得多。代表性方法参见文献《Fast motion deblurring》,ACM Transactions on Graphics(SIGGRAPH ASIA),2009,vol.28,no.5,article no.145,《Psf estimation using sharp edge prediction》,IEEE Conf.Computer Vision and PatternRecognition,2008,pp.1–8。但是,当前的最大后验估计方法往往先要利用平滑滤波器与冲击滤波器预测图像边缘信息,然后再利用预测边缘图像估计模糊核,经过这两个过程的反复迭代最终估计获得模糊核,缺乏严格的最优化理论支撑,所以无法保证解的(局部)最优性。
发明内容
本发明的目的在于解决现有技术方法在理论和实现两个方面的不足,提出一种基于L0稀疏先验的最大后验模糊核估计新方法。
本发明技术方案是:首先,引入基于L0范数的显著边缘稀疏先验,利用迭代硬阈值收缩法实现显著边缘特征的隐性自动预测;其次,引入基于L0范数的相机抖动模糊核稀疏先验,利用迭代重新加权最小二乘法实现模糊核的快速估计;最后,利用基于超拉普拉斯先验的图像非盲去模糊方法,实现高质量图像去模糊。
1、本发明具体实施步骤:
(1)利用水平方向和垂直方向的一阶导数算子▽h=[1,-1;0,0],▽v=[1,0;-1,0],得到相机抖动模糊图像y的梯度图像
(2)令待估计模糊核k的尺寸为Z×Z,为提高模糊核估计方法的收敛性,采用多尺度实现方式迭代估计模糊核;
(3)令初始模糊核k(1)的尺寸为3×3,且k(1)=[0,0,0;1,1,1;0,0,0]3,根据(3.1)-(3.12)确定尺度总数scales_num及各尺度下模糊核k(s)的尺寸ksize(s)×ksize(s):
(3.1) s=1;
(3.2) tmp=3;
(3.3) while(tmp<Z)
(3.4) ksize(s)=tmp;
(3.5) s=s+1;
(3.6) tmp=ceil(tmp*sqrt(2));
(3.7) if(mod(tmp,2)==0)
(3.8) tmp=tmp+1;
(3.9) end;
(3.10) end;
(3.11) ksize(s)=Z;
(3.12) scales_num=s;
(4.1) r(s)=floor(size(yd,1)*ksize(s)/Z);
(4.2) c(s)=floor(size(yd,2)*ksize(s)/Z);
(4.3)
(5)令显著边缘图像与模糊核在各尺度下的循环迭代次数为out_iter_num,循环迭代初始次数为out_iter=1,初始尺度为s=1,显著边缘L0稀疏先验的正则化参数为η,模糊核L0稀疏先验的正则化参数为γ,Kold对应模糊核k(1)的二维矩阵,对应模糊插值图像的一维向量,kold对应模糊核k(1)的一维向量;
(6.1)
(6.2)
(6.3)τnew=τold2;
(6.4)
其中,硬阈值算子ΦHARD(·,·)定义如下:
(7.1)估计
(7.2)将knew投影到约束集
(7.3)更新
(7.4)kold=knew;
(8)更新out_iter:out_iter=out_iter+1;
(9)如果out_iter<out_iter_num,转到(6),否则转到(10);
(10)更新第s+1个尺度下模糊核k的初始值:k(s+1)=imresize(k(s),[ksize(s),ksize(s)],'bilinear');
(11)更新s:s=s+1;
(12)如果s<scales_num,返回(6),否则转到(13);
有益效果:
(1)发明方法的模糊核估计是个严格意义上的稀疏最优化问题;
(2)发明方法的实现方式简单,无需平滑滤波、冲击滤波等任何迭代预处理;
(3)发明方法的计算复杂度低;
(4)发明方法的模糊核估计准确度高;
估计出相机抖动模糊核之后,利用基于超拉普拉斯先验的图像非盲去模糊方法,从而获得更高质量的去模糊图像。
2、本发明模型推导过程:
不失一般性,相机抖动模糊可利用如下卷积型观察模型进行表述
y=Kf+n
其中,y,f,以及n分别是y,f,以及n的向量表达形式,K是模糊核k的矩阵表达形式。采取分而治之的策略解决相机抖动去模糊问题,分成两大步骤:(1)模糊核估计;(2)非盲图像恢复;
模糊核估计是在图像梯度域实现的,为此,建立如下梯度域观察模型:
其中, 类似的,上述梯度域卷积型观察模型可以写成如下矩阵-向量表达形式
yd=Kfd+nd=Fdk+nd
其中,yd,fd,nd,k分别是yd,fd,nd,k的向量表达形式,Fd是fd的矩阵表达形式。
图像中的显著边缘是精确估计模糊核的重要所在。为了便于自动预测显著边缘,提高模糊核估计的准确性,首先提出基于L0范数的显著边缘先验及其优化模型:
根据迭代硬阈值收缩方法,上述最优化问题可以下步骤(1)-(3)进行迭代求解:
(2)τnew=τold2,τold对应当前的调整参数,τnew代表更新的调整参数;
(3)τold=τnew;
当显著边缘给定时,利用相机抖动模糊核的稀疏特性,从而能够建立基于L0范数的模糊核先验及其优化模型:
其中,Ξ(k)=k0是基于L0范数的模糊核先验,kest是对应k的最优化结果,γ是正则化参数;
根据迭代重新加权最小二乘方法,首先上述将最优化问题转化为如下近似形式:
其中,W=diag({1(kl+0.0001)2},然后利用以下步骤(4)-(6)进行迭代求解:
(5)进一步将knew投影到约束集
(6)更新
附图说明
附图1.发明方法流程图;
附图2.仿真实验的真实图像Lena及模糊核;
附图3.仿真实验的真实图像Cameraman及模糊核;
附图4.仿真实验的真实图像House及模糊核;
附图5.仿真实验的真实图像Boat及模糊核;
附图6.对应图像Lena的相机抖动模糊图;
附图7.基于规范化稀疏度量盲去模糊方法的恢复图像及估计模糊核;
附图8.基于L0稀疏先验的盲去模糊方法的恢复图像及估计模糊核;
附图9.基于L0稀疏先验的盲去模糊方法的显著边缘图像。
具体实施方式
(1)利用水平方向和垂直方向的一阶导数算子▽h=[1,-1;0,0],▽v=[1,0;-1,0],利用MATLAB函数conv2获得相机抖动模糊图像y的梯度图像
yh=conv2(y,▽h,'valid');
yv=conv2(y,▽v,'valid');
(2)令待估计模糊核k的尺寸为Z×Z,为提高模糊核估计方法的收敛性,采用多尺度实现方式迭代估计模糊核;
(3)令初始模糊核k(1)=[0,0,0;1,1,1;0,0,0]3(尺寸为3×3),并且利用以下MATLAB代码确定尺度总
数scales_num及各尺度下模糊核k(s)的尺寸ksize(s)×ksize(s):
s=1;
tmp=3;
while(tmp<Z)
ksize(s)=tmp;
s=s+1;
tmp=ceil(tmp*sqrt(2));
if(mod(tmp,2)==0)
tmp=tmp+1;
end;
end;
ksize(s)=Z;
scales_num=s;
r(s)=floor(size(yd,1)*ksize(s)/Z);
c(s)=floor(size(yd,2)*ksize(s)/Z);
(5)令显著边缘图像与模糊核在各尺度下的循环迭代次数为out_iter_num,循环迭代初始次数为out_iter=1,初始尺度为s=1,显著边缘L0稀疏先验的正则化参数为η,模糊核L0稀疏先验的正则化参数为γ,Kold对应模糊核k(1)的二维矩阵,kold对应模糊核k(1)的一维向量,对应模糊插值图像的一维向量;
(6)令τold=1,利用迭代硬阈值收缩法估计第s个尺度下的显著边缘(迭代硬阈值收缩法具体参见文献《Iterative hard thresholding for compressed sensing》,(T.Blumensath,M.E.Davies),Appliedand Computational Harmonic Analysis,2009,vol.27,no.3,pp.265–274)。具体利用以下MATLAB代码进行m次循环估计:
τnew=τold/2;
τold=τnew;
(7)令权矩阵利用迭代重新加权最小二乘方法估计第s个尺度下的模糊核k(s)=knew。迭代重新加权最小二乘方法可参见文献《Iteratively reweighted least squares minimization forsparse recovery》(I.Daubechies,et al.),Communications on Pure and Applied Mathematics,2010,vol.63,pp.1–38。具体利用以下步骤进行n次循环估计:
(7.1)利用MATLAB函数pcg(共轭梯度法)求解如下方程组,估计knew:
共轭梯度法可参见参考文献《Convex optimization》(S.Boyd,L.Vandenberghe),Cambridge University Press,2004。
knew(knew<0)=0;
sumk=sum(knew(:));
knew=knew./sumk;
(7.3)更新
(7.4)kold=knew;
(8)更新out_iter:out_iter=out_iter+1;
(9)如果out_iter<out_iter_num,转到(6),否则转到(10);
(10)更新第s+1个尺度下模糊核k的初始值:k(s+1)=imresize(k(s),[ksize(s),ksize(s)],'bilinear');
(11)更新s:s=s+1;
(12)如果s<scales_num,返回(6),否则转到(13);
(14)利用基于超拉普拉斯先验的图像非盲去模糊方法恢复出去模糊图像具体参见文献《Fast imagedeconvolution using hyper-laplacian priors》(D.Krishnan and R.Fergus),Advances in Neural InformationProcessing Systems,2009,vol.22,pp.1033-1041.
(15)为了验证基于L0稀疏先验的盲相机抖动去模糊方法的有效性,设计4组仿真实验,每组实验均提供了真实图像及真实模糊核(附图2、附图3、附图4、附图5),并与最大似然估计盲去模糊方法(简写为MaxL)以及基于规范化稀疏度量的盲去模糊方法(简写为NormSpar)进行比较,MaxL可调用MATLAB图像处理工具箱中的函数deconvblind进行实现,NormSpar具体参见文献《Blind deconvolution using a normalizedsparsity measure》(D.Krishnan,et al.),IEEE Conf.Computer Vision and Pattern Recognition,2011:233–240。为了衡量不同方法的去模糊效果,采用峰值信噪比(简称PSNR)和结构相似性度量(简称SSIM)这两个客观标准进行度量。所有方法的参数均采用人工方式选取,以取得最佳PSNR、SSIM、以及视觉效果为准。此外,具体的实验环境包括:软件MATLAB版本号为v7.0,电脑的配置为Pentium(R)Core-Duo T4200CPU(2GHZ),2GB RAM,Microsoft Windows XP(version2002,Service Pack3).
(17)SSIM的具体定义与实现参见文献《Image quality assessment:from error measurement to structuralsimilarity》(Z.Wang,et al.),IEEE Transactios on Image Processing,2004,vol.13,no.4,pp.600–612.
(18)表1提供了上述3种盲去模糊方法在4组仿真实验中的PSNR和SSIM。由表1可见,本发明在所有中均取得了最高PSNR,特别是对于图像Lena、House、Boat,本发明的优势明显;从SSIM的角度看,本发明在对应图像Lena、House、Boat的3组实验中的去模糊效果优于MaxL和NormS,对图像Cameraman本发明取得了与NormS基本相似的去模糊效果,且两者均优于MaxL。
表1.各种盲去模糊方法在4组仿真实验中的PSNR和SSIM
(19)从视觉感知的角度看,本发明也取得了最优的去模糊效果。附图6给出了对应图像Lena的相机抖动模糊图;附图7给出了基于规范化稀疏度量盲去模糊方法的恢复图像及估计模糊核;附图8给出了本发明基于L0稀疏先验的盲去模糊方法的恢复图像及估计模糊核。容易看出,本发明估计的模糊核更接近真实模糊核,且恢复的图像具有更清晰的边缘和较少的振铃效应;而基于规范化稀疏度量盲去模糊方法的恢复图像与估计模糊核与真实情况相去甚远。此外,附图9给出了基于L0稀疏先验的盲去模糊方法的显著边缘图像。
以上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术人员应了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内,本发明要求保护范围由所附的权利要求书其等效物界定。
Claims (3)
1.一种基于L0稀疏先验的盲相机抖动去模糊方法,其特征在于,首先,引入基于L0范数的显著边缘稀疏先验,利用迭代硬阈值收缩法实现显著边缘特征的隐性自动预测;其次,引入基于L0范数的相机抖动模糊核稀疏先验,利用迭代重新加权最小二乘法实现模糊核的快速估计;最后,利用基于超拉普拉斯先验的图像非盲去模糊方法,实现高质量图像去模糊。
2.根据权利要求1所述的,基于L0稀疏先验的盲相机抖动去模糊方法,具体实施步骤:
(1)利用水平方向和垂直方向的一阶导数算子▽h=[1,-1;0,0],▽v=[1,0;-1,0],得到相机抖动模糊图像y的梯度图像
(2)令待估计模糊核k的尺寸为Z×Z,为提高模糊核估计方法的收敛性,采用多尺度实现方式迭代估计模糊核;
(3)令初始模糊核k(1)的尺寸为3×3,且k(1)=[0,0,0;1,1,1;0,0,0]3,利用以下MATLAB代码(3.1)-(3.12)确定尺度总数scales_num及各尺度下模糊核k(s)的尺寸ksize(s)×ksize(s):
(3.1) s=1;
(3.2) tmp=3;
(3.3) while(tmp<Z)
(3.4) ksize(s)=tmp;
(3.5) s=s+1;
(3.6) tmp=ceil(tmp*sqrt(2));
(3.7) if(mod(tmp,2)==0)
(3.8) tmp=tmp+1;
(3.9) end;
(3.10) end;
(3.11) ksize(s)=Z;
(3.12) scales_num=s;
(4.1) r(s)=floor(size(yd,1)*ksize(s)/Z);
(4.2) c(s)=floor(size(yd,2)*ksize(s)/Z);
(4.3)
(5)令显著边缘图像与模糊核在各尺度下的循环迭代次数为out_iter_num,循环迭代初始次数为out_iter=1,初始尺度为s=1,显著边缘L0稀疏先验的正则化参数为η,模糊核L0稀疏先验的正则化参数为γ,Kold对应模糊核k(1)的二维矩阵,对应模糊插值图像的一维向量,kold对应模糊核k(1)的一维向量;
(6.1)
(6.2)
(6.3)τnew=τold2;
(6.4)
其中,硬阈值算子ΦHARD(·,·)定义如下:
(7)令权矩阵利用迭代重新加权最小二乘方法估计第s个尺度下的模糊核k(s)=knew,具体而言利用步骤(7.1)-(7.4)进行n次循环估计:
(7.1)估计
(7.2)将knew投影到约束集
(7.3)更新
(7.4)kold=knew;
(8)更新out_iter:out_iter=out_iter+1;
(9)如果out_iter<out_iter_num,转到(6),否则转到(10);
(10)更新第s+1个尺度下模糊核k的初始值:k(s+1)=imresize(k(s),[ksize(s),ksize(s)],'bilinear');
(11)更新s:s=s+1;
(12)如果s<scales_num,返回(6),否则转到(13);
3.根据权利要求1所述的,基于L0稀疏先验的盲相机抖动去模糊方法,模型推导过程:不失一般性,相机抖动模糊可利用如下卷积型观察模型进行表述
y=Kf+n
其中,y,f,以及n分别是y,f,以及n的向量表达形式,K是模糊核k的矩阵表达形式。采取分而治之的策略解决相机抖动去模糊问题,分成两大步骤:(1)模糊核估计;(2)非盲图像恢复;
模糊核估计是在图像梯度域实现的,为此,建立如下梯度域观察模型:
其中, 类似的,上述梯度域卷积型观察模型可以写成如下矩阵-向量表达形式
yd=Kfd+nd=Fdk+nd
其中,yd,fd,nd,k分别是yd,fd,nd,k的向量表达形式,Fd是fd的矩阵表达形式;
图像中的显著边缘是精确估计模糊核的重要所在。为了便于自动预测显著边缘,提高模糊核估计的准确性,首先提出基于L0范数的显著边缘先验及其优化模型:
根据迭代硬阈值收缩方法,上述最优化问题可以下步骤(1)-(3)进行迭代求解:
(2)τnew=τold2,τold对应当前的调整参数,τnew代表更新的调整参数;
(3)τold=τnew;
当显著边缘给定时,利用相机抖动模糊核的稀疏特性,从而能够建立基于L0范数的模糊核先验及其优化模型:
根据迭代重新加权最小二乘方法,首先上述将最优化问题转化为如下近似形式:
其中,W=diag({1(kl+0.0001)2},然后利用以下步骤(4)-(6)进行迭代求解:
(6)更新
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310362933.0A CN103413277B (zh) | 2013-08-19 | 2013-08-19 | 基于l0稀疏先验的盲相机抖动去模糊方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310362933.0A CN103413277B (zh) | 2013-08-19 | 2013-08-19 | 基于l0稀疏先验的盲相机抖动去模糊方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103413277A true CN103413277A (zh) | 2013-11-27 |
CN103413277B CN103413277B (zh) | 2017-04-05 |
Family
ID=49606282
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310362933.0A Active CN103413277B (zh) | 2013-08-19 | 2013-08-19 | 基于l0稀疏先验的盲相机抖动去模糊方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103413277B (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103761710A (zh) * | 2014-01-08 | 2014-04-30 | 西安电子科技大学 | 基于边缘自适应的高效图像盲去模糊方法 |
CN104036473A (zh) * | 2014-05-30 | 2014-09-10 | 南京邮电大学 | 基于分裂Bregman迭代的快速鲁棒图像运动去模糊方法 |
CN104605853A (zh) * | 2015-01-21 | 2015-05-13 | 北京师范大学 | 一种基于脑功能特征的实时神经解码系统 |
CN105493140A (zh) * | 2015-05-15 | 2016-04-13 | 北京大学深圳研究生院 | 图像去模糊方法及系统 |
CN105809642A (zh) * | 2016-03-11 | 2016-07-27 | 中山大学 | 一种基于l0正则化的自然图像盲去运动模糊的方法 |
CN105957025A (zh) * | 2016-04-21 | 2016-09-21 | 天津大学 | 基于稀疏表示的非一致图像盲复原方法 |
CN106097256A (zh) * | 2016-05-31 | 2016-11-09 | 南京邮电大学 | 一种基于图像盲去模糊的视频图像模糊度检测方法 |
CN106530251A (zh) * | 2016-11-03 | 2017-03-22 | 北京航星机器制造有限公司 | 一种盲卷积运动模糊图像复原方法 |
CN107146202A (zh) * | 2017-03-17 | 2017-09-08 | 中山大学 | 基于l0正则化和模糊核后处理的图像盲去模糊的方法 |
CN107688198A (zh) * | 2017-07-25 | 2018-02-13 | 中国科学院地质与地球物理研究所 | 地震数据重构方法及装置、计算机可读存储介质、设备 |
CN113989128A (zh) * | 2020-07-31 | 2022-01-28 | 山东财经大学 | 基于改进的l0梯度最小化模型的图像平滑方法及系统 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070009169A1 (en) * | 2005-07-08 | 2007-01-11 | Bhattacharjya Anoop K | Constrained image deblurring for imaging devices with motion sensing |
CN102073993A (zh) * | 2010-12-29 | 2011-05-25 | 清华大学 | 一种基于摄像机自标定的抖动视频去模糊方法和装置 |
CN102208100A (zh) * | 2011-05-31 | 2011-10-05 | 重庆大学 | 基于Split Bregman 迭代的全变差正则化图像盲复原方法 |
CN102354395A (zh) * | 2011-09-22 | 2012-02-15 | 西北工业大学 | 基于稀疏表示的模糊图像盲复原方法 |
CN102708550A (zh) * | 2012-05-17 | 2012-10-03 | 浙江大学 | 一种基于自然图像统计特性的盲去模糊算法 |
CN102800055A (zh) * | 2012-06-28 | 2012-11-28 | 西安电子科技大学 | 图像盲去模糊的低秩分解方法 |
-
2013
- 2013-08-19 CN CN201310362933.0A patent/CN103413277B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070009169A1 (en) * | 2005-07-08 | 2007-01-11 | Bhattacharjya Anoop K | Constrained image deblurring for imaging devices with motion sensing |
CN102073993A (zh) * | 2010-12-29 | 2011-05-25 | 清华大学 | 一种基于摄像机自标定的抖动视频去模糊方法和装置 |
CN102208100A (zh) * | 2011-05-31 | 2011-10-05 | 重庆大学 | 基于Split Bregman 迭代的全变差正则化图像盲复原方法 |
CN102354395A (zh) * | 2011-09-22 | 2012-02-15 | 西北工业大学 | 基于稀疏表示的模糊图像盲复原方法 |
CN102708550A (zh) * | 2012-05-17 | 2012-10-03 | 浙江大学 | 一种基于自然图像统计特性的盲去模糊算法 |
CN102800055A (zh) * | 2012-06-28 | 2012-11-28 | 西安电子科技大学 | 图像盲去模糊的低秩分解方法 |
Non-Patent Citations (4)
Title |
---|
LI XU 等: "Unnatural L0 Sparse Representation for Natural Image Deblurring", 《CVPR 2013》, 28 June 2013 (2013-06-28) * |
NEEL JOSHI 等: "Image Deblurring and Denoising using Color Priors", 《CVPR 2009》, 25 June 2009 (2009-06-25) * |
QI SHAN 等: "High-quality Motion Deblurring from a Single Image", 《ACM TRANSACTIONS ON GRAPHICS》, vol. 27, no. 3, 31 August 2008 (2008-08-31), XP 055073383, DOI: doi:10.1145/1360612.1360672 * |
王爱齐 等: "图像去模糊的l0范数最小化方法", 《信号处理》, vol. 28, no. 11, 30 November 2012 (2012-11-30) * |
Cited By (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103761710B (zh) * | 2014-01-08 | 2016-07-06 | 西安电子科技大学 | 基于边缘自适应的高效图像盲去模糊方法 |
CN103761710A (zh) * | 2014-01-08 | 2014-04-30 | 西安电子科技大学 | 基于边缘自适应的高效图像盲去模糊方法 |
CN104036473B (zh) * | 2014-05-30 | 2018-01-05 | 南京邮电大学 | 基于分裂Bregman迭代的快速鲁棒图像运动去模糊方法 |
CN104036473A (zh) * | 2014-05-30 | 2014-09-10 | 南京邮电大学 | 基于分裂Bregman迭代的快速鲁棒图像运动去模糊方法 |
CN104605853A (zh) * | 2015-01-21 | 2015-05-13 | 北京师范大学 | 一种基于脑功能特征的实时神经解码系统 |
CN105493140A (zh) * | 2015-05-15 | 2016-04-13 | 北京大学深圳研究生院 | 图像去模糊方法及系统 |
WO2016183716A1 (zh) * | 2015-05-15 | 2016-11-24 | 北京大学深圳研究生院 | 图像去模糊方法及系统 |
CN105493140B (zh) * | 2015-05-15 | 2018-02-02 | 北京大学深圳研究生院 | 图像去模糊方法及系统 |
CN105809642A (zh) * | 2016-03-11 | 2016-07-27 | 中山大学 | 一种基于l0正则化的自然图像盲去运动模糊的方法 |
CN105809642B (zh) * | 2016-03-11 | 2018-10-09 | 中山大学 | 一种基于l0正则化的自然图像盲去运动模糊的方法 |
CN105957025A (zh) * | 2016-04-21 | 2016-09-21 | 天津大学 | 基于稀疏表示的非一致图像盲复原方法 |
CN106097256B (zh) * | 2016-05-31 | 2019-01-01 | 南京邮电大学 | 一种基于图像盲去模糊的视频图像模糊度检测方法 |
CN106097256A (zh) * | 2016-05-31 | 2016-11-09 | 南京邮电大学 | 一种基于图像盲去模糊的视频图像模糊度检测方法 |
CN106530251A (zh) * | 2016-11-03 | 2017-03-22 | 北京航星机器制造有限公司 | 一种盲卷积运动模糊图像复原方法 |
CN106530251B (zh) * | 2016-11-03 | 2019-05-17 | 北京航星机器制造有限公司 | 一种盲卷积运动模糊图像复原方法 |
CN107146202A (zh) * | 2017-03-17 | 2017-09-08 | 中山大学 | 基于l0正则化和模糊核后处理的图像盲去模糊的方法 |
CN107146202B (zh) * | 2017-03-17 | 2020-05-19 | 中山大学 | 基于l0正则化和模糊核后处理的图像盲去模糊的方法 |
CN107688198A (zh) * | 2017-07-25 | 2018-02-13 | 中国科学院地质与地球物理研究所 | 地震数据重构方法及装置、计算机可读存储介质、设备 |
CN113989128A (zh) * | 2020-07-31 | 2022-01-28 | 山东财经大学 | 基于改进的l0梯度最小化模型的图像平滑方法及系统 |
CN113989128B (zh) * | 2020-07-31 | 2022-07-26 | 山东财经大学 | 基于改进的l0梯度最小化模型的图像平滑方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN103413277B (zh) | 2017-04-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103413277A (zh) | 基于l0稀疏先验的盲相机抖动去模糊方法 | |
US9262815B2 (en) | Algorithm for minimizing latent sharp image cost function and point spread function cost function with a spatial mask in a regularization term | |
CN108475415B (zh) | 用于图像处理的方法和系统 | |
CN104599242B (zh) | 使用多尺度非局部正则的模糊核估计方法 | |
Ren et al. | Single image super-resolution via adaptive high-dimensional non-local total variation and adaptive geometric feature | |
US9589328B2 (en) | Globally dominant point spread function estimation | |
US9245328B2 (en) | Algorithm for minimizing latent sharp image cost function and point spread function with a spatial mask in a fidelity term | |
CN111861925B (zh) | 一种基于注意力机制与门控循环单元的图像去雨方法 | |
CN106920220A (zh) | 基于暗原色和交替方向乘子法优化的湍流图像盲复原方法 | |
CN107133923B (zh) | 一种基于自适应梯度稀疏模型的模糊图像非盲去模糊方法 | |
CN102254309A (zh) | 一种基于近红外图像的运动模糊图像去模糊方法和装置 | |
CN102576454A (zh) | 利用空间图像先验的图像去模糊法 | |
US20150294186A1 (en) | Point spread function classification using structural properties | |
CN103761710A (zh) | 基于边缘自适应的高效图像盲去模糊方法 | |
CN105913392A (zh) | 复杂环境下退化图像综合质量提升方法 | |
CN110796616A (zh) | 基于分数阶微分算子的l0范数约束和自适应加权梯度的湍流退化图像恢复方法 | |
CN112634163A (zh) | 基于改进型循环生成对抗网络去图像运动模糊方法 | |
CN104282003A (zh) | 一种基于梯度筛选的数字模糊图像盲复原方法 | |
CN115345791A (zh) | 一种基于注意力机制残差网络模型的红外图像去模糊算法 | |
WO2022233252A1 (zh) | 图像处理方法、装置、计算机设备和存储介质 | |
CN115526777A (zh) | 一种盲超分网络建立方法、盲超分方法及存储介质 | |
CN107742278A (zh) | 结合l0范数和空间尺度信息的运动模糊图像盲复原方法 | |
CN104966269A (zh) | 一种多帧超分辨率成像的装置及方法 | |
KR101921608B1 (ko) | 깊이 정보 생성 장치 및 방법 | |
CN114331913A (zh) | 基于残差注意力块的运动模糊图像复原方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
EE01 | Entry into force of recordation of patent licensing contract | ||
EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20131127 Assignee: Nanjing causal Artificial Intelligence Research Institute Co., Ltd Assignor: Nanjing Post & Telecommunication Univ. Contract record no.: X2019320000168 Denomination of invention: Blind camera shake deblurring method based on L0 sparse prior Granted publication date: 20170405 License type: Common License Record date: 20191028 |