CN106530251B - 一种盲卷积运动模糊图像复原方法 - Google Patents
一种盲卷积运动模糊图像复原方法 Download PDFInfo
- Publication number
- CN106530251B CN106530251B CN201610973795.3A CN201610973795A CN106530251B CN 106530251 B CN106530251 B CN 106530251B CN 201610973795 A CN201610973795 A CN 201610973795A CN 106530251 B CN106530251 B CN 106530251B
- Authority
- CN
- China
- Prior art keywords
- image
- spread function
- point spread
- iteration
- observed
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 44
- 239000000654 additive Substances 0.000 claims description 3
- 230000000996 additive effect Effects 0.000 claims description 3
- 238000012804 iterative process Methods 0.000 claims description 3
- 241000208340 Araliaceae Species 0.000 claims 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims 1
- 235000003140 Panax quinquefolius Nutrition 0.000 claims 1
- 235000008434 ginseng Nutrition 0.000 claims 1
- 238000011084 recovery Methods 0.000 abstract description 11
- 230000000694 effects Effects 0.000 abstract description 10
- 238000004422 calculation algorithm Methods 0.000 abstract description 6
- 238000001914 filtration Methods 0.000 abstract description 4
- 238000007689 inspection Methods 0.000 description 4
- 238000005096 rolling process Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 241000406668 Loxodonta cyclotis Species 0.000 description 1
- 238000007476 Maximum Likelihood Methods 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000012152 algorithmic method Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 210000004209 hair Anatomy 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000011017 operating method Methods 0.000 description 1
- 238000013341 scale-up 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)
Abstract
本发明一种盲卷积运动模糊图像复原方法。该方法包括如下步骤:S1:采集得到观测图像,使用微分滤波器对观测图像进行导数滤波,生成高频退化图像;S2:更新清晰图像;S3:采用迭代最小二乘法求解下式更新求解点扩散函数;S4:若点扩散函数的尺寸小于设定值,增大点扩散函数的尺寸,返回S2。否则,中止迭代。S5:根据S4得到的最佳点扩散函数和观测图像通过现有的非盲卷积方法,求解得到复原图像。本发明求解的点扩散函数稳定,从而抑制了噪声对恢复图像的影响,使得算法的模糊恢复效果显著、自动化程度高操作简单、执行速度快。
Description
技术领域
本发明涉及图像中的信号处理技术,特别是一种盲卷积运动模糊图像复原方法。
背景技术
在用摄像机获取景物图象时,如果在相机曝光期间景物和摄像机之间有相对运动,那么往往会使得到的照片变得模糊。这类图象模糊称运动模糊,通过建立运动图象的复原模型来解决图象的复原问题是研究解决运动模糊的主要手段。如果考虑噪声的影响,运动模糊图象的退化模型可以描述为一个退化函数和一个加性噪声项n(x,y),处理一幅输入图象f(x,y)产生一幅退化图象g(x,y)。
其中,h(x,y)称为模糊核或点扩散函数,表示卷积,f(x,y)表示原始清晰图象,g(x,y)表示观察到的退化图象。
由于空间域的卷积等同于频率域的乘积,所以式(1)的频率域描述为:
G(u,v)=H(u,v)F(u,v)+N(u,v) (2)
式(2)中的大写字母项是式(1)中相应项的傅里叶变换。
运动模糊图象复原技术可以划分为两组:确定性的和随机性的。确定性的方法对于带有很小噪声且退化函数已知的图象有效。原始图象从退化图象的逆变换得到。随机性的技术根据特定的随机准则,即最小二乘方法,找到最优的复原。逆滤波、有约束的最小二乘法、维纳滤波和Richardson-Lucy迭代去卷积等恢复方法是图像复原中最基础最经典的方法。在这里我们回顾一下有约束的最小二乘法与Richardson-Lucy迭代去卷积。
应用有约束的最小二乘方恢复方法时,只需有关噪声均值和方差的知识就可对每幅给定的图象给出最佳恢复结果。有约束最小二乘方法在无噪声或者噪声很小的情况下恢复效果比较理想,对于含有一定强度噪声的情况下,恢复效果也不令人满意。
Richardson-Lucy算法是目前应用很广泛的一种图象迭代恢复方法。Richardson-Lucy算法能够按照泊松噪声统计标准求出与给定点扩散函数卷积后,最有可能成为输入模糊图象的图象。当点扩散函数已知但图象噪声信息未知时,也可以使用这种恢复方法进行有效的恢复。
但是Richardson-Lucy算法方法存在一些较严重的缺陷,问题一是噪声放大问题,这也是这类方法(利用最大似然性求解)的通病。问题二是对于恢复图象中的不同部分分别执行多少迭代才合适的问题。因为图象中信噪比高的部分可能需要数百次迭代才能获得满意的结果;而另一些光滑的对象可能只需很少次数即可达到满意的结果。Richardson-Lucy算法随着迭代次数的增加计算时间也大幅度增加,不利于运动模糊的实时恢复。
发明内容
本发明解决的技术问题是:克服现有技术的不足,提供了一种盲卷积运动模糊图像复原方法,通过引入最小化方案,将非凸问题转化为求解一系列L1范数问题,适用于由不同模糊类型引起的运动模糊图像复原。
本发明的技术方案是:一种盲卷积运动模糊图像复原方法,步骤如下:
(1)采集得到观测图像g(x,y),使用微分滤波器和对观测图像g(x,y)进行导数滤波,生成高频退化图像
(2)更新清晰图像q:
设置外迭代i=0,1..M-1;内迭代j=0,1...N-1;M、N均为正整数;λ为设定的正则化参数;t为设定的收缩阈值;qi为第i次外迭代后的图像;利用下式
经过两层迭代后,调整正则化参数λ'=λ||qi||2,得到更新后的清晰图像q;其中为外迭代i次、内迭代j次后的清晰图像;
(3)采用迭代最小二乘法求解下式更新求解点扩散函数k:
在迭代过程中,通过使用共轭梯度迭代,根据前一次的点扩散函数k更新计算新的约束权重ψ;根据最小二乘迭代法求出最终满足约束条件的点扩散函数k;
(4)若点扩散函数k的尺寸小于设定值kmax,增大点扩散函数k的尺寸,返回步骤(2)交替更新清晰图像q和点扩散函数k;否则,中止迭代并进入步骤(5);
(5)根据步骤(4)得到的最佳点扩散函数k(x,y)和观测图像g(x,y)通过现有的非盲卷积方法,求解得到复原图像f(x,y);其在频域上的求解公式如下:
F(u,v)=(G(u,v)-N(u,v))/(K(u,v)+e)
其中,N(u,v)表示加性噪声,e为常数。
本发明与现有技术相比的优点在于:
1)本发明通过引入尺度正则化弥补了模糊图像中高频部分的衰减,大大提高了点扩散函数估计的稳定性,从而提高了不同模糊类型形成的运动模糊图像恢复的自适应能力,从而获得更好的图像恢复效果。
2)本发明通过引入最小二乘法方案,将非凸问题转化为求解一系列L1范数问题,抑制了噪声对恢复图像的影响,使得算法的模糊恢复效果显著、自动化程度高操作简单、执行速度快。
附图说明
图1为本发明实施例系统的结构示意图;
图2为实施例中采集的观测图像;
图3为采用有约束的最小二乘法对图2观测图像处理得到的清晰图像;
图4为采用Richardson—Lucy方法对图2观测图像处理得到的清晰图像;
图5为采用本发明方法对图2观测图像处理得到的清晰图像。
具体实施方式
采用本发明方法对小型X射线扫描检查系统采集的图像进行复原。该系统结构示意如图1所示,它包括以下几个部分:
X光机:主要由控制器和X射线发生器组成,用于产生连续低能(几十到几百kV)X射线。
机械传送装置:包括传送带、电机和系统骨架、箱体。
探测与数据采集分系统:用探测器阵列和通讯模块对X射线产生的高低能量进行同步测量,并与运行检查分系统进行高速数据传递。
运行检查分系统:包括主机、高分辨率显示器。执行接受用户指令,控制系统运行和接受探测器数据,进行多功能图像处理功能。
电气控制分系统:执行安全联锁、与控制台交互信息、通道货物指示和传送带电机控制的操作。
在小型x射线扫描检查系统中,探测与数据采集分系统采集到被检物的图像后传输到运行检查分系统,在这里进行图像的复原操作。
其中,运动模糊图像复原操作步骤如下:
(1)输入系统采集到的观测图像g(x,y),使用离散滤波器和对观测图像g(x,y)进行导数滤波,生成一幅高频退化图像
(2)由高频退化图像p求出更新评估图像
设外迭代i=0,1..M-1;内迭代j=0,1...N-1;初始化M、N为正整数2;初始化正则化参数λ=20;初始化收缩阈值t=0.001;qi为第i次外迭代后的图像;利用下式
经过两层迭代后,调整正则化参数λ'=λ||qi||2,得到更新后的清晰图像;其中为外迭代i次、内迭代j次后的图像;k为点扩散函数的估计。
(3)采用迭代最小二乘法求解下式更新求解点扩散函数k:
在迭代过程中,通过使用共轭梯度迭代,根据前一次的点扩散函数k更新计算新的约束权重ψ;初始化化(h为点扩散函数k的尺寸,初始化设置模糊核k尺寸为3×3),根据最小二乘迭代法求出满足约束条件的点扩散函数k。详细求解方法参见A.Levin等人发表的《Image and depth from a conventional camera with a coded aperture》,SIGGRAPH,26(3):70,2007
(4)若点扩散函数k的尺寸小于设定值kmax=27×27,设置增长倍数设为按照增长倍数扩大点扩散函数k的尺寸,返回步骤(2)交替更新图像q和点扩散函数k。否则,中止迭代。
(5)根据步骤(4)得到的最佳点扩散函数k(x,y)和观测图像g(x,y)通过现有的非盲卷积方法如Richardson—Lucy算法等,便可求解得到复原图像f(x,y)。本发明选择一种具有快速、强鲁棒性、模糊核误差小的非盲解卷积方法,其在频域上的求解公式如下:
F(u,v)=(G(u,v)-N(u,v))/(K(u,v)+e) (5)
如图2所示,该图为实施例中采集的观测图像;使用有约束的最小二乘法对观测图像进行去模糊处理的效果图如图3所示,由于本组实验的观测图像噪声很小,我们可以看出使用有约束的最小二乘法得到的复原图像效果还算比较理想;图4为使用Richardson—Lucy方法对观测图像进行去模糊处理得到的复原图像,从图中可以看出由于对于未对观测图象中的不同部分分别执行不同迭代次数,使得处理后的复原图像出现了震铃效应和边缘锐化;如图5所示,使用本发明方法进行去模糊操作得到的复原图像是三种方法中复原效果最好的,同时它避免了采用Richardson—Lucy方法时出现的严重的震铃效应。
Claims (1)
1.一种盲卷积运动模糊图像复原方法,其特征在于步骤如下:
(1)采集得到观测图像g(x,y),使用微分滤波器和对观测图像g(x,y)进行导数滤波,生成高频退化图像
(2)更新清晰图像q:
设置外迭代i=0,1..M-1;内迭代j=0,1...N-1;M、N均为正整数;λ为设定的正则化参数;t为设定的收缩阈值;qi为第i次外迭代后的图像;利用下式
经过两层迭代后,调整正则化参数λ'=λ||qi||2,得到更新后的清晰图像q;其中为外迭代i次、内迭代j次后的清晰图像;
(3)采用迭代最小二乘法求解下式更新求解点扩散函数k:
在迭代过程中,通过使用共轭梯度迭代,根据前一次的点扩散函数k更新计算新的约束权重ψ;根据最小二乘迭代法求出最终满足约束条件的点扩散函数k;
(4)若点扩散函数k的尺寸小于设定值kmax,增大点扩散函数k的尺寸,返回步骤(2)交替更新清晰图像q和点扩散函数k;否则,中止迭代并进入步骤(5);
(5)根据步骤(4)得到的最佳点扩散函数k(x,y)和观测图像g(x,y)通过现有的非盲卷积方法,求解得到复原图像f(x,y);其在频域上的求解公式如下:
F(u,v)=(G(u,v)-N(u,v))/(K(u,v)+e)
其中,N(u,v)表示加性噪声,e为常数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610973795.3A CN106530251B (zh) | 2016-11-03 | 2016-11-03 | 一种盲卷积运动模糊图像复原方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610973795.3A CN106530251B (zh) | 2016-11-03 | 2016-11-03 | 一种盲卷积运动模糊图像复原方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106530251A CN106530251A (zh) | 2017-03-22 |
CN106530251B true CN106530251B (zh) | 2019-05-17 |
Family
ID=58349754
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610973795.3A Active CN106530251B (zh) | 2016-11-03 | 2016-11-03 | 一种盲卷积运动模糊图像复原方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106530251B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107730466A (zh) * | 2017-10-11 | 2018-02-23 | 北京工业大学 | 基于结构自相似性和低秩矩阵表示的图像盲复原方法 |
CN108537746B (zh) * | 2018-03-21 | 2021-09-21 | 华南理工大学 | 一种基于深度卷积网络的模糊可变图像盲复原方法 |
CN108694715A (zh) * | 2018-05-15 | 2018-10-23 | 清华大学 | 基于卷积稀疏编码的单相机rgb-nir成像系统 |
CN110533617B (zh) * | 2019-08-30 | 2022-05-27 | Oppo广东移动通信有限公司 | 图像处理方法及装置、存储介质 |
CN110852973B (zh) * | 2019-11-12 | 2022-09-23 | 华中科技大学 | 一种脉冲噪声模糊图像非线性复原方法及系统 |
CN111476722B (zh) * | 2020-03-12 | 2023-04-28 | 哈尔滨工业大学(深圳)(哈尔滨工业大学深圳科技创新研究院) | 基于点扩散函数的图像复原方法、装置及其相关设备 |
CN112116541B (zh) * | 2020-09-24 | 2024-05-14 | 南京航空航天大学 | 基于梯度l0范数和总变分正则化约束的模糊图像复原方法 |
CN112116544B (zh) * | 2020-10-08 | 2022-04-22 | 东南数字经济发展研究院 | 一种抵抗图像平滑滤波的预处理方法 |
CN113256565A (zh) * | 2021-04-29 | 2021-08-13 | 中冶华天工程技术有限公司 | 运动模糊图像智能复原方法 |
CN113628200A (zh) * | 2021-08-19 | 2021-11-09 | 国网上海市电力公司 | 一种基于盲去卷积的输电线路图片去模糊方法 |
CN117115025B (zh) * | 2023-08-25 | 2024-04-05 | 中国人民解放军国防科技大学 | 一种基于多帧盲解卷积的图像快速复原方法、设备及介质 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103413277A (zh) * | 2013-08-19 | 2013-11-27 | 南京邮电大学 | 基于l0稀疏先验的盲相机抖动去模糊方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013148142A1 (en) * | 2012-03-29 | 2013-10-03 | Nikon Corporation | Algorithm for minimizing latent sharp image cost function and point spread function cost function with a spatial mask in a regularization term |
US9692939B2 (en) * | 2013-05-29 | 2017-06-27 | Yeda Research And Development Co. Ltd. | Device, system, and method of blind deblurring and blind super-resolution utilizing internal patch recurrence |
-
2016
- 2016-11-03 CN CN201610973795.3A patent/CN106530251B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103413277A (zh) * | 2013-08-19 | 2013-11-27 | 南京邮电大学 | 基于l0稀疏先验的盲相机抖动去模糊方法 |
Non-Patent Citations (2)
Title |
---|
凹凸范数比值正则化的快速图像盲去模糊;余义斌 等;《电子学报》;20160531;第44卷(第5期);第3节 |
基于高频图像L1范数的模糊图像的复原算法;张寒冰;《工业控制计算机》;20160531;第29卷(第5期);第1.1节 |
Also Published As
Publication number | Publication date |
---|---|
CN106530251A (zh) | 2017-03-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106530251B (zh) | 一种盲卷积运动模糊图像复原方法 | |
CN106846268B (zh) | 一种高斯-脉冲混合图像噪声去除方法 | |
CN106920220B (zh) | 基于暗原色和交替方向乘子法优化的湍流图像盲复原方法 | |
US10489892B2 (en) | Method for signal processing | |
CN104282003B (zh) | 一种基于梯度筛选的数字模糊图像盲复原方法 | |
Xia et al. | A discrete-time learning algorithm for image restoration using a novel L2-norm noise constrained estimation | |
CN104091314A (zh) | 基于边缘预测和稀疏比值正则约束的湍流退化图像盲复原方法 | |
CN1904941A (zh) | 一种图像处理中的去模糊方法 | |
CN112991199B (zh) | 一种基于残差密集网络的图像高低频分解噪声去除方法 | |
Trambadia et al. | Design and analysis of an image restoration using wiener filter with a quality based hybrid algorithms | |
Garvanova et al. | Quality Assessment оf Defocused Image Recovery Algorithms | |
Liu et al. | Total variation with overlapping group sparsity and Lp quasinorm for infrared image deblurring under salt-and-pepper noise | |
CN108376393A (zh) | 一种面向高速直线运动对象的模糊图像盲复原方法 | |
CN103784164B (zh) | 超声信号的预处理方法及系统 | |
CN104933713B (zh) | 一种利用边缘分析的图像mtf估计方法 | |
CN106504213A (zh) | 一种非盲卷积运动模糊图像复原方法 | |
CN106488312B (zh) | 视频源的检测方法及装置 | |
CN110121016B (zh) | 一种基于双重曝光先验的视频去模糊方法及装置 | |
Zhu et al. | Image Restoration Based on Wiener Filter and Constrained Least Square Filter | |
Yang et al. | Guided filter based edge-preserving image non-blind deconvolution | |
Su et al. | Defocused image restoration using RBF network and iterative Wiener filter in wavelet domain | |
Hall et al. | Foreground and background reconstruction in Poisson video | |
CN111275692B (zh) | 基于生成对抗网络的红外小目标检测方法 | |
CN106023106A (zh) | 一种自适应参数更新的贝叶斯图像反卷积恢复方法 | |
Patil et al. | Quality Improvements of Camera Captured Pictures using Blind and Non-blind Deconvolution Algorithms |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 |