gaborfilter函数源代码里面有三个程序

gaborfilter函数源代码（里面有三个程序），对研究gaborfilter有很大帮助！

%%%%%%%VERSION 1

%The Gabor filter is basically a Gaussian （with variances sx and sy along x and y-axes respectively）
%modulated by a complex sinusoid （with centre frequencies U and V along x and y-axes respectively） 
%described by the following equation
%%
%               1                -1     x  ^    y  ^
%%% G（xy） = ---------- * exp （[----{（----） 2+（----） 2}+2*pi*i*（Ux+Vy）]）
%            2*pi*sx*sy           2    sx       sy

%% Describtion :

%% I : Input image
%% Sx & Sy : Variances along x and y-axes respectively
%% U & V : Centre frequencies  along x and y-axes respectively

%% G : The output filter as described above
%% gabout : The output filtered image

%%  Author : Ahmad poursaberi  e-mail : a.poursaberi@ece.ut.ac.ir
%%          Faulty of Engineering Electrical&Computer DepartmentTehran
%%          UniversityIranJune 2004

function [Ggabout] = gaborfilter（ISxSyUV）;

if isa（I‘double‘）~=1 
    I = double（I）;
end

for x = -fix（Sx）:fix（Sx）
    for y = -fix（Sy）:fix（Sy）
        G（fix（Sx）+x+1fix（Sy）+y+1） = （1/（2*pi*Sx*Sy））*exp（-.5*（（x/Sx）^2+（y/Sy）^2）+2*pi*i*（U*x+V*y））;
    end
end

Imgabout = conv2（Idouble（imag（G））‘same‘）;
Regabout = conv2（Idouble（real（G））‘same‘）;

%gabout = uint8（sqrt（Imgabout.*Imgabout + Regabout.*Regabout））;
gabout = sqrt（Imgabout.*Imgabout + Regabout.*Regabout）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1860  2004-06-14 09:39  gaborfilter2.m

     文件       1365  2010-03-28 22:54  gaborfilter.m

     文件       1553  2010-03-29 14:58  gaborfilter1.m

----------- ---------  ---------- -----  ----

                 4778                    3