Gabor滤波

Halcom

Gabor滤波:
Sx = 0.5;  % x方向的差异系数
Sy = 0.5;  % y方向的差异系数
U = 1.0;   % x方向的中心频率
V = 1.0;   % y方向的中心频率
[G,gabout] = gabor_filter(im,Sx,Sy,U,V);

1. function [G,gabout] = gabor_filter(I,Sx,Sy,U,V)
   
2. % 函数输入：
   
3. %          I : 输入的二维图像
   
4. %          Sx & Sy : 方差 along x and y-axes respectively
   
5. %          U & V : 中心频率  along x and y-axes respectively
   
6. % 函数输出：
   
7. %          G : G（x，y）
   
8. %          gabout : Gabor滤波图像
   
9. % G(x,y)表达式如下：
   
10. %               1                -1     x  ^    y  ^
    
11. %%% G(x,y) = ---------- * exp ([----{(----) 2+(----) 2}+2*pi*i*(Ux+Vy)])
    
12. %            2*pi*sx*sy           2    sx       sy
    
13. 
    
14. if isa(I,'double')~=1
    
15.     I = double(I);
    
16. end
    
17. 
    
18. for x = -fix(Sx):fix(Sx)
    
19.     for y = -fix(Sy):fix(Sy)
    
20.         G(fix(Sx)+x+1,fix(Sy)+y+1) = (1/(2*pi*Sx*Sy))*exp(-.5*((x/Sx)^2+(y/Sy)^2)+2*pi*i*(U*x+V*y));
    
21.     end
    
22. end
    
23. 
    
24. Imgabout = conv2(I,double(imag(G)),'same');  % 卷积
    
25. Regabout = conv2(I,double(real(G)),'same');  % 卷积
    
26. 
    
27. gabout = uint8(sqrt(Imgabout.*Imgabout + Regabout.*Regabout)); % 输出滤波图像

复制代码

参考：【1】http://halcom.cn/forum.php?mod=viewthread&tid=737