用逆滤波和维纳滤波进行图像复原.pdfVIP

用逆滤波和维纳滤波进行图像复原 在图像的获取、传输以及记录保存过程中，由于各种因素，如成像设备与目标 物体的相对运动，大气的湍流效应，光学系统的相差，成像系统的非线性畸变，环 境的随机噪声等原因都会使图像产生一定程度的退化，图像退化的典型表现是图像 出现模糊、失真，出现附加噪声等。由于图像的退化，使得最终获取的图像不再是 原始图像，图像效果明显变差。为此，要较好地显示原始图像，必须对退化后的图 像进行处理，恢复出真实的原始图像，这一过程就称为图像复原。 图像复原技术是图像处理领域一类非常重要的处理技术，主要目的就是消除或 减轻在图像获取及传输过程中造成的图像质量下降即退化现象，恢复图像的本来面 目。 图像复原的过程是首先利用退化现象的某种先验知识，建立退化现象的数学模 型，然后再根据退化模型进行反向的推演运算，以恢复原来的景物图像。 一、 实 验目的 1 了解图像复原模型 2 了解逆滤波复原和维纳滤波复原 3 掌握维纳滤波复原、逆滤波的 Matlab 实现 二、实验原理 1、逆滤波复原 gxy,fxy,如果退化图像为，原始图像为，在不考虑噪声的情况下，其，，，， 退化模型可用下式表示 ，,，, gxyfxydd,,,,,,,,,,,,,，，，，，，,,,,,, (12-25) 由傅立叶变换的卷积定理可知有下式成立 GuvHuvFuv,,,, ，，，，，， (12-26) Guv,Huv,Fuv,gxy,式中，、、分别是退化图像、点扩散函数，，，，，，，， hxy,fxy,、原始图像的傅立叶变换。所以 ，，，， ，，Guv,，，,,11fxyFFuvF,,,,，， ，，，，,,，，Huv,，，，， (12-27) 由此可见，如果已知退化图像的傅立叶变换和系统冲激响应函数(“滤被”传 递函数)，则可以求得原图像的傅立叶变换，经傅立叶反变换就可以求得原始 fxy,Guv,Huv,图像，其中除以起到了反向滤波的作用。这就是逆滤波 复，，，，，， 原的基本原理 在有噪声的情况下，逆滤波原理可写成如下形式 GuvNuv,,，，，， Fuv,,,，，HuvHuv,,，，，， (12-28) 式中，Nuv,是噪声 nxy,的傅立叶变换 ，，，， 2、维纳滤波复原 ˆfxy,维纳滤波就是最小二乘滤波，它是使原始图像 fxy,与其恢复图 像，，，，之间的均方误差最小的复原方法。对图像进行维纳滤波主要是为了消除 图像中存 在的噪声，对于线性空间不变系统，获得的信号为 ，,，, gxyfhxyddnxy,,,,,,,，,,,,,,，，，，，，，，,,,,,, (12-29) ˆfxy,为了去掉 gxy,中的噪声，设计一个滤波器mxy,，其滤波器输出 为，，，，，，，即 ，,，,ˆ ,,,,,,,,,fxygmxydd,,,，，，，，，,,,,,, (12-30) 使得均方误差式 22ˆ，，eEfxyfxy,,min,, ，，，，,,,,，， (12-31) ˆfxy,gxy,fxy,成立，其中称为给定时的最小二乘估计值 ，，，，，， Suv,fxy,Rxy,Suv, 设为的相关函数的傅立叶变换，分别，，，，，，，，ffn nxy,Rxy,Huv,hxy,为的相关函数的傅立叶变换，为冲激响应函数 的，，，，，，，，n Suv,Suv,fxy,nxy,傅立叶变换，有时也把和分别称为和的功率谱 密，，，，，，，，fn mxy,度，则滤波器的频域表达式为 ，， 2Huv,，，1 Muv,,，，Suv,2，，Huv,，，nHuv,，，，Suv,，，f (12-32) 于是，维纳滤波复原的原理可表示为 ，，2,,Huv,，，1,,ˆFuvGuv,,, ，，，，,,Suv,2Huv,，，，， nHuv,，，，,,Suv,，，f，， (12-33) Suv,，，n 对于维纳滤波，由上式可知，当 Huv,0,时，由于存在项，所以，， Suv,，，fHuv,不会出现被 0 除的情形，同时分子中含有 Huv,项，在 Huv,0, 处，，，，，，， 1Huv,0,。当SuvSuv,,时，