在MATLAB下的gist特征提取及傅里叶变换

内含完全代码以及函数，测试有效，拥有完整备注

% function G = createGabor（or n）

% G = createGabor（numberOfOrientationsPerScale n）;
%
% Precomputes filter transfer functions. All computations are done on the
% Fourier domain. 
%
% If you call this function without output arguments it will show the
% tiling of the Fourier domain.
%
% Input
%     numberOfOrientationsPerScale = vector that contains the number of
%                                orientations at each scale （from HF to BF）
%     n = imagesize （square images）
%
% output
%     G = transfer functions for a jet of gabor filters

or = [8 8 4];
n = 128;

Nscales = length（or）;
Nfilters = sum（or）;

l=0;
for i=1:Nscales
    for j=1:or（i）
        l=l+1;
        param（l:）=[.35 .3/（1.85^（i-1）） 16*or（i）^2/32^2 pi/（or（i））*（j-1）];
    end
end

% Frequencies:
[fx fy] = meshgrid（-n/2:n/2-1）;
fr = fftshift（sqrt（fx.^2+fy.^2））;
t = fftshift（angle（fx+sqrt（-1）*fy））;

% Transfer functions:
G=zeros（[n n Nfilters]）;
for i=1:Nfilters
    par=param（i:）;
    tr=t+param（i4）; 
    tr=tr+2*pi*（tr<-pi）-2*pi*（tr>pi）;

    G（::i）=exp（-10*param（i1）*（fr/n/param（i2）-1）.^2-2*param（i3）*pi*tr.^2）;
end


if nargout == 0
    figure
    for i=1:Nfilters
        max（max（G（::i）））
        contour（fftshift（G（::i））[1 .7 .6]‘r‘）;
        hold on
        drawnow
    end
    axis（‘on‘）
    axis（‘square‘）
    axis（‘ij‘）
   
end
G


 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2017-03-20 19:31  Gist\
     文件        1399  2017-03-18 15:02  Gist\createGabor.m
     文件         734  2017-03-18 13:37  Gist\demo.m
     文件         392  2010-09-25 15:27  Gist\downN.m
     文件        1456  2017-03-18 14:21  Gist\gistGabor.asv
     文件        1578  2017-03-18 14:35  Gist\gistGabor.m
     文件        1050  2007-02-10 17:07  Gist\prefilt.m
     文件       14278  2001-06-24 16:32  Gist\test.jpg