圖像頻譜圖畫圖——matlab

I =imread('C:\Users\wangd\Desktop\in000155.jpg');  %讀入原圖像文件 I1 = rgb2gray(I); subplot(1,2,1);imshow(I1);          %顯示原圖像 fftI1=fft2(I1);                       %二維離散傅立葉變換 sfftI1=fftshift(fftI1);              %直流分量移到頻譜中心 RR1=real(sfftI1);                    %取傅立葉變換的實部 II1=imag(sfftI1);                    %取傅立葉變換的虛部 A1=sqrt(RR1.^2+II1.^2);             %計算頻譜幅值 A1=(A1-min(min(A1)))/(max(max(A1))-min(min(A1)))*225;%歸一化 subplot(1,2,2);imshow(A1);          %顯示原圖像的頻譜