برنامه فیلتر کردن حوزه فرکانس

در برنامه زیر فیلتر کردن حوزه فرکانس تصویر ورودی به وسیله یک فیلتر پایین گذر گوسی انجام شده است. در مورد سایر فیلترهای حوزه فرکانس نیز روند کار به همین شکل است.

clearvars;
close all;
img = double(rgb2gray(imread('FIG2_3.jpg')));
figure;imshow(uint8(img));
title('input image');

[M,N] = size(img);
imgz = [img,zeros(M,N)];
imgz = [imgz;zeros(M,2*N)];
[X,Y] = meshgrid(0:2*N-1,0:2*M-1);
imgz = imgz.*(-1).^(X+Y);

IMGZ = fft2(imgz);
figure;imshow(20*log10(abs(IMGZ)),[]);
title('input image in freq. domain');

D0 = 50;
D = sqrt((X-N).^2+(Y-M).^2);
GLPF = exp(-D.^2/2/D0^2);
figure;imshow(20*log10(GLPF),[]);
title('GLPF in freq. domain');

IMGO = GLPF.*IMGZ;
figure;imshow(20*log10(abs(IMGO)),[]);
title('output image in freq. domain');

imgo = real(ifft2(IMGO));
imgo = imgo.*(-1).^(X+Y);
imgo = imgo(1:M,1:N);
figure;imshow(uint8(imgo));
title('output image');

برنامه فیلتر کردن حوزه فرکانس