function?[G1gabout1]?=?gaborfilter（ISxSyftheta）;

if?isa（I‘double‘）~=1?
???I=double（I）;
end
k=6*Sx/f;
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）;
???????G1（fix（Sx）+x+1fix（Sy）+y+1）=（1/（2*pi*Sx*Sy））?*?exp（-.5*（（x/Sx）^2+（y/Sy）^2））*sin（2*pi*xPrime/k）;
???????G2（fix（Sx）+x+1fix（Sy）+y+1）=（1/（2*pi*Sx*Sy））?*?exp（-.5*（（x/Sx）^2+（y/Sy）^2））*cos（2*pi*xPrime/k）;
???end
end

GG1=conv2（Idouble（G1）‘same‘）;
GG2=conv2（Idouble（G2）‘same‘）;


gabout1=sqrt（GG1.*GG1+GG2.*GG2）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????目錄??????????0??2010-06-23?16:37??gabor

?????文件????????424??2010-05-17?21:59??gabor\gaborfilter1.m

?????文件????????576??2010-06-21?15:24??gabor\Gabor.m

-----------?---------??----------?-----??----

?????????????????1000????????????????????3