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

发表于 分类于 ， 阅读次数： Changyan：
本文字数： 1.9k 阅读时长 ≈ 2 分钟
  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('中值滤波图');

img