非高斯二维Gabor滤波器图像纹理提取及识别方法
【技术领域】
[0001] 本发明设及一种非高斯二维G油or滤波器图像纹理提取及识别方法，属于模式识 别技术领域。
【背景技术】
[0002] 二维G油or滤波器是一种有效的图像纹理特征提取方法。关于二维G油or滤波器 的研究主要集中在二维G油or滤波器参数选择，二维G油or滤波器的快速计算和二维G油or 滤波器的各种各样应用方面。
[0003] 传统的二维G油or滤波器可W看作是在频域由复正弦函数调制高斯核函数而成， 故传统的二维G油or滤波器变换到频域仍为高斯型，采用传统二维G油or滤波器对图像进 行滤波处理后，只能提取图像纹理的高斯频率信息，图像中处于非高斯频带部分的图像纹 理信息就被丢失。非高斯二维G油or滤波器在传统二维高斯滤波器的基础上引入一个新的 参数r控制非高斯二维G油or滤波器的形状。非高斯二维G油or滤波器可W看作是在频域 由复正弦函数调制非高斯核函数而成，是传统二维G油or滤波器的一种推广，能提取图像 纹理的非高斯频率信息。
[0004] 人脸识别方面的实验表明通过选择适当的参数r，非高斯二维G油or滤波器在人 脸图像纹理特征提取及识别方面要大大优于传统的二维G油or滤波器。
【发明内容】
[0005] 本发明提供了一种非高斯二维G油or滤波器图像纹理提取及识别方法，W用于解 决非高斯图像纹理特征提取及识别问题。针对局部传统二维Gabor滤波器仅能提取图像纹 理高斯频率信息的缺陷，本发明提出的基于非高斯二维G油or滤波器的图像纹理特征提取 方法能有效提取图像的非高斯纹理信息。
[0006] 本发明非高斯二维G油or滤波器图像纹理提取及识别方法是该样实现的；首先构 造出非高斯二维G油or滤波器，然后采用非高斯二维G油or滤波器对图像进行滤波处理，计 算滤波后的图像的幅度矩阵，接着将此幅度矩阵进行分块，计算各块的直方图特征向量，然 后连接各块的直方图特征向量，得到原始图像的直方图特征向量。最后把该特征向量送入 基于卡方统计量的最近邻分类器进行分类识别。
[0007] 所述非高斯二维G油or滤波器图像纹理提取及识别方法的具体步骤如下：
[000引 Stepl、首先构造非高斯二维G油or滤波器；非高斯二维G油or滤波器构造为；
[0010]其中 口"，kv=kmax/fV，q>,。二孔片/4,kmax是最大频率，一般kmax=n/2， / = ^/^，尺度参数V取值为V= 0,. . .，4，方向参数y取值为y= 0,. . .，3，z= (X，y)和 z|| = (x2+y2)，x，y表示传统二维G油or滤波器的坐标，i表示虚单位，r是一个引进的参 数，r〉l，用于控制非高斯二维G油or滤波器的形状；图2是本发明中不同参数r下非高斯二 维滤波器示意图。对于不同参数下的每一组滤波器，水平方向是方向变化，垂直方向是尺度 变化。图2表明，随着r增大，滤波器掩模非零部分缩小，故更能提取更加精细的局部纹理 f目息。
[0011] 其中，传统二维G油or滤波器的定义如下；
[0012]
[0013] Step2、采用非高斯二维G油or滤波器对一幅图像进行滤波处理；由于非高斯二维 G油or滤波器是一个多尺度多方向的二维滤波器，则在不同尺度和不同方向下对图像进行 滤波处理；
[0014]首先假定图像为IW(X，y)，非高斯二维G油or滤波器为NGu,V, Y(X，y)，那么具 有参数V、y和r的非高斯二维Gabor滤波器对图像I(X，y)进行滤波的计算式为；
[0015] 其中，X，y表示图像的坐标，*表示卷积操作，巧表示坐标X，y处滤波后的 图像，是一个复数，整个滤波后的图像变成一个复数矩阵，表示为的大小与原 图像相同；
[0016]Step3、求复数矩阵巧(心的模，记为诚y作，心），姑5巧(心)是一个实数矩阵，即幅 度矩阵，然后对矩阵进行分块：把均匀分成4X4的无重叠方块，一共 16 块，表示为/,">(i= 0, 1，2,…，15);
[0017]St巧4、提取各分块（i=0, 1，2,…，15)的直方图特征向量，表示为邱 (i= 0, 1，2,…，15)，然后依次连接所有各块的直方图特征向量，得到原始图像的多直方图 特征向量为
[0018] steps、把原始图像的多直方图特征向量送入基于卡方统计量的最近邻分类器进 行分类，来识别原始人脸图像的身份。
[0019] 所述步骤steps中，基于卡方统计量的最近邻分类器进行分类时，先 计算卡方统计量距离；设定两幅人脸图像I?和IW的直方图特征向量分别为；
，则此两 个向量之间的距离，即卡方统计量距离采用如下公式计算：
[0020]
[002U 其中/%G,w?和J分别表示纹理特征向量w和户\6。w的第 i个元素，K'表示纹理向量的长度，eps为一固定值，为Matl油中最小的正数。
[0022] 本发明的有益效果是：
[0023] 1、本发明构造的非高斯二维G油or滤波器是传统二维G油or滤波器的一种扩展， 传统二维G油or滤波器仅仅是非高斯二维G油or滤波器的一种特例；
[0024] 2、本发明构造的非高斯二维G油or滤波器可W看作是在频域由复正弦函数调制 非高斯核函数而成，故能提取图像纹理的非高斯频率信息，而传统的二维Gabor滤波器可 W看作是在频域由复正弦函数调制高斯核函数而成，故只能提取图像纹理的高斯频率信 息，在采用传统二维Gabor滤波器对图像进行滤波后，处于非高斯频带部分的图像纹理信 息就被丢失；
[0025] 3、本发明对采用非高斯二维G油or滤波器滤波后图像的幅度矩阵进行分块，计算 各分块的直方图特征，连接各分块的直方图特征得到原始图像总的直方图特征，具有很强 的分类能力；
[0026] 4、在匹配识别阶段，本发明采用卡方统计量烟lisquarestatistic)作为两个纹 理特征向量之间的距离度量，采用最近邻分类器进行分类，算法简单、计算方便，可W做到 实时的图像匹配识别。
【附图说明】
[0027] 图1是本发明中采用非高斯二维G油or滤波器提取图像纹理特征的步骤示意框 图；
[002引图2是本发明中不同参数r下非高斯二维滤波器示意图，对于不同参数下的每一 组滤波器，水平方向是方向变化，垂直方向是尺度变化。
[0029] 图3是本发明所用theextendedYaleB人脸数据库光照子集中一个人的64张样 本图像；
[0030] 图4是本发明实施例2中不同参数下本方法人脸识别率。
【具体实施方式】
[0031] 实施例1;如图1-4所示，非高斯二维G油or滤波器图像纹理提取及识别方法，首 先构造出非高斯二维G油or滤波器，然后采用非高斯二维G油or滤波器对图像进行滤波处 理，计算滤波后的图像的幅度矩阵，接着将此幅度矩阵进行分块，计算各块的直方图特征向 量，然后连接各块的直方图特征向量，得到原始图像的直方图特征向量。最后把该特征向量 送入基于卡方统计量的最近邻分类器进行分类识别。
[0032] 所述非高斯二维G油or滤波器图像纹理提取及识别方法的具体步骤如下：
[0033]Stepl、首先构造非高斯二维G油or滤波器；非高斯二维G油or滤波器构造为；
[0034]
[0035]其中 口"，kv=kmayfV，CP/,=7l///4,kmax是最大频率，一般kmax=n/2， /二>/5，尺度参数V取值为V = 0,. . .，4，方向参数y取值为y = 0,. . .，3，z = (X，y)和 IIzII= (x2+y2)，x,y表示传统二维G油or滤波器的坐标，i表示虚单位，r是一个引进的参 数，r〉l，用于控制非高斯二维G油or滤波器的形状；