第五章_图像复原

ˆ ( x, y) min {g (s, t )} f
( s ,t )S xy


寻找图像中最暗点 消除“盐”噪声
32
5.3.2 统计排序滤波器
最大值滤波器可以去除”胡椒”噪声,但 会从黑色物体边缘移除一些黑色像素. 最小值滤波器可以去除”盐”噪声,但会 从亮色物体边缘移除一些白色像素.
图5.8(a)
5
图像恢复准则：要用某一客观标准来度量，则为某种准则下的 最优估计。 ◘图像恢复与图像增强的异同点 相同点：图像增强与图像恢复都是改善给定图像的质量。 不同点： （1）图像恢复是利用退化过程的先验知识，来建立图像的退 化模型，再采用与退化相反的过程来恢复图像，而图像增强一 般无需对图像降质过程建立模型。 （2）图像恢复是针对图像整体，以改善图像的整体质量。而 图像增强是针对图像的局部，以改善图像的局部特性，如图像 的平滑和锐化。
5.3 仅有噪声存在时的空域滤波复原

当一幅图像中唯一存在的退化是噪声 时,(5.1.1)式和(5.1.2)式变成: g(x,y) = f(x,y) +η(x,y) (5.3-1) 和 G(u,v) = F(u,v) + N(u,v) (5.3-2)
当只有加性噪声存在时，使用空域滤波法。 此时图像增强和复原基本上是等同的。除 了介绍一些新的滤波器外，执行所有滤波 器的机制与3.5节讨论的完全一样。
若b a, 灰度值b将显示为一个亮点, a的值将显示为一个暗点. 若Pa或Pb为零, 则脉冲噪声称为单极脉冲. 若Pa或Pb均不为零, 尤其是近似相等时, 脉冲噪声值类似于随机 分布在图像上的胡椒和盐粉细粒.
16
5.2 噪声模型
例5.1：样本噪声图 像和它们的直方图
17
高斯
瑞利
伽马
指数
均匀
椒盐

g ( s, t )Q1 g ( s, t ) Q

其中Q为滤波器的阶数。该滤波器适合减少甚至消 除椒盐噪声。但注意它不能同时消除“椒”和“盐” 噪 当Q值为正数时,滤波器用于消除"椒"噪声； 声。 当Q值为负数时,滤波器用于消除"盐"噪声；
当Q=0时,反调和均值滤波器退化为算术均值滤波器； 当Q=-1时,反调和均值滤波器退化为谐波均值滤波器.
图像复原
代数恢复方法
●无约束最小二乘法 ●有约束最小二乘法 －能量约束 －平均约束 －均方误差最小 （维纳滤波） －最大熵约束
10
频域恢复
●逆滤波 ●维纳滤波 ●功率谱均衡 ●几何平均
其它
● 卡尔曼滤波 ● 运动模糊恢复 ● 几何畸变消除
5.1 图像退化/复原过程的模型


退化过程：退化函数 H 和加性噪声η(x, y)联合作 用于原始图像 f(x, y)产生退化图像 g(x, y)。 复原过程：给定退化图像g(x, y) ，以及关于退化 函数 H 和加性噪声η(x, y)的一些知识，获得原始 图像f(x, y)的一个估计 fˆ ( x, y ) 。
这两个公式是本章大部分内容的基础。
12
5.2 噪声模型

数字图像的噪声主要来源于图像的获取(数 字化过程)和传输过程。
噪声的空域和频域特性： 空域：假设噪声值与空间坐标和像素值无关 频域：考察噪声的频谱，如白噪声的频谱是 一个常量。
13
5.2 噪声模型
高斯噪声
高斯随机变量z的PDF为：
1 ( z u ) 2 / 2 2 p( z ) e 2
2
图像复原与图像增强的区别



图像复原方法通常都包含有一个评判准则， 以产生所期望结果的最佳估计。 对比而言，图像增强是一个启发式过程，通 常利用人类视觉系统的生理特点去操作图像。 例如：


对比度拉伸是一种增强技术 用去模糊函数去除图像模糊是一种图像恢复技术
3
4
◘基本概念
图像降质：宇航、卫星、航空测绘、遥感、天文学中所得照
25
5.3.1 均值滤波器
结论： 算术和几何均值滤波器(尤其后者)适合处 理高斯噪声。 反谐波均值滤波器适合脉冲噪声，但缺点 是必须知道噪声是暗噪声还是亮噪声，以 便正确选择Q的符号，否则后果严重。且 其不能同时消除“椒”和“盐”噪声。
26
5.3.1 均值滤波器
习题5.1-5.5：下图中的白条为7个像素宽210个像素高，白条 之间的间隔为17个像素。分别用算术均值滤波器、几何均值 滤波器、反谐波均值滤波器(Q = 1和Q = -1两种情况)对该图 进行处理，结果会是怎样？
23
5.3.1 均值滤波器
(a) 以0.1的概率被”胡椒” 噪声污染的图像 (b) 以0.1的概率被”盐” 噪声污染的图像 (c) 用3×3大小、阶数为 1.5的反谐波滤波器滤波 的结果 (d) 用Q=-1.5滤波(b)的结果
反谐波均值滤 波器对处理椒 盐噪声的效果
24
5.3.1 均值滤波器
在反谐波滤波中错误地选择符号的结果 (a) 原图像 (b) 用3×3 的大小和Q＝－1.5的反谐波滤波器滤波的结果 (c) 用Q=1.5滤波的结果
8
图像退化



在宇航、遥感、天文和侦察照片中，退化可 能是由大气扰动、光学系统的象差、相机和 对象之间的相对运动等引起的； 电子显微图片中的退化，常常是由电子透镜 的球面象差引起的； 医学射线图片的退化，是由射线图像系统本 身特性所导致的低分辨率和低灰对比度引起 的。
9
第五章 图像复原
图像复原所包括的主要内容如图所示
退化函数 H 复原滤波
退化 图5.1 图像退化/复原过程的模型
复原
11
5.1 图像退化/复原过程的模型


如果H是一个线性、与位臵无关的过程，那么退化 图像在空间域可表示为： g(x, y) = h(x, y) * f(x, y) +η(x, y) (5.1-1) 其中h(x, y)是退化函数的空域表示，*表示卷积。 等价的频域表示为： G(u, v) = H(u, v)F(u, v) + N(u, v) (5.1-2) 大写字母表示的是(5.1-1)中对应项的傅里叶变换。
1 p( z ) b a 0
a zb 其他
概率密度的均值和方差由下式给定 : ab 2 2 2 (b a ) = 12

Hale Waihona Puke Baidu
15
5.2 噪声模型
脉冲(椒盐)噪声
(双极)脉冲噪声的PDF为：
Pa p ( z ) Pb 0
za z b 其他
1 mn d
( s ,t )S xy
g ( s, t )
r
d的取值范围：[0, mn - 1] 当d = 0时，退变为算术均值滤波器 当d = mn - 1时，退变为中值滤波器 作用：适合高斯和椒盐混合噪声。 Why?
35
5.3.2 统计排序滤波器
(a) 由加性均匀噪声污染的图像 均值为0,方差为800的均匀噪声 (b) 图(a)加上椒盐噪声污染的图像 Pa=Pb=0.1的椒盐噪声 (c) 对(b)算术均值滤波处理结果 (d) 对(b)几何均值滤波器处理结果 (e) 对(b)中值滤波器处理结果 (f) 对(b) 进行d=5的阿尔法修正均值滤 波处理结果 对高噪声级别图像，宜采用较大的滤波器尺寸 由于脉冲噪声的存在,算术均值 滤波器和几何均值滤波器没有 起到良好作用. 中值滤波器和阿尔法修正滤波 器效果更好,阿尔法修正最好. 当d较大时，阿尔法修正滤波器 性能与中值滤波器相近，但保持 有一定平滑能力
图5.8(b) (a)用大小为3×3的最 大值滤波器对图5.8(a) 滤波的结果 (b)用同样大小的最小 值滤波器对图5.8(b)滤 波的结果
33
5.3.2 统计排序滤波器

中点滤波器 当前像素位臵的新灰度值为邻域内像素最 大值和最小值的平均值：
1 ˆ f ( x, y ) max {g ( s, t )} min {g ( s, t )} ( s ,t )S xy 2 ( s ,t )S xy
“中点”含义：邻域像素灰度范围的中点 作用：最适合均值为0的高斯噪声和均匀分布噪声
34
5.3.2 统计排序滤波器

阿尔法修剪均值滤波器
从邻域像素值序列中删去最高的d / 2个像素值和 最低的d / 2个像素值，然后对剩下的mn - d个像 素求均值，得到当前像素的新像素值：
ˆ ( x, y ) f
数字图像处理
Digital Image Processing
第五章 图像复原
图像复原与图像增强的区别



和图像增强一样，图像复原的最终目标是 改善图像的质量。尽管二者有交叉，但图 像增强很大程度上是一个主观的过程，而 图像复原多半是客观的过程。 图像复原试图利用退化现象的先验知识来 重建或恢复图像。 图像复原一般先对退化过程建模，然后应 用其相反的过程来恢复出原始图像。
( s ,t )S xy
尤其适合于脉冲噪声(即冲击噪声或椒盐噪 声)的处理(无论单极或双极)
29
5.3.2 统计排序滤波器
对噪声图像多次应用中值滤波器 (a)由概率Pa=Pb=0.1的椒盐 噪声污染的图像 (b) 用尺寸为3×3的中值滤波 器处理的结果 (c) 用该滤波器处理(b)的结果 (d) 用相同的滤波器处理(c)的 结果 经过多次处理,逐渐消除 噪声；但多次应用中值滤 波器,会使图像模糊
18

5.3.1 均值滤波器

算术均值滤波器
最简单的均值滤波器。令Sxy表示中心在点(x,y)、窗 口尺寸为m×n的矩形子图坐标集合，g(x,y)为污染 图像。则复原图像 fˆ 在点(x,y)处的值为区域Sxy内像 素的算术平均值：
1 ˆ f ( x, y ) g ( s, t ) mn ( s ,t )S xy
片，由于大气湍流，光学系统的相差及摄像机与物体之间的相
对运动等，会使图像降质。 退化图像：由于各种原因，使得原清晰图像变模糊，或者原 图像没有达到应有的质量而形成的降质图像。 图像恢复（复原）： 使退化图像恢复本来面目。
图像恢复过程及其关键：根据图像降质过程的某些先验知识，
建立“退化（降质）模型”，运用和退化相反的过程，将退化 图像恢复。
其中z表示灰度值, 表示z的平均值或期望值, 表示z的标准差. 标准差的平方 2称为z的方差. 当z服从高斯分布时, 其值70%落在[( ), ( )]范围内, 且有95%落在[( 2 ), ( 2 )]范围内.
14
5.2 噪声模型
均匀分布噪声
均匀分布噪声的PDF为：
27
5.3.2 统计排序滤波器
回顾：什么是统计排序滤波器？
本节介绍四类统计排序滤波器： 中值滤波器 最大和最小值滤波器 中点滤波器 阿尔法修剪均值滤波器
28
5.3.2 统计排序滤波器

中值滤波器 当前像素位臵的新灰度值为邻域中像素的 灰度中值：
ˆ ( x, y) median{g (s, t )} f
ˆ ( x, y ) f
mn
( s ,t )S xy

1 g ( s, t )
谐和均值滤波器对“盐”噪声效果较好，但不适合 处 理“椒”噪声。对其他噪声如高斯噪声也有较好的 效 果。
22
5.3.1 均值滤波器

反谐波均值滤波器
ˆ ( x, y) ( s ,t )Sxy f
( s ,t )S xy
这一操作可用所有系数为1/mn的卷积掩膜来实现。 如前所述，该滤波器是以模糊为代价来达到降噪的 目的。
19
5.3.1 均值滤波器

几何均值滤波器
ˆ ( x, y ) f g ( s, t ) ( s ,t )S xy
1 mn
几何均值滤波器的平滑效果与算术均值滤波器相 当，但损失的图像细节较少。 缺点：计算复杂度高
30
5.3.2 统计排序滤波器

最大值滤波器 当前像素位臵的新灰度值为邻域中像素的 灰度最大值：
ˆ ( x, y) max {g (s, t )} f
( s ,t )S xy


寻找图像中最亮点 消除“胡椒”噪声
31
5.3.2 统计排序滤波器

最小值滤波器 当前像素位臵的新灰度值为邻域中像素的 灰度最小值：
6
（3）图像恢复主要是利用图像退化过程来恢复图 像的本来面目，它是一个客观过程，最终的结果 必须要有一个客观的评价准则。而图像增强主要
是用各种技术来改善图像的视觉效果，以适应人
的心理、生理需要，而不考虑处理后图像是否与
原图像相符，也就很少涉及统一的客观评价准则。
7
图像退化



在一个图像系统中存在着许多退化源。一些 退化因素只影响一幅图像中某些个别点的灰 度，称为点退化。 而另外一些退化因素则可以使一幅图像的一 个空间区域变得模糊，称为空间退化。 此外还有数字化器、显示器、时间、彩色、 以及化学作用引起的退化。
20
5.3.1 均值滤波器
(a) 电路板的X射线图像 (b) 由加性高斯噪声污染 的图像 (c) 用3×3算术均值滤波器 滤波的结果 (d) 用3×3的几何均值滤波 器滤波的结果
算术均值和几何 均值都能衰减噪 声,但比较而言, 几何均值滤波器 模糊程度较弱.
21
5.3.1 均值滤波器

谐波均值滤波器