【图像处理】Gabor滤波器

Gabor的核函数参考的wiki

使用实数Real的公式计算核函数代码：

Mat getGaborFilter(float lambda, float theta, 
	float sigma2,float gamma, 
	float psi = 0.0f){
		if(abs(lambda-0.0f)<1e-6){
			lambda = 1.0f;
		} 
		float sigma_x = sigma2;
		float sigma_y = sigma2/(gamma*gamma);
		int nstds = 3;
		float sqrt_sigma_x = sqrt(sigma_x);
		float sqrt_sigma_y = sqrt(sigma_y);
		int xmax = max(abs(nstds*sqrt_sigma_x*cos(theta)),abs(nstds*sqrt_sigma_y*sin(theta)));
		int ymax = max(abs(nstds*sqrt_sigma_x*sin(theta)),abs(nstds*sqrt_sigma_y*cos(theta)));
		int half_filter_size = xmax>ymax ? xmax:ymax;
		int filter_size = 2*half_filter_size+1;
		Mat gaber = Mat::zeros(filter_size,filter_size,CV_32F);
		for(int i=0;i<