数字图像处理实验六噪声添加和空域滤波

1. 进入 Matlab 7.11 集成开发环境，并打开脚本编辑器。 2. 在脚本编辑器中编写一段程序，要求： (1) 读入存放在 eight.tif 中的原始图像； (2) 调用 imnoise 函数向该图像中添加 2% 的椒盐噪声 ( 参数为“ salt & pepper ” ) ； (3) 在同一窗口中，显示原始图像和加噪后的图像； (4) 将编写的程序保存为 Example3_2_5AddSaltPepper1.m 。 3. 修改上述程序，将添加的噪声改为均值为 0 、方差为 0.001 的高斯噪声 ( 参数 为“ gaussian ” ) ，并将编写的程序另存为 Example3_2_5AddSaltPepper2.m 。 4. 运行这两个程序，修改噪声类型和参数大小，比较噪声对原始图像的影响。 5. 在脚本编辑器中编写一段程序，要求： (1) 读入存放在 eight.tif 中的原始图像； (2) 分别向该图像中添加椒盐噪声； (3) 分别对加噪图像进行 3× 3 平均值滤波和中值滤波，并在同一窗口中， 显示原始图像、加噪后的图像和两种滤波结果图像； (4) 将编写的程序保存为 Example3_2_6MedianFilter1.m 。 6. 修改上述程序，将添加的噪声改为均值为 0 、方差为 0.001 的高斯噪声 ( 参数 为“ gaussian ” ) ，并将编写的程序另存为 Example3_2_6MedianFilter2.m 。 7. 运行这两个程序，修改噪声类型和参数大小，比较不同滤波方法的滤波效果

%%Example3_2_5AddSaltPepper1
clear;
clc;
TIF=imread('eight.tif');%%读入TIF图像%%该图像中添加2%的椒盐噪声
J=imnoise(TIF,'salt & pepper',0.02);
H.color=[1 1 1];%%设置背景为白色，打开原图像
figure(H);
subplot(121);
imshow(J);
title('0.02加噪图');
subplot(122);
imshow(TIF);
title('原图');

%%Example3_2_5AddSaltPepper2
clear;
clc;
TIF=imread('eight.tif');%%读入TIF图像%%将添加的噪声改为均值为0方差为0.001的高斯噪声
J=imnoise(TIF,'gaussian',0,0.01);
H.color=[1 1 1];%%设置背景为白色，打开原图像
figure(H);
subplot(121);
imshow(J);
title('均值为 0、方差为 0.001 的高斯噪声图');
subplot(122);
imshow(TIF);
title('原图');

%%Example3_2_6MedianFilter1
clear;
clc;
TIF=imread('eight.tif');%%读入TIF图像%%该图像中添加2%的椒盐噪声
J=imnoise(TIF,'salt & pepper',0.02);
J=rgb2gray(J);%%行3×3平均值滤波
Y=ones(3,3)/9;
B_1=imfilter(J,Y);%%中值滤波
B_2=medfilt2(J);%%对比
H.color=[1 1 1];%%设置背景为白色，打开原图像
figure(H);
subplot(221);
imshow(TIF);
title('原图');
subplot(222);
imshow(J);
title('0.02加噪图');
subplot(223);
imshow(B_1);
title('3×3平均值滤波图');
subplot(224);
imshow(B_2);
title('中值滤波图');

%%Example3_2_6MedianFilter2
clear;
clc;
TIF=imread('eight.tif');%%读入TIF图像%%将添加的噪声改为均值为0方差为0.001的高斯噪声
J=imnoise(TIF,'gaussian',0,0.01);
J=rgb2gray(J);%%行3×3平均值滤波
Y=ones(3,3)/9;
B_1=imfilter(J,Y);%%中值滤波
B_2=medfilt2(J);%%对比
H.color=[1 1 1];%%设置背景为白色，打开原图像
figure(H);
subplot(221);
imshow(TIF);
title('原图');
subplot(222);
imshow(J);
title('均值为 0、方差为 0.001 的高斯噪声图');
subplot(223);
imshow(B_1);
title('3×3平均值滤波图');
subplot(224);
imshow(B_2);
title('中值滤波图');

本文来自互联网用户投稿，文章观点仅代表作者本人，不代表本站立场，不承担相关法律责任。如若转载，请注明出处。 如若内容造成侵权/违法违规/事实不符，请点击【内容举报】进行投诉反馈！