gabor 滤波器 函数:
function [GaborReal, GaborImg] = MakeGaborKernal(GaborH, GaborW, U, V, Kmax,f,sigma )
% function [GaborReal, GaborImg] = MakeGaborKernal[GaborH, GaborW, U, V]
% 用以生成 Gabor 核
% GaborReal: 核实部 GaborImg: 虚部
% GaborH,GaborW: Gabor窗口 高宽.
% U,V: 方向 大小
% ||Ku,v||^2
% G(Z) = ---------------- exp(-||Ku,v||^2 * Z^2)/(2*sigma*sigma)(exp(i*Ku,v*Z)-exp(-sigma*sigma/2))
% sigma*sigma
HarfH = fix(GaborH/2);
HarfW = fix(GaborW/2);
Qu = pi*U/8;
sqsigma = sigma*sigma;
% Kv = 2.5*pi*(2^(-(V+2)/2));
Kv = Kmax/(f^V);
postmean = exp(-sqsigma/2);
for j = -HarfH : HarfH
for i = -HarfW : HarfW
tmp1 = exp(-(Kv*Kv*(j*j+i*i)/(2*sqsigma)));
tmp2 = cos(Kv*cos(Qu)*i+Kv*sin(Qu)*j) - postmean;
% tmp3 = sin(Kv*cos(Qu)*i+Kv*sin(Qu)*j) - exp(-sqsigma/2);
tmp3 = sin(Kv*cos(Qu)*i+Kv*sin(Qu)*j);