close all;
clear all;
clc;
load gindb   %取點的數(shù)據(jù)庫（每圖40點）
load imgdb  %圖片數(shù)據(jù)庫（20張）

% Creating 40 Gabor Filters
%G = cell(5,8);
Kmax = pi/2   
f = sqrt(2)
sig = 2*pi
m=1:20  %循環(huán)運算20張圖片
n=1:40  %循環(huán)運算每張圖片上的40個點
x=B(n,1,m)
y=B(n,2,m)  %x與y為取的點，也用循環(huán)實現(xiàn)全部運算
%Z(:,:,n,m)=[x;y]
for s = 0:4
for t = 0:7
K = Kmax/f^s * exp(i*t*pi/8);
Kreal = real(K);
Kimag = imag(K);
NK = Kreal^2+Kimag^2;
G=(sig^(-2)*K^2)*exp((-1)*NK*(x^2+y^2)/(2*sig^2))*(exp(i*[Kreal Kimag]*[x;y])-exp(-(sig^2)/2));
W=zeros(1,1)
W(1)=G
%W(2)=imag(G{t,s,n,m})
C=zeros(1,1)
C(1)=A(round(x),round(y)) 
R=conv(C,W)  
%用J存儲幅值,以備計算
%J{:,:,n,m}=sqrt(real(R)^2+imag(R)^2)  %用一個四維細胞數(shù)組存儲計算結(jié)果
end
end
save gabor J




Kmax = pi/2   
f = sqrt(2)
sig = 2*pi
%Z(:,:,n,m)=[x;y]
s = 0:4
t = 0:7
K = Kmax/f^s * exp(i*t*pi/8);