西交大数字图像处理第四次作业

《西交大数字图像处理第四次作业》由会员分享，可在线阅读，更多相关《西交大数字图像处理第四次作业（21页珍藏版）》请在金锄头文库上搜索。

1、 数字图像处理第四次作业数字图像处理第四次作业姓名：班级：学号：提交日期：2015 年 3 月 31 日1、空域低通滤波器：分别用高斯滤波器和中值滤波器去平滑测试图像、空域低通滤波器：分别用高斯滤波器和中值滤波器去平滑测试图像 test1 和和 test2，模板，模板 大小分别是大小分别是 3x3，5x5，7x7； （1）问题分析：问题分析：空域滤波是直接对图像的数据做空间变换达到滤波的目的。它是一种邻域 运算，其机理就是在待处理的图像中逐点地移动模板，滤波器在该点地响应通过事先定义 的滤波器系数与滤波模板扫过区域的相应像素值的关系来计算。如果输出像素是输入像素 邻域像素的线性组合则称为线性滤

2、波（例如最常见的均值滤波和高斯滤波） ，否则为非线性 滤波（中值滤波、边缘保持滤波等） 。 空域滤波器从处理效果上可以平滑空间滤波器和锐化空间滤波器：平滑空间滤波器用 于模糊处理和减小噪声，经常在图像的预处理中使用；锐化空间滤波器主要用于突出图像 中的细节或者增强被模糊了的细节。模板在源图像中移动的过程中，当模板的一条边与图像轮廓重合后，模板中心继续向 图像边缘靠近，那么模板的某一行或列就会处于图像平面之外，此时最简单的方法就是将 模板中心点的移动范围限制在距离图像边缘不小于（n1）/2 个像素处，单处理后的图像 比原始图像稍小。如果要处理整幅图像，可以在图像轮廓边缘时用全部包含于图像中的模

3、板部分来滤波所有图像，或者在图像边缘以外再补上一行和一列灰度为零的像素点（或者 将边缘复制补在图像之外） 。 中值滤波器的设计：中值滤波器的设计：中值滤波器是一种非线性统计滤波器，它的响应基于图像滤波器包围的图像区域中像 素的排序，然后由统计排序的中间值代替中心像素的值。它比小尺寸的线性平滑滤波器的 模糊程度明显要低，对处理脉冲噪声（椒盐噪声）非常有效。中值滤波器的主要功能是使 拥有不同灰度的点看起来更接近于它的邻近值，去除那些相对于其邻域像素更亮或更暗， 并且其区域小于滤波器区域一半的孤立像素集。在一维的情况下，中值滤波器是一个含有奇数个像素的窗口。在处理之后，位于窗口 正中的像素的灰度值，

4、用窗口内各像素灰度值的中值代替。例如若窗口长度为 5，窗口中 像素的灰度值为 80、90、200、110、120，则中值为 110，因为按小到大（或大到小）排序 后，第三位的值是 110。于是原理的窗口正中的灰度值 200 就由 110 取代。如果 200 是一 个噪声的尖峰，则将被滤除。然而，如果它是一个信号，则滤波后就被消除，降低了分辨 率。因此中值滤波在某些情况下抑制噪声，而在另一些情况下却会抑制信号。 将中值滤波推广到二维的情况。二维窗口的形式可以是正方形、近似圆形的或十字形 等。本次作业使用正方形模板进行滤波，它的中心一般位于被处理点上。窗口的大小对滤 波效果影响较大。（2）处理结果

5、：处理结果：test1.pgm 中值滤波：中 中 中 中 中 中 中 中 test1.pgm(3x3)(中 中 medfilt2中 中 中中 中 中 中 中 中 中 中 test1.pgm(5x5)(中 中 medfilt2中 中 中中 中 中 中 中 中 中 中 test1.pgm(7x7)(中 中 medfilt2中 中 中test2.tif 中值滤波：中 中 中 中 中 中 test2.tif中 3x3中中 中 中 中 中 中 test2.tif中 5x5中中 中 中 中 中 中 test2.tif中 7x7中中 中 中 中 中 中 中 中 test2.tif(3x3)(中 中 medf

6、ilt2中 中 中中 中 中 中 中 中 中 中 test2.tif(5x5)(中 中 medfilt2中 中 中中 中 中 中 中 中 中 中 test2.tif(7x7)(中 中 medfilt2中 中 中（3）结果分析及总结：）结果分析及总结： 中值滤波器的主要功能是使拥有不同灰度的点更接近于它的相邻点。事实上，使用mm 中值滤波器来去除那些相对于其邻域像素更亮或更暗并且其区域小于 m2/2（滤波 器区域的一半）的孤立像素族。在这种情况下， “去除“的意思是强制为邻域的中值灰度。 较大的族所受到的影响明显较小。从处理后的图像看，图像的平滑效果较为明显，且受窗 口的影响，窗口越大，平滑效果

7、越明显，图像细节越模糊，尤其是 test2 中人脸图像的眼睛 部分，随着滤波器模板的增大，可以明显的感受到图像模糊的效果。另外，对比编写的程 序的处理结果和 medfilt2 函数的处理结果，可以看出处理效果基本一致。 2、利用固定方差、利用固定方差 sigma=1.5 产生高斯滤波器产生高斯滤波器. 附件有产生高斯滤波器的方法。附件有产生高斯滤波器的方法。 （1）问题分析：）问题分析： 高斯滤波器的设计：高斯滤波器的设计： 高斯滤波是一种根据高斯函数的形状来选择模板权值的线性平滑滤波方法。高斯平滑滤波器对去除服从正态分布的噪声是很有效果的。一维零均值高斯函数为 。其中，高斯分布参数决定了高斯

8、滤波器的宽度。对图像来说，常用二维零均值离散高斯函数作平 滑滤波器，函数表达式如下：（1）高斯函数具有 5 个重要性质：二维高斯函数具有旋转对称性，即滤波器在各个方向上的平滑程度是相同的。一般来说 一幅图像的边缘方向是不知道的。因此，在滤波之前是无法确定一个方向比另一个方向上 要更多的平滑的。旋转对称性意味着高斯滤波器在后续的图像处理中不会偏向任一方向。高斯函数是单值函数。高斯滤波器用像素邻域的加权均值来代替该点的像素值，而每一 邻域像素点的权值是随着该点与中心点距离单调递减的。这一性质是很重要的，因为边缘 是一种图像局部特征。如果平滑运算对离算子中心很远的像素点仍然有很大的作用，则平 滑运算

9、会使图像失真。高斯函数的傅立叶变换频谱是单瓣的。这一性质是高斯函数傅立叶变换等于高斯函数本 身这一事实的直接推论。图像常被不希望的高频信号所污染，而所希望的图像特征，既含 有低频分量，又含有高频分量。高斯函数傅立叶变换的单瓣意味着平滑图像不会被不需要 的高频信号所污染，同时保留了大部分所需要的信号。高斯滤波器的宽度(决定着平滑程度)是由参数 表征的，而且 和平滑程度的关系是 非常简单的。 越大，高斯滤波器的频带就越宽，平滑程度就越好。通过调节平滑程度参 数 ，可在图像特征分量模糊(过平滑)与平滑图像中由于噪声和细纹理所引起的过多的不 希望突变量(欠平滑)之间取得折衷。由于高斯函数的可分离性，大

10、高斯滤波器可以有效实现。通过二维高斯函数的卷积可以 分两步来进行，首先将图像与一维高斯函数进行卷积，然后将卷积的结果与方向垂直的相 同一维高斯函数进行卷积。因此，二维高斯滤波的计算量随滤波模板宽度成线性增长而不 是成平方增长。这些性质使得它在早期的图像处理中特别有用，表明高斯平滑滤波器无论 在空间域还是在频率域都是十分有效的低通滤波器。 （2）处理结果：）处理结果：test1.pgm 高斯滤波：中 中 中 中 中 中 中 中 test1.pgm(3x3)(中 中 fspecial中 中 )中 中 中 中 中 中 中 中 test1.pgm(5x5)(中 中 fspecial中 中 )中 中 中

11、 中 中 中 中 中 test1.pgm(7x7)(中 中 fspecial中 中 )test2.tif 高斯滤波：中 中 中 test2.tif5010015020025030035040045050050100150200250300350400450500中 中 中 中 中 中 中 中 test2.tif中 3x3中5010015020025030035040045050050100150200250300350400450500中 中 中 中 中 中 中 中 test2.tif中 5x5中5010015020025030035040045050050100150200250300350

12、400450500中 中 中 中 中 中 中 中 test2.tif中 7x7中5010015020025030035040045050050100150200250300350400450500中 中 中 中 中 中 中 中 test2.tif(3x3)(中 中 fspecial中 中 )中 中 中 中 中 中 中 中 test2.tif(5x5)(中 中 fspecial中 中 )中 中 中 中 中 中 中 中 test2.tif(7x7)(中 中 fspecial中 中 )（3）结果分析及总结：）结果分析及总结：高斯滤波是一种线性平滑滤波，适用于消除高斯噪声，广泛应用于图像处理的减噪过

13、程。通俗的讲，高斯滤波就是对整幅图像进行加权平均的过程，每一个像素点的值，都由 其本身和邻域内的其他像素值经过加权平均后得到。高斯滤波的具体操作是：用一个模板 （或称卷积、掩模）扫描图像中的每一个像素，用模板确定的邻域内像素的加权平均灰度 值去替代模板中心像素点的值。从处理后的图像看，图像的平滑效果较为明显，且受窗口 的影响，窗口越大，平滑效果越明显，图像细节越模糊。将编写的程序处理结果和调用 MATLAB 中函数处理结果对比，发现处理结果基本相同。 3、利用高通滤波器滤波测试图像、利用高通滤波器滤波测试图像 test3,4：包括：包括 unsharp masking, Sobel edge

14、detector, and Laplace edge detection；Canny algorithm。 （1）问题分析：问题分析： 锐化滤波能减弱或消除图像中的低频率分量，但不影响高频率分量。因为低频分量对 应图像中灰度值缓慢变化的区域，因而与图像的整体特性，如整体对比度和平均灰度值等 有关。锐化滤波将这些分量滤去可使图像反差增加，边缘明显。在实际应用中，锐化滤波 可用于增强被模糊的细节或者低对比度图像的目标边缘。图像锐化的主要目的有两个：一是增强图像边缘，使模糊的图像变得更加清晰，颜色 变得鲜明突出，图像的质量有所改善，产生更适合人眼观察和识别的图像；二是希望经过 锐化处理后，目标物体的

15、边缘鲜明，以便于提取目标的边缘、对图像进行分割、目标区域识别、区域形状提取等，为进一步的图像理解与分析奠定基础。由于锐化使噪声受到比信号还要强的增强，所以要求锐化处理的图像有较高的信噪比； 否则，锐化后图像的信噪比更低。 反锐化掩膜图像增强（反锐化掩膜图像增强（unsharpunsharp maskingmasking） 图像的反锐化掩蔽算法可以表示为：(1)其中 fs(x,y)表示经过反锐化掩蔽得到的锐化图像，是 f(x,y)的模糊形式。反锐化掩蔽进一步的普遍形式称为高提升滤波。在图像中任何一点(x,y)处，高提升滤波后的图像可hbf定义如下：(2)其中 A1，与前式一样，是f的模糊形式，此

16、式也可以写成：(3)结合式(1)，可以得到：(4)这一表达式可计算高提升滤波图像。如果选择拉普拉斯变换，式(4)变成：(5)高提升滤波处理可以通过任何一个图 1 所示的掩模得以实现。当 A=1 时，高提升滤波处理就是标准的拉普拉斯变换。随着 A 超过 l 不断增大，锐化处理的效果越来越不明显。最终，当 A 足够大时，高提升图像将近似等于经常数调制的图像。图 1 高频提升滤波技术可以用其中一种掩膜来实现（）1A 索贝尔边缘检测（索贝尔边缘检测（SobelSobel edgeedge detectordetector）索贝尔算子（Sobel operater）主要用作边缘检测，在技术上，它是一离散性差分算子， 用来运算图像亮度函数的灰度之近似值。在图像的任何一点使用此算子，将会产生对应的 灰度矢量或其法矢量。 该算子包含两组 3x3 的矩阵，分别为横向及纵向，将之与图像作平面卷积，即可分别 得出横向及纵向的亮度差分近似值。如果以 A 代表原始图像，Gx 及 Gy 分别代表经横向及 纵向边缘检测