数字图像处理之一傅立叶变换

clear all;
imgrgb = imread('D:/car.jpg');
imshow( imgrgb);
X = rgb2gray(imgrgb);
A=fft2(double(X));        %傅里叶变换
ifcoef=ifft2(A);
temp1=log(1+abs(A));
spectrum=fftshift(A);      %使图像对称
temp2=log(1+abs(spectrum));
figure
subplot(2,2,1) ,imshow(X),title('Source image');
subplot(2,2,2) ,imshow(temp1,[]),title('FFT image');
subplot(2,2,3) ,imshow(temp2,[]),title('  Shift FFT image');
subplot(2,2,4), imshow(ifcoef,[]),title('  IFFT image');