基于MATLAB图像复原论文.doc

河南师范大学XXXX毕业设计河南师范大学XXXX学校本科毕业设计 学号： 0808224041 基于MATLAB的离焦模糊图像复原学院名称： 计算机与信息技术学院 专业名称： 通信工程 年级班别： 2008级31班 姓 名： 宋亚威 指导教师： 范海 菊 2012年5月II目 录目 录1前 言31 图像退化/复原处理的模型32 噪声的特征43噪声的分类54直接逆滤波65维纳滤波76 约束的最小二乘方滤波97 使用Lucy-Richardson算法的迭代非线性复原118 盲去卷积13总 结15参考文献16致 谢17基于MATLAB的离焦模糊图像复原摘要 摘 要 图像在获取、传输和存储过程中由于会受多种原因到如模糊、失真、噪声等原因的影响，这些原因会造成使图像的质量的下降1。因此，我们必须需要采取一定的方法尽可能地减少或消除图像质量的下降，恢复图像的本来面目，这就是称为图像复原，也称为图像恢复。在对目前有关数字图像复原技术的文献进行理解和综合的基础上通过阅读图像复原技术相关资料，本文主要对探讨了维纳(Wiener)滤波、约束最小二乘滤波算法、Lucy-Richardson算法和盲解卷积算法进行了探讨，并对其使用相关的工具箱函数deconvwnr函数、deconvreg函数、deconvlucy函数、deconvblind函数的用法进行简单介绍仿真。另外本文对上述算法进行了仿真实现，并分析对另外本文对上述算法进行了仿真实现，并分析了四种算法的实验结果。比了四种算法的实验结果。 关键词 图像复原；维纳滤波恢复；约束最小二乘滤波恢复；Lucy-Richardson恢复；盲解卷积恢复Based on the MATLAB of defocus blurred image restoration Abstract Images in access, transfer and storage process, due to various reasons such as fuzzy, distortion, the influence of noise and so on, the quality of image must be taken downImage in the acquisition, transmission and storage process will be subject to such as blurring, distortion, noise and other reasons, these reasons will make the image quality degradation., so, the method must as far as possible to reduce or eliminate image quality declines and restore images, this is image restoration.Therefore, we needed to take a certain amount of ways to reduce or eliminate image quality to fall, to restore the image of self, this is known as image restoration. For the current figures in the literature image restoration technology on the basis of understanding and comprehensive,By reading the image restoration technology related data. this paper focuses on the Wiener (Wiener) filtering, constrained least squares filtering algorithm, Lucy - Richardson algorithm and blind solution convolution algorithm,This paper mainly discusses the Wiener filter, constrained least squares filtering algorithm, Lucy-Richardson algorithm and blind deconvolution algorithm， and its relevant toolbox deconvwnr function, deconvreg function,deconvlucy function, deconvblind function have been given a brief introduction of usage. and the deconvwnr function，the deconvreg function ，the deconvlucy function and the deconvblind function are used for emulation.In addition, this paper also simulates the above algorithms and analyzes the experimental results of the four algorithm.This article on the above algorithm to simulation and experimental result analysis of four kinds of algorithms.录Keywords image restoration; Wiener filtering restore; Constrained least squares filtering restore; Lucy-Richardson recovery; Blind solution convolution recovery目录1 图像退化/复原处理的模型12 噪声的特征23噪声的分类34直接逆滤波45维纳滤波56 约束的最小二乘方滤波77 使用Lucy-Richardson算法的迭代非线性复原98 盲去卷积11总 结12参考文献13致 谢14前 言在实际的日常生活中，人们要接触很多图像，画面。而在景物成像这个过程里可能会出现模糊、失真或混入噪声，最终导致图像质量下降，这种现象称为图像“退化”。因此我们可以采取一些技术手段来尽量减少甚至消除图像质量的下降，还原图像的本来面目，即在预定义的意义上改善给定的图像，这就是图像复原。尽管图像增强和图像复原之间有重叠部分，但前者主要是主观处理，而图像复原大部分是客观处理。复原通过使用退化现象的先验知识试图重建或恢复一副退化的图像。因此，复原技术趋向于将退化模型化并用相反的处理来恢复原图像，即考虑用模糊函数来消除图像的模糊。引起图像模糊有多种多样的原因，举例来说有运动引起的，高斯噪声引起的，斑点噪声引起的，椒盐噪声引起的等等。 本文主要研究离焦模糊图像的复原，离焦模糊图像是指在拍摄时景物与相机的相对运动引起的离焦 ,或是成像区域内不同深度的对象所引起不同程度的离焦 ,还有由于在成像区域中存在不同深度的对象会使自动调焦系统引起混淆而导致拍摄的相片离焦等。因此本文研究使用MATLAB把退化现象模型化，并对几种常用的滤波方法用MATLAB进行了仿真实现，为人们在不同的应用场合及不同的图像数据条件下选择不同的复原算法提供了一定的依据1 图像退化/复原处理的模型如图1所示， 退化 复原 图1本文中用退化函数把退化过程模型化，它和加性噪声项一起，作用于输入图像，产生一副退化的图像： (1.1)给定、一些关于退化函数H的知识以及一些关于加性噪声的知识，复原的目标就是得到原图像的一个估计。我们要是这个估计尽可能地接近原始的输入图像。通常，我们对和知道越多，就越接近。若H是线性的、空间不变的过程，则退化图像在空间域通过下式给出： (1.2)其中，是退化函数的空间表示，且空间域的卷积和频域的乘法组成了一个傅立叶变换对，所以可以用等价的频域表示写出恰面的模型： (1.3)其中，用大写字母表示的项是卷积方程式中相应项的傅立叶变换。退化函数有时称为光学传递函数（DTF）。在空间域，成为点扩散函数（PSF）。对于任何种类的输入，让作用于光源的一个店来得到退化的特征，点扩散函数就是来源此的一个名词。由于退化是线性的，所以空间不变的退化函数可以被模型化为卷积，同样地，复原处理有时也称为反卷积。2 噪声的特征对图像信号来说，可将黑白图像看作二维亮度分布，则噪声可看作是对亮度的干扰，用来表示。噪声是随机的，因而需要用随机过程来描述，即要求知道其分布函数和概率密度函数。在许多情况下，这些函数很难测定和描述，甚至不能得到，所以常用统计特征来描述噪声，如均值、方差和相关函数等。图像噪声是图像在摄取或传输时所受的随机信号干扰，是图像中各种妨碍人们对其信息接受的因素。很多时候将图像噪声看成是多维随机过程，因而描述噪声的方法完全可以借用随机过程的描述，即用其概率分布函数和概率密度分布函数。 在许多情况下，这些函数很难测定和描述，甚至不能得到，所以常用统计特征来描述噪声，如均值、方差和相关函数等。均方值，描述噪声的总功率： 方差，描述噪声的交流功率：均值的平方，表示噪声的直流功率： 3噪声的分类1、按干扰源分类 图像噪声按照其干扰源可以分为外部噪声和内部噪声。 外部噪声是指从处理系统以外来的影响，例如天线干扰或电磁波从电源线窜入系统的噪声。 内部噪声主要有四种基本形式： (1)由光和电的基本性质引起：如电流可看作电子或空穴的运动，这些粒子的运动将产生随机散粒噪声；导体中电子流动的热噪声；光量子运动的光量子噪声等。 (2)机械运动产生的噪声：接头震动造成的电流不稳；磁头或磁带、磁盘抖动造成的噪声等 。(3)元器件噪声：如光学底片的颗粒噪声；磁带、磁盘缺陷噪声；光盘的疵点噪声等。 (4)系统内部电路噪声：如CRT的偏转电路二次发射电子等噪声。 外部噪声，即指系统外部干扰以电磁波或经电源串进系统内部而引起的噪声。如电气设备，天体放电现象等引起的噪声。 内部噪声：一般又可分为以下四种： （1） 由光和电的基本性质所引起的噪声。如电流的产生是由电子或空穴粒子的集合，定向运动所形成。因这些粒子运动的随机性而形成的散粒噪声；导体中自由电子的无规则热运动所形成的热噪声；根据光的粒子性，图像是由光量子所传输，而光量子密度随时间和空间变化所形成的光量子噪声等。（2） 电器的机械运动产生的噪声。如各种接头因抖动引起电流变化所产生的噪声；磁头、磁带等抖动或一起的抖动等。 （3） 器材材料本身引起的噪声。如正片和负片的表面颗粒性和磁带磁盘表面缺陷所产生的噪声。随着材料科学的发展，这些噪声有望不断减少，但在目前来讲，还是不可避免的。 （4） 系统内部设备电路所引起的噪声。如电源引入的交流噪声；偏转系统和箝位电路所引起的噪声等。2、按按对信号的影响分类噪声与信号的关系分 按照噪声对信号的影响这里我们可以将噪声可分为加性噪声模型和乘性噪声模型两大类。设为信号，为噪声，影响信号后的输出为。 (1)加性噪声 加性噪声 (3.1)加性嗓声和图像信号强度是不相关的，如运算放大器，又如图像在传输过程中引进“信道噪声”电视摄像机扫描图像的噪声的，这类带有噪声的图像可看成为理想无噪声图像与噪声之和。形成的波形是噪声和信号的叠加，其特点是和信号无关。如一般的电子线性放大器，不论输入信号的大小，其输出总是与噪声相叠加的。 (2)乘性噪声 (3.2) 其输出是两部分的叠加，第二个噪声项受信号的影响越大则第二项越大，即噪声项受信号的调制。如光量子噪声、底片颗粒噪声都随信号增大而增大。乘性噪声模型和分析计算都比较复杂，通常当信号变化很小时，第二项近似不变，此时可用加性噪声模型来处理。通常总是假定信号和噪声是相互独立的。乘性嗓声和图像信号是相关的，往往随图像信号的变化而变化，如飞点扫描图像中的嗓声、电视扫描光栅、胶片颗粒造成等，由于载送每一个象素信息的载体的变化而产生的噪声受信息本身调制。在某些情况下，如信号变化很小，噪声也不大。为了分析处理方便，常常将乘性噪声近似认为是加性噪声，而且总是假定信号和噪声是互相统计独立。3、 按统计特征分类按概率密度函数分 按照噪声幅度分布的统计特性来看又可分为以下几种：这是比较重要的，主要因为引入数学模型，这就有助于运用数学手段去除噪声。 (1)白噪声（White Noise）：具有常量的功率谱。白噪声的一个特例是高斯噪声（Gaussian Noise）。在空间域和频域中，由于高斯噪声在数学上的易处理性，这种噪声(也称为正态噪声)模型经常被用在实践中。事实上，这种易处理性非常方便，使高斯模型经常适用于临街情况下。它的直方图曲线服从一维高斯型分布： (3.3)(2)椒盐噪声（Pepper Noise）：一种在图像中产生黑色、白色点的脉冲噪声。椒盐噪声是由图像传感器，传输信道，解码处理等产生的黑白相间的亮暗点噪声，往往由图像切割引起。椒盐噪声是指两种噪声，一种是盐噪声（salt noise），另一种是胡椒噪声（pepper noise）。盐=白色，椒=黑色。前者是高灰度噪声，后者属于低灰度噪声。一般两种噪声同时出现，呈现在图像上就是黑白杂点。该噪声在图像中较为明显，对图像分割、边缘检测、特征提取等后续处理具有严重的破坏性。(3)冲击噪声（Impulsive Noise）：指一幅图像被个别噪声像素破坏，而且这些噪声像素的亮度与其领域的亮度明显不同。 冲击噪声呈突发状，常由外界因素引起；其噪声幅度可能相当大，无法靠提高信噪比来避免，是传输中的主要差错。(4)量化噪声（Quatization Noise）：是指在量化级别不同时出现的噪声。例如,将图像的亮度级别减少一半的时候会出现伪轮廓。4直接逆滤波 用于复原一副退化图像的最简单方法是构成一个如下形式的估计： (4.1) 然后采用的傅里叶逆变换来得到图像的相应估计。这种方法成为逆滤波。由前文中讨论的，模型，我们可以将估计表示为 (4.2) 这个简单的表达式告诉我们，及时准确的知道了，也不能恢复，因为噪声分量是一个随机函数，他的傅里叶变换是未知的。另外，在实际中，有许多为零的情况也是个问题。即使项也可以忽略，用为零的值来除它也将会控制复原估计。采用逆滤波时，典型的方法是形成比率，然后为了得到逆，将频率的范围限制在接近原点的频率。概念是中的零不太可能在接近原点的地方出现，因为变换的数值通常是该区域中的最高值。5维纳滤波维纳滤波（N.Wiener最先在1942年提出的方法）是一种最早也最为人们所熟知的线性图像复原方法。维纳滤波器寻找一个使统计误差函数 (5.1)最小的估计。其中，是期望值操作符，是未退化的图像。该表达式在频域可表示为 (5.2)式中，表示退化函数 表示的复共轭 表示噪声的功率谱 表示未退化图像的功率谱比率称为噪信功率比。我们看到，若对于和的所有相关值，噪声功率谱为零，则这个比率就为零，且维纳滤波器就成为前文讨论的逆滤波器。我们感兴趣的两个量平均噪声功率和平均图像功率，分别定义为 (5.3) (5.4)其中，和分别表示图像和噪声数组的垂直和水平大小。这些量都是标准量，且它们的比率是一个标量，有时用来代替函数，以便产生一个常量数组。在这种情况下，即使真实的比率未知，交互式地变化常量并观察复原的结果的实验就变成了一件简单的事。当然，假设函数为常量是一种粗糙的的近似。在前述滤波器方程中，用一个常量数组来代替就产生所谓的参数维纳滤波器。在IPT中，维纳滤波是使用函数deconvwnr来实现的，函数deconvwnr有三种可能的语法形式。在所有的这些形式中，g代表退化图像，fr是复原图像。第一种语法形式fr=deconvwnr（g，PSF）假设信噪比为零。从而，维纳滤波器的这种形式就是前文中的逆滤波器。语法 fr=deconvwnr（g，PSF，NSPR）假设噪信功率比已知，或是个常量或是个数组；函数接受其中的任何一个。这是用于实现参数维纳滤波器的语法，在这种情况下，NSPR可以是一个交互的标量输入。最后，语法 fr=deconvwnr（g，PSF，NACORR,FAVORR）假设噪声和未退化图像的自相关函数NACORR和FAVORR是已知的。注意deconvwnr的这种形式使用和的自相关来代替这些函数的功率谱。从相关理论我们可知 (5.5)其中，“”表示相关操作，表示傅立叶变换。这个表达式说明了对于deconvwnr的使用，通过计算功率谱的傅里叶变换，可以得到自相关函数。噪声的自相关有类似的注释。若复原图像呈现出由算法中使用的离散傅立叶变换所引入的振铃，则它往往会提醒我们在调用J=edgetaper(I,PSF)该函数利用点扩散函数PSF模糊了输入图像I的边缘。这样，输出图像J就是图像I和I的模糊版本的加权和。这个由PSF的自相关函数所决定的加权数组在它的中心区域取J等于I，而在接近边缘的地带等于I的模糊版本。使用deconvwnr函数复原模糊噪声图像 f=checkerboard(60); PSF=fspecial(motion,30,100); gb=imfilter(f,PSF,circular); noise=imnoise(zeros(size(f),gaussian,0,0.001); g=gb+noise; fr1=deconvwnr(g,PSF); Sn=abs(fft2(noise).2; nA=sum(Sn(:)/prod(size(noise); Sf=abs(fft2(f).2; fA=sum(Sf(:)/prod(size(f); R=nA/fA;deconvwnr函数复原模糊噪声图像： fr2=deconvwnr(g,PSF,R); NCORR=fftshift(real(ifft2(Sn); ICORR=fftshift(real(ifft2(Sf);f=imread(cameraman.jpg); LEN=30; THETA=45; PSF=fspecial(motion,LEN,THETA); MF=imfilter(f,PSF,circular,conv); wnr=deconvwnr(MF,PSF); subplot(2,2,1);imshow(f);title(原始图像);subplot(2,2,2);imshow(MF);title(模糊后的图像);subplot(2,2,3);imshow(wnr);title(恢复后的图像); fr3=deconvwnr(g,PSF,NCORR,ICORR); 5.1 模糊的噪声图像 5.2逆滤波的结果 5.3 使用常数比率的维纳滤波的结果 5.4使用自相关函数的维纳滤波的结果经过仿真，如图所示，我们得到的结果虽然仍有一些噪声存在，但是已经和原图很接近了。因为原图像和噪声函数都是已知的，所以可以正确地估算参量，并且图5.4便是在这种情况下能够由wiener反卷积所得到的最佳结果。在实践中，当这些量之一（或更多）未知时，挑战便是在试验中智能地选择所用的函数，知道获得可接受的结果为止。6 约束的最小二乘方滤波还有一种容易实现的线性复原的方法称为约束最小二乘方滤波，约束复原除要求了解关于退化系统的传递函数之外，还需要知道某些噪声的统计特性或噪声与图像的某些相关情况。在最小二乘约束复原中，要设法寻找一个最优估计，使得形式为的函数最小化.。在此准则下，可把图像的复原问题看作对求式( 6.1) 目标泛函的最小值。复原寻找最优估计，使为最小，且服从。因此准则函数可以写成： (6.1)其中：Q为的线性算子，对应一高通卷积滤波算子,表示对作某些线性操作的矩阵，通常选择拉普拉斯算子，如Laplacian算子，则且；为拉格朗日乘子（Langrange multiple）。为了使式（6.1）最小，对上式求导，并令为0就可以得到最小二乘解：推导出复原式为 （） ， (6.2)对应的频域表示为： (6.23) 这就式（6.3）构成了所谓约束最小二乘复原滤波算法，显然该算法无须获知原图像的统计值，便可以有效的实施最优估计，这点与维纳滤波明显不同。 约束最小二乘滤波恢复函数deconvreg的调用格式：deconvreg(I，PSF，NP，LRANGE，REGOP)，其中，I表示输入像，PSF表示点扩散函数，NP、LRANGE（ 输入） 和 REGO是可选参数，分别表示图像的噪声强度、拉氏算子的搜索范围和约束算子，同时，该函数也可以在指定的范围内搜索最优的拉氏算子。下面通过程序事例说明约束最小二乘滤波恢复图像：I=imread (lena.jpg); PSF=fspecial(gaussian,10,4); Blurred=imfilter(I,PSF,conv); V=.03;BN=imnoise(Blurred,gaussian,0,V); NP=V*prod(size(I); I=imread (cameraman.jpg); PSF=fspecial(gaussian,10,4); Blurred=imfilter(I,PSF,conv); V=.03;BN=imnoise(Blurred,gaussian,0,V); NP=V*prod(size(I); reg LAGRA=deconvreg(BN,PSF,NP); Edged=edgetaper(BN,PSF); reg2=deconvreg(Edged,PSF,NP/1.2); reg3=deconvreg(Edged,PSF,LAGRA); figuresubplot(2,3,1);imshow (I);title(原始图像);subplot(2,3,2);imshow (BN);title(加入高斯噪声的图像);subplot(2,3,3);imshow (reg);title(恢复后的图像);subplot(2,3,4);imshow(reg2);title(振铃抑制图像);subplot(2,3,5);imshow(reg3);title(拉格朗日算子恢复图像);reg LAGRA=deconvreg(BN,PSF,NP); Edged=edgetaper(BN,PSF); reg2=deconvreg(Edged,PSF,NP/1.2); reg3=deconvreg(Edged,PSF,LAGRA); figuresubplot(2,3,1);imshow (I);title(原始图像);subplot(2,3,2);imshow (BN);title(加入高斯噪声的图像);subplot(2,3,3);imshow (reg);title(恢复后的图像);subplot(2,3,4);imshow(reg2);title(振铃抑制图像);subplot(2,3,5);imshow(reg3);title(拉格朗日算子恢复图像);图6约束最小二乘恢复对比图 利用振铃抑制恢复图像是3几种中恢复效果最好的，其他几种方法也可以恢复但是比较模糊，效果不是很明显。7 使用Lucy-Richardson算法的迭代非线性复原 Lucy-Richardson(LR)算法假设图像服从Poission 分布，采用最大似然法进行估计，是一种基于贝叶斯分析的迭代算法。Lucy-Richardson算法是目前世界上应用最广泛的函数恢复技术之一，它是一种迭代非线性复原算法，能够按照泊松噪声统计标准求出与给定的点扩散函数PSF (Point Spread Function)卷积后最有可能成为输入模糊图像的图像。当PSF已知但图像噪声信息未知时，也可以使用这个函数进行有效的恢复。其最优估计以最大似然准则作为标准，即要使概率密度函数最大，推导出的迭代式为： (7.1) (7.1)其中，和分别为卷积运算和相关运算；为迭代次数，可以令进行迭代，可以证明，当噪声可以忽略当不断增大时会依概率收敛于，从而恢复出原始图像。 当噪声不可忽略时，把式(1.2)代入式(7.1)可得到： (7.2) 从上式可看出，若噪声不可忽略，则以上过程的收敛性将难以保证，即 LR 存在放大噪声的缺陷。因此，处理噪声项是 LR 算法应用于低信噪比图像复原的关键。 是未退化图像的估计，这个算法的迭代本质是显而易见的。它的非线性本质是载方程的右边用来除产生的。就像大多数非线性方法一样，关于什么时候停止L-R算法通常是很难回答。通常，接下来的处理途径是对于给定的应用，在我们获得满意的结果时，观察输出并终止算法。在IPT中，L-R算法是由名为deconvlucy的函数完成的，该函数的语法为 fr=deconvlucy(g,PSF,NUMIT,DAMPAR,WEIGHT)其中，fr代表复原的图像，g代表退化的图像，PSF是点扩散函数，NUMIT为迭代次数（默认为10次），DAMPAR和WEIGHT定义如下。DAMPAR是一个标量，它指定了结果图像与原图像g之间的偏离阈值。当像素偏离原值的范围在DAMPAR之内时，就不用再迭代。这既抑制了这些像素上的噪声，又保存了必要的图像细节。默认值为0（无衰减）。WEIGHT是一个与g同样大小的数组，它为每一个像素分配一个权重来反映其重量。当用一个指定的PSF来模拟模糊时，WEIGHT可以从计算像素中剔除那些来自图像边界的像素点，因此，PSF造成的模糊是不同的。若PSF的大小为，则在WEIGHT中用到的零边界的宽度是ceil（n/2）。默认值是同输入图像g同等大小的一个单位数组。若复原图像呈现出由算法中所用的离散傅里叶变换所引入的振铃，则在调用函数deconvlucy之前，要利用函数edgetaper。下面通过程序事例说明 Lucy-Richardson算法恢复图像：I=imread(lajiao000.jpg); PSF=fspecial(gaussian,5,5) ;Blurred=imfilter(I,PSF,symmetric,conv); V=.003; BN=imnoise(Blurred,gaussian,0,V); luc=deconvlucy(BN,PSF,5); figuresubplot(2,2,1);imshow(I);title(原始图像);subplot(2,2,2);imshow (Blurred);title(模糊后的图像);subplot(2,2,3);imshow (BN);title(加噪后的图像);subplot(2,2,4);imshow (luc);title(恢复后的图像); 图7.1 Lucy-Richardson恢复对比图用Lucy-Richardson算法可以较好的恢复图像9。但由于迭代产生的噪声痕迹是最大化可能性数据逼近法的常见问题，在低信噪比条件下，恢复图像可能会出现一些斑点，这些斑点并不代表图像的真实结构，只不过是恢复图像过于逼近噪声所产生的结果。另外本方法存在一些较严重的缺陷，问题一是噪声放大问题，这也是这类方法（利用最大似然性求解）的通病。问题二是对于恢复图像中的不同部分，分别执行多少迭代才合适的问题。因为图像中信噪比高的部分可能需要数百次迭代才能获得满意的结果；而另一些光滑的对象可能只需很少次数即可达到满意的结果，所以适当选择迭代次数对图像恢复也很重要。这两个问题若得不到解决，将会对最终结果产生不利影响从图中我们可以看到，。经多次迭代，尤其是在低信噪比情况下，重建图像可能会出现一些斑点，这些斑点并不代表图像的真实结构，是输出图像过于逼近噪声所产生的结果。8 盲去卷积 在图像复原过程中，最困呐的问题之一是，如何向前面所讨论的那样获得复原算法中使用的PSF的恰当估计。通常图像恢复方法均在成像系统的点扩展函数PSF已知下进行, 实际上它通常是未知的. 在 PSF未知的情况下, 盲去卷积是实现图像恢复的有效方法。因此，把那些不以PSF知识为基础的图像复原方法统称为盲去卷积算法。 在过去的20年里，一种盲去卷积的方法已经受到了人们的极大重视，它是以最大似然估计(MLE)为基础，即一种用被随机噪声所干扰的量进行估计的最优化策略。简要的说，关于MLE方法的一种解释就是将图像数据看成随机量，它们与另外一族可能的随机量之间有着某种似然性。似然函数用、和来加以表达，然后，问题就变成了寻求最大似然函数。在盲去卷积中，最优化问题规定的约束条件并假定收敛时通过迭代来求解，得到的最大和就是还原的图像和PSF。工具箱通过函数deconvblind来执行盲去卷积，它有如下语法： f,PSFe=deconvblind(g,INITPSF)其中，g代表退化函数，INITPSF是点扩散函数的出事估计。PSFe是这个函数最终计算到的估计值，fr是利用估计的PSF复原的图像。 用来去的复原图像的算法是L-R迭代复原算法。PSF估计受其初始推测尺寸的巨大影响，而很少受其值的影响。若复原图像呈现出由算法中使用的离散傅里叶变换所引入的振铃，则我们在调用函数deconvblind值钱，通常要使用函数edgetaper。下面我们使用函数deconvblind估计PSF： f=checkerboard(8); imshow(pixeldup(f,8); PSF=fspecial(gaussian,7,10); figure,imshow(pixeldup(PSF,73),); SD=0.01; g=imnoise(imfilter(f,PSF),gaussian,0,SD2); INITPSF=ones(size(PSF); NUMIT=5; LIM=ceil(size(PSF,1)/2); WEIGHT=zeros(size(g); WEIGHT(LIM+1:end-LIM,LIM+1:end-LIM)=1; DAMPAR=10*SD; fr,PSFe=deconvblind(g,INITPSF,NUMIT,DAMPAR,WEIGHT); imshow(pixeldup(fr,73),);I=imread(lajiao.jpg); PSF=fspecial(motion,10,30); Blurred=imfilter(I,PSF,circ,conv) ;INITPSF=ones(size(PSF); J P=deconvblind (Blurred,INITPSF,20); figuresubplot(2,2,1);imshow (I);title(原始图像);subplot(2,2,2);imshow (Blurred);title(模糊后的图像)subplot(2,2,3);imshow (J);title(初步恢复后的图像); 图8.1原图像 图8.2退化图像原始PSF 图8.3进行5次迭代估计的PSF 图8.4进行10次迭代 图8.5进行20次迭代上面利用盲去卷积进行了仿真，分别利用5次，10次，20次的迭代后所得到的PSF，显然相比之下，迭代次数越多越接近原始的PSF。该算法优点是,同时恢复了图像和点扩张函数,在对失真情况毫无先验知识的情况下, 仍能实现对模糊图像的恢复操作。利用 MATLAB实现的图像恢复, 并对恢复图像的失真情况做了改善。在进行图像恢复时,重建 PSF,对图像进行重建, 得到恢复的图像。总 结本文介绍了图像退化的原因并且简要介绍了当前主流的图像复原方法，并通过对各种复原方法的仿真，了解了各种方法的优劣性，为我们在实际生活提供依据。但是无论是哪一种方法都有所局限性，我们应该努力致力于研究新型的优秀的图像复原方法，来获得更好的图像复原效果。同时，我们知道总而言之，图像恢复的算法很多，无论何种算法，原则上，他们都依据或取得相关信息（包括关于退化系统、原图像、噪声等的确定性信息和统计信息）才能有效的实施。算法利用的信息越多信息的准确性越高，则复原图像的质量就越高。而且采用 MATLAB实现图像恢复，通过几条简单的MATLAB命令就可完成一大串高级计算机语言才能完成的任务,简捷明快。大多数图像处理模型是可以通过使用MATLAB的基本函数通过编程实现的。参考文献1张燕妮.，离焦模糊复原方法研究D硕士研究生论文，大连：大连理工大学，2005.32徐飞.施晓红 . MATLAB应用图像处理M西安电子科技大学出版社，2002.3路敬