Изображение текстуры с использованием фильтра Габора - PullRequest
Новая
       24

Изображение текстуры с использованием фильтра Габора

4 голосов
/ 16 сентября 2010

У меня есть следующий фильтр Габора для извлечения функции текстуры изображения .. 

a=imread('image0001.jpg');
a=double(a); 
a=a-mean(a(:)); 
[r,c,l]=size(a); 

 K=5;  S=6;
 Uh=0.4;
 Ul=0.05;
 alpha=(Uh/Ul)^(1/(S-1));  
 sigmau=(alpha-1)*Uh/((alpha+1)*sqrt(2*log(2))); 
 sigmav=tan(pi/(2*K))*(Uh-2*log(2)*((sigmau^2)/Uh))/sqrt((2*log(2))-(((2*log(2))^2)*(sigmau^2)/(Uh^2)));
 sigmax=1/(2*pi*sigmau);
 sigmay=1/(2*pi*sigmav);
 b=fft2(a);
 [e d]=size(b);
  i=1;
 G=zeros(r,c,S*K);
 IZ=zeros(r,c,S*K);
for m=0:S-1
    for n=0:K-1 
       fprintf(1,'.');
       for x=-r/2+1:r/2;
          for y=-c/2+1:c/2;
             xdash=(alpha^(-m))*((x)*cos(n*pi/K)+(y)*sin(n*pi/K));
             ydash=(alpha^(-m))*((y)*cos(n*pi/K)-(x)*sin(n*pi/K));
             g(r/2+x,r/2+y)=(alpha^(-m))*((1/(2*pi*sigmax*sigmay))*exp(-0.5*(((xdash^2)/(sigmax^2))+((ydash^2)/(sigmay^2)))+0.8i*pi*xdash));
          end
       end
       [rr cc]=size(g); 
       G(:,:,i)=g;
       h=fft2(g); 
       z=b.*h;
       iz=ifft2(z);
       IZ(:,:,i)=iz;
       FeatureVector(i)=mean(abs(iz(:)));
       i=i+1;
    end
 end
 fprintf(1,'\n');

%%%%%%%%%

Когда я запускаю этот код, я получаю эту ошибку:

Ошибка при использовании ==> разРазмеры матрицы должны совпадать.Ошибка в ==> ComputeGaborFeatures4 при 37 z = b. * H;

Пожалуйста, если кто-нибудь может помочь мне решить эту ошибку, или любой может дать мне еще один простой фильтр Габора?

Ответы [ 2 ]

0 голосов
/ 27 февраля 2013

ошибка может изменить g (r / 2 + x, r / 2 + y) на g (r / 2 + x, c / 2 + y), значение

0 голосов
/ 10 февраля 2012

Ошибка связана с вызовом умножения массива (. *) С b и h неравного размера, поскольку rr не равно r, а cc не равно c.

Либо вы хотите использовать умножение матриц (*), либо вам нужно сделать g и такой же размер перед вызовом fft2.

...