图像复原及重建

1、东北大学研究生考试试卷评分考试科目：图像处理分析与识别课程编号：y2014600802阅 卷 人: 魏 颖、陈东岳考试日期： 2014年11月29日姓 名： 傅仁德学 号： 1401029 注 意 事 项1. 考前研究生将上述项目填写清楚.2. 字迹要清楚,保持卷面清洁.3. 交卷时请将本试卷和题签一起上交.4. 课程考试后二周内授课教师完成评卷工作,公共课成绩单与试卷交研究生院培养办公室,专业课成绩单与试卷交各学院,各学院把成绩单交研究生院培养办公室.东北大学研究生院培养办公室运动模糊图像复原的算法分析与实现摘要本文通过建立图像降质的数学模型，详细分析了图像在匀速直线运动时及加性高斯白噪声作

2、用下的退化过程，提出了运动模糊图像复原的三种算法，即逆滤波、维纳滤波和约束最小二乘方滤波。并分别用MATLAB编程实现其算法的可行性。实验结果表明，逆滤波产生了一幅不大适合的图像，其噪声非常的明显，视觉效果不仅没有改善，反而使图像更加模糊。维纳滤波的结果也不好，不过就恢复图像的细节（即高频分量）而言，它比逆滤波效果显然要好。约束最小二乘方滤波很好地恢复了图像中的文字，不过图像背景有不少失真。就总体效果而言，该算法比前两者要好得多。关键字：退化 运动 模糊 复原 滤波目 录1 图像退化与复原11.1 图像降质的数学模型11.2 匀速直线运动模糊的退化模型22 运动模糊图象的复原方法及原理42.1

3、 逆滤波复原原理42.2 维纳滤波复原原理52.3 有约束最小二乘复原原理63 运动模糊图像复原的实现73.1 运动模糊图像复原的MATLAB实现73.2 复原结果比较8实验小结9参考文献101 图像退化与复原数字图像在获取的过程中，由于光学系统的像差、光学成像衍射、成像系统的非线性畸变、摄影胶片的感光的非线性、成像过程的相对运动、大气的湍流效应、环境随机噪声等原因， 图像会产生一定程度的退化。因此，必须采取一定的方法尽可能地减少或消除图像质量的下降，恢复图像的本来面目， 这就是图像复原，也称为图像恢复。图像复原与图像增强有类似的地方，都是为了改善图像。但是它们又有着明显的不同。图像复原是试图

4、利用退化过程的先验知识使已退化的图像恢复本来面目，即根据退化的原因， 分析引起退化的环境因素，建立相应的数学模型，并沿着使图像降质的逆过程恢复图像。从图像质量评价的角度来看，图像复原就是提高图像的可理解性。而图像增强的目的是提高视感质量，图像增强的过程基本上是一个探索的过程，它利用人的心理状态和视觉系统去控制图像质量，直到人们的视觉系统满意为止。图像复原是利用退化现象的某种先验知识，建立退化现象的数学模型，再根据模型进行反向的推演运算，以恢复原来的景物图像。因而，图像复原可以理解为图像降质过程的反向过程。建立图像复原的反向过程的数学模型，就是图像复原的主要任务。经过反向过程的数学模型的运算，要

5、想恢复全真的景物图像比较困难。所以，图像复原本身往往需要有一个质量标准，即衡量接近全真景物图像的程度，或者说，对原图像的估计是否到达最佳的程度。由于引起退化的因素众多而且性质不同，为了描述图像退化过程所建立的数学模型往往多种多样，而恢复的质量标准也往往存在差异性，因此图像复原是一个复杂的数学过程，图像复原的方法、技术也各不相同。1.1 图像降质的数学模型图像复原处理的关键问题在于建立退化模型。输入图像经过某个退化系统后输出的是一幅退化的图像。为了讨论方便，把噪声引起的退化即噪声对图像的影响一般作为加性噪声考虑，这也与许多实际应用情况一致，如图像数字化时的量化噪声、随机噪声等就可以作为加性噪声，

6、即使不是加性噪声而是乘性噪声，也可以用对数方式将其转化为相加形式。原始图像经过一个退化算子或退化系统的作用，再和噪声进行叠加，形成退化后的图像。图1-1表示退化过程的输入和输出的关系，其中概括了退化系统的物理过程，就是所要寻找的退化数学模型。图1-1 图像的退化模型数字图像的图像恢复问题可看作是：根据退化图像和退化算子的形式沿着反向过程去求解原始图像，或者说是逆向地寻找原始图像的最佳近似估计。 图像退化的过程可以用数学表达式写成如下的形式： (1-1) 在这里，是一种统计性质的信息。在实际应用中，往往假设噪声是白噪声，即它的频谱密度为常数，并且与图像不相关。在图像复原处理中，尽管非线性、时变和

