G(u,v) =F(u,v)+N(u, v)
⑶
实验名称：图像退化与复原 实验目的
1. 了解光电图像的退化原因；
2. 掌握和理解基本的噪声模型，并能对图像进行加噪处理；
3. 了解点扩展函数(PSF)与光学传递函数(OTF)的关系，熟悉几种经典的退化模 型的
模拟试验和OTF 估计方法；
4. 熟悉和掌握几种经典的图像复原方法及其基本原理；
5. 能熟练利用MATLAB 或C/C++工具进行图像的各种退化处理， 并能编程实现 退化
图像的复原。
三. 实验原理
光电成像系统出现图像退化的过程是复杂多变的，为了研究的需要，通常情 况下都把退化简化为化为一个线性移不变过程，见下图
1所示。
障质过稈 | 屯原
图1光电图像退化与复原原理图
因此，在空域中退化过程可以表示如下：
g (x,y) = f (x,y) * h(x,y) + h(x,y)
(1)
只有加性噪声不存在情况下，退化过程可以模型化如下表达式：
g(x,y) = f (x,y) + h(x,y)
(2)
其频域表达式为
:
针对这种退化图像的复原，除了周期噪声以外，通常都可以采用空间域滤波 的方法进行图像复原，此时图像复原与图像增强几乎是没有区别的。常见的空间 域滤波方法有均值滤波器和统计排序滤波器。
当退化图像存在线性移不变退化时，
图像的复原不能采用简单空间域滤波器
来实现，要实现线性移不变退化图像的复原， 必须知道退化系统的退化函数，即 点扩展函数h(x,y)。在点扩展函数已知的情况下，常见图像复原方法有逆滤波 和维纳滤波两种。
在考虑噪声的情况下，逆滤波的原理可以表示如下：
通常情况下，N (u,v)是未知的，因此即使知道退化模型也不能复原图像 此外，当H (u,v )的任何元素为零或者值很小时，N (u,v )/H (u,v )的比值决定 着复原的结果，从而导致图像复原结果出现畸变。对于这种情况， 通常采用限制
滤波频率使其难以接近原点值，从而减少遇到零值的可能性。
维纳滤波则克服了逆滤波的缺点，其数学模型表示如下：
然而，为退化图像的功率谱很少是已知的，因此常常用下面表达式近似:
因此，本实验的内容就是利用上述经典图像复原的原理，对降质退化图像进 行复原。
四. 实验步骤
本次实验主要包括光电图像的退化模型和复原方法实现两大部分内容。
(一)图像的退化图像 1、大气湍流的建模
°
F(u,v) = G(u,v) U
F(u,v) =
G(u,v) H(u,v) F(u,v) +
N(u,v)
H(u,v)
°
犏 F (u,v)=犏 J _________
(u,v) H (u,v) H *(u,v)2
+ S h (u,v)/S f (u,v)
G(u,v)
1) 湍流引起图像退化的光学传递函数(OTF 生成。已知湍流退化模型的 OTF 表达式如下：
H (u,v) = exp[-k(u 2+ v 2)5/6 ]
(7)
其中，k 为一个常数，反映了大气湍流的严重程度。
(u,v)分别代表了(x,y )
方向的频率坐标。为了生成中心化的 OTF 可以考虑将式(7)改写为：
H (u,v)= exp[-k ((u- M /2)2 + (v- N/2)2)5"]
(8)
其中，M ，N 为图像的长和宽。
2) 读入一幅灰度图像，设定式(8)中k= 0.0025,进行退化试验。分别显示 原始图
像、退化模型和退化图像。
3) 设定k= 0.001、0.00025重复上一步的试验。
2、运动模糊的图像退化试验
1)
匀速直线运动引起图像退化的光学传递函数(
OTF )生成。已知相机匀
速直线运动的OTF 表达式如下:
对原图形进行灰度处理
由FP 与H 进行相关处
理,
读入原始图像
V
显示原图像
将上述结果进行 fft 处理
得到FP
设计退化湍流模型为 H 显示传递函数
显示退化图像
结束
图2大气湍流的退化过程
结束
图3运动模糊的图像退化
（二）图像复原试验 1、逆滤波
1） 根据试验（一）设计一幅退化图像（包括噪声污染+模糊退化两部分）， 其中模
糊退化可选高斯模糊、大气湍流模糊或运动模糊（
方向可任意指定，如
10度、20度、45度等），噪声模型可自行设定。
H (u,v) =
sin[p(ua+ vb)]e - jp(ua+vb)
p(ua + vb)
(9)
其中，T 为相机曝光时间，a ,b 分别表示（x,y ）方向的速度；（u,v ）分别对
应（x, y ）方向的频率坐标。
2）读入一幅灰度图像，设定式（9）中T = 1.0, a=b=0.1,编写MATLAB 代
码进行模糊退化试验。要求分别显示原始图像、退化模型和退化图像。
3）设定不同的值，a ，b 值，重复上一步的试验。
4）利用数字显微镜或其它图像采集设备，移动物体过程中，采集图像。
读入原始图像
对原图形进行灰度处理
显示原图像
将上述结果进行fft 处理
得到FP
设计运动模糊模型为 H 显示传递函数
由FP 与H 进行相关处
理,
显示退化图像
