% 画像を読み込んで表示
imagedata = imread('201107301600.jpg');
figure;imgesc(imagedata);colormap(gray(256));

% フーリエ変換
[h w] = size(imagedata);
fftdata = fft2(double(imagedata), h, w);

% 振幅
abs_ex = log(abs(fftshift(fftdata)));
fftdata_abs = abs_ex/max(abs_ex(:));
figure;imagesc(fftdata_abs);colormap(gray(256));

% 位相
fftdata_ang = abs(ifft2(double(angle(fftdata)), h, w));
figure;imagesc(fftdata_abs);colormap(gray(256));

%逆フーリエ変換
ifftdata = ifft2(double(fftdata), h, w);
figure;imagesc(ifftdata);colormap(gray(256));
