www.pudn.com > GaitRecognition.rar > MakeGaborKernal.m, change:2010-07-08,size:1037b


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); 
        
        GaborReal(j+HarfH+1, i+HarfW+1) = Kv*Kv*tmp1*tmp2/sqsigma; 
        GaborImg(j+HarfH+1, i+HarfW+1) = Kv*Kv*tmp1*tmp3/sqsigma; 
    end 
end