Filtracja

close all;
clear all;

lena = imread('lena3.bmp');
lenaF = fft2(lena);
lenaO = fftshift(lenaF);
lenaAb = abs(lenaO);
lenaAb = log10(lenaAb+1);
lenaFz = angle(lenaF.*(lenaAb > 0.0001));

figure(1);
subplot(1,3,1)
imshow(lena,[]);
title('Oryginał');
subplot(1,3,2);
imshow(lenaAb,[]);
title('Amplituda');
subplot(1,3,3);
imshow(lenaFz,[]);
title('Faza');

figure(2);
[f1,f2] = freqspace(512, 'meshgrid');
Hd = ones(512);
r = sqrt(f1.^2 + f2.^2);
Hd((r>0.1)) = 0;
Hd2 = ones(512);
Hd2((r<0.1)) = 0;
Hd3 = ones(512);
Hd3(r>0.1 & r<0.4) = 0;
colormap(jet(64));
mesh(f1,f2,Hd);

lenaMn = lenaO.* Hd;
lenaO2 = ifftshift(lenaMn);
lenaF2 = ifft2(lenaO2);

lenaMn2 = lenaO.* Hd2;
lenaO3 = ifftshift(lenaMn2);
lenaF3 = ifft2(lenaO3);

lenaMn3 = lenaO.* Hd3;
lenaO4 = ifftshift(lenaMn3);
lenaF4 = ifft2(lenaO4);
%rodzielam wyswietlanie filtrów i obrazów ze wzgledu na słaba widocznosc
%róznic w obrazach w przypadku gdy są pokazane w małej rozdzielczosci.
figure(3);
subplot(1,3,1);
imshow(lenaMn);
title('Filtr dolnoprzepustowy');
subplot(1,3,2);
imshow(lenaMn2);
title('Filtr górnoprzepustowy');
subplot(1,3,3);
imshow(lenaMn3);
title('Filtr pasmowo przepustowy');

figure(4);
subplot(1,3,1);
imshow(lenaF2,[]);
title('r > 0.1');
subplot(1,3,2);
imshow(lenaF3,[]);
title('r < 0.1');
subplot(1,3,3);
imshow(lenaF4,[]);
title('r > 0.1 i r < 0.4');

h = fwind1(Hd, hamming(21));
h2 = fwind1(Hd2, hamming(21));
[H1 f1 f2] = freqz2(h,512,512);
[H2 f1 f2] = freqz2(h2,512,512);
figure(5);
mesh(f1,f2,H1);

lenaMn4 = lenaO.* H1;
lenaO5 = ifftshift(lenaMn4);
lenaF5 = ifft2(lenaO5);

lenaMn5 = lenaO.* H2;
lenaO6 = ifftshift(lenaMn5);
lenaF6 = ifft2(lenaO6);

figure(6);
subplot(1,2,1);
imshow(lenaF5,[]);
title('Dolnoprzepustowa');
subplot(1,2,2);
imshow(lenaF6,[]);
title('Górnoprzepustowa');