基于MATLAB中fft2函数的二维图像数据的频谱生成代码

如何利用MATLAB中的fft2函数对二维的图像数据进行频域变换

clc

clear all

close all

I = imread('testwomen.jpg');    %读取源图像

size(I) %显示图像采样值

figure(1);

imshow(I);

title('the original photo')

%---------------------------

I1 = I(:,:,1);

fftI = fft2(I1);

A1 = abs(fftI);         %取模值

B1 = (A1 - min(min(A1))) / (max(max(A1)) - min(min(A1))) * 255;  %归一化显示

figure(2)

imshow(B1);

title('the frequency spectrum without shiftfft()')

%---------------------------

sfftI = fftshift(fftI);

A = abs(sfftI);

B = (A - min(min(A))) / (max(max(A)) - min(min(A))) * 255;

figure(3)

%imshow(B,[0,255]);  

imshow(B);
title('the frequency spectrum with shiftfft()')







黑白可以用数值表示，即信号值的大小。上图是二维图像信号的频谱，横轴代表x方向的频率，纵轴代表y方向的频率，黑白表示不同频率分量振幅的强弱。

频谱的中心代表了低频信号的振幅，频谱远离中心位置的灰度代表了高频成分振幅的大小。