fft2_example.m

matlab7.0从入门到精通(光盘内容)

%fft2_example.m
%分析图像的频谱
img = imread('coins.png');                        %读图像文件
f_img = fft2(double(img));                                %2维fft变换
f_img = fftshift(f_img);                          %把直流成分移动频谱的中间
imshow(img);                                      %显示图像
figure;
f_img_abs = abs(f_img);                           %得到频谱幅度
f_img_abs = (f_img_abs-min(min(f_img_abs)))./ ... %把频谱幅度变换到[0 255]范围内
    (max(max(f_img_abs))-min(min(f_img_abs)))*255;
imshow(f_img_abs);                                %显示频谱幅度
title('图像的幅频分布');
figure;
f_img_angle = angle(f_img);                       %得到频谱的相位
f_img_angle = (f_img_angle-min(min(f_img_angle)))... %把频谱相位变换到[0 255]范围内
    ./(max(max(f_img_angle))-min(min(f_img_angle)))*255;
imshow(f_img_angle);                              %显示频谱相位
title('图像的相频分布');