毕业设计说明书：基于运动图像复原的维纳滤波器设计.doc

1、兰州理工大学毕业设计ILANZHOUUNIVERSITYOFTECHNOLOGY毕业设计说明书题目基于运动图像复原的维纳滤波器设计学生姓名学号专业班级通信工程指导教师学院计算机与通信学院答辩日期2013年6月13日兰州理工大学毕业设计II基于运动图像复原的维纳滤波器设计THEDESIGNOFTHEWIENERFILTERFORIMAGERESTORATIONBASEDON兰州理工大学毕业设计III摘要运动模糊图像成像过程中可能会出现模糊、失真或混入噪声，最终导致图像质量下降。这种质量的下降会造成图像中的目标很难识别或者图像中的特征无法提取，必须对其进行恢复,维纳滤波是一种常见的图像复原方法。本

2、设计主要对维纳滤波的基本原理进行研究，并结合MATLAB中的函数，设计相应的维纳滤波器，对运动模糊图像和它的加噪图像进行复原。之后，对逆滤波和维纳滤波进行图像复原仿真实验，并对比它们的复原效果。关键词维纳滤波；图像恢复；退化模型兰州理工大学毕业设计IVABSTRACTMOTIONBLURINIMAGINGPROCESSMAYAPPEARFUZZY,DISTORTEDORMIXEDWITHNOISE,RESULTINGINTHEDECREASEOFIMAGEQUALITYTHEDROPINQUALITYWILLCAUSETHATTHEIMAGEOFTHETARGETCANNOTBEEXTRACT

3、EDANDAREDIFFICULTTOIDENTIFY,SOITMUSTBERESTOREDWIENERFILTERINGISACOMMONMETHODFORIMAGERESTORATIONSTUDYONTHEDESIGNOFTHEMAINPRINCIPLEOFTHEWIENERFILTER,ANDCOMBINEDWITHTHEFUNCTIONOFMATLAB,DESIGNTHECORRESPONDINGWIENERFILTER,TOTHERESTORATIONOFMOTIONBLURREDIMAGEANDNOISEIMAGEITTHEN,ONTHEINVERSEFILTERANDWIENER

4、FILTERFORIMAGERESTORATIONSIMULATIONEXPERIMENT,ANDCOMPAREDTHEIRRESTORATIONEFFECTKEYWORDSWIENERFILTERIMAGERESTORATIONDEGRADEDIMAGE兰州理工大学毕业设计V目录第1章绪论111图像复原的背景及意义112图像复原方法213维纳滤波简介2第2章图像基本退化模型及恢复421图像噪声422图象退化模型5221退化模型5222连续函数退化模型7223离散函数退化模型9224匀速直线运动图像的退化模型1223图像的恢复方法14231逆滤波复原法14232约束最小平方复原法15233维纳

5、滤波复原法16第3章维纳滤波实现退化图像的复原1831维纳滤波的基本原理18311维纳滤波概述18312运动模糊参数的确定19313维纳霍夫（WIENERHOPF）方程2332维纳滤波仿真实现24321K参数对运动模糊图像复原的影响24322图像的恢复效果对比26总结33参考文献33附录一外文文献翻译35附录二源程序清单54致谢62兰州理工大学毕业设计1第1章绪论在实际的日常生活中，人们要接触很多匀速运动图像，画面，而在景物成像这个过程里可能会出现模糊、失真或混入噪声，最终导致图像质量下降，这种现象称为图像“退化”。因此我们可以采取一些技术手段来尽量减少甚至消除图像质量的下降，还原图像的本来面

6、目。这就是图像复原。引起图像模糊有很多种的原因，举例来说有运动引起的，高斯噪声引起的，斑点噪声引起的，椒盐噪声引起的等等。图像复原的算法数字图像复原问题实际上是在一定的准则下，采用数学最优化方法从退化的图像去推测原图像的估计问题。不同的准则及不同的数学最优化方法就形成了各种各样的算法。常见的复原方法有，逆滤波复原算法，维纳滤波复原算法，盲卷积滤波复原算法，约束最小二乘滤波复原算法等等。图像复原是图像处理中的重要技术，图像复原可以在某种意义上对图像进行改进，即可以改善图像的视觉效果，又能够便于后续处理。其中维纳滤波是最典型的一种，20世纪40年代，维纳奠定了最佳滤波器研究的基础。即假定输入时有用

