圖像二維離散傅里葉變換、幅度譜、相位譜

clear, clc
I = imread('...');F = fftshift(fft2(I));                     % 對圖像進行二維 DFT(fft2)，并移至中心位置
magn = log(abs(F));                        % 加 log 是便于顯示，縮小值域
phase = log(angle(F)*180/pi);              % 轉換為度數subplot(2, 2, 1), imshow(I, []), title('原圖像')
subplot(2, 2, 2), imshow(log(F), []), title('圖像頻譜')
subplot(2, 2, 3), imshow(magn, []), title('圖像幅度譜')
subplot(2, 2, 4), imshow(phase, []), title('圖像相位譜')

• abs(F) == sqrt(real(F).^2 + imag(F).^2)
• angle(F) == log(imag(F)) == atan2(imag(F), real(F))