Gabor特征与cvGabor源码解读

一、简介

        在图像方面，Gabor filter是用于边缘检测的线性滤波器。Gabor filters的频率与方向表示与人的视觉系统非常相似，所以其能很好的用于纹理表示和区分。在空间域，一个2D Gabor filter是一个正弦平面波与高斯核函数调制而成。Gabor滤波器是自相似的，也就是说，所有Gabor滤波器都可以从一个母小波经过膨胀和旋转产生。实际应用中，Gabor滤波器可以在频域的不同尺度，不同方向上提取相关特征。

       Gabor变换是一种加窗短时Fourier变换（Window Fourier transform or Short Time Fourier Transform）。Fourier变换是整体上将信号分解为不同的频率分量（任何信号都可分解为复正弦信号之和），对确定性信号及平稳信号使用。其缺点为缺乏时间的局部性信息，并且对时变信号、非平稳信号的分析存在严重不足，(1)无法告知某些频率成分发生在哪些时间内；(2)无法告知某个时刻信号频谱的分布情况。

二、小波变换

         小波变换与Fourier变换相比，是一个时间和频域的局域变换因而能有效地从信号中提取信息，通过伸缩和平移等运算功能对函数或信号进行多尺度细化分析（Multiscale Analysis），解决了Fourier变换不能解决的许多困难问题。小波变换具有良好的时频局部化、尺度变换和方向特征，是分析纹理的有力工具。其优点在于：

        (1)小波分解可以覆盖整个频域（提供了一个数学上完备的描述)

        (2)小波变换通过选取合适的滤波器，可以极大的减小或去除所提取得不同特征之间的相关性

        (3)小波变换具有“变焦”特性，在低频段可用高频率分辨率 和低时间分辨率（宽分析窗口），在高频段，可用低频率分辨率和高时间分辨率（窄分析窗口）

        (4)小波变换实现上有快速算法（Mallat小波分解算法）

三、Gabor函数

        通过频率参数和高斯函数参数的选取，Gabor变换可以选取很多纹理特征，但是Gabor是非正交的，不同特征分量之间有冗余。

        Gabor是有高斯核函数与复正弦函数调制而成，如图所示。

        图(a)为偏移x轴30°的正弦函数，图(b)为高斯核，图(c)为对应的Gabor filter。可以看出正弦函数是如何在空间上具有局部性的。

        二维Gabor函数的第一种表示形式：

 

             为正弦函数的波长，核函数方向，为相位偏移，为高斯标准差，为x、y两个方向的纵横比（指定了Gabor函数的椭圆率）。

       二维Gabor函数的第二种形式：

       v的取值决定了Gabor滤波的波长，u的取值表示Gabor核函数的方向，K表示总的方向数。参数决定了高斯窗口的大小，这里取。下面注释的cvgabor代码就是基于第二种形式的。

      

//cvgabor.h
/***************************************************************************
 *   Copyright (C) 2006 by Mian Zhou   *
 *   M.Zhou@reading.ac.uk   *
 *                                                                         *
 *   This program is free software; you can redistribute it and/or modify  *
 *   it under the terms of the GNU General Public License as published by  *
 *   the Free Software Foundation; either version 2 of the License, or     *
 *   (at your option) any later version.                                   *
 *                                                                         *
 *   This program is distributed in the hope that it will be useful,       *
 *   but WITHOUT ANY WARRANTY; without even the implied warranty of        *
 *   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the         *
 *   GNU General Public License for more details.                          *
 *                                                                         *
 *