7、信号和噪声信号的合成，并且它们都是广义平稳过程和他们的二阶统计特性都已知。维纳根据最小均方准则（即滤波器的输出信号与需要信号的均方值最小），求得了最佳线性滤波器的参数，这种滤波器被称为维纳滤波器。MATLAB是一款主要用于数值计算和图像处理的工具软件。由于它采用了矩阵的形式存贮数据，因此在图像处理领域能够发挥速度快，效率高的优点。它包含了许多功能强大的工具箱，借助于这些工具箱，用户可以非常方便地进行图像分析和处理工作。此外，和其它软件比较，由于MATLAB对于图像处理的针对性，它还具有代码简洁的优势。正是基于上述情况，本文采用了MATLAB来实现文中提到的算法，并且取得了不错的效果。11图像复

8、原的背景及意义图像复原就是研究如何从所得的变质图像中复原出真实图像，或说是研究如何从获得的信息中反演出有关真实目标的信息。造成图像变质或者说使图像模糊的原因很多，如果是因为在摄像时相机和被摄景物之间有相对运动而造成的图像模糊则称为运动模糊。所得到图像中的景物往往会模糊不清，我们称之为运动模糊图像。运动模糊图像在日常生活中普遍存在，给人们的实际生活带来了很多不便。近年来，在数字图像处理领域，关于运动模糊图像的复原处理成为了国内外研究的热点问题之一，也出现了一些行之有效的算法和方法。这些算法和方法在不同的情况下，具有不同的复原效果。因为这些算法都是其作者在假定的前提条件下提出的，而实际上的模糊图像

9、，并不是兰州理工大学毕业设计2一定能够满足这些算法前提，或者只满足其部分前提。作为一具实用的图像复原系统，就得提供多种复原算法，使用户可以根据情况来选择最适当的算法以得到最好的复原效果。图像复原关键是要知道图像退化的过程，即要知道图像退化后的图像进行复原处理非常具有现实意义。图像复原的目的就是根据图像退化的先验知识，找到一种相应的反过程的方法来处理图像，从而尽量得到原来图像的质量，以满足人类视觉系统的要求，以便观赏、识别或者其它应用的需要。12图像复原方法图像复原技术在实际生活中有着很广泛的应用。图像复原算法有线性和非线性两类。常用的几种图像复原方法，如维纳滤波法、正则滤波法、LR算法、盲去卷

10、积等，它们都有自己的特点，也都能满足一定条件下对退化图像的处理。1）维纳滤波法1维纳滤波法是由WIENER首先提出的，应用于一维信号处理，取得了很好的效果。之后，维纳滤波法被用于二维信号处理，也取得了不错的效果，尤其在图像复原领域由于维纳滤波计算量小，复原效果好，从而得到了广泛的应用和发展。2）正则滤波法另一个容易实现线性复原的方法称为约束的最小二乘方滤波，在IPT中称为正则滤波，并且通过函数DECONVREG来实现。3）LUCYRICHARDSON算法LR算法是一种迭代非线性复原算法，它是从最大似然公式印出来的，图像用泊松分布加以模型化的。4）盲去卷积在图像复原过程中，最困难的问题之一是，如

11、何获得PSF的恰当估计。那些不以PSF为基础的图像复原方法统称为盲去卷积。它以MLE为基础的，即一种用被随机噪声所干扰的量进行估计的最优化策略。13维纳滤波简介维纳滤波器（WIENERFILTER）是由数学家维纳（RORBERTWIENER）提出的一种以最小平方为最优准则的线性滤波器。在一定的约束条件下，其输出与一给定函数（通常称为期望输出）的差的平方达到最小，通过数学运算最终可变为一个托布利兹方程的求解问题。维纳滤波器又被称为最小二乘滤波器或最小平方滤波器，目前是基本的滤波方法之一。维纳滤波是利用兰州理工大学毕业设计3平稳随机过程的相关特性和频谱特性对混有噪声的信号进行滤波的方法。维纳滤波，

12、从连续的（或离散的）输入数据中滤除噪声和干扰以提取有用的信息的过程称为滤波，滤波器研究的一个基本课题就是如何设计和制造最佳的或最优的滤波器。所谓最佳滤波器是指能够根据某一最佳准则进行滤波的滤波器。20世纪40年代，维纳奠定了关于最佳滤波器研究的基础，即假定线性滤波器的输入为有用信号和噪声之和，两者均为广义平稳过程且知它们的二阶统计特性，维纳根据最小均方误差准则（滤波器的输出信号与需要信号之差的均方值最小），求的了最佳线性滤波器的参数，这种滤波器称为维纳滤波器。在维纳研究的基础上，人们还根据最大输出信噪比准则、统计检测准则以及其他最佳准则求得的最佳线性滤波器。实际上，在一定条件下，这些最佳滤波器

