بسم الله الرحمن الرحیم

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

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');