Gaborfilter matlab

Gabor FILTER MATLAB IMAGE PROCESSING

% GABORFILTER Bi-dimensional Gabor filter with DC component compensation.
%   [GGABOUT]=GABORFILTER（ISFWP） filters the input image I with the 2D
%   Gabor filter G described by the parameters S F W and P to create the 
%   output filtered image GABOUT.
%   This version of the 2D Gabor filter is basically a bi-dimensional 
%   Gaussian function centered at origin （00） with variance S modulated by
%   a complex sinusoid with polar frequency （FW） and phase P described by 
%   the following equation:
%
%   G（xySFWP）=k*Gaussian（xyS）*（Sinusoid（xyFWP）-DC（FSP））
%   where:                        
%       Gaussian（xyS）=exp（-pi*S^2*（x^2+y^2））
%       Sinusoid（xyFWP）=exp（j*（2*pi*F*（x*cos（W）+y*sin（W））+P）））
%       DC（FSP）=exp（-pi*（F/S）^2+j*P）
%
%   PS: The term DC（FSP） compensates the inherent DC component produced 
%   by the Gaussian envelop as shown by Movellan in [1].
%
%	Tips:
%       1） To get the real part and the imaginary part of the complex 
%       filter output use real（gabout） and imag（gabout） respectively;
% 
%       2） To get the magnitude and the phase of the complex filter output 
%       use abs（gabout） and angle（gabout） respectively.

%	Author: Stiven Schwanz Dias e-mail: stivendias@gmail.com
%           Cognition Science Group Informatic Department
%           University of Esp韗ito Santo Brazil January 2007.
%
%   References:
%       [1] Movellan J. R. - Tutorial on Gabor Filters. Tech. rep. 2002.

function [GGABOUT]=gaborfilter（ISFWP）;

if isa（I‘double‘）~=1
    I=double（I）;
end
size=fix（1.5/S）; % exp（-1.5^2*pi） < 0.1%
%k=2*pi*S^2;
%F=S^2/sqrt（2*pi）;
k=1;
for x=-size:size
    for y=-size:size
        G（size+x+1size+y+1）=k*exp（-pi*S^2*（x*x+y*y））*...
            （exp（j*（2*pi*F*（x*cos（W）+y*sin（W））+P））-exp（-pi*（F/S）^2+j*P））;
    end
end

GABOUT=conv2（Idouble（G）‘same‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        2093  2006-12-30 02:14  colormaps.mat
     文件        1903  2007-01-02 19:33  gaborfilter.m
     文件         857  2008-05-31 16:15  gabortest.m
     文件       69214  2007-05-22 20:18  lena.jpg
     文件        1340  2009-05-19 23:41  license.txt