clear all;
close all;
f = imread('cameraman.tif');
[M,N] = size(f);
f = imnoise(f,'salt & pepper',0.25);

imwrite(f,'w:\stuff\meinard\bilder\noisy\median_filter\original_image.tif','TIFF');

% Spektrum des Bildes berechnen und anzeigen
F = fft2(f);
Fdisplay = fftshift(abs(F));
Fdisplay = log(1+Fdisplay);
Fdisplay = Fdisplay/max(max(Fdisplay));
imwrite(Fdisplay,'w:\stuff\meinard\bilder\noisy\median_filter\original_image_log_spectrum.tif','TIFF');

add_imagecol = [padarray(f,[4 4],255,'both'); padarray(256*Fdisplay,[4 4],255,'both')];
multi = [add_imagecol];

for filtersize = [3 5 9]
str_filtersize = num2str(filtersize);

% Bild Median-filtern und abspeichern
fk = medfilt2(f,[filtersize filtersize]);
imwrite(fk,strcat('w:\stuff\meinard\bilder\noisy\median_filter\filtered_image_',str_filtersize,'.tif'),'TIFF');

% Bild mittels FFT und punktweiser Multiplikation filtern, Spektrum
% abspeichern
FKdisplay = fftshift(log(1+abs(fft(fk))));
MFK = max(max(FKdisplay));
imwrite(FKdisplay/MFK,strcat('w:\stuff\meinard\bilder\noisy\median_filter\filtered_image_log_spectrum_',str_filtersize,'.tif'),'TIFF');
add_imagecol = [padarray(fk,[4 4],255,'both'); padarray(256*FKdisplay/MFK,[4 4],255,'both')];
multi = [multi, add_imagecol];
imwrite(multi,strcat('w:\stuff\meinard\bilder\noisy\median_filter\collection.tif'),'TIFF');
end;