1、数字图像处理学 第6章 图像复原(第三讲),6.5 中值滤波,对受到噪声污染的退化图像的复原可以采用线性滤波方法来处理，有许多情况下是很有效的。但是多数线性滤波具有低通特性，在去除噪声的同时也使图像的边缘变得模糊了。中值滤波方法在某些条件下可以作到既去除噪声又保护了图像边缘的较满意的复原。,中值滤波是一种去除噪声的非线性处理方法。它是由图基（Turky）在1971年提出的。开始，中值滤波用于时间序列分析，后来被用于图像处理，在去噪复原中得到了较好的效果。,6.5.1 中值滤波的基本原理,6.5.2 加权的中值滤波,中值滤波的基本原理是把数字图像或数字序列中一点的值用该点的一个邻域中各点值的中值
2、代替。中值的定义如下：,一组数 x1,x2,x3 xn,把个数按值的大小顺序排列于下,(6144),y称为序列 x1,x2,x3 xn 的中值。例如有一序列为（80,90,200,110,120），这个序列的中值为110。,把一个点的特定长度或形状的邻域称作窗口。在一维情形下，中值滤波器是一个含有奇数个像素的滑动窗口。窗口正中间那个像素的值用窗口内各像素值的中值代替。,其中,(6145),设输入序列为，I为自然数集合或子集，窗口长度为n。则滤波器输出为,例如，有一输入序列如下,0000002222222233333330023455555000在此序列中前面的8是脉冲噪声，中间一段是一种寄生振
3、荡，后面是希望保留的斜坡和跳变。在此采用长度为3的窗口，得到的结果为,显然，经中值滤波后，脉冲噪声8被滤除了，振荡被平滑掉了，斜坡和阶跃部分被保存了下来。,yi=0000002222222233333330023455555000,中值滤波的运算方法可以在有限程度上作些分析。例如常数 K 与序列 f(i)相乘的中值有如下关系存在,(6146),而常数 K 与序列 f(i)相加的中值有如下关系,(6147),对几种基本信号进行中值滤波的例子如图64所示。图中(a)是阶跃信号，经中值滤波后仍然保持了阶跃部分；图(b)原始信号是斜坡，滤波后也保持了其形状；图(c)的原始信号是单脉冲信号，经滤波后消去
4、了这个脉冲；,图(d)中的原始信号是双脉冲，经中值滤波后也被消去了；图(e)的原始信号是三脉冲，滤波后对其没有影响；图(f)的原始信号是三角形，滤波后虽然有少许变形，但也还基本保持了原来的形状。,图64 对几种基本信号中值滤波的结果的例子,中值滤波的概念很容易推广到二维，此时可以利用某种形式的二维窗口。设 表示数字图像各点的灰度值，滤波窗口为A的二维中值滤波可定义为,(6148),二维中值滤波的窗口可以取方形，也可以取近似圆形或十字形。,图65是二维中值滤波的实例。图中(a)是原始图像，图(b)是混有高斯白噪声的图像，(c)是33窗口中值滤波结果图像，(d)是55窗口中值滤波结果图像，(e)是
5、33窗口均值滤波结果图像，(f)是55窗口均值滤波结果图像，(g)是加有椒盐噪声的图像，(h)是33窗口中值滤波结果图像，(i)是55窗口中值滤波结果图像，(j)是33窗口均值滤波结果图像，(k)是采用55窗口均值结果图像。,二维中值滤波及均值滤波实例,(a)原像,(b)加有高斯白噪声图像,(c)中值滤波图像,(d)均值滤波图像,(e)加有椒盐噪声图像,(f)中值滤波图像,(g)均值滤波图像,6.5.1 中值滤波的基本原理,6.5.2 加权的中值滤波,以上讨论中的中值滤波，窗口内各点对输出的作用是相同的。如果希望强调中间点或距中间点最近的几个点的作用，可以采用加权中值滤波法。加权中值滤波的基本
6、原理是改变窗口中变量的个数，可以使一个以上的变量等于同一点的值，然后对扩张后的数字集求中值。,(1)一维加权的中值滤波以窗口为3的一维加权中值滤波为例，表示如下,由公式（6149）可见，在窗口内，中间点取奇数，两边点取对称数，也就是位于窗口中间的像素重复两次，位于窗口边缘的两个像素重复一次，形成新的序列，然后对新的序列在施以常规中值滤波处理。,（2）二维的加权中值滤波 二维加权中值滤波与一维情况类似。如果适当地选取窗口内各点的权重，加权中值滤波比简单中值滤波能更好地从受噪声污染的图像中恢复出阶跃边缘以及其他细节。二维加权中值滤波以33窗口为例，表示如下,原始窗口为：,加权后的中值滤波如下式所示
7、：,(6150),即中间的点取三个值（重复两次），上、下、左、右的点各取两个（重复一次），对角线上的点取一个（不重复）。,加权中值滤波与普通中值滤波有时会有不同的效果。例如，对于普通中值滤波有 y=Med(1 1 1 1 5 5 1 5 5)=1；而加权后的中值滤波为 y=weighted_Med(1 1 1 1 5 5 1 5 5)=5。加权中值滤波保持了方块角上的一点的值。,中值滤波可有效地去除脉冲型噪声，而且对图像的边缘有较好的保护。但是它也有其固有的缺陷，如果使用不当，会损失许多图像细节。例如，采用33窗口对图66(a)所示的原始图像滤波。滤波结果如图(b)所示，其结果不但削去了方块的
8、4个角，而且把中间的小方块也滤掉了。因此，中值滤波在选择窗口时要考虑其形状及等效带宽，以避免滤波处理造成的信息损失。,图66 中值滤波的实例一,图67是中值滤波的另一实例。图(a)是一条细线条图像，经33窗口滤波后，图像中的细线条完全滤掉了，如图(b)所示。,图67 中值滤波的实例二,以上两例可以直观地看到，中值滤波对图像中的细节处很不理想，所以，中值滤波对所谓的椒盐噪声(pepper salt Noise)的滤除非常有效，但是它对点、线等细节较多的图像却不太适用。,在本章给出一个中值滤波的计算机程序（见附录六），以供读者参考。,在图64中，为了比较中值滤波的效果，也给出了均值滤波的处理结果。
9、均值滤波的滤波过程也是使一个窗口在图像（或序列）上滑动，窗中心位置的值用窗内各点值的平均值来代替。以二维均值滤波为例，它的定义如下：,设 xij 表示数字图像各像素的灰度值，A为一个33的窗口，则二维均值滤波的定义为,(6151),一般均值滤波的边缘保护特性不如中值滤波。附录（七）是一个二维均值滤波的程序，以便于读者参照比较。,6.6 几种其他空间复原技术,前边讨论了几种基本的图像复原技术。除此之外，尚有一些其他的空间图像复原方法，本节将对这些方法作一些简要的讨论。,6.6.1 几何畸变校正,6.6.2 盲目图像复原,在图像的获取或显示过程中往往会产生几何失真。例如成像系统有一定的几何非线性。
10、这主要是由于视像管摄像机及阴极射线管显示器的扫描偏转系统有一定的非线性，因此会造成如图68所示的枕形失真或桶形失真。图(a)为原始图像，图(b)和图(c)为失真图像。,除此之外还有由于斜视角度获得的图像的透视失真。另外，由卫星摄取的地球表面的图像往往覆盖较大的面积，由于地球表面呈球形，这样摄取的平面图像也将会有较大的几何失真。对于这些图像必须加以校正，以免影响分析精度。,图68 几何畸变,正常图像,枕形失真,桶形失真,由成像系统引起的几何畸变的校正有两种方法。一种是预畸变法，这种方法是采用与畸变相反的非线性扫描偏转法，用来抵消预计的图像畸变；,另一种是所谓的后验校正方法。这种方法是用多项式曲线
11、在水平和垂直方向去拟合每一畸变的网线，然后求得反变化的校正函数。用这个校正函数即可校正畸变的图像。图像的空间几何畸变及其校正过程如图69所示。,图69 空间几何畸变及校正的概念,任意几何失真都可由非失真坐标系 变换到失真坐标系 的方程来定义。方程的一般形式为,(6152),(6153),在透视畸变的情况下，变换是线性的，即,设 是无失真的原始图像，是 畸变的结果，这一失真的过程是已知的并且用函数 和 定义。于是有,(6154),这说明在图像中本来应该出现在像素 上的灰度值由于失真实际上却出现在 上了。这种失真的复原问题实际上是映射变换问题。在给定了，的情况下，其复原处理可如下进行：,(1)对于
12、 中的每一点，找出在 中相应的位置。由于 和 不一定是整数，所以通常(,)不会与 中的任何点重合。,(2)找出 中与(,)最靠近的点，并且令=，也就是把 点的灰度值赋于。如此逐点作下去，直到整个图像，则几何畸变得到校正。,(3)如果不采用(2)中的灰度值的代换方法也可以采用内插法。这 种 方 法 是 假 定(,)点 找 到 后，在 中找出包围着(,)的四个邻近的数字点，并且有：,f(x,y)中点(x0,y0)的灰度值由 中四个点的灰度值间的某种内插法来确定。,(6155),在以上方法的几何校正处理中，如果(,)处在图像 之外，则不能确定其灰度值，而且校正后的图像多半不能保持其原来的矩形形状。,
