添加高斯噪声和椒盐噪声,并用理想低通滤波器进行去噪
t=imread('C:\Users\Administrator\Desktop\123.jpg');
subplot(3,2,1)
imshow(t),title('原图');
picture_gaussian_noise=imnoise(t,'gaussian',0,0.02);
subplot(3,2,2)
imshow(picture_gaussian_noise),title('添加均值为0,方差为0.02的高斯噪声');
picture_salt_noise=imnoise(t,'salt & pepper',0.1);
subplot(3,2,3)
imshow(picture_salt_noise),title('加入噪声密度:0.1的椒盐噪声');
s=fftshift(fft2(picture_salt_noise));
[a,b]=size(s);
a0=round(a/2);
b0=round(b/2);
d=50;
for i=1:a
for j=1:b
distance=sqrt((i-a0)^2+(j-b0)^2);
if distance<=d
h=1;
else
h=0;
end
s(i,j)=h*s(i,j);
end
end
s=uint8(real(ifft2(ifftshift(s))));
subplot(3,2,4),imshow(s);
title('高斯噪声理想低通滤波所得图像');
subplot(3,2,5),imshow(s);
title('椒盐噪声理想低通滤波所得图像