这段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)