tuxiangchuli.rar_给图像加噪声

I = imread('Jellyfish.jpg'); I = im2double(I); figure, imshow(I); title('原始图像'); %% 加随机噪声 I1 = I + rand(size(I)); I1 = mat2gray(I1); figure, subplot(2,3,1); imshow(I1,[]); title('随机噪声图像'); %% 加颗粒噪声图像 I2 = imnoise(I, 'salt & pepper',0.02); subplot(2,3,2); imshow(I2,[]); title('颗粒噪声图像'); %% 加高斯噪声 I3 = imnoise(I, 'gaussian', 0.02); subplot(2,3,3); imshow(I3,[]); title('高斯噪声图像'); I=rgb2gray(I); %% 加随机噪声 I1 = I + rand(size(I)); I1 = mat2gray(I1); %% 加颗粒噪声图像 I2 = imnoise(I, 'salt & pepper',0.02); %% 加高斯噪声 I3 = imnoise(I, 'gaussian', 0.02); %% 设计低通滤波器 [M,N] = size(I); D0 = 600; m = fix(M/2); n = fix(N/2); for i = 1:M for j = 1:N D = sqrt(i -m)^2 + (j - n)^2; if D <= D0; H(i,j) = 1; else H(i,j) = 0; end % H(i,j) = 1/(1+(D/D0)^(2*nn)); end end %% 图像降噪 %% F1 = fft2(I1); F1 = fftshift(F1); result1 = F1 .* H; result1 = ifftshift(result1); J1 = ifft2(result1); J1 = mat2gray((abs(real(J1)))); subplot(2,3,4); imshow(J1,[]); title('降噪后'); %% F2 = fft2(I2); F2 = fftshift(F2); result2 = F2 .* H; result2 = ifftshift(result2); J2 = ifft2(result2); J2 = mat2gray((abs(real(J2)))); subplot(2,3,5); imshow(J2,[]); title('降噪后'); %% F3 = fft2(I3); F3 = fftshift(F3); result3 = F3 .* H; result3 = ifftshift(result3); J3 = ifft2(result3); J3 = mat2gray((abs(real(J3)))); subplot(2,3,6); imshow(J3,[]); title('降噪后');