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

【问题标题】：How to show 40 gabor filter in matlab如何在matlab中显示40个gabor滤波器
【发布时间】：2013-04-27 17:38:40
【问题描述】：

谁能帮助我如何在 matlab 中显示 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(G{m,n,i,j},[]);

结果：

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

【问题讨论】：

我认为您需要发布更多代码以供我们帮助。例如，您显示了gaborfilter 函数，但没有显示gaborfilter1 函数。他们是一样的吗？我假设您在循环中运行gaborfilter1，您可以发布该代码吗？我还怀疑您的 Gabor 过滤器公式是错误的。西格玛在哪里？你打算f做什么？

标签： matlab image-processing filter

【解决方案1】：

使用以下函数。我希望这是有用的。
-------------------------------------------------- --------------

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

【讨论】：

【解决方案2】：

用复数部分替换“cos”组件->complex(0, (2*pi*f*xprime)) 还可以将方程乘以(1/sqrt(2*Sy*Sx)) 的比例因子。

【讨论】：