% demonstration file for seminar  16/Nov/00
load clown
x=ind2gray(X,map);
imshow(x,64), title('original Image');
disp('press a key')
pause
[nr,nc]=size(x);
f2=fftshift(fft2(x));
imagesc(log(abs(f2)));colorbar,colormap(jet), title('Fourier transform of Clown image');
disp('press a key')
pause
[x1,y1,z1]=peaks(20);
surf(x1,y1,z1,del2(z1));
colormap(jet(64)); colorbar('horiz');
disp('press a key')
pause
colormap(gray);
f=ifft2(f2);
imshow(f,64), title('Inverse Fourier transform of Clown image');
disp('press a key')
pause