ff.m

Index to gray conversion in matlab

close all
clear all
clc

% Load a gray scale image
I=imread('pout.tif');
%

% Show the image
figure(1);imshow(I,[]);
%

% Take the 2D Fast Fourier Transform
FI = fft2(I);
%

% Shift the DC component to the center 
SFI = fftshift(FI);
%

% Use logarithmic compression of dynamic range
% and display the image

figure(2);imshow(log(1+abs(SFI)),[]);