Improved 2D Gabor Filter

Bi-dimensional Gabor filter with DC component compensation.

This version of the 2D Gabor filter [1] is basically a bi-dimensional Gaussian function centered at origin (0,0) with variance S modulated by a complex sinusoid with polar frequency (F,W) and phase P described by the following equation:

G(x,y,S,F,W,P)=k*Gaussian(x,y,S)*(Sinusoid(x,y,F,W,P)-DC(F,S,P)),
where:
Gaussian(x,y,S)=exp(-pi*S^2*(x^2+y^2))
Sinusoid(x,y,F,W,P)=exp(j*(2*pi*F*(x*cos(W)+y*sin(W))+P)))
DC(F,S,P)=exp(-pi*(F/S)^2+j*P)

The following reference will be helpful to understand how setting the Gabor parameters:
[1] Movellan, J. R. - Tutorial on Gabor Filters. Tech. rep., 2002.