how to apply low pass and high pass filters using fft2 and fftshift?

19 vues (au cours des 30 derniers jours)
Hadi Ghahremannezhad
Hadi Ghahremannezhad le 23 Sep 2020
I am trying to apply low-pass and high-pass filters on a jpg image.
The filters are calculated by:
function [cL,cH] = getfilters(radius)
[x,y] = meshgrid(-128:127,-128:127);
z = sqrt(x.^2+y.^2);
cL = z < radius;
cH = ~cL;
end
so to filter the image, I have:
img = imread('a jpg image')
ft = fftshift(fft2(img));
[cL, cH] = getfilters(40);
l_ft = ft .* cL;
h_ft = ft .* cH;
low_filtered_image = ifft2(ifftshift(l_ft), 'symmetric');
high_filtered_image = ifft2(ifftshift(h_ft), 'symmetric');
To display the filtered inverse Fourier transform:
low_f = uint8(abs(low_filtered_image));
high_f = uint8(abs(high_filtered_image));
figure, imshow(low_f); figure, imshow(high_f);
Which parts are wrong in this code?

Réponses (1)

Soroush Vahidi
Soroush Vahidi le 11 Fév 2022
It is exactly my question too.

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by