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
生成的滤波器组图片：
