function fout = gaborFilter(gamma,theta,lambda,sigma)
    sigma_x = sigma;
    sigma_y = sigma/gamma;

    sz=fix(4*max(sigma_y,sigma_x));

    [x y]=meshgrid(-sz:sz,sz:-1:-sz);

    % Rotation 
    x_theta=x*cos(theta)+y*sin(theta);
    y_theta=-x*sin(theta)+y*cos(theta);

    fout=filterNorm(exp(-0.5*(x_theta.^2/sigma_x^2+y_theta.^2/sigma_y^2)).*exp(2*i*pi/lambda*x_theta));