matlab中提供滤波器的种类有,数字图像处理实验报告2
数字图像处理实验报告
3. 傅立叶变换
1) 读出woman.tif这幅图像，对其进行快速傅立叶变换，分别显示其幅度图像和相位 图像。
仅对相位部分进行傅立叶反变换后查看结果图像。
2) 仅对幅度部分进行傅立叶反变换后查看结果图像。 3) 将图像的傅立叶变换F置为其共轭后进行反变换，比较新生成图像与原始图像的差 异。
4. 平滑频域滤波
1) 设计理想低通滤波器、巴特沃斯低通滤波器和高斯低通滤波器，截至频率自选，分 别给出各种滤波器的透视图。
2) 读出test_pattern.tif这幅图像，分别采用理想低通滤波器、巴特沃斯低通滤波器和高 斯低通滤波器对其进行滤波(截至频率自选)，再做反变换，观察不同的截止频率下采用不同低通滤波器得到的图像与原图像的区别，特别注意振铃效应。(提示:1)在频率域滤波同样要注意到填充问题；2)注意到(-1)x+y;)
5. 锐化频域滤波
1) 设计理想高通滤波器、巴特沃斯高通滤波器和高斯高通滤波器，截至频率自选，分 别给出各种滤波器的透视图。
2) 读出test_pattern.tif这幅图像，分别采用理想高通滤波器、巴特沃斯高通滤波器和高 斯高通滤波器对其进行滤波(截至频率自选)，再做反变换，观察不同的截止频率下采用不同高通滤波器得到的图像与原图像的区别。
三、实验原理与算法分析
1. 空间域增强
空间域滤波是在图像空间中借助模板对图像进行领域操作，处理图像每一个像素的取值都是根据模板对输入像素相应领域内的像素值进行计算得到的。空域滤波基本上是让图像在频域空间内某个范围的分量受到抑制，同时保证其他分量不变，从而改变输出图像的频率分布，达到增强图像的目的。
空域滤波一般分为线性滤波和非线性滤波两类。线性滤波器的设计常基于对傅立叶变换 的分析，非线性空域滤波器则一般直接对领域进行操作。各种空域滤波器根据功能主要分为 平滑滤波器和锐化滤波器。平滑可用低通来实现，平滑的目的可分为两类：一类是模糊，目 的是在提取较大的目标前去除太小的细节或将目标内的小肩端连接起来；另一类是消除噪 声。锐化可用高通滤波来实现，锐化的目的是为了增强被模糊的细节。结合这两种分类方法， 可将空间滤波增强分为四类：
线性平滑滤波器(低通) 非线性平滑滤波器(低通) 线性锐化滤波器(高通) 非线性锐化滤波器(高通)
空间滤波器都是基于模板卷积，其主要工作步骤是：
1) 将模板在图中移动，并将模板中心与图中某个像素位置重合； 2) 将模板上的系数与模板下对应的像素相乘； 3) 将所有乘积相加；
4) 将和(模板的输出响应)赋给图中对应模板中心位置的像素。
2. 平滑滤波器
第 12 页 共 38 页
12
数字图像处理实验报告
1) 线性平滑滤波器
线性低通平滑滤波器也称为均值滤波器，这种滤波器的所有系数都是正数，对3×3 的模板来说，最简单的是取所有系数为1，为了保持输出图像任然在原来图像的灰度值范围 内，模板与象素邻域的乘积都要除以9。
MATLAB 提供了fspecial 函数生成滤波时所用的模板，并提供filter2 函数用指定的滤 波器模板对图像进行运算。函数fspecial 的语法格式为： h=fspecial(type);
h=fspecial(type,parameters);
其中参数type 指定滤波器的种类，parameters 是与滤波器种类有关的具体参数。
表2.1 MATLAB中预定义的滤波器种类
MATLAB 提供了一个函数imnoise 来给图像增添噪声，其语法格式为：
J=imnoise(I,type);
J=imnoise(I,type,parameters);
参数type 指定噪声的种类，parameters 是与噪声种类有关的具体参数。参数的种类见 表2.2。
第 13 页 共 38 页 13
数字图像处理实验报告
表2.2 噪声种类及参数说明
2) 非线性平滑滤波器
中值滤波器是一种常用的非线性平滑滤波器，其滤波原理与均值滤波器方法类似，但计 算的非加权求和，而是把领域中的图像的象素按灰度级进行排序，然后选择改组的中间值作为输出象素值。
MATLAB 提供了medfilt2 函数来实现中值滤波，其语法格式为：
B=medfilt2(A,[m n]); B=medfilt2(A);
其中，A 是原图象，B 是中值滤波后输出的图像。[m n]指定滤波模板的大小，默认模 板为3×3。
3. 锐化滤波器
图像平滑往往使图像中的边界、轮廓变得模糊，为了减少这类不利效果的影响，需要利 用图像锐化技术，使图像的边缘变得清晰。
1) 线性锐化滤波器
线性高通滤波器是最常用的线性锐化滤波器。这种滤波器的中心系数都是正的，而周围 的系数都是负的，所有的系数之和为0。
对3×3 的模板来说，典型的系数取值为：
[-1 -1 -1; -1 8 -1; -1 -1 -1]
事实上这是拉普拉斯算子。语句h=-fspecial(‘laplacian’,0.5)得到的拉普拉斯算子为：
h =-0.3333 -0.3333 -0.3333
-0.3333 2.6667 -0.3333 -0.3333 -0.3333 -0.3333 2) 非线性锐化滤波
邻域平均可以模糊图像，因为平均对应积分，所以利用微分可以锐化图像。图像处理中 最常用的微分方法是利用梯度。常用的空域非线性锐化滤波微分算子有sobel 算子、prewitt 算子、log 算子等。
4. 频域增强
频域增强是利用图像变换方法将原来的图像空间中的图像以某种形式转换到其他空间 中，然后利用该空间的特有性质方便地进行图像处理，最后再转换回原来的图像空间中，从 而得到处理后的图像。
频域增强的主要步骤是：
第 14 页 共 38 页
14
数字图像处理实验报告
选择变换方法，将输入图像变换到频域空间。
在频域空间中，根据处理目的设计一个转移函数，并进行处理。 将所得结果用反变换得到增强的图像。
常用的频域增强方法有低通滤波和高通滤波。
5. 低通滤波
图像的能量大部分集中在幅度谱的低频和中频部分，而图像的边缘和噪声对应于高频部 分。因此能降低高频成分幅度的滤波器就能减弱噪声的影响。由卷积定理，在频域实现低通 滤波的数学表达式：
G(u,v) =H(u,v)F(u,v)
1) 理想低通滤波器(ILPF)
2) 巴特沃斯低通滤波器(BLPF)
3) 指数型低通滤波器(ELPF)
6. 高通滤波
由于图像中的细节部分与其高频分量相对应，所以高通滤波可以对图像进行锐化处理。高通滤波与低通滤波相反，它是高频分量顺利通过，使低频分量受到削弱。高通滤波器和低通滤波器相似，其转移函数分别为:
1) 理想高通滤波器(IHPF)
2) 巴特沃斯高通滤波器(BLPF)
3) 指数型高通滤波器(ELPF)
图像经过高通滤波处理后，会丢失许多低频信息，所以图像的平滑区基本上会消失。所 以，可以采用高频加强滤波来弥补。高频加强滤波就是在设计滤波传递函数时，加上一 个大于0小于1的常数c，即：
H′(u,v) =H(u,v)+c
第 15 页 共 38 页 15