FFTAnalysis

clf;

clear all;

b = imread('image_name');

c = cast(b,'int8');

[x,y]=size(c);

for k = 1:x

for l = 1:y

if (c(k,l) > 100)

c(k,l)=1;

else

c(k,l)=0;

endif

endfor

endfor

a = fftshift(fft2(c(:,:,1)));

[x1,y1]=size(a);

surf(1:y1,1:x1,abs(a),'edgecolor','none')