图像复原方法综述

----骤。数字图像处理就是将图像信号转换成数字格式,并利用计算机进行加工和处理的过程。模型并对原始图像进行合理估计。本文主要介绍了图像退化的原因、图像复原技术的分类和目前常用的几种图像复原方法,性和有效性。关键词:图像退化、图像复原、维纳滤波、正则滤波、LR算法、盲区卷积、2、图像复原概述感光胶卷的非线性及胶片颗粒噪声以及电视摄像扫描的非线性等所引起的几何失真,都难免起的质量下降为图像退化。出现模糊、失真,出现附加噪声等。由于图像的退化，在图图像复原技术是图像处理领域中一类非常重要的处理技术,与图像增强等其他基本图像处理技术类似，也是以获取视觉质量某种程度的改善为目的,所不同的是图像复原过程实际复原的处理过程就是对退化图像品质的提升,并通过图像品质的提升来达到图像在视觉上的由于引起图像退化的因素众多,且性质各不相同，目前没有统一的复原方法,众多研究人员根据不同的应用物理环境，采用了不同的退化模型、处理技巧和估计准则,从而得到了不图像复原算法是整个技术的核心部分。目前,国内在这方面的研究才刚刚起步,而国外却----数字图像复原技术最早则是从对天文观测图像的后期处理中逐步发展起来的。其中一个成功例子是NASA的喷气推进实验室在1964年用计算机处理有关月球的照片。照片是在空间飞行器上用电视摄像机拍摄的,图像的复原包括消除干扰和噪声,校正几何失真和对比度损片是在相机移动过程中拍摄的，图像复原的主要目的就是消除移动造成的失真[2]。早期的复原方法有：非邻域滤波法,最近邻域滤波法以及效果较好的维纳滤波和最小二乘滤波等。随着数字信号处理和图像处理的发展,新的复原算法不断出现，在应用中可以根据具体情况加以选择。目前国内外图像复原技术的研究和应用主要集中于诸如空间探索、天文观测、物质研究、遥感遥测、军事科学、生物科学、医学影象、交通监控、刑事侦察等领域。如生物方面,主要是用于生物活体细胞内部组织的三维再现和重构，通过复原荧光显微镜所采集的细胞内部逐层切片图，来重现细胞内部构成;医学方面，如对肿瘤周围组织进行显微观察，以获取肿瘤安全切缘与癌肿原发部位之间关系的定量数据；天文方面,如采用迭代盲反卷积进行气动光学效应图像复原研究等。图像复原问题的有效性关键之一取决于描述图像退化过程模型的精确性。要建立图像的退化模型，则首先必须了解、分析图像退化的机理并用数学模型表现出来。在实际的图像处理过程中,图像均需以数字离散函数表示,所以必须将退化模型离散化[3]。-w-wABCDABABeeeeij=0--如果用矩阵表示上式，则可写为:现实中造成图像降质的种类很多，常见的图像退化模型及点扩展函数有如下情景[15]：)线性移动降质在拍照时，成像系统与目标之间有相对直线移动会造成图像的降质。水平方向线性移述：|l0其他式中,d是降质函数的长度。在应用中如果线性移动降质函数不在水平方向，则可类似质函数。(2)散焦降质当镜头散焦时，光学系统造成的图像降质相应的点扩展函数是一个均匀分布的圆形光斑。此时,降质函数可表示为:|l0其他式中，R是散焦半径。(3)高斯(Gauss)降质Gauss降质函数是许多光学测量系统和成像系统最常见的降质函数。对于这些系统,决T(6)----4、几种较经典的复原方法介绍图像复原算法有线性和非线性两类。线性算法通过对图像进行逆滤波来实现反卷积,这法保证图像的非负性。而非线性方法通过连续的迭代过程不断提高复原质量，直到满足预先设定的终止条件，结果往往令人满意。但是迭代程序导致计算量很大,图像复原时耗较长,有时甚至需要几个小时。所以实际应用中还需要对两种处理方法综合考虑，进行选择[4]。(1)维纳滤波法Wiener号处理,取得了很好的效果。之后，维纳滤波法被用于二维信号处理，也取得了不错的效果,尤其在图像复原领域,由于维纳滤波计算量小,复原效果好,从而得到了广泛的应用和发展。维纳滤波器寻找一个使统计误差函数(7)fE符，f是未退化的图像。该表达式在频域可表示为nnnfnn------(2)正则滤波法另一个容易实现线性复原的方法称为约束的最小二乘方滤波,在IPT中称为正则滤在最小二乘复原处理中,常常需要附加某种约束条件。例如令Q为f的线性算子，222g2gHf=n2的最小化问题，这种有附加条件的极值问题可以用拉格朗日乘数法寻找一个f，使下述准则函数为最小:2222----(3)Lucy-Richardson算法L-R算法是一种迭代非线性复原算法，它是从最大似然公式印出来的，图像用泊松分布加以模型化的。当下面这个迭代收敛时模型的最大似然函数就可以得到一个令人满意的方程:kIPTLRdeconvlucy------(4)盲去卷积在图像复原过程中，最困难的问题之一是,如何获得PSF的恰当估计。那些不以PSF为基础的图像复原方法统称为盲区卷积。它以MLE为基础的，即一种用被随机噪声所干扰的量进行估计的最优化策略。工具箱nd模拟实验结果如下:图像分别迭代5次，10次，20次的结果。------[1]罗军辉冯平等。MATLAB7．0在图像处理中的应用[M]，北京:机械工业出版社[2]刘维一,于德月,王肇圻等。用迭代法消除数字图像放大后的模糊[J]。光电子.13(4)：398~400。[3]陶洪。数字共焦显微技术及其图像复原算法研究四川大学硕士学位论文[4]冈萨雷斯等著。数字图像处理(MATLAB版)[M]，电子工业出版社2009。----程序附录清单(1)维纳滤波I＝imread(＇C:\DocumentsandSettings\Administrator\桌面\原始图.jpg');ubplot(231);Imshow(I)；I＝rgb2gray(I);imshow(I)；BlurredNoisy=im2uint8(Blurred);NPOW=sum(NP(：)/prod(size(noise)));NCORR＝fftshift(real(ifftn(NP)))；IP=abs(fftn(I)).^2；IPOW=sum(IP(：)/prod(size(noise)))；ICORR=fftshift(real(ifftn(IP)));NSR=NPOW/IPOW;subplot(233);imshow(BlurredNoisy,[])；title('A=BlurredandNoisy');subplot(234);imshow(deconvwnr(BlurredNoisy,PSF),[]);title('deconbwnr(A,PSF，NSR)＇);RRtitle(＇deconbwnr(A，PSF，NCORR,ICORR)')；subplot(236);imshow(deconvwnr(BlurredNoisy,PSF,NPOW,ICORR1),[]);title('deconbwnr(A，PSF,NPOW,ICORR_1_D)');closeallI=imread(＇C:\DocumentsandSettings\Administrator\桌面\原始图．jpg');imshow(I);imshow(I);--BlurredNoisy＝imnoise(imfilter(I，PSF),'gaussian＇,0,V);NOISEPOWER=V*prod(size(I));[JLAGRA]=deconvreg(BlurredNoisy,PSF,NOISEPOWER)；subplot(233);imshow(BlurredNoisy);title('A=BlurredandNoisy');subplot(234);imshow(J);title('[JLAGRA]＝deconvreg(A,PSF,NP)');subplot(235)；imshow(deconvreg(BlurredNoisy,PSF,[]，LAGRA/10))；title('deconvreg(A,PSF,[]，10*LAGRA');I＝imread('C:\DocumentsandSettings\Administrator\桌面\原始图.jpg')；imshow(I);I=rgb2gray(I)；imshow(I);BlurredNoisy=imnoise(imfilter(I,PSF),'gaussian＇，0,V);BlurredNoisy=double(BlurredNoisy);WT(5：end－4,5：end-4)=1;J1＝deconvlucy(BlurredNoisy,PSF);imshow(mat2gray(BlurredNoisy));title('A=BlurredandNoisy')；imshow(mat2gray(J1));subplot(235)；imshow(mat2gray(J2));subplot(236);imshow(mat2gray(J3));title(＇deconvlucy(A,PSF，NI，DP，WT)');--(4)盲去卷积clearallI=imread('C:\DocumentsandSettings\Administrator\桌面\原始图．jpg');imshow(I);imshow(I);PSF＝fspecial(＇gaussian',7,10)；WT=zeros(size(I))；[JP]=deconvblind(BlurredNoisy,INITPSF,5，10*sqrt(V),WT,FUN,0);[KP]=deconvblind(BlurredNoisy