7、空间变化的系统模型更具有普遍性和准确性，更与复杂的退化环境相接近，但它给实际处理工作带来了巨大的困难，常常找不到解或者很难用计算机来处理。因此，在图像复原处理中， 往往用线性系统和空间不变系统模型来加以近似。这种近似的优点使得线性系统中的许多理论可直接用于解决图像复原问题，同时又不失可用性。1.2 匀速直线运动模糊的退化模型在所有的运动模糊中，由匀速直线运动造成图象模糊的复原问题更具有一般性和普遍意义。因为变速的、非直线运动在某些条件下可以被分解为分段匀速直线运动。本节只讨论由水平匀速直线运动而产生的运动模糊。假设图象有一个平面运动，令和分别为在和方向上运动的变化分量，表示运动的时间。记录介质

8、的总曝光量是在快门打开后到关闭这段时间的积分。则模糊后的图象为： (1-2)式中为模糊后的图象。以上就是由于目标与摄像机相对运动造成的图象模糊的连续函数模型。如果模糊图象是由景物在方向上作匀速直线运动造成的，则模糊后图象任意点的值为： (1-3)式中是景物在方向上的运动分量，若图象总的位移量为，总的时间为，则运动的速率为。则上式变为： (1-4)以上讨论的是连续图象，对于离散图象来说，对上式进行离散化得： （1-5）其中为照片上景物移动的像素个数的整数近似值。是每个像素对模糊产生影响的时间因子。由此可知，运动模糊图象的像素值是原图象相应像素值与其时间的乘积的累加。从物理现象上看，运动模糊图象实

9、际上就是同一景物图象经过一系列的距离延迟后再叠加，最终形成的图象。如果要由一幅清晰图象模拟出水平匀速运动模糊图象，可按下式进行： （1-6）这样可以理解此运动模糊与时间无关，而只与运动模糊的距离有关，在这种条件下，使实验得到简化。因为对一幅实际的运动模糊图象，由于摄像机不同，很难知道其曝光时间和景物运动速度。我们也可用卷积的方法模拟出水平方向匀速运动模糊。其过程可表示为： (1-7)其中 (1-8)称为模糊算子或点扩散函数，表示卷积，表示原始(清晰)图象，表示观察到的退化图象。如果考虑噪声的影响,运动模糊图象的退化模型可以描述为一个退化函数和一个加性噪声项，处理一幅输入图象产生一幅退化图象。

10、(1-9)由于空间域的卷积等同于频率域的乘积，所以式(2-9)的频率域描述为： (1-10)式(1-9)中的大写字母项是式(1-10)中相应项的傅里叶变换。2 运动模糊图象的复原方法及原理运动模糊图象的恢复是图象复原的主要分支之一，它的恢复算法有很多种。有些算法虽然有很好的恢复效果，但算法复杂，恢复时间比较长（如最大熵法）。有些算法虽然计算速度较快，但恢复效果不尽人意（如空间域逆向恢复）。下面介绍逆滤波、维纳滤波和有约束最小二乘滤波三种恢复方法的原理。2.1 逆滤波复原原理 在六十年代中期，逆滤波（去卷积）开始被广泛地应用于数字图象复原。Nathan用二维去卷积方法来处理由漫游者、探索者等外星

11、探索发射得到的图象。由于和噪声相比，信号的频谱随着频率升高下降较快，因此高频部分主要是噪声。Nathan采用的是限定逆滤波传递函数最大值的方法。在同一时期，Harris采用PSF的解析模型对望远镜图象总由于大气扰动造成的模糊进行了逆滤波处理，Mcglamery则采用由实验确定的PSF来对大气扰动图象进行逆滤波。从此以后，逆滤波就成了模糊图象复原的一种标准技术。恢复退化图象最简单的方法是直接逆滤波。在该方法中，用退化函数除退化图象的傅里叶变换来计算原始图象的傅里叶变换估计，由式2-1可以得到逆滤波退化公式： (2-1)这个公式说明逆滤波对于没有被噪声污染的图象很有效，这里不考虑在空间的某些位置上

12、当接近0时可能遇到的计算问题，幸运的是忽略这些点在恢复结果中并不会产生可感觉到的影响。但是，如果出现噪声就会引起几个问题：第一，对于幅值比较小的频率处噪声的影响可能变得显著起来。这种状况通常对于高频。在实际中，通常幅值衰减得比快得多，因此噪声的影响可能支配整个复原结果。将复原限定在足够大得原点处得一个小邻域中，可以克服这个问题。第二个问题针对噪声本身的频谱，我们通常没有充分的有关噪声的信息来足够好地确定。为了克服接近0所引起的问题，在分母中加入一个小的常数，将式(2-1)修改为： （2-2）2.2 维纳滤波复原原理在大部分图象中，邻近的像素是高度相关的，而距离较远的像素其相关性较弱。由此，我们

