实验4 线性空间滤波器操作与非线性空间滤波器操作

• 1.对fatBaby.bmp分别添进行均值滤波和高斯滤波（使用fspecial函数定义模板，分别都使用33,55及7*7模板），并比较高斯滤波与均值滤波的效果。

I = imread('miaoa.jpg');
w1 = fspecial('average',3);
w2 = fspecial('gaussian',3);
g1 = imfilter(I,w1,'replicate');
g2 = imfilter(I,w2,'replicate');
subplot(3,2,1); imshow(g1);title('3*3均值滤波');
subplot(3,2,2); imshow(g2);title('3*3高斯滤波');

w1 = fspecial('average',5);
w2 = fspecial('gaussian',5);
g1 = imfilter(I,w1,'replicate');
g2 = imfilter(I,w2,'replicate');
subplot(3,2,3); imshow(g1);title('5*5均值滤波');
subplot(3,2,4); imshow(g2);title('5*5高斯滤波');


w1 = fspecial('average',7);
w2 = fspecial('gaussian',7);
g1 = imfilter(I,w1,'replicate');
g2 = imfilter(I,w2,'replicate');
subplot(3,2,5); imshow(g1);title('7*7均值滤波');
subplot(3,2,6); imshow(g2);title('7*7高斯滤波');

• 2. 对fatBaby.bmp使用定义矩阵的方式定义模板（3*3模板），再进行均值滤波和高斯滤波，并比较高斯滤波与均值滤波的效果。

f=imread('miaoa.jpg');
subplot(331);
imshow(f);
title('原图')
subplot(332);
w=[1,1,1;1,1,1;1,1,1;]/9;
g1=imfilter(f,w);
imshow(g1);
title('均值3*3')
subplot(333);
w = [1,2,1;2,4,2;1,2,1;]/9;
g2=imfilter(f,w);
imshow(g2);
title('高斯3*3');

• 3.读入一幅自己拍摄的图像，转化为灰度图，添加椒盐噪声后，分别进行均值滤波，高斯滤波与中值滤波操作（使用fspecial函数定义模板，medfilt2函数实现中值滤波），比较均值滤波、高斯滤波与中值滤波的去噪声效果。

A = imread('miaoa.jpg');
f = rgb2gray(A);
subplot(232);
g = imnoise(f, 'salt & pepper', 0.2);
imshow(g);
title('添加椒盐噪声');
subplot(234)
w1 = fspecial('average');
g1 = imfilter(g, w1, 'replicate');
imshow(g1);
title('均值滤波');
subplot(235)
w2 = fspecial('gaussian');
g2 = imfilter(g, w2, 'replicate');
imshow(g2);
title('高斯滤波');
subplot(236)
g3 = medfilt2(g);
imshow(g3);
title('中值滤波');

• 4.读入一幅自己拍摄的图像，转化为灰度图，添加椒盐噪声后，分别进行最小值滤波、最大值滤波及中值滤波，并比较运行效果（使用ordfilt2函数）。

f = imread('miaoa.jpg');
f = rgb2gray(f);
f = imnoise(f,'salt & pepper',0.02);
f1 = ordfilt2(f,5,ones(3,3));
subplot(131);
imshow(f1);
title('最小值滤波');
 
subplot(132);
g2=ordfilt2(f,3*3,ones(3,3));
imshow(g2);
title('最大值滤波');
 
subplot(133);
g3=medfilt2(f);
imshow(g3);
title('中值滤波');