matlab实现gabor filter (2)

这段gabor filter生成了5个尺度，8个方向的gabor filter 滤波器组，生成的滤波器组还是挺漂亮的：

 

 

function [Gs] = mygabor;
% 生成的gabor filter 还是不错的
% Img = imread('1.tiff');

sigma = 2*pi;
sigma2 = sigma^2;

GaborZ = 71;

n=1;
figure;
for v=0:1:4
    for u=0:1:7
        j = u+8*v;
    n
        Kv = pi*2^(-(v+2)/2);
        faiu = pi * u/8;
%         Kj = Kv * exp( i * faiu );
        Kj = [Kv *cos(faiu) Kv *sin(faiu)];
        K2 = norm(Kj');
        K2 = K2.^2;
        Gab1 = (K2 /(sigma2));
        for zx = -GaborZ:GaborZ-1
            for zy = -GaborZ:GaborZ-1
                x = [zx zy];
                x=x';
                Gab2 = exp(-K2 * (zx^2 + zy^2)/(2*sigma2));
                Gab3 = (exp(i * Kj * x) - exp(-(sigma2)/2));
                Gr(zx+GaborZ+1,zy+GaborZ+1) = real(Gab1 * Gab2 * Gab3);
%                 Gab(zx+GaborZ+1,zy+GaborZ+1) = norm(Gab1 * Gab2 * Gab3);
            end
        end
        subplot(5,8,n),imshow(Gr,[]);
%         figure;imshow(Gab,[]);
        n=n+1; 
    end
%     figure;
end

生成的滤波器组图片：