main.m

来自「一个用matlab写的对图象处理方面的函数」· M 代码 · 共 57 行

clear;
url='boy.bmp';%输入图象路径
[A,map]=imread(url);
source=ind2gray(A,map);
%在灰度图像中加入一定量的高斯噪声和椒盐噪声
gaussian=imnoise(source,'gaussian',0,0.02);%叠加均值为0，方差为0.02的高斯噪声
saltpepper=imnoise(source,'salt & pepper',0.04);%叠加密度为0.04的椒盐噪声
figure(1);
subplot(211);
imshow(source);
title('原图像');
subplot(223);
imshow(gaussian);
title('添加高斯噪声后的图像');
subplot(224);
imshow(saltpepper);
title('添加椒盐噪声后的图像');

%对加入噪声的两幅图象进行8邻域的均值滤波
H =  fspecial('average',[3 3]);
gaussian1=imfilter(gaussian,H);
saltpepper1=imfilter(saltpepper,H);
figure(2);
subplot(211);
imshow(source);
title('原图像');
subplot(223);
gaussian1=uint8(gaussian1);
imshow(gaussian1);
title('均值滤波含高斯噪声的图像');
subplot(224);
saltpepper1=uint8(saltpepper1);
imshow(saltpepper1);
title('均值滤波含椒盐噪声的图像');

%对加入噪声的两幅图象进行中值滤波
gaussian2=medfilt2(gaussian);
saltpepper2=medfilt2(saltpepper);
figure(3);
subplot(211);
imshow(source);
title('原图像');
subplot(223);
imshow(gaussian2);
title('中值滤波含高斯噪声的图像');
subplot(224);
imshow(saltpepper2);
title('中值滤波含椒盐噪声的图像');

%PSNR
valueg=psnr(source,gaussian);
valueg1=psnr(source,gaussian1);
valueg2=psnr(source,gaussian2);

values=psnr(source,saltpepper);
values1=psnr(source,saltpepper1);
values2=psnr(source,saltpepper2);