教学课件第4章 图像增强（第4－2讲）（研究生学位课）.ppt
教学课件第4章 图像增强（第42讲）（研究生学位课）数字图像处理学数字图像处理学第第4章章 图像增强图像增强（第二讲）（第二讲）阮秋琦教授阮秋琦教授4.1.1 直方图 4.1.2 直方图修改技术的基础 4.1.3 直方图均衡化处理 4.1.4 直方图规定化处理 4.1.5 4.1.5 图像对比度处理图像对比度处理 4.1.5 图图像像对对比度比度处处理理 由于图像的亮度范围不足或非线性会使图像的对比度不甚理想，可用像素幅值重新分配的方法来改善图像对比度。扩大图像的亮度范围可以用线性映射的方法，这种方法如图310所示。由图可以看出原图像的范围较小，经映射后的图像亮度范围展宽了。图4-10 数字图像对比度增强 在这种转换中，设计转换函数应考虑到灰度量化问题，如果原始图像的灰度级为 k 级，映射后输出图像的灰度级仍然是 k 级，这样由于输出图像的灰度范围加大了。在对比度处理法中，根据不同的目的可以设计出不同的转换函数。图4-11是线性转换函数，这种函数将图像在整个灰度范围内作线性映射。图4-11 图像灰度的线性映射变换 另外一种映射转换函数如图4-12所示。这种转换是将图像中两个极端的灰度值加以限幅，这种限幅的比例也是可以选择的。图4-12 限幅的线性映射变换 除此之外，为了不同的目的还有其他一些类型的转换函数。这些转换函数的形式如图4-13(a)、(b)、（c）所示。图4-13 其他一些转换函数 灰度变换的效果如图4-14(a)(b)所示，其中（a）是原像，(b)是处理后的图像。图 4-14 灰度变换处理效果 灰度反转的转换函数是把图像的低亮度区域转到较高的亮度区，而高亮度区转换为低亮度区，其效果如图4-15所示，其中(a)是原像，(b)是处理后的图像。图 4-15 灰度反转处理效果 锯齿形转换可以把几段较窄的输入灰度区间都扩展到整个输出灰度范围内，这种处理可以把灰度变化较平缓的区域也较鲜明地显示出来。其效果如图4-16所示，其中（a）是原像，（b）是处理后的图像，这里选 n=2 。图 4-16 锯尺状变换函数处理效果（n=2）开窗式转换的目的是只对部分输入灰度区间进行转换，通过窗口位置的选择可以观察某些灰度区间的灰度分布，并且对这一区域的灰度进行映射变换。当然，图313只是举出几种常用的转换函数的形状。根据不同的需要还可以设计出更多的转换函数，其基本原理都是一样的，只不过处理效果不同罢了。经开窗式转换函数处理的图像效果如图4-17所示，图(a)是原像(b)是处理后的图像。图 4-17 经开窗变换函数处理的效果 利用直方图修正技术增强图像简便而有效。直方图均衡化处理可大大改善图像灰度的动态范围，利用直方图规定化方法能得到更加符合需要的结果，通过对比度转换函数的正确设计可以方便灵活地改善图像。因此，这种方法在数字图像处理中得到广泛应用。(a)为原始图像(b)为均衡化处 理后的图像(c)为原始图像(d)为均衡化处 理后的图像图4-184.2 4.2 图像平滑化处理图像平滑化处理 一幅图像可能存在着各种寄生效应。这些寄生效应可能在传输中产生，也可能在量化等处理过程中产生。一个较好的平滑方法应该是既能消掉这些寄生效应又不使图像的边缘轮廓和线条变模糊。这就是研究图像平滑化处理要追求的主要目标。图像平滑化处理方法有空域法和频域法两大类。主要有:邻域平均法，低通滤波法，多图像平均法等等。4.2.1 邻邻域平均法域平均法 邻域平均法是简单的空域处理方法。这种方法的基本思想是用几个像素灰度的平均值来代替每个像素的灰度。假定有一幅NN个像素的图像 f(x,y)，平滑处理后得到一幅图像 g(x,y)。g(x,y)由下式决定 式中，S 是(x,y)点邻域中点的坐标的集合，但其中不包括(x,y)点，是集合内坐标点的总数。(4-(4-2 21)1)式(4-21)说明，平滑化的图像 g(x,y)中的每个像素的灰度值均由包含在(x,y)的预定邻域中的f(x,y)的几个像素的灰度值的平均值来决定。例如，可以以点(x,y)为中心，取单位距离构成一个邻域，其中点的坐标集合为：图4-19(a)的方法是一个点的邻域，定义为以该点为中心的一个圆的内部或边界上的点的集合。图中像素间的距离为x ，选取x 为半径作圆，那么，点 R 的灰度值就是圆周上四个像素灰度值的平均值。图4-19 在数字图像中选取邻域的方法 四邻域：八邻域：处理结果表明，上述选择邻域的方法对抑制噪声是有效的，但是随着邻域的加大，图像的模糊程度也愈加严重。为克服这一缺点，可以采用阈值法减少由于邻域平均所产生的模糊效应。其基本方法由下式决定：（4-224-22）式中 T 就是规定的非负的阈值。这个表达式的物理概念是：当一些点和它的邻域内的点的灰度的平均值的差不超过规定的阈值 T 时，就仍然保留其原灰度值不变，如果大于阈值 T 时就用它们的平均值来代替该点的灰度值。这样就可以大大减少模糊的程度。实现方法：以图4-19（a)和(b)作模板，扫过全部图像，即可完成平滑处理。边缘处理：1）、在原图像上补上行和列，在处理；2）、处理后重复一下边缘行或列的结果。4.2.2 4.2.2 低通滤波法低通滤波法 这种方法是一种频域处理法。在分析图像信号的频率特性时，一幅图像的边缘、跳跃部分以及颗粒噪声代表图像信号的高频分量，而大面积的背景区则代表图像信号的低频分量。用滤波的方法滤除其高频部分就能去掉噪声，使图像得到平滑。由卷积定理可知 (4-23)(4-23)其中:F(u,v)是含有噪声的图像的傅立叶变换;G(u,v)是平滑处理后的图像之傅立叶变换;H(u,v)是传递函数。选择传递函数 H(u,v)，利用 H(u,v)使 F(u,v)的高频分量得到衰减，得到 G(u,v),后再经反傅立叶变换就可以得到所希望的平滑图像 g(x,y)了。根据前面的分析，显然 H(u,v)应该具有低通滤波特性，所以这种方法叫低通滤波法平滑化处理。低通滤波平滑化处理流程如图4-20所示。FFTH(u,v)IFFT图4-20 线性滤波器处理框图 常用的低通滤波器有如下几种：理想低通滤波器 布特沃斯（Butterworth）低通滤波器 指数低通滤波器 梯形低通滤波器 一个理想的二维低通滤波器的传递函数由下式表示：1.理想低通滤波器理想低通滤波器 (4(4-24)24)式中 D0 是一个规定的非负的量，叫做理想低通滤波器的截止频率。D(u,v)是从频率颊的原点到(u,v)点的距离，即(4-25(4-25)H(u,v)对(u,v)来说是一幅三维图形。H(u,v)的剖面图如图4-21所示。将剖面图绕纵轴旋转360就可以得到整个滤波器的传递函数。所谓理想低通滤波器是指以截频 D0 为半径的圆内的所有频率都能无损地通过，而在截频之外的频率分量完全被衰减。理想低通滤波器可以用计算机模拟实现，但是却不能用电子元器件来实现。图4-21 理想低通滤波器剖面图 理想低通滤波器平滑处理的概念是清晰的，但在处理过程中会产生较严重的模糊和振铃现象。这种现象正是由于傅立叶变换的性质决定的。因为滤波过程是由式（4-23）描述的，由卷积定理可知在空域中则是一种卷积关系，即：(4-26)(4-26)既然 H(u,v)是理想的矩形特性，那么它的反变换 h(x,y)的特性必然会产生无限的振铃特性。经与 f(x,y)卷积后则给 g(x,y)带来模糊和振铃现象，D0 越小这种现象越严重，当然，其平滑效果也就较差。这是理想低通不可克服的弱点。一个 n 阶布特沃斯低通滤波器的传递函数由下式表示 2.2.布特沃斯（布特沃斯（ButterworthButterworth）低通滤波器）低通滤波器 (4-27)(4-27)式中 D0 为截止频率，D(u,v)的值由下式决定(4-28)(4-28)布特沃斯低通滤波器又称最大平坦滤波器。它的通带与阻带之间没有明显的不连续性。也就是说，在通带和阻带之间有一个平滑的过渡带。通常把 H(u,v)下降到某一值的那一点定为截止频率 D0。在式（4-29）中是把 H(u,v)下降到原来值的1/2 时的 D(u,v)定为截频点 D0。一般情况下常常采用下降到 H(u,v)最大值的 那一点为截止频点。这样，式(4-27)可修改为式(4-29)的形式(4-29)(4-29)布特沃斯低通滤波器 H(u,v)的剖面图如图4-22所示。与理想低通滤波器的处理结果相比，经布特沃斯滤波器处理过的图像模糊程度会大大减少。因为它的 H(u,v)不是陡峭的截止特性，它的尾部会包含有大量的高频成分。图4-22 布特沃斯低通滤波器传递函数特性 布特沃斯低通滤波器的特点：1）、由于有平缓的过渡带，图像将不会有振铃现象。2）、模糊程度大大减小。3.指数低通指数低通滤滤波器波器 在图像处理中常用的另一种平滑滤波器是指数低通滤波器。它的传递函数如下式表示(4-30)(4-30)式中 D0 为截频，D(u,v)由下式决定(4-31)(4-31)(4-32)(4-32)(4-33)(4-33)由于指数低通滤波器有更快的衰减率，所以，经指数低通滤波的图像比布特沃斯低通滤波器处理的图像稍模糊一些。由于指数低通滤波器的传递函数也有较平滑的过渡带，所以图像中也没有振铃现象。图4-23 指数低通滤波器传递函数径向剖面图 4.4.梯形低通滤波器梯形低通滤波器 梯形低通滤波器传递函数的形状介于理想低通滤波器和具有平滑过渡带的低通滤波器之间。它的传递函数由下式表示(4-34)(4-34)其中 D(u,v)=u2+v21/2，在规定D0 和 D1 时要满足 D0 D1 的条件。一般为了方便，把传递函数的第一个转折点 D0 定义为截止频率；第二个变量 D1 可以任意选取只要 D1 大于 D0 就可以。梯形低通滤波器传递函数的剖面如图4-24所示。图4-24 梯形低通滤波器传递函数剖面图 由于梯形滤波器的传递函数特性介于理想低通滤波器和具有平滑过渡带滤波器之间，所以其处理效果也介于其两者中间。梯形滤波法的结果有一定的振铃现象。用低通滤波器进行平滑处理可以使噪声伪轮廓等寄生效应减低到不显眼的程度，但是由于低通滤波器对噪声等寄生成分滤除的同时，对有用高频成分也滤除，因此，这种去噪的美化处理是以牺牲清晰度为代价而换取的。如果一幅图像包含有加性噪声，这些噪声对于每个坐标点是不相关的，并且其平均值为零，在这种情况下就可能采用多图像平均法来达到去掉噪声的目的。设 g(x,y)为有噪声图像,n(x,y)为噪声，f(x,y)为原始图像，可用下式表示：4.2.3 4.2.3 多图像平均法多图像平均法多图像平均法是把一系列有噪声的图像gj(x,y)迭加起来，然后再取平均值以达到平滑的目的。(4-35)具体做法如下：取 M 幅内容相同但含有不同噪声的图像，将它们迭加起来，然后作平均计算，如下式所示(4-36)(4-36)由此得出(4-37)(4-37)(4-38)(4-38)(4-39)(4-39)由上二式可见，M 增加则像素值的方差就减小，这说明由于平均的结果使得由噪声造成的像素灰度值的偏差变小。当作平均处理的噪声图像数目增加时，其统计平均值就越接近原始无噪声图像。这种方法在实际应用中的最大困难在于把多幅图像配准起来，以便使相应的像素能正确地对应排列。图4-25 示出了图像平滑处理的效果，其中（a）是待处理图像，（b）是处理后的图像。图 4-25 图像平滑处理效果
