filtracja

clear all;
close all;

img = imread('lena.bmp');

img_f= fft2(img)
img_s = fftshift(img_f)
img_abs = abs(img_s)
img_abs = log10(img_abs+1)
img_faza = angle(img_f.*(img_abs>0.0001))

figure(1)

subplot(1,3,1)
imshow(img,[])

subplot(1,3,2)
imshow(img_abs,[])

subplot(1,3,3)
imshow(img_faza,[])
%------------------------------------
[f1,f2] = freqspace(512,'meshgrid');
Hd = ones(512);
r = sqrt(f1.^2 + f2.^2);
Hd((r>0.1)) = 0;

colormap(jet(64))

img=img_f.*Hd;
img_fi = ifft2(img);
img_i = ifftshift(img_fi);

figure(2)
subplot(1,3,1)
imshow(img,[])

subplot(1,3,2)
imshow(img_i,[])
%-----------------------------------
[f1,f2] = freqspace(512,'meshgrid');
Hd = ones(512);
r = sqrt(f1.^2 + f2.^2);
Hd((r<0.1)) = 0;

colormap(jet(64))

img=img_f.*Hd;
img_fi = ifft2(img);
img_i = ifftshift(img_fi);

subplot(1,3,3)
imshow(img_i,[])
%------------------------------------
[f1,f2] = freqspace(512,'meshgrid');
Hd = ones(512);
r = sqrt(f1.^2 + f2.^2);
Hd((r<0.1)) = 0;

colormap(jet(64))

img_i = ifftshift(img);
img_fi = ifft2(img_i);

subplot(1,3,3)
imshow(img_fi,[])
%------------------------------------
h = fwind1(Hd,hamming(21));
[H f1 f2] = freqz2(h,512,512);
mesh(f1,f2,H);