Gabor wavelets transform 【Gabor 小波变换】
18:22:14
| 一：定义 There is considerable evidence (reviewed in MacLennan 1991) that images in primary visual cortex (V1) are represented in terms of Gabor wavelets, that is, hierarchically arranged, Gaussian-modulated sinusoids (equivalent to the pure states of quantum mechanics). The Gabor-wavelet transform of a two-dimensional visual field generates a four-dimensional field: two of the dimensions are spatial, the other two represent spatial frequency and orientation. To represent this four-dimensional field in two-dimensional cortex, it is necessary to ``slice'' the field, which gives rise to the columns and stripes of striate cortex. The representation is nearly optimal, as defined by the Gabor Uncertainty Principle (a generalization of the Heisenberg Uncertainty Principle to information representation and transmission). Time-varying two-dimensional visual images may be viewed as three-dimensional functions of space-time, and it is possible that time-varying images are represented in vision areas by a three-dimensional Gabor-wavelet transform, which generates a time-varying five-dimensional field (representing two spatial dimensions, spatial frequency, spatial orientation and temporal frequency). The effect is to represent the ``optic flow'' of images in terms spatially fixed, oriented grating patches with moving gratings. (See MacLennan 1991 for more details.) Finally, Pribram provides evidence that Gabor representations are also used for controlling the generation of motor fields (see citations in MacLennan 1997, p. 64). Gabor 滤波器，
其中， 在上式 ，λ代表余弦函数波长，θ代表Gabor函数平行条纹法线的方向【θ epresents the orientation of the normal to the parallel stripes】，ψ为相位补偿，σ为高斯包络的sigma，γ为空间相位比和Gabor函数中的椭圆率 以下是MATLAB中实现： function gb=gabor_fn(sigma,theta,lambda,psi,gamma) ============================================================== 另外一个MATLAB实现函数： function [gab]=gaborcreate(M,N) |
Gabor滤波实现(1)
2007年12月23日 星期日 10:04
| 最近忙着论文，需要Gabor滤波代码，可是网上总找不到合适的代码，于是就自己编了一个，不当之处请指点。参考论文为 L. Wiskott，J. M. Fellous，N. Kruger，C. v. d.Malsburg. Face Recognition by Elastic Bunch Graph Matching，IEEE Trans. On PAMI，Vol.19，No.7，pp775-779，1997 首先实现滤波器： function [bank] = do_createfilterbank(imsize,varargin) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% conf = struct(..., conf = do_getargm(conf,varargin); bank = cell(1,conf.freqnum*conf.orientnum); for ind = 1:length(bank) function [filter,Kv,Phiu] = do_gabor (imsize,nu,mu,Kmax,f,sigma,orientnum) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% rows = imsize(1); Kv = Kmax/f^nu; F1 = (Kv ^ 2)/ (sigma^2) * exp(-Kv^2 * abs(X.^2 + Y.^2) / (2*sigma^2)) ; |
非常感谢
Gabor 滤波实现(1)已经创建了Gabor滤波组，现在可以使用该滤波组对图像进行转换，得到振幅和相位。
function [result] = do_filterwithbank(im,bank)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% 函数实现：对图像使用Gabor滤波组进行转换换
%
%%% 参数：
% im — 被转换的图像
% bank — 由函数do_createfilterbank得到的滤波组
%
%%% 返回：
% result — 图像被转换后的结果
% .amplitudes — 不同像素点的振幅向量
% .phases — 不同像素点的相位向量
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
[N1 N2] = size(im);
N3 = length(bank);
phases = zeros(N1,N2,N3);
amplitudes = zeros(N1,N2,N3);
imagefft = fft2(im);
for ind = 1:N3
fprintf('正在处理滤波 %d /n',ind);
temp = ifft2(imagefft .* bank{ind}.filter);
phases(:,:,ind) = angle(temp);
amplitudes(:,:,ind) = abs(temp);
end
result.phases = phases;
result.amplitudes = amplitudes;
整个程序可以如下使用。 im = imread('image.jpg'); im = rgb2gray(im); bank = do_createfilterbank(size(im)); result = do_filterwithbank(im,bank);