Spectre de puissance d'une image

J'ai commencé (un petit projet) afin de calculer le spectre de puissance d'une image dans le domaine fréquentiel.

Donc, ce que j'ai jusqu'à maintenant est la suivante:

%//close all; clear all; %//not generally appreciated
img   = imread('ajw_pic.jpg','jpg'); % it is a color image
img = rgb2gray(img); %//change to gray
psd = 10*log10(abs(fftshift(fft2(img))).^2 );
figure(2); clf
mesh(psd)

Jusqu'à présent, il semble bon; je obtenir le tracé de maillage qui ressemble à des spectres, je vois dans divers travaux académiques.

Cependant, ce que je cherche, c'est un graphique de l'intrigue de ce spectre de puissance en fonction de leur fréquence, et je ne suis pas entièrement sûr de savoir comment obtenir cette fréquence de vecteur. Je pourrais faire par exemple:

N=400;        %//the image is 400 x 400
f=-N/2:N/2-1; %//possible frequencies?

mais je ne suis pas convaincu que c'est tout à fait correct que cela donne lieu à des fréquences négatives.

Je serais vraiment reconnaissant si quelqu'un pouvait me diriger dans la bonne direction pour le tracé de la fréquence journal contre le spectre de puissance.

OriginalL'auteur AJW | 2013-12-09