clc;clear;close all;
%% I : Input image
I = imread('barbara.png');
%% Sx & Sy : Variances along x and y-axes respectively
Sx=1;Sy=8;
%% f : The frequency of the sinusoidal function
f=4;
%% theta : The orientation of Gabor filter
theta=-pi/4;
[G,gabout] = gaborfilter1(I,Sx,Sy,f,theta);
figure,mesh(real(G));
figure,imshow(uint8(gabout)); 
theta=0;
[G,gabout] = gaborfilter1(I,Sx,Sy,f,theta);
figure,mesh(real(G));
figure,imshow(uint8(gabout)); 
theta=pi/4;
[G,gabout] = gaborfilter1(I,Sx,Sy,f,theta);
figure,mesh(real(G));
figure,imshow(uint8(gabout)); 
%figure,imshow(uint8(real(G))); 
%figure,imshow(uint8(imag(G))); 
%figure(1);imagesc(real(G));colormap(gray);
%figure(2);imagesc(imag(G));colormap(gray);
figure,mesh(real(G));
%figure,mesh(imag(G));
figure,imshow(uint8(gabout)); 

U=2;V=8;
[G,gabout] = gaborfilter(I,Sx,Sy,U,V);
figure,mesh(real(G));
figure,mesh(imag(G));
figure,imshow(uint8(gabout)); 

[G1,G2,gabout1,gabout2] = gaborfilter2(I,Sx,Sy,f,theta);
figure,mesh(G1);
figure,mesh(G2);
figure,imshow(uint8(gabout1)); 
figure,imshow(uint8(gabout2));