I = imread('cameraman.tif');
[G,gabout] = gaborfilter1(I,2,4,16,pi/3);
figure,imshow(uint8(gabout));