数字图像的预处理——从轮胎到硬币到火星。

I = imread('eight.tif');
J = imnoise(I, 'salt & pepper', 0.09);
K = medfilt2(J);
figure, imshow(I); figure, imshow(J); figure, imshow(K);

I = imread('eight.tif');
J = imnoise(I, 'salt & pepper', 0.09);
K = medfilt2(J);
figure, imshow(I); figure, imshow(J); figure, imshow(K);

I = imread('eight.tif');
J = imnoise(I, 'salt & pepper', 0.09);
K = medfilt2(J);
figure, imshow(I); figure, imshow(J); figure, imshow(K);
figure, imshow(I-K); figure, imshow((I-K).^(2));

RGB = imread('saturn.png');
I = rgb2gray(RGB);
I2G = imnoise(I,'gaussian', 0.02);
I2 = imnoise(I,'salt & pepper', 0.02);
I20 = imnoise(I,'salt & pepper', 0.2);
figure, imshow(RGB);
figure, imshow(I);
figure, imshow(I2G);
figure, imshow(I2);
figure, imshow(I20);

RGB = imread('saturn.png');
I = rgb2gray(RGB);
J = imnoise(I,'gaussian', 0, 0.005);
K = wiener2(J, [10, 10]);
H = fspecial('disk', 10);
blurred = imfilter(J, H, 'replicate');
figure, imshow(RGB);
figure, imshow(I);
figure, imshow(J);
figure, imshow(H);
figure, imshow(blurred);
figure, imshow(K);

RGB = imread('saturn.png');
I = rgb2gray(RGB);
J = imnoise(I,'gaussian', 0, 0.005);
K = wiener2(J, [10, 10]);
H = fspecial('disk', 10);
blurred = imfilter(J, H, 'replicate');
figure, imshow(RGB);
figure, imshow(I);
figure, imshow(J);
figure, imshow(H);
figure, imshow(blurred);
figure, imshow(K);

RGB = imread('saturn.png');
I = rgb2gray(RGB);
J = imnoise(I,'gaussian', 0, 0.005);
K = wiener2(J, [10, 10]);
H = fspecial('disk', 10);
blurred = imfilter(J, H, 'replicate');
figure, imshow(RGB);
figure, imshow(I);
figure, imshow(J);
figure, imshow(H);
figure, imshow(blurred);
figure, imshow(K);