Gabor小波变换

Matlab Gabor小波变换 可以直接运行 只需换图片

% function [GGABOUT]=gaborfilter（ISFWP）;
I=imread（‘columbia.bmp‘）
if isa（I‘double‘）~=1
    I=double（I）;
end
size=fix（1.5/S）; % exp（-1.5^2*pi） < 0.1%
%k=2*pi*S^2;
%F=S^2/sqrt（2*pi）;
k=1;
for x=-size:size
    for y=-size:size
        G（size+x+1size+y+1）=k*exp（-pi*S^2*（x*x+y*y））*...
            （exp（j*（2*pi*F*（x*cos（W）+y*sin（W））+P））-exp（-pi*（F/S）^2+j*P））;
    end
end

GABOUT=conv2（Idouble（G）‘same‘）; 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件     231480  2009-05-09 22:19  gabor\columbia.bmp

     文件        425  2009-11-04 11:20  gabor\gabor.m

     目录          0  2009-11-05 14:04  gabor

----------- ---------  ---------- -----  ----

               231905                    3