Conception d'un banc de filtres Gabor dans le domaine de Fourier

Question

Je suis en train de concevoir un banc de filtres de gabor dans le domaine de Fourier en utilisant le code suivant dans Matlab:

for i=1:u 
for j = 1:v 
    theta =theta_list(j); 
    w_1_x=(centerFreq(1,i)*((m+1)/2)*cos(theta)); % size of filter and image is m-by-n 
    w_1_y=(centerFreq(2,i)*((n+1)/2)*sin(theta)); 
    w_2_x=(-centerFreq(1,i)*((m+1)/2)*cos(theta)); 
    w_2_y=(-centerFreq(2,i)*((n+1)/2)*sin(theta)); 
    sigma=sqrt((w_1_x)^2+(w_1_y)^2); 
    sigma_H=sigma/(2*pi); 
    gF=fspecial('gaussian', [m,n], sigma_H); 
    gFilter1=circshift(gF, [round(w_1_x), round(w_1_y)]); 
    gFilter2=circshift(gF, [round(w_2_x), round(w_2_y)]); 
    gaborArray1{i,j} = 0.5*gFilter1+0.5*gFilter2; 
    gaborArray2{i,j} = 0.5*1i*gFilter1-0.5*1i*gFilter2; 
    gaborFilterBank{i,j}=gaborArray1{i,j}+1i.*gaborArray2{i,j}; 
end 
end 

Cependant, quand je transformerai le filtre au domaine de l'image en utilisant le code suivant et visualiser la carte de phase du filtre:

filter=ifft2(ifftshift(gaborFilterBank{i,j})); 
phase=zeros(m, n); 
for p=1:m 
for q=1:n 
    phase(p,q)=atan(imag(filter(p,q)/real(filter(p,q))); 
    end 
end 
figure, imshow(phase, []); 

j'obtenir des motifs étranges au centre du filtre (comme représenté dans la figure ci-jointe). Voulez-vous s'il vous plaît laissez-moi savoir ce qui ne va pas avec mon code?

Answer 1

Je trouve le problème. Au lieu de multiplier les coordonnées du centre par [(m + 1)/2, (n + 1)/2], j'aurais dû les multiplier par [m, n]. Plus précisément, au lieu de:

w_1_x=(centerFreq(1,i)*((m+1)/2)*cos(theta)); % size of filter and image is m-by-n 
w_1_y=(centerFreq(2,i)*((n+1)/2)*sin(theta)); 
w_2_x=(-centerFreq(1,i)*((m+1)/2)*cos(theta)); 
w_2_y=(-centerFreq(2,i)*((n+1)/2)*sin(theta)); 

Il devrait être:

w_1_x=(centerFreq(1,i)*m*cos(theta)); % size of filter and image is m-by-n 
w_1_y=(centerFreq(2,i)*n*sin(theta)); 
w_2_x=(-centerFreq(1,i)*m*cos(theta)); 
w_2_y=(-centerFreq(2,i)*n*sin(theta));