如何在matlab中显示40 gabor滤波器

问题描述：

有人可以帮助我如何在matlab中显示gabor滤波器，我可以显示它但它不是我想要的。这是我的代码：如何在matlab中显示40 gabor滤波器

[Gf,gabout] = gaborfilter1(B,sx,sy,f,theta(j)); 


G{m,n,i,j} = Gf;

，这是Gabor滤波器类：

function [Gf,gabout] = gaborfilter(I,Sx,Sy,f,theta); 

if isa(I,'double')~=1 
    I = double(I); 
end 

for x = -fix(Sx):fix(Sx) 
    for y = -fix(Sy):fix(Sy) 
     xPrime = x * cos(theta) + y * sin(theta); 
     yPrime = y * cos(theta) - x * sin(theta); 
     Gf(fix(Sx)+x+1,fix(Sy)+y+1) = exp(-.5*((xPrime/Sx)^2+(yPrime/Sy)^2))*cos(2*pi*f*xPrime); 
    end 
end 

Imgabout = conv2(I,double(imag(Gf)),'same'); 
Regabout = conv2(I,double(real(Gf)),'same'); 

gabout = sqrt(Imgabout.*Imgabout + Regabout.*Regabout);

然后，我imshow与此代码：

imshow(G{m,n,i,j},[]);

，结果：

enter image description here

但我想要这个结果，有人可以帮我解决这个问题吗？

gabor filter

我认为您需要发布更多代码才能提供帮助。例如，您显示'gaborfilter'功能，但不显示'gaborfilter1'功能。他们是一样的吗？我假设你正在循环中运行'gaborfilter1'，你可以发布该代码吗？我也怀疑你的Gabor滤波器的公式是错误的。西格玛在哪里？你打算做什么？ – Tokkot 2013-10-03 17:04:46

答

使用下面的功能。我希望这是有用的。
------------------------------------------------ ----------------

gfs = GaborFilter(51,0.45,0.05,6,4); 
n=0; 
for s=1:6 
    for d=1:4 
     n=n+1; 
     subplot(6,4,n) 
     imshow(real(squeeze(gfs(s,d,:,:))),[]) 
    end 
end

Sample Image
------------------------ ----------------------------------------

function gfs = GaborFilter(winLen,uh,ul,S,D) 

% gfs(SCALE, DIRECTION, :, :) 

winLen = winLen + mod(winLen, 2) -1; 
x0 = (winLen + 1)/2; 
y0 = x0; 

if S==1 
    a = 1; 
    su = uh/sqrt(log(4)); 
    sv = su; 
else 
    a = (uh/ul)^(1/(S-1)); 
    su = (a-1)*uh/((a+1)*sqrt(log(4))); 
    if D==1 
     tang = 1; 
    else 
     tang = tan(pi/(2*D)); 
    end 
    sv = tang * (uh - log(4)*su^2/uh)/sqrt(log(4) - (log(4)*su/uh)^2); 
end 

sx = 1/(2*pi*su); 
sy = 1/(2*pi*sv); 
coef = 1/(2*pi*sx*sy); 

gfs = zeros(S, D, winLen, winLen); 

for d = 1:D 
    theta = (d-1)*pi/D; 
    for s = 1:S 
     scale = a^(-(s-1)); 
     gab = zeros(winLen); 
     for x = 1:winLen 
      for y = 1:winLen 
       X = scale * ((x-x0)*cos(theta) + (y-y0)*sin(theta)); 
       Y = scale * (-(x-x0)*sin(theta) + (y-y0)*cos(theta)); 
       gab(x, y) = -0.5 * ((X/sx).^2 + (Y/sy).^2) + (2*pi*1j*uh)*X ; 
      end 
     end 
     gfs(s, d, :, :) = scale * coef * exp(gab); 
    end 
end

答

将“ cos“分量 - >complex(0, (2*pi*f*xprime)) ans还将公式乘以比例因子(1/sqrt(2*Sy*Sx))。

如何在matlab中显示40 gabor滤波器

相关推荐