13、可以认为典型图象的自相关函数通常随着与原点的距离增加下降。由于图象的功率谱是图象本身自相关函数的傅里叶变换，我们可以认为的功率谱随着频域的升高而下降。一般地，噪声源往往具有平坦的功率谱，即使不是如此，其随着频率的升而下降的趋势也要比典型图象的功率谱慢得多。因此，可以料想功率谱的低频分以信号为主，然而高频部分则主要被噪声所占据。由于逆滤波滤波器的幅值常随着频率的升高而升高，因此会增强高频部分的噪声。为了克服以上缺点，出了采用最小均方误差的方法（维纳滤波）进行模糊图象恢复。维纳(wiener)滤波可以归于反卷积(或反转滤波)算法一类，它是由Wiener首提出的，并应用于一维信号，并取得很好的效果。

14、以后算法又被引入二维信号理，也取得相当满意的效果，尤其在图象复原领域，由于维纳滤波器的复原效良好，计算量较低，并且抗噪性能优良，因而在图象复原领域得到了广泛的应用并不断得到改进发展，许多高效的复原算法都是以此为基础形成的。如果取 (2-3)和分别是图象和噪声的自相关矩阵。即，并且都是正定对称矩阵，则有 (2-4)的模方最小，实际上就意味着使噪声和信号的比对复原图象影响最小。因为图象和噪声的相关矩阵都是把图象当作随机过程来研究，从而描述其统计特性的量，在这里最小二乘方的最佳已经演变成均方误差最小准则下的最佳。同样根据式(2-4)可求得频域维纳滤波公式如下 (2-5)时，为标准维纳滤波器；时，为含

15、参维纳滤波器。若没有噪声时,即,维纳滤波器则退化成理想反滤波器。实际应用中必须调节以满足式(2-4)。因为,实际很难求得因此，可以用一个比值代替两者之比，从而得到简化的维纳滤波公式: (2-6)2.3 有约束最小二乘复原原理 由于大多数图象恢复问题都不具有唯一解，或者说恢复具有病态特征。为了克服这一问题，通常需要在恢复过程中对运算施加某种约束。设对图象施加某一线性运算，求在约束条件 （2-7）下，使为最小的作为原图的最佳估计。利用拉格朗日乘数法，先构造一辅助函数： （2-8）令可得： （2-9）解之得： （2-10） 式中。把式(2-10)代入式(2-7)中可以证明，是的单调递增函数。因此可以

16、用迭代法求出满足约束条件(2-7)式的待定系数，首先任取一个,代入(2-10)，把求得的再代入式(2-7)，若结果大于时，便减少；反之增大，再重复上述过程，直到约束条件式(2-11)被满足为止（实际求解时，只要能之差小于某一给定值就可以了）。把求得的代入，便最后求得最佳估计。 我们可以直接从空间域的有约束最小二乘方恢复式(2-10)得到它的频域解: (2-11) 应用有约束最小二乘方恢复方法时，只需有关噪声均值和方差的知识就可对每幅给定的图象给出最佳恢复结果。3 运动模糊图像复原的实现 数字图像在获取的过程中，由于光学系统的像差、光学成像衍射、成像系统的非线性畸变、摄影胶片的感光的非线性、成像

17、过程的相对运动、大气的湍流效应、环境随机噪声等原因，图像会产生一定程度的退化.因此，必须采取一定的方法尽可能地减少或消除图像质量的下降，恢复图像的本来面目，这就是图像复原，也称为图像恢复。图像复原是试图利用退化过程的先验知识使已退化的图像恢复本来面目，即根据退化的原因，分析引起退化的环境因素，建立相应的数学模型，并沿着使图像降质的逆过程恢复图像。前面讲述的三种运动模糊图像复原方法均可在matlab环境下实现，下面将对每一种方法进行举例。3.1 运动模糊图像复原的MATLAB实现 I=imread(C:UsersAdministratorDesktopFig0526(a).tif); figur

18、e(1);imshow(I,); title(原图像); PSF=fspecial(motion,30,45); MF=imfilter(I,PSF,circular); noise=imnoise(zeros(size(I),gaussian,0,0.0001); MFN=imadd(MF,im2uint8(noise); figure(2);imshow(MFN,); title(运动模糊图像); NSR=sum(noise(:).2)/sum(MFN(:).2); figure(3); imshow(deconvwnr(MFN,PSF),); title(逆滤波复原); figure(4); imshow(deconvwnr(MFN,PSF,NSR),); title(维纳滤波复原); NP=0.002*prod(size(I); reg1 LAGRA=deconvreg(MFN,PSF,NP/300); figure(5);imshow(reg1); title(最小二乘滤波复原);下面是程序输出结果： 3.2 复原结果比较实验结果表明，逆滤波产生了一幅不