function [G,gabout] = gaborfilter2(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); G(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(G)),'same'); Regabout = conv2(I,double(real(G)),'same');
gabout = sqrt(Imgabout.*Imgabout + Regabout.*Regabout);
0 Comments
Sign in to comment.