图像退化与复原培训讲学.docx

《图像退化与复原培训讲学.docx》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、图像退化与复原精品文档一．实验名称：图像退化与复原二．实验目的1．了解光电图像的退化原因；2．掌握和理解基本的噪声模型，并能对图像进行加噪处理；3．了解点扩展函数(PSF)与光学传递函数(OTF)的关系，熟悉几种经典的退化模型的模拟试验和OTF估计方法；4．熟悉和掌握几种经典的图像复原方法及其基本原理；5．能熟练利用MATLAB或C/C++工具进行图像的各种退化处理，并能编程实现退化图像的复原。三．实验原理图1光电图像退化与复原原理图光电成像系统出现图像退化的过程是复杂多变的，为了研究的需要，通常情况下都把退化简化为化为一个线性移不变过程，见下图1

2、所示。因此，在空域中退化过程可以表示如下：(1)收集于网络，如有侵权请联系管理员删除精品文档只有加性噪声不存在情况下，退化过程可以模型化如下表达式：(2)其频域表达式为：(3)针对这种退化图像的复原，除了周期噪声以外，通常都可以采用空间域滤波的方法进行图像复原，此时图像复原与图像增强几乎是没有区别的。常见的空间域滤波方法有均值滤波器和统计排序滤波器。当退化图像存在线性移不变退化时，图像的复原不能采用简单空间域滤波器来实现，要实现线性移不变退化图像的复原，必须知道退化系统的退化函数，即点扩展函数。在点扩展函数已知的情况下，常见图像复原方法有逆滤波和维

3、纳滤波两种。在考虑噪声的情况下，逆滤波的原理可以表示如下：(4)通常情况下，是未知的，因此即使知道退化模型也不能复原图像。此外，当的任何元素为零或者值很小时，的比值决定着复原的结果，从而导致图像复原结果出现畸变。对于这种情况，通常采用限制滤波频率使其难以接近原点值，从而减少遇到零值的可能性。维纳滤波则克服了逆滤波的缺点，其数学模型表示如下：(5)然而，为退化图像的功率谱很少是已知的，因此常常用下面表达式近似：收集于网络，如有侵权请联系管理员删除精品文档(6)因此，本实验的内容就是利用上述经典图像复原的原理，对降质退化图像进行复原。一．实验步骤本次实

4、验主要包括光电图像的退化模型和复原方法实现两大部分内容。(一)图像的退化图像1、大气湍流的建模1)湍流引起图像退化的光学传递函数（OTF）生成。已知湍流退化模型的OTF表达式如下：(7)其中，为一个常数，反映了大气湍流的严重程度。分别代表了方向的频率坐标。为了生成中心化的OTF，可以考虑将式(7)改写为：(8)其中，，为图像的长和宽。2)读入一幅灰度图像，设定式(8)中，进行退化试验。分别显示原始图像、退化模型和退化图像。3)设定重复上一步的试验。收集于网络，如有侵权请联系管理员删除精品文档图2大气湍流的退化过程1、运动模糊的图像退化试验1)匀速直

5、线运动引起图像退化的光学传递函数（OTF）生成。已知相机匀速直线运动的OTF表达式如下：(9)其中，为相机曝光时间，,分别表示方向的速度；分别对应方向的频率坐标。2)读入一幅灰度图像，设定式（9）中T=1.0，a=b=0.1，编写MATLAB代码进行模糊退化试验。要求分别显示原始图像、退化模型和退化图像。3)设定不同的值，，值，重复上一步的试验。4)利用数字显微镜或其它图像采集设备，移动物体过程中，采集图像。收集于网络，如有侵权请联系管理员删除精品文档图3运动模糊的图像退化(一)图像复原试验1、逆滤波1)根据试验（一）设计一幅退化图像（包括噪声污染

6、+模糊退化两部分），其中模糊退化可选高斯模糊、大气湍流模糊或运动模糊（方向可任意指定，如10度、20度、45度等），噪声模型可自行设定。2)利用MATLAB编程实现利用全逆滤波方法对退化图像的复原。要求在同一个窗口下显示原始退化图像、复原结果及复原结果与理想图像的差值图共3个图，并对复原结果进行必要的分析。逆滤波复原公式如下：(10)收集于网络，如有侵权请联系管理员删除精品文档其中，为退化图像的傅立叶变换，为退化系统的光学传递函数（OTF）。1)伪逆滤波：为了防止逆滤波中过小，使得复原后的图像数据过大和放大噪声，可采用频谱半径（阈值）限制下的逆滤波

7、方法，即(11)其中，为中心化频谱中某点到原点（零频）的距离或半径。另一种替代方法是直接限制的值，即(12)其中，σ为一个阈值，用于限制频谱的幅度值。这种方法被称为伪逆滤波。实验要求利用式(11)方式的伪逆滤波重复实验步骤内容2）所涉及的图像。1、Wiener滤波1)针对以上逆滤波设计的退化图，编程实现利用Wiener滤波对其进行复原。滤波原理如下：(13)其中，为退化图像的傅立叶变换，为退化系统的光学传递函数（OTF），为一个与信噪比有关的调节因子。要求在同一个窗口下显示理想图像（退化前）、退化图像、复原结果等共3个图，并对复原结果进行必要的分析

8、。2)改变值，重复试验内容1）。收集于网络，如有侵权请联系管理员删除精品文档以上应根据原理自行编写代码，不允许直接调用MA