%画像読み込み
imagedata=imread('test.jpg');
figure;imagesc(imagedata);colormap(gray(256)); %画像表示

%フーリエ変換
fftdata = fft2(imagedata);

%振幅
fftdata_abs = uint8(fftshift(abs(fftdata))/10);
figure;imagesc(fftdata_abs);colormap(gray(256)); %画像表示

%位相
fftdata_ang = abs(ifft2(double(angle(fftdata))));
figure;imagesc(fftdata_ang);colormap(gray(256)); %画像表示

%逆フーリエ変換
ifftdata = ifft2(fftdata);
figure;imagesc(ifftdata);colormap(gray(256)); %画像表示
