%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]范圍內(nèi)
    (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]范圍內(nèi)
    ./(max(max(f_img_angle))-min(min(f_img_angle)))*255;
imshow(f_img_angle);                              %顯示頻譜相位
title('圖像的相頻分布');