matlab图像的空域滤波.docVIP

实验四 图像的空域滤波 一、实验目的 1、理解图像滤波的基本定义及目的； 2、掌握空域滤波的基本原理及方法； 3、掌握用MATLAB语言进行图像的空域滤波的方法。 1、读出“girl.bmp”这幅图像。给读出的图像加入高斯噪声。 2、分别采用不同大小的模板对加有噪声的图像进行均值滤波，用一个图像处理对话框显示原图像、加有噪声的图像及均值滤波的图像。比较结果。 二、实验原理 均值滤波 2、中值滤波 MATLAB图像处理工具箱提供了基于卷积的图像滤波函数filter2。filter2的语法格式为： Y = filter2(h, X) 其中Y = filter2(h,X)返回图像X经算子h滤波后的结果，默认返回图像Y与输入图像X大小相同。 其实filter2和conv2是等价的。MATLAB在计算filter2时先将卷积核旋转180度，再调用conv2函数进行计算。 fspecial函数用于创建预定义的滤波算子，其语法格式为： h = fspecial(type) h = fspecial(type, parameters) 参数type指定算子类型，parameters指定相应的参数，具体格式为： type=average，为均值滤波，参数parameters为n，代表模版尺寸，用向量表示，默认值为[3,3]。 type= gaussian，为高斯低通滤波器，参数parameters有两个，n表示模版尺寸，默认值为[3,3]，sigma表示滤波器的标准差，单位为像素，默认值为0.5。 type= laplacian，为拉普拉斯算子，参数parameters为alpha，用于控制拉普拉斯算子的形状，取值范围为[0,1]，默认值为0.2。 type= log，为拉普拉斯高斯算子，参数parameters有两个，n表示模版尺寸，默认值为[3,3]，sigma为滤波器的标准差，单位为像素，默认值为0.5 type= prewitt，为prewitt算子，用于边缘增强，无参数。 type= sobel，为著名的sobel算子，用于边缘提取，无参数。 type= unsharp，为对比度增强滤波器，参数alpha用于控制滤波器的形状，范围为[0,1]，默认值为0.2。 三、实验要求 1、读出“girl.bmp”这幅图像。给读出的图像加入高斯噪声。 分别采用不同大小的模板对加有噪声的图像进行均值滤波，用一个图像处理对话框显示原图像、加有噪声的图像及均值滤波的图像。比较结果。 A=imread(C:\Documents and Settings\Administrator\桌面\数字图像实验实验\girl.bmp, bmp); B=imnoise(A, gaussian, 0.02); subplot(2,2,1), imshow(A), title(原图像); subplot(2,2,2), imshow(B), title(加噪声图像); C=fspecial(average, 3); D=filter2(C, B); E=fspecial(average, 5); F=filter2(E, B); subplot(2,2,3), imshow(D,[]), title(3乘3模板); subplot(2,2,4), imshow(F,[]), title(5乘5模板); 分别采用不同大小的模板对加有噪声的图像进行中值滤波，用一个图像处理对话框显示原图像、加有噪声的图像及中值滤波的图像。比较结果。 A=imread(C:\Documents and Settings\Administrator\桌面\数字图像实验实验\girl.bmp, bmp); B=imnoise(A, gaussian, 0.02); subplot(2,2,1), imshow(A), title(原图像); subplot(2,2,2), imshow(B), title(加噪声图像); C=medfilt2(B, [3,3]); D=medfilt2(B, [5,5]); subplot(2,2,3), imshow(C), title(3乘3模板); subplot(2,2,4), imshow(D), title(5乘5模板); 采用大小相同的模板对加有噪声的图像分别进行均值滤波和中值滤波，用一个图像处理对话框显示原图像、加有噪声的图像、均值滤波的图像及中值滤波的图像。比较结果。 A=imread(C:\Documents and Settings\Administrator\桌面\数字图像实验实验\girl.bmp, bmp); B=imnoise(A, gaussian, 0.02); subplot(2,2,1), imshow(A), title(原