# How to get location of central frequency and poles in FFT

4 vues (au cours des 30 derniers jours)
Anurag Gaur le 16 Avr 2022
Commenté : Anurag Gaur le 21 Avr 2022
Thanks in advance to help me out,
I have an image, which transformed into frequency domain using 'Fast Fourier Transformation (fft2)' and then processed using a low pas filter 'gaussian' (part of my code is as follows)
image=imread('FFT_image.jpg'); % my_image is only taken as example (RGB)
subcrop=rgb2gray(imcrop(image,[560 87 500 480])); % region of interest has been selected and converted into gray scale
FFT=fft2(subcrop); % image transformed into frequency domain using in-built command fft2 (for 2D)
FFTS= fftsfift(FFT); % frequency then shifted to center
FFT2=(FFTS-min(min(FFTS)))./(max(max(FFTS))).*255; % I have used this command to display the (image file attached)
figure(1)
imshow(FFT2) % image file is attached, named 'Shifted FFT frequency domain'
TS=circshift(FFTS,66); % shifting the location of the pole (first order) to central frequency (zero order) manually
TS2= (TS-min(min(TS)))./(max(max(TS))).*255;
figure (2)
imshow(TS2) % image is attached 'Poles_shifted_FFT'
Now, how can I locate the position of central frequency (red arrow) and the poles (blue arrow) using any command or inbuilt function. I have to interchange the postion of these two and then convolute with a gaussian low pass filter. Currently I am doing this task manually using the cursor to locate these two points and the get the difference. Then using 'circshift' and put the value by the amount poles needs to be shifted at location of central frequency.
The final output of my image after using the low pass filter is a 'double complex' type. I would like to add false/pseudo color for better understanding. I tried using 'rgb2ind' but these ways are not working. Kindly help me for this also.
(both the images are attached)
##### 0 commentairesAfficher -2 commentaires plus anciensMasquer -2 commentaires plus anciens

Connectez-vous pour commenter.

### Réponses (1)

yanqi liu le 19 Avr 2022
yes，sir，may be use Butterworth band stop filter to location the bright area
##### 1 commentaireAfficher -1 commentaires plus anciensMasquer -1 commentaires plus anciens
Anurag Gaur le 21 Avr 2022
Is it possible to get in terms of X and Y?

Connectez-vous pour commenter.

### Catégories

En savoir plus sur Morphological Operations dans Help Center et File Exchange

### Community Treasure Hunt

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

Start Hunting!

Translated by