??通过点扩散函数 PSF 进行图像复原的实验添加适当的噪声。
??fspecial函数中’motion’指的是运动模糊算子。
% checkerboard产生测试板图像第一个参数是每个正方形一边的像素数,第二个參数行数第三为列数(缺省则等于行数)
??周期噪声本身表现为类似冲击的串,这种串在傅里叶频谱中通常是可见的滤除这些成分嘚主要途径是通过陷波滤波。n阶巴特沃兹陷波滤波器的传递函数如下:
正则滤波(约束的最小二乘方滤波)
使用L-R算法的迭代非线性复原
??前媔所讨论的图像复原方法都是线性的而LR算法是非线性方法,又称Lucy-Richardson算法
??不以PSF知识为基础的图像复原方法统称为盲去卷积算法。
摘要: 摘 要: 提出了一种改进的運动模糊复原算法首先利用改进的高通滤波(方向微分)快速算法判定运动模糊方向,然后利用边缘检测算法计算模糊尺度,从而确定点扩散函数最后通过维纳滤波算法对模糊图像进行恢复。针对K值估计问题提出了一种新的K值自动估计算法。实验结果表明该方法鲁棒性强、实时性好,图像复原效果好 关键词: 图像复原; 点扩散函数; 维纳滤波; K值估计 &nb
摘 要: 提出了一种改进的运动模糊复原算法,首先利鼡改进的高通滤波(方向微分)快速算法判定运动模糊方向然后利用边缘检测算法计算模糊尺度,从而确定点扩散函数,最后通过维纳滤波算法对模糊图像进行恢复针对K值估计问题,提出了一种新的K值自动估计算法实验结果表明,该方法鲁棒性强、实时性好图像复原效果恏。
关键词: 图像复原; 点扩散函数; 维纳滤波; K值估计
运动模糊图像复原是图像恢复的重要课题之一随着科学技术的不断进步和发展,运动模糊图像复原在天文、军事、道路交通[1]、医学图像及侦破等各个领域中的应用越来越多要求也越来越高,从事运动模糊图像的恢複研究具有重要的现实意义
目前,对于不同原因造成的模糊图像的复原虽然有广泛的研究不少学者也已经提出了一些不同的具体算法,但这些方法一般仅适用于已知点扩散函数的参数情况而对于运动模糊图像,点扩散函数通常是未知的如果知道模糊的方向和尺度,僦可以据此构造点扩散函数用前人提出的方法进行恢复。可见寻求一种行之有效的方法来自动鉴别运动模糊参数并实现运动模糊图像嘚复原是一个迫切并且必要的任务。
对运动模糊图像的盲复原关键是精确确定点扩散函数(即模糊长度和模糊方向)并寻求一种好的複原方法。本文对模糊图像进行方向性的高通滤波(方向微分)来判定运动模糊方向并针对原算法运算速度慢的缺点进行改进,采用由粗到细的检测方法从而大大提高了检测速度;然后利用边缘检测算法,计算模糊的尺度;最后采用维纳滤波算法对模糊图像进行恢复並提出一种新的K值估计算法,将K值自动估计与传统的K取经验值相结合进行两次维纳滤波,从而得到一个较为精确的K值
本算法可以较恏地实现图像复原,特别是针对一些模糊程度较强、人的肉眼根本不能辨别的模糊图像也能较好地恢复具有很强的实用价值。
1 图像复原嘚原理 图像复原处理建立在图像退化的数学模型基础上这个退化数学模型应该能够反映图像退化的原因。由于图像的退化因素较多洇此图像处理过程中把退化原因作为线性系统退化的一个因素来对待,从而建立系统退化模型来近似描述图像函数的退化模型如图1所示。若H是线性的、空间不变的过程则退化图像在空间域的表达式为:
可见,用维纳滤波实现运动模糊图像的复原实质上就是求取退化函数H(u,v)囷K值(信噪比)的过程其取值的准确性将直接影响复原效果。下面介绍H(u,v)和K值的求取
由前面分析可知,要实现图像复原首先要建立正确的退囮模型图像复原处理可看成是一个估计过程,如果已经给出了退化图像g(x,y)并估计出点扩散函数PSF的参数(运动模糊方向和运动模糊长度)则任意方向的匀速直线运动[4]模糊图像的点扩散函数也就可以根据任意方向匀速直线运动模糊的退化模型而唯一确定,进而由最小二乘准则使用參变维纳滤波来近似复原出f(x,y)因此,运动模糊图像的复原首先要确定运动模糊图像的退化模型,再根据退化模型和原始图像估计出PSF的参數,最后由相应的复原方法进行完成图像的复原工作
2.2 模糊方向鉴别 运动方向的提取是运动模糊图像确定点扩散函数的关键因素之一。准确地估计出运动模糊方向就可以通过图像旋转,将运动模糊方向旋转到水平轴方向对应的运动模糊点扩散函数也随之变为一维的,运动模糊点扩散函数的估计及图像复原就由二维问题转化为一维问题从而降低运动模糊图像复原的难度。
将原始图像看成是自相关及其功率谱是各向同性的一阶马尔科夫过程,運动模糊降低了运动方向上图像的高频成分而对于其他方向上高频成分影响较少,且随着方向偏离越大影响越小若对运动模糊图像进荇方向性的高通滤波(方向微分),则当滤波方向为运动方向时由于此方向对应的高频成分最小,因此高通滤波后模糊图像能量损失最大,得箌的微分图像对应的灰度值之和必然最小因此方向微分后的图像的灰度值之和最小的方向即为运动模糊方向。方向微分示意如图2所示
妀进的模糊方向求解方法采用由粗到细的方法逐步减小感兴趣的角度区域,而将大量的非目标角度以较大的步长跳过从而简化了运算,夶大提高了检测速度
2.3 基于边缘检测的运动模糊长度的检测 模糊长度是点扩散函数的另一个重要参数,本文将边缘检测应用于运动模糊长喥的鉴别
为了使算法具有更强的鲁棒性在此用梯度变化来近似表礻图像的边缘,进而求出模糊长度将Prewitt边缘检测子扩展为n×n(n为奇数)的矩阵, 在此用5×5的矩阵即n=5。
运动模糊图像的边缘延伸到周围像素甴于运动模糊部分总是发生在模糊区域的两端,且模糊具有对称性,因此对模糊图像进行边缘检测,分别求取两端的边缘宽度,取其平均值为运動模糊长度在此用迭代算法求取两端的边缘宽度,具体步骤如下:
(2)对边缘图像进行逐行扫描记录每行中连续边缘的边缘宽度,并对整個边缘图像中的边缘宽度进行统计,做出统计直方图设定边缘初始宽度阈值T等于边缘宽度统计平均值以及迭代次数为n;
(3)将每一行中连續边缘的边缘宽度Tij与上述阈值进行比较,将符合0.8T≤Tij≤1.2T的边缘宽度保留下来;
(4)对每一行中最大的两个边缘宽度求和,并对所有行的和值进行統计求出统计平均值;
(5)用步骤(4)中得到的值的1/2更新边缘宽度阈值;
(6)返回步骤(2)重复以上步骤,直到宽度值收敛或达到迭代次数
3 改进的K徝自动估计 用维纳滤波进行图像复原的关键是确定点扩散函数和K值。前文已经求出模糊角度和模糊方向将其带入公式(5)即可得到点扩散函數h(x,y)。因此接下来的主要的任务是求K值。
4 实验研究 以运动模糊車牌[6]的复原为例在摄像头拍摄时,因为车辆高速运动加之曝光时间通常较短,得到的图像往往非常模糊这就要首先做图像复原处理,最后得到一幅大致上比较清晰的图像以供后续处理
从结果可以看出利用本文提出的改进的K值自动估计算法与传统的算法相比,复原效果有明显的改善。该算法对于模糊程度未知、造成模糊的原因未知的模糊图像的复原问题具有重要的实际应用价值
本文對运动模糊图像复原进行了研究,对运动模糊图像的退化模型、点扩散函数的求取以及维纳滤波K值的自动估计等几个问题给出了相应的改進算法以运动模糊车牌复原为例进行了实验验证,并与传统的K取经验值进行了比较实验结果表明,本算法能够较好地估计运动模糊参數通过对方向微分算法的改进提高了运算速度,利用改进的K值自动估计对图像复原效果与传统算法相比较复原效果具有明显的改善
图像的老化可以视为以下这样嘚一个过程。一个是退化函数的影响(致使图片模糊褪色等),一个可加性噪声的影响
简单的来考虑,加法的逆运算是减法乘法的逆运算的除法,微分的逆运算是积分(严密一点说是不定积分)那么就可以得到一个简单的结论了,要出去卷积的话肯定需要用到卷积的逆运算。卷积的逆运算是---------反卷积额,好像是一个理所应当的名字 我们建立了一个关于卷积的直观认识,将信号反转与滤波器系数求积囷那么,反卷积是一种什么样的运算呢或者具体的来讲,反卷积的空间运算表现形式是什么样的这样的考虑其实是多余的,或者说不用考虑的那么复杂。
综上所述,其实逆滤波的问题点有两个:、
从表示式上可以看出,这个模型是不会有0值的不过,这个模型与低通滤波器佷相似阻带的值都是极小的。这可能会使得图像的直接逆滤波失败这个之后再说。
这个模型其实在Photoshop上也有一个同名的滤镜详细的推倒我就不做了,这个模型的表达式如下所示
注意,运动模糊后的图像的尺寸会变化如果还是按照原图截取,会造成图像成分的损失茬复原图像时候效果不是太好,而且不知道导致效果不好的原因是由于成分的缺失,还是噪声的干扰所以,这里我适当的扩展了图像嘚尺寸以保留图像的所有成分
下面是实验用图像。图像的噪声选用的是高斯噪声均值为0,方差为0.08退化函数则选用先前叙述的两种,┅个个大气湍流模型一个是运动模糊。
所谓直接逆滤波就是不管噪声的影响,直接进行逆滤波的方法
对于大气湍流模型而言,直接逆滤波会得到很不理想的结果下面是直接逆滤波的实验结果。
实验结果完全没有任何价值观察其频谱,频谱的四角很亮而原本频谱朂亮的直流分量都看不到了。所以这里做一个限制处理。也就是仅仅只处理靠近直流分量的部分,其他的不做处理然后处理完的结果,过一个10阶巴特沃斯低通滤波器可以得到如下结果。
这样的话只需要调整限制半径,可以得到一个比之前较好的结果当然,这招茬运动模糊的图片面前就略显得无力了。结果我就不贴了
维纳滤波器的推导,其实是个很复杂的过程这里就不推导了,直接看结果可以得到一些有用的结论。
下面的维纳滤波器的实验结果:
其实这个方法是一个很好的想法,将图像的能量作为评价图像平滑程度的度量尽可能的将其平滑。设噪声的能量是一个定值使用拉普拉斯未定系数法,将其进行迭代然后解开。 这种方法有了很多的变种包括很著名的TV(Total Variation,全变分)模型这个我之后的博文会讲到。
在这里使用本方法的目的是,减尐噪音对于逆滤波的影响表达式如下所示。
这个滤波器可以消除很严重的噪声,并且复原图像将实验用图像的噪声的方差提升到0.2,洅进行滤波可以得到如下结果。