13、与维纳滤波器是等价的。设维纳滤波器的输入为含噪声的随机信号。期望输出与实际输出之间的差值为误差，对该误差求均方，即为均方误差。因此均方误差越小，噪声滤除效果就越好。为使均方误差响应。如果能够满足维纳霍夫方程，就可使维纳滤波器达到最佳。根据维纳霍夫方程，最佳维纳滤波器的冲激响应，完全由输入自相关函数以及输入与期望输出的互相关函数所决定。基本维纳滤波就是用来解决从噪声中提取信号问题的一种过滤方法。它基于平稳随机过程模型，且假设退化模型为线性空间不变系统的。实际上这种线性滤波问题，可看成是一种估计问题或一种线性估计问题。基本的维纳滤波是根据全部过去和当前的观察数据来估计信号的当前值，它的解是以均方误

14、差最小条件下所得到的系统的传递函数或单位样本响应的形式给出的，因此更常称这种系统为最佳线性滤波器。设计维纳滤波器的过程就是在寻求在最小均方误差下滤波器的单位样本响应或传递函数的表达式，其实质是在解维纳霍夫方程。兰州理工大学毕业设计4第2章图像基本退化模型及恢复21图像噪声噪声可以理解为妨碍人们感觉器官对所接收的信源信息理解的因素。而图像中各种妨碍人们对其信息接受的因素即可称为图像噪声，噪声在理论上可以定义为不可预测，只能用概率统计方法来认识的随机误差，因此将图像噪声看成是多维随机过程是合适的，因而描述噪声的方法完全可以借用随机过程的描述，即用其概率分布函数和概率密度分布函数15。设图像信号对黑

15、白图像可看作是二维亮度分布了,FXY，则噪声可看作是对亮度的干扰，可用,YXN来表示。噪声是随机的，在许多情况下这些很难测出或描述，甚至不可能得到，因而需用随机过程来描述，即要求知道其分布函数和密度函数，所以常用统计特征来描述噪声，如均值、方差、相关函数等。描述噪声的总功率,2YXNE方差，描述噪声的交流功率2,YXNEYXNE均值的平均，表示噪声的直流功率2,YXNE图像噪声可分为外部噪声和内部噪声。L外部噪声从处理系统以外来的影响，如天线的干扰或电磁波从电源线窜入系统的噪声。2内部噪声有四种基本形式由光和电的基本性质引起如电流可看作电子或空穴运动，这些粒子运动产生随机散粒噪声；导体中电子流

16、动的热噪声；光量子运动的光量子噪声等。机械运动产生韵噪声接头振动使电流不稳，磁头或磁带、磁盘抖动等。元器件噪声如光学底片的颗粒噪声，磁带、磁盘缺陷噪声，光盘的疵点噪声等。系统的内部电路噪声如CRT的偏转电路二次发射电子等噪声。从噪声的分类来看是多种多样的，但从统计的观点来看，凡是统计特征不随时间变化的称作平稳噪声，统计特征随时间变化的称作非平稳噪声。从噪声的幅运动模糊图像的恢复与处理度分布的统计特征来看，其密度函数有高斯型、瑞利型，分别称为高斯噪声和瑞利噪声。高斯噪声的概率密度函数为（21）222/21ZEZP（21）兰州理工大学毕业设计5式（21）中Z表示灰度级，表示Z的平均值或期望值，表示

17、Z的标准差。标准差的平方2称为Z的方差。当Z服从上式的分布时，其值有70落,范围内，且有95落在2,2范围内。瑞利噪声的概率密度函数为22AZAZAZBZPBAZ02/2（22）其中均值和方差分别为4/BA442B按噪声对信号的影响可分为加性噪声模型和乘性噪声模型两大类。设,YXF为信号，,YXN外为噪声，影响信号后的输出为,YXG。L加法性噪声,YNNYXFYXG（23）形成波形是噪声和信号的叠加，其特点是,YXN对和信号无关，如一般的电子线性放大器，不论输入信号大小，其输出总是与噪声相叠加。2乘法性噪声,1,YXNYXFYXFYXNYXFYXG（24）其输出是两部分的叠加，第二个噪声项信号

18、受,YXF的影响。,YXF越大，则第二项越大，即噪声项受信号的调制。如光电子噪声、底片颗粒噪声都随信号增大而增大。乘法性噪声模型和分析计算都比较复杂，通常信号变化很小时，第二项近似不变，此时可以用加法性噪声模型来处理。通常总是假定信号和噪声是相互独立的。22图象退化模型221退化模型要进行图像恢复，必须弄清楚退化现象有关的某些知识，用相反的过程去掉它，这就要兰州理工大学毕业设计6了解、分析图像退化的机理，建立起退化图像的数学模型15。一些退化因素只影响一幅图像中某些个别点的灰度，而另外一些退化因素则可以使一幅图像中的一个空间区域变得模糊起来。前者称为点退化，后者称为空间退化。在一个图像系统中存

19、在着许多退化源，其机理比较复杂，因此要提供一个完善的数学模型是比较复杂和困难的。但是在通常遇到的很多实例中，我们将退化原因作为线性系统退化的一个因素来对待，从而建立系统退化模型来近似描述图像函数的退化。如图21所示，这是一种简单的通用图像退化模型，输入图像,YXF经过一个退化系统或退化算子H后产生的退化图像,YXG，我们可以表示为（25）的形式。,YXNYXFHYXG（25）式中H为退化系统图21图像退化模型如果暂不考虑加性噪声。,YXN的影响，即令。0,YXN，则有（26）,YXFHYXG（26）设K，1K，2K为常数，,11YXFHYXG，,22YXFHYXG，则退化系统H具有如下性质L齐

20、次性,YXKGYXFKHYXKFH）（27）即系统对常数与任意图像乘积的响应等于常数与该图像的响应的乘积。2叠加性,2121YXFHYXFHYXFYXFH,21YXGYXG（28）H（X，Y）F,YXF,YXG,YXN兰州理工大学毕业设计7即系统对两幅图像之和的响应等于它对两个输入图像的响应之和。3线性同时具有齐次性与叠加性的系统就称为线性系统。线性系统有式（29）,22112211YXFHKYXFHKYXFKYXFKH,2211YXGKYXGK（29）不满足齐次性或叠加性的系统就是非线性系统。显然，线性系统为求解多个激励情况下的响应带来很大方便。4位置空间不变性，有式（210）,BYAXFH

21、BYAXG（210）式中的A和B占分别是空间位置的位移量。这就说明了图像上任何一点通过该系统的响应只取决于在该点的灰度值，而与该点的坐标位置无关由上述基本定义可知，如果系统具有式210的关系，那么系统就是线性空间不变的系统。在图像恢复处理中，尽管非线性和空间变化的系统模型具有普遍性和准确性。但是，它却给处理工作带来巨大的困难，通常没有解或者很难用计算机来处理15。因此在图像恢复处理中，往往用线性和空间不变性的系统模型加以近似。这种近似的优点是可直接利用线性系统中的许多理论与方法来解决图像恢复问题。所以图像恢复处理中主要采用线性的、空间不变的恢复技术。222连续函数退化模型空间坐标位置和景物明暗

22、程度均为连续变化的图像，称为连续图像。在图像线性运算的分析中，常常用到点源的概念。事实上，一幅图像可以看成由无穷多极小的像素所组成，每一个像素都可以作为一个点源。在数学上，点源可以用狄拉克石函数来表示，二维占函数可定义为式（211）000000,1,0XXYYDXDYXXYYXXYY其它211如果二维单位冲激信号沿X轴和Y轴分别有位移0X和0Y，则如式（212）,10,0,0XYDXDYXYXY其它212兰州理工大学毕业设计8,YX具有取样特性。由式211和212很容易得（213）Y,FXDXDYYY,XXYFX,0000（213）此外，任意二维信号,YXF与,YX卷积的结果就是该二维信号本身

23、，即（214）,YXFYXYXF（214）而任意二维信号,YXF与,00YYXX卷积的结果就是该二维信号产生相应位移后的结果（215）,0000YYXXFYYXXYXF215由二维卷积定义，有（216）,YXYXFYXF216考虑退化模型中韵H是线性空间不变系统，因此，根据线性系统理论，系统H的性能就可以由其单位冲撤响应,YXH来表征，即（217）,YXHYXH217而线性空间不变系统H对任意输入信号,YXF的响应则为该信号与系统的单位冲激响应的卷积为（218）FX,Y,DDYXHF218在不考虑加性噪声的情况下，上述退化模型的响应为（219）,DDYXHFYXFHYXG219由于系统H是空间

24、不变的，则它对移位信号的响应为（220）,0000YYXXGYXHYYXXF220在有加性噪声的情况下，上述线性退化模型可以表示为（221）,YXNYXFHYXG兰州理工大学毕业设计9,YXNDDYXHF221简记为（222）,YXNYXHYXFYXG222在上述情况中，都假设噪声与图像中的位置无关。式219和式222都是连续图像的退化模型。由此可见，如果把降质过程看成为一个线性空间不变系统，那么，在不考虑噪声影响时，系统输出的退化图像,YXG应为输入原始图像,YXF和引起系统退化图像的点扩散函数,YXH的卷积。因此，系统输出或影像被其输入景物和点扩散函数唯一确定。显然，系统的点扩散函数是描述

25、图像系统特性的重要函数。223离散函数退化模型为了用数字计算机对图像进行处理，首先必须把连续图像函数,YXF进行空间的和幅值的离散化处理空间连续坐标,YX的离散化，称为圈像的采祥，幅值,YXF的离散化称为灰度级的整量。将这两种离散化和在一起，称为图像的数字化。如图22所示，连续的模拟图像经过离散化处理后变成计算机能够辨识的点阵图像，称为数字图像。严格的数字图像是一个经过等距离矩形网格采样，对幅度进行等间隔量化的二维函数。将一幅图像进行数字化的过程就是在计算机内生成一个二维矩阵的过程15。图22离散退化模型数字图像可以由以下三种途径得到1将传统的可见光图像经过数字化处理转换为数字图像，例如将一幅

26、照片通过扫描仪输入到计算机中，扫描的过程实质上就是一个数字化的过程。2应用各种光电转换设备直接得到数字图像，例如卫星上搭载的推帚式扫描仪和光机扫描仪可以直接获取地表甚至地下物体的图像并实时存入存储器中。H采样C,YXF,YXN,21NNG兰州理工大学毕业设计103直接由二维离散数学函数生成数字图像无论哪种方式，最终得到的数字图像都是一个二维矩阵。对于一幅连续图像,YXF，若X，Y方向的相等采样间隔分别为XY，并均取N点，则数字图像,JIF。可用如下矩阵表示（223）,JIF1,11,10,11,11,10,11,01,00,0NNFNFNFNFFFNFFF223图像像素矩阵的产生，为图像处理提

27、供了一种新的途径，对于许多图像的处理，都可以转化为对矩阵的分析，从而使问题变得准确、简便、易行。数字图像处理实质就是对二维矩阵的处理，是将一幅图像变为另一幅经过修改的图像，是将一个二维矩阵变为另一个二维矩阵的过程。首先讨论一维的情况，然后再推广至二维情况。假设对两个函数XF和XH进行均匀采样，其结果放到尺寸为A和B的两个数组中，X的取值范围是0,1,2,1A；对XH，X的取值范围是0,1,2,1B。我们可以利用离散卷积来计算XG。为了避免卷积的各个周期重叠，并将函数用零扩展补齐。用XFE和XHE来表示扩展后的函数，则有224和（225）XFE1010MXAAXXF224XHE1010MXBBX

28、XH225则它们的卷积为226XGE10MMEEMXHMF226因为XFE和XHE的周期为M，XGE的周期也为M。引入矩阵表示法，则式（226）可写为227HFG227其中兰州理工大学毕业设计11110MGGGGEEE228110MFFFFEEE229021201110EEEEEEEEEHMHMHMHHHMHHHH）（230根据XHE的周期性可知，XHEMXHE，所以上式又可以写成（231）021201110EEEEEEEEEHMHMHHHHHMHHH）（231H是个循环矩阵，即每行最后一项等于下一行的最前一项，最后一行最后一项等于第一行最前一项。将一维结果推广到二维，可首先做成大小NM的周期

29、延拓图像，即,YXFE1101010,NYBMXABYAXYXF，232,YXHE1101010,NYDMXCDYCXYXH，233这样延拓后，,YXFE和,YXHE分别为二维周期函数。它们在X和Y方向上的周期分别为M和N。于是得到二维退化模型为一个二维卷积形式234,YXGE1010,MMEENNNYMXHNMF234如果考虑噪声将噪声项加上，上式可写成为（235）兰州理工大学毕业设计12,YXGE1010,MMEEENNYXNNYMXHNMF235同样，可以用矩阵来表示（236）021201110HHHHHHHHHNHFGMMM110MNFFFEEE110MNNNNEEE236其中每个IH

30、是由扩展函数气,YXHE的第I行而来，即（237）0,2,1,2,0,1,1,1,0,IHNIHNIHIHIHIHIHNIHIHHEEEEEEEEEI）（237237是一个循环矩阵。因为H中的每块是循环标注的，所以H是块循环矩阵。224匀速直线运动图像的退化模型在所有的运动模糊中，由匀速直线运动造成图象模糊的复原问题更具有一般性和普遍意义。因为变速的、非直线运动在某些条件下可以被分解为分段匀速直线运动。假设图象有一个平面运动,令XT和YT分别为在X和Y方向上运动的变化分量，T表示运动的时间。记录介质的总曝光量是在快门打开后到关闭这段时间的积分。则模糊后的图象为238DTTYY,TXXFYGX,

31、T000（238）式（238）中GX,Y为模糊后的图象。以上就是由于目标与摄像机相对运动造成的图象模糊的连续函数模型。如果模糊图象是由景物在X方向上作匀速直线运动造成的，则模糊后图象任意点的值为239DTYTXXFGTYX00,（239）式（239）中TX0是景物在X方向上的运动分量，若图象总的位移量为A，总的时间为兰州理工大学毕业设计13T，则运动的速率为TX0AT/T。则上式变为240DTYTATXFYXGT0,（240）以上讨论的是连续图象，对于离散图象来说，对上式进行离散化得241TYTATXFYXGLI10,（241）其中L为照片上景物移动的像素个数的整数近似值。是每个像素对模糊产生

32、影响的时间因子2。由此可知，运动模糊图象的像素值是原图象相应像素值与其时间的乘积的累加。从物理现象上看，运动模糊图象实际上就是同一景物图象经过一系列的距离延迟后再叠加，最终形成的图象。如果要由一幅清晰图象模拟出水平匀速运动模糊图象，可按式242进行10,1,LIYXFLYXG（242）这样可以理解此运动模糊与时间无关，而只与运动模糊的距离有关，在这种条件下，使实验得到简化。因为对一幅实际的运动模糊图象，由于摄像机不同，很难知道其曝光时间和景物运动速度。我们也可用卷积的方法模拟出水平方向匀速运动模糊。其过程可表示为243,YXHYXFYXG243其中其它0101,LXLYXH244HX,Y称为模

33、糊算子或点扩散函数，表示卷积，,YXF表示原始清晰图象，,YXG表示观察到的退化图象。如果考虑噪声的影响，运动模糊图象的退化模型可以描述为一个退化函数和一个加性噪声项,YXN，处理一幅输入图象,YXF产生一幅退化图象,YXG。,YXNYXHYXFYX245由于空间域的卷积等同于频率域的乘积，所以式245的频率域描述为246,VUNVUFVUHVUG246式246中的大写字母项是式245中相应项的傅里叶变换。兰州理工大学毕业设计1423图像的恢复方法231逆滤波复原法对于线性移不变系统而言247,YXG,YXF,YXNYXH247上式两边进行傅里叶变换得248,VUG,VUH,VUF,VUN24

34、8式中,VUG，,VUF，,VUH和,VUN分别是,YXG,YXF,YXH和,YXN的二维傅里叶变换。通常在无噪声的理想情况下，上式可简化249,VUG,VUF,VUH则,VUF,VUG/,VUH249,1VU称为逆滤波器。对式（249）再进行傅里叶反变换可得到,YXF。但实际上碰到的问题都是有噪声，因而只能求,VUF的估计值250,VUF,VUF,VUHVUN250然后再作傅里叶逆变换得251,YXFYXF,1VUHVUNDUDVEVYUXJ2251这就是逆滤波复原的基本原理。其复原过程可归纳如下对退化图像,YXG作二维离散傅里叶变换，得到,VUG计算系统点扩散函数,YXH的二维傅里叶变换，

35、得到,VUH。（这一步值得注意的是，通常,YXH的尺寸小于,YXG的尺寸。为了消除混叠效应引起的误差，需要把,YXH的尺寸延拓。计算,VUF的傅里叶变换，求得,YXF。逆滤波复原法的缺陷，0,VUH无确定,VUF，0,VUH放大噪声。若噪声为零，则采用逆滤波恢复法能完全再现原图像。若噪声存在，而且,VUH很小或为零时，则噪声被放大。这意味着退化图像中小噪声的干扰在,VUH较小时，会对逆滤波恢复的图像产生很大的影响，有可能使恢复的图像和,YXF相差很大，甚至面目全非。逆滤波复原法解决方法兰州理工大学毕业设计15解决该病态问题的唯一方法就是避开,VUH的零点即小数值的,VUH两种途径一是在0,VU

36、H及其附近，认为地仔细设置的值。二是使,VUH具有低通滤波性质。,1VUH202220220U1DVUDVUVH，252232约束最小平方复原法约束最小平方复原是一种以平滑度为基础的图像复原方法。如前所述，在进行图像恢复计算时，由于退化算子矩阵H的病态性质，多数在零点附近数值起伏过大，使得复原后的图像产生了多余的噪声和边缘。约束最小平方复原仍然是以最小二乘方滤波复原公式为基础，通过选择合理的Q，并优化2QF，从而去掉被恢复图像的这种尖锐部分，即增加图像的平滑性。我们知道，图像增强的拉普拉斯算子,它具有突出边缘的作用，则恢复了图像的平滑性，因此，在作图像恢复时可将其作为约束。现在的问题是如何将其

37、表示成2QF的形式，以便使用式253。在离散情况下，拉普拉斯算子可用式（253）的差分运算实现,2222YYXFXYXF（,41,1,1,1YXFYXFYXFYXFYXF（253）利用,YXF与式（254）的模板算子进行卷积可实现式（253）的运算,YXP010141010254在离散卷积的过程中，可利用延伸,YXF和,YXP来避免交叠误差。延伸后的函数为,YXPE。建立分块循环矩阵，将平滑准则表示为矩阵形式（255）C032130122101CCCCCCCCCCCCMMMM255兰州理工大学毕业设计16式255中每个子矩阵JC1,1,0MJ是,YXPE的第J行组成的NN循环矩阵。即JC如式（

38、256）表示JC0,2,1,2,0,1,1,1,0,JPNJPNJPJPJPJPJPNJPJPEEEEEEEEE）256根据循环矩阵的对角化可知，可利用前述的矩阵W进行对角化，即（257）E1WCW257式中，E为对角矩阵，其元素为（258）KIKIKKODNNKPIKE0,258则TFTCCF2QF，两边同乘以1W，得（259）1WF1EEDDGWD1259式中，D为D的共轭矩阵。所以有（260）,24222VUGVUPNVUHNVUHNVUF260式中，1,1,0UNV，而且2,VUH,VUHVUH。本滤波器也称为最小平方滤波器。233维纳滤波复原法维纳滤波法是由WIENER首先提出的，应

39、用于一维信号处理，取得了很好的效果。之后，维纳滤波法被用于二维信号处理，也取得了不错的效果，尤其在图像复原领域，由于维纳滤波计算量小，复原效果好，从而得到了广泛的应用和发展8。维纳滤波器寻找一个使统计误差函数达到最小的准则函数来实现图像复原的。22FFEE261兰州理工大学毕业设计17式（261）中，E表示数学期望。设FR和RN分别是F和N的自相关矩阵，定义如（262）FRETFF262RNTENN263根据上述定义可知，FR和RN均为实对称矩阵。在大多数实际图像中，相近像素点是高度相关的，而距离教远的像素点的相关性则相对较弱。通常情况下，无论是F还是N，其元素之间的相关不会延伸到2030个像

40、素的距离之外，因此。一般来说，自相关矩阵FR和RN在主对角线附近有一个非零元素区域，而矩阵的右上角和左上角的区域内将接近零值。如果像素之间的相关是像素距离的函数，而不是像素位置的函数，则可将FR和RN近似分为线循环矩阵。因而，用循环矩阵的对角化，可写成如（264）形式1FRWAW2641NRWBW265W为MNMN矩阵，包含MM个NN子矩阵。以WI，M表示W的I和M列分块矩阵，则（266）2,WEJIMMIMWN266其中，I，M0，1,2，M1，WN是NN矩阵，以W（K，N）表示K行N列元素，则有2,WJKNMKNENK，N0,1,2，M1267矩阵A，B的元素分别为矩阵FR和RN中的自相关

41、元素的傅里叶变换，这些自相关的傅里叶变换分别定义为,EFXY和,NXYE的谱密度,RUVF和,PUVN。则（268）111111TTTFHHQQHGWDDWWABWWDWG268兰州理工大学毕业设计18第3章维纳滤波实现退化图像的复原31维纳滤波的基本原理311维纳滤波概述维纳WIENER滤波是用来解决从噪声中提取信号问题的一种滤波的方法。实际上这种线性滤波问题，可以看成是一种估计问题或一种线性估计问题。一个线性系统，如果它的单位样本响应为NH，当输入一个随机信号NX，且（31）NXNSNV31其中NS表示信号，NV表示噪声，则输出NY为（32）NYMMNXMH32我们希望NX通过线性系统NH

42、后得到的NY尽量接近于NS,因此称NY为NS的估计值，用NS表示，即（33）NYNS33图31维纳滤波器的输入一输出关系如图31所示。这个线性系统NH称为对于NS的一种估计器。实际上，式33的卷积形式可以理解为从当前和过去的观察值NX来估计信号的当前值NS。因此，用NH进行过滤的问题可以看成是一个估计问题。由于我们现在涉及的信号是随机信号，所以这样一种过滤问题实际上是一种统计估计问题。一般，从当前的和过去的观察值估计当前的信号值NYNS称为过滤或滤波从过去的观察值，估计当前的或将来的信号值NYNNS0N称为预测或外推从过去的观察值，估计过去的信号值NYNNS1N称为平滑或内插。因此维纳过滤又常

43、常被称为最佳线性过滤与预测或线性最优估计。这里所谓最佳与最优是以最小均方误差为准则的。这里只讨论过滤与预测问题。H（N）NXNSNVNSNY兰州理工大学毕业设计19如果我们以与S分别表示信号的真值与估计值，而用NE表示它们之间的误差，即（34）NENSNS34显然，NE可能是正的，也可能是负的，并且它是一个随机变量。因此，用它的均方值来表达误差是合理的，所谓均方误差最小即它的平方的统计平均值最小如（35）MIN2NEEMIN2NSNSE35采用最小均方误差准则作为最佳过滤准则的原因还在于它的理论分析比较简单，不要求对概率的描述。并且在这种准则下导出的最佳线性系统对其它很广泛一类准则而言也是最佳

44、的。312运动模糊参数的确定1算法理论分析假设快门的开启和关闭所用时间非常短，那么光学成像过程不会受到运动的干扰，图像也不会出现运动模糊退化现象。如果设T为曝光时间，则运动模糊退化模型为（36）000,TGXYFXXYYDTNXY（36）式（36）中G（X，Y）表示模糊退化图像，F（X，Y）表示原始图像。N（X，Y）表示噪声。首先考虑没有噪声的情况，对式36进行傅里叶变换得（37）0020,TJUXTVYTGUVFUVEDT（37）0020,TJUXTVYTHUVEDT（38）式（38）中HU，V表示退化图像的点扩散函数PSF10。假设当前图像做匀速直线运动，匀速直线运动模糊退化函数由式38变

45、换为（39）,SINJUAVBTHUVUAVBEUAVB（39）由于图像在PC机上存储为离散形式，需要将上述传递函数表示为离散表达式，设图像尺寸为MN，由二维离散傅里叶变换的公式得（310）兰州理工大学毕业设计20,SINUAVBJMNTUAVBHUVEUAVBMNMN（310）其中，U取值为0到M1，V取值为0到N1如（311）,VUHVUFVUG（311）当N为其它整数值时，H（U，V）0，从而G（U，V）0因此，G（U，V）的图像在非零整数的线上显示为黑色条纹黑色表示最小灰度，白色表示最大灰度。如果M，N为素数，虽然U，V在各自取值范围内无法为非零正整数，但对于一般图像其频谱图依然会呈现

46、规则的明暗条纹状12。这是由于SIN为周期函数，它在自己的前后两个半周期内呈现明显的递减和递增特性，从而也形成规则的明暗条纹。容易证明，退化图像频谱中条纹倾斜角度即为直线斜率所对应角度，可用公式表示为（312）TANANABM312TANTAN2NAM313默认图像频谱暗条纹方向与运动模糊的方向相垂直，由式313可以看出，仅当NM时，条纹角度与模糊角度是垂直的，但当所处理图片长和宽不相等时，简单认为模糊角度和条纹倾斜角度垂直是不准确的。而如果对图片进行不当的裁剪会破坏原始图像信息，尤其对于抓拍到的高速车辆图像，其背景静止而只有车辆运动，原始像素信息会的到较好的保留，如强行将图片修剪为正方形会对模糊参数的检测带来不利影响。而根据式313，对任意尺寸的图像，一旦检测出退化图像频谱条纹角度，就可以有效的确定运动模糊角度12。对图像频谱处理过程中，通常将图像通过循环移位方式把U0，V0移到中心位置，由点到直线的距离公式，中心点0,0到直线1AMBN的距离D为（314）221DABMN314由对称性，图像中心两个暗条纹之间的间距D2D，设图像的模糊长度为L，则SIN,COSALBL。令MN，得（315）222COSSINMDL（315）兰州理工大学毕业设计21由式3