《数字图象处理-第6章图像复原.ppt》由会员分享,可在线阅读,更多相关《数字图象处理-第6章图像复原.ppt(47页珍藏版)》请在金锄头文库上搜索。
1、 第第6 6章章 图图 像像 复复 原原Hf(x,y)g(x,y)n(x,y)图图61 图像退化模型图像退化模型输入图像输入图像系统系统加性噪声加性噪声退化图像退化图像6.1 6.1 退化模型退化模型 系统的分类方法很多系统的分类方法很多, ,例如例如: : 线性系统和非线性系统;线性系统和非线性系统; 时变系统和非时变系统;时变系统和非时变系统; 集中参数系统和分布参数系统;集中参数系统和分布参数系统; 连续系统和离散系统等等。连续系统和离散系统等等。 线性系统就是具有均匀性和相加性的系统。对于图线性系统就是具有均匀性和相加性的系统。对于图6161所所示的系统来说,可表示成下式示的系统来说,
2、可表示成下式 (61)(61) 6.1.1 6.1.1 系统系统 的基本定义的基本定义 (62)(62) 如果暂不考虑加性噪声如果暂不考虑加性噪声 n ( (x, y) ) 的影响,而令的影响,而令 n (x, y)=0 时,则时,则 如果输入信号为如果输入信号为 , ,对应的输出信号为,对应的输出信号为 , ,通过系统后有下式成立,通过系统后有下式成立 (63)(63) 如果一个系统的参数不随时间变化,即称为时不变系统或非时变如果一个系统的参数不随时间变化,即称为时不变系统或非时变系统。否则,就称该系统为时变系统。与此概念相对应,对于二维函系统。否则,就称该系统为时变系统。与此概念相对应,对
3、于二维函数来说,如果数来说,如果 (65)(65) 在图像复原处理中,往往用线性和空间不变性的系统模型加以近似。在图像复原处理中,往往用线性和空间不变性的系统模型加以近似。这种近似的优点是使线性系统理论中的许多理论可直接用来解决图像这种近似的优点是使线性系统理论中的许多理论可直接用来解决图像复原问题,所以图像复原处理特别是数字图像复原处理主要采用线性复原问题,所以图像复原处理特别是数字图像复原处理主要采用线性的,空间不变的复原技术的,空间不变的复原技术。 在线性系统理论中,曾定义了单位冲激信号在线性系统理论中,曾定义了单位冲激信号 。它是一个振幅。它是一个振幅在原点之外所有时刻为零,在原点处振
4、幅为无限大、宽度无限小,在原点之外所有时刻为零,在原点处振幅为无限大、宽度无限小,面积为的窄脉冲。其时域表达式为面积为的窄脉冲。其时域表达式为 (66)(66) 6.1.2 6.1.2 连续函数退化模型连续函数退化模型 根据根据 的关系的关系, , 如果令如果令 , , 则有下式成立则有下式成立 由于由于 是线性算子,所以是线性算子,所以 (612)(612) 令令 则则 (613)(613) 如果如果 f(x) ,h(x) 都是具有周期为都是具有周期为 N 的序列,那么,它们的时的序列,那么,它们的时域离散卷积可定义为下式之形式。域离散卷积可定义为下式之形式。 (617)(617) 显然,显
5、然, 也是具有周期也是具有周期 N 的序列。周期卷积可用常规卷的序列。周期卷积可用常规卷 积法计算也可用卷积定理进行快速卷积计算。积法计算也可用卷积定理进行快速卷积计算。 6.1.3 6.1.3 离散的退化模型离散的退化模型 逆滤波复原法也叫做反向滤波法。基本原理如下:逆滤波复原法也叫做反向滤波法。基本原理如下: 如果退化图像为如果退化图像为 g( (x,y) ),原始图像为原始图像为 f( (x,y) ),在不考虑噪声的情况下,在不考虑噪声的情况下,其退化模型用下式表示其退化模型用下式表示 (645)(645) 这显然是一卷积表达式。由傅立叶变换的卷积定理可知有下式成立这显然是一卷积表达式。
6、由傅立叶变换的卷积定理可知有下式成立 6.2 6.2 逆滤波逆滤波 6.2.1 6.2.1 逆滤波的基本原理逆滤波的基本原理 式中式中G( (u,v),),H( (u,v),),F( (u,v) ) 分别是退化图像分别是退化图像 g( (x,y) ),点扩散函数点扩散函数 h( (x,y) ),原始图像原始图像 f( (x,y) ) 的傅立叶变换。的傅立叶变换。 (6(646)46) 由式由式(646)(646),可得,可得 (647)(647) (648)(648) 图图 62 62 实际的逆滤波处理框图实际的逆滤波处理框图 图图6262的模型包括了退化和恢复运算。退化和恢复总的传递函数可的
7、模型包括了退化和恢复运算。退化和恢复总的传递函数可用用 来表示。此时有来表示。此时有 (651) (651) 式中式中 是是 的估计值,的估计值, 是是 的傅立的傅立叶变换。叶变换。 叫做输入传递函数,叫做输入传递函数, 叫做处理传递函数,叫做处理传递函数, 叫做输出传递函数。叫做输出传递函数。 6.3 6.3 中值滤波中值滤波 对受到噪声污染的退化图像的复原可以采用线性滤波方法来处理,有许多对受到噪声污染的退化图像的复原可以采用线性滤波方法来处理,有许多情况下是很有效的。但是多数线性滤波具有低通特性,在去除噪声的同时情况下是很有效的。但是多数线性滤波具有低通特性,在去除噪声的同时也使图像的边
8、缘变得模糊了。中值滤波方法在某些条件下可以作到既去除也使图像的边缘变得模糊了。中值滤波方法在某些条件下可以作到既去除噪声又保护了图像边缘的较满意的复原。噪声又保护了图像边缘的较满意的复原。 中值滤波的基本原理是把数字图像或数字序列中一点的值用该点中值滤波的基本原理是把数字图像或数字序列中一点的值用该点的一个邻域中各点值的中值代替。中值的定义如下:的一个邻域中各点值的中值代替。中值的定义如下: 一组数一组数 x x1 1, x, x2 2, x, x3 3 x xn n , , 把个数按值的大小顺序排列于下把个数按值的大小顺序排列于下 6.5.1 6.5.1 中值滤波的基本原理中值滤波的基本原理
9、 (6144)(6144) 例如,有一输入序列如下例如,有一输入序列如下 在此采用长度为在此采用长度为3 3的窗口,得到的结果为的窗口,得到的结果为: : 显然,经中值滤波后,脉冲噪声显然,经中值滤波后,脉冲噪声8 8被滤除了,振荡被平滑掉了,斜坡被滤除了,振荡被平滑掉了,斜坡和阶跃部分被保存了下来。和阶跃部分被保存了下来。 图图64 64 对几种基本信号中值滤波的结果的例子对几种基本信号中值滤波的结果的例子 二维中值滤波及均值滤波实例二维中值滤波及均值滤波实例(a) 原像(b)加有高斯白噪声图像(c)中值滤波图像(d) 均值滤波图像(e) 加有椒盐噪声图像(f) 中值滤波图像(g) 均值滤波
10、图像以上讨论中的中值滤波,窗口内各点对输出的作用是相同的。如果希以上讨论中的中值滤波,窗口内各点对输出的作用是相同的。如果希望强调中间点或距中间点最近的几个点的作用,可以采用加权中值滤望强调中间点或距中间点最近的几个点的作用,可以采用加权中值滤波法。加权中值滤波的基本原理是改变窗口中变量的个数,可以使一波法。加权中值滤波的基本原理是改变窗口中变量的个数,可以使一个以上的变量等于同一点的值,然后对扩张后的数字集求中值。个以上的变量等于同一点的值,然后对扩张后的数字集求中值。6.5.2 6.5.2 加权的中值滤波加权的中值滤波 (1) (1) 一维加权的中值滤波一维加权的中值滤波以窗口为以窗口为3
11、 3的一维加权中值滤波为例,表示如下的一维加权中值滤波为例,表示如下 (2 2)二维的加权中值滤波)二维的加权中值滤波 二维加权中值滤波以二维加权中值滤波以3333窗口为例,表示如下窗口为例,表示如下 原始窗口为:原始窗口为: 加权后的中值滤波如下式所示:加权后的中值滤波如下式所示: (6150)(6150) 即中间的点取三个值(重复两次),上、下、左、右的点各取两个即中间的点取三个值(重复两次),上、下、左、右的点各取两个(重复一次),对角线上的点取一个(不重复)。(重复一次),对角线上的点取一个(不重复)。 图图66 66 中值滤波的实例一中值滤波的实例一 图图6767是中值滤波的另一实例
12、。图是中值滤波的另一实例。图(a)(a)是一条细线条图像,经是一条细线条图像,经3333窗口滤波后,图像中的细线条完全滤掉了,如图窗口滤波后,图像中的细线条完全滤掉了,如图(b)(b)所示。所示。图图67 67 中值滤波的实例二中值滤波的实例二 6.6.1 6.6.1 几何畸变校正几何畸变校正 6.6 6.6 几种其他空间复原技术几种其他空间复原技术 图图68 68 几何畸变几何畸变 正常图像正常图像枕形失真枕形失真桶形失真桶形失真l一种是一种是预畸变法预畸变法,这种方法是采用与畸变相反的非线性扫描,这种方法是采用与畸变相反的非线性扫描偏转法,用来抵消预计的图像畸变;偏转法,用来抵消预计的图像
13、畸变;l另一种是所谓的另一种是所谓的后验校正方法后验校正方法。这种方法是用多项式曲线在。这种方法是用多项式曲线在水平和垂直方向去拟合每一畸变的网线,然后求得反变化的校水平和垂直方向去拟合每一畸变的网线,然后求得反变化的校正函数。用这个校正函数即可校正畸变的图像。正函数。用这个校正函数即可校正畸变的图像。图图69 69 空间几何畸变及校正的概念空间几何畸变及校正的概念 多数的图像复原技术都是以图像退化的某种先验知识为基础,也就多数的图像复原技术都是以图像退化的某种先验知识为基础,也就是假定系统的脉冲响应是已知的。但是,在许多情况下难以确定退化是假定系统的脉冲响应是已知的。但是,在许多情况下难以确
14、定退化的点扩散函数。在这种情况下,必须从观察图像中以某种方式抽出退的点扩散函数。在这种情况下,必须从观察图像中以某种方式抽出退化信息,从而找出图像复原方法。化信息,从而找出图像复原方法。6.6.2 6.6.2 盲目图像复原盲目图像复原对具有加法性噪声的模糊图像作盲目图像复原的方法有两种,就是对具有加法性噪声的模糊图像作盲目图像复原的方法有两种,就是直直接测量法接测量法和和间接估计法间接估计法。直接测量法盲目图像复原通常要测量图像的模糊脉冲响应和噪声功率直接测量法盲目图像复原通常要测量图像的模糊脉冲响应和噪声功率谱或协方差函数。谱或协方差函数。图像复原的新方法图像复原的新方法 (Partial differential equation) (PDE) (Partial differential equation) (PDE) 偏微分方程图像修复方法偏微分方程图像修复方法