周波数領域で画像のパワー スペクトルを計算する (小さなプロジェクト) を開始しました。
だから、私が今まで持っているものは次のとおりです。
%// 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)
ここまでは良さそうです。さまざまな学術論文で見られるスペクトルに似たメッシュ プロットが得られます。
ただし、私が探しているのは、このパワー スペクトルと周波数のグラフ プロットであり、この周波数ベクトルを取得する方法が完全にはわかりません。たとえば、次のことができます。
N=400; %// the image is 400 x 400
f=-N/2:N/2-1; %// possible frequencies?
しかし、これが負の周波数を発生させるため、これが完全に正しいとは確信していません.
ログ周波数とパワースペクトルをプロットする正しい方向に誰かが私を向けることができれば、本当に感謝しています。