数字图像处理实验报告3.docx

数字图像处理与分析

实验报告

学院：班级：姓名：学号：

实验三空域滤波

一、 实验目的

了解空域滤波的方法

掌握几种模版的基本原理

二、 实验内容

使用函数fspecial()生成几种特定的模板

使用函数imfilter()配合模板对图象数据进行二维卷积

比较各种滤波器的效果

三、 实验步骤

线性平滑（低通）滤波器。

用h=fspecial(‘average’) 得到的h为3×3的邻域平均模板，然后用h进行平滑处理。相关程序如下：

x=imread(C:\DocumentsandSettings\Administrator\桌面\1.TIF);h=fspecial(average)

y=imfilter(x,h);

得到相应的结果为：

改变模板大小重试

改变fspecial()的参数为高斯函数，再进行同样的平滑处理，观察其结果相关程序如下：

x=imread(C:\DocumentsandSettings\Administrator\桌面\1.TIF);h=fspecial(gaussian)

y=imfilter(x,h);

得到相应的结果为:

改变模板大小为[7,7]、[9,9]、[11,11]，观察实验得到的结构有何变化相关程序分别如下：

x=imread(C:\DocumentsandSettings\Administrator\桌面\1.TIF);h=fspecial(gaussian,[7,7])

y=imfilter(x,h);

x=imread(C:\DocumentsandSettings\Administrator\桌面\1.TIF);h=fspecial(gaussian,[9,9])

y=imfilter(x,h);

x=imread(C:\DocumentsandSettings\Administrator\桌面\1.TIF);h=fspecial(gaussian,[11,11])

y=imfilter(x,h);

得到的结果分别为：

用mesh函数直观观察模板相关程序如下：

h1=fspecial(gaussian,100,3);[xy]=meshgrid(1:100);mesh(x,y,h1);h2=fspecial(gaussian,100,10);figure

mesh(x,y,h2)

得到相应的图像如下：

图1 图2

进行中值滤波：对图象添加胡椒盐噪声（saltpepper）,然后分别用前面的均值滤波器和中值滤波器medfilt2()进行去噪声处理，并给出比较结果

相关程序如下：

x=imread(C:\DocumentsandSettings\Administrator\桌面\1.TIF);h=imnoise(x,saltpepper,0.02);

y=medfilt2(rgb2gray(h));subplot(1,2,1)

imshow(x)subplot(1,2,2)imshow(y)

得到的图像如下：

5. 线性锐化滤波器的模板为h=[ -1 -1 -1;-1 8 -1;-1 -1 -1;]，用该模板对图象进行锐化处理

相关程序如下：

x=imread(C:\DocumentsandSettings\Administrator\桌面\1.TIF);h=[ -1 -1 -1; -1 8 -1;-1 -1 -1 ]

y=imfilter(x,h);subplot(1,2,1)imshow(x)subplot(1,2,2)imshow(y)

得到相应的图像如下：

fspecial(laplacian)函数得到锐化模板代码：

x=imread(C:\DocumentsandSettings\Administrator\桌面\1.TIF);h=fspecial(laplacian)

y=imfilter(x,h);subplot(1,2,1)imshow(x)subplot(1,2,2)imshow(y)

处理结果：

6. 索贝尔（sobel）算子 h1=[-1 0 1;-2 0 2;-1 0 1;]

h2=[ 1 2 1;0 0 0;-1 -2-1;]

分别使用以上算子对图象进行处理可以得到图象的边缘

相关程序如下：

x=imread(C:\DocumentsandSettings\Administrator\桌面\1.TIF);h1=[-1 0 1;-2 0 2;-1 0 1;]

y=imfilter(x,h);subplot(1,2,1)imshow(x)subplot(1,2,2)imsh