[发明专利]一种基于核聚类的虹膜分类方法

权利要求书

1.一种基于核聚类的虹膜分类方法，包括以下两个过程：

一、针对现存虹膜数据库建立聚类的过程，具体包括以下步骤：

步骤1：设现存虹膜数据库中虹膜样本原始图像的数量为N，对现存虹膜数据库的所有 虹膜样本原始图像进行编号；

步骤2：提取虹膜数据库中一张经编号的虹膜样本原始图像I；

步骤3：将步骤2中的虹膜样本原始图像I进行去噪、定位和归一化处理，得到其大小为 64×256像素的虹膜归一化图像A；

步骤4：将步骤3所得的虹膜归一化图像A进行光照补偿和直方图拉伸后，等分成大小 为32×256像素的上下两层图像：上层图像A_up和下层图像A_down；

步骤5：对上层图像A_up采用Daubechies-4小波进行三层分解，删除三层小波分解过程 中得到的三个高频通道下的系数，保留其余7个通道并计算每个通道的平均能量E^up和标准 差STD^up；

然后将各个通道上的平均能量E^up和标准差STD^up按照LL³，LH³，HL³，LH²，HL²，LH¹，HL¹的顺序排列，得到上层图像A_up的分类特征向量Feature_A_up：

Feature_A_up=[E1up,STD1up,E2up,STD2up,E3up,STD3up,E4up,STD4up,E5up,STD5up,E6up,STD6up,E7up,STD7up];]]>

步骤6：对下层图像A_down采用与步骤5相同的处理方法，得到下层图像A_down的分 类特征向量Feature_A_down：

Feature_A_down=]]>

[E1down,STD1down,E2down,STD2down,E3down,STD3down,E4down,STD4down,E5down,STD5down,E6down,STD6down,E7down,STD7down];]]>

步骤7：将上层图像A_up的分类特征向量Feature_A_up与下层图像A_down的分类特征 向量Feature_A_down一一对应相除，得到虹膜样本原始图像I内外半圆在小波系数统计量上 的比值，作为虹膜归一化图像的特征向量：

Feature_A=]]>

[E1upE1down,STD1upSTD1down,E2upE2down,STD2upSTD2down,E3upE3down,STD3upSTD3down,E4upE4down,STD4upSTD4down,E5upE5down,STD5upSTD5down,E6upE6down,STD6upSTD6down,E7upE7down,STD7upSTD7down]]]>

步骤8：反复执行步骤2至步骤7，对虹膜数据库中的所有编号的虹膜样本原始图像，计 算其虹膜归一化图像A的特征向量，得到N个虹膜归一化图像的特征向量，并组成一个N×14 的虹膜归一化图像特征向量矩阵X；

步骤9：定义一个N×2的聚类矩阵B，聚类矩阵B的行号与虹膜数据库中虹膜样本原始 图像的编号一一对应，其中每一行中的两个元素中一个为0，另一个为1，“1”表示对应行号 的虹膜样本原始图像属于“1”所在列所对应的聚类；

步骤10：定义对角矩阵D=1/Σi=1NBi,1001/Σi=1NBi,2,]]>其对角元素依次是聚类矩阵B的第一、 二列所有元素之和的倒数；计算质心矩阵C＝X′BD，其中X′为虹膜归一化图像特征向量矩阵 X的转置矩阵，质心矩阵C的第k列C_k，k＝1，2，表示第k类聚类的质心；

步骤11：计算虹膜数据库中所有虹膜样本原始图像的归一化图像的特征向量Feature_A 与两聚类质心C_k的距离K(φ(x)，C_k)；其中φ(x)表示归一化图像的特征向量Feature_A；x表示 虹膜样本原始图像的编号，即聚类矩阵B的行号；然后比较κ(φ(x)，C₁)与κ(φ(x)，C₂)的大小， 把编号为x的虹膜样本原始图像分到距离聚类的质心最小的那个聚类，即对应修改聚类矩阵 B的x行中“1”的位置；得到更新后的聚类矩阵B′；

步骤12：比较聚类矩阵B′与聚类矩阵B，若二者对应元素完全相同，则执行步骤13；若 二者对应元素不完全相同，则令B＝B′并返回步骤10；

步骤13：聚类矩阵B′的行号与虹膜数据库中虹膜样本原始图像的编号一一对应，其中每 一行中的元素“1”所在列的列号代表对应行号的虹膜样本原始图像所属的聚类，这样就把所 述虹膜数据库的虹膜样本原始图像分成了两个聚类；

通过步骤1至步骤12，就能够得到所述虹膜数据库的第一级分类的两个聚类；若将第一 级分类的两个聚类分别视为两个新的虹膜数据库，再分别采取步骤1至步骤13的操作，就能 够得到所述虹膜数据库的第二级分类的四个聚类；以此类推，则能够得到所述虹膜数据库的 第n级分类的2ⁿ个聚类；

二、待检测虹膜样本原始图像的归类过程，具体包括以下步骤：

步骤14：对于待检测虹膜样本原始图像，进行去噪、定位和归一化处理，得到其大小为 64×256像素的虹膜归一化图像A；

步骤15：将步骤14所得的虹膜归一化图像A进行光照补偿和直方图拉伸后，等分成大 小为32×256像素的上下两层图像：上层图像A_up和下层图像A_down；

步骤16：对上层图像A_up采用Daubechies-4小波进行三层分解，删除三层小波分解过 程中得到的三个高频通道下的系数，保留其余7个通道并计算每个通道的平均能量E^up和标 准差STD^up；

然后将各个通道上的平均能量E^up和标准差STD^up按照LL³，LH³，HL³，LH²，HL²，LH¹，HL¹的顺序排列，得到上层图像A_up的分类特征向量Feature_A_up：

Feature_A_up=[E1up,STD1up,E2up,STD2up,E3up,STD3up,E4up,STD4up,E5up,STD5up,E6up,STD6up,E7up,STD7up];]]>

步骤17：对下层图像A_down采用与步骤16相同的处理方法，得到下层图像A_down 的分类特征向量Feature_A_down：

Feature_A_down=]]>

[E1down,STD1down,E2down,STD2down,E3down,STD3down,E4down,STD4down,E5down,STD5down,E6down,STD6down,E7down,STD7down];]]>

步骤18：将上层图像A_up的分类特征向量Feature_A_up与下层图像A_down的分类特 征向量Feature_A_down一一对应相除，得到待检测虹膜样本原始图像内外半圆在小波系数统 计量上的比值，作为虹膜归一化图像的特征向量：

Feature_A=]]>

[E1upE1down,STD1upSTD1down,E2upE2down,STD2upSTD2down,E3upE3down,STD3upSTD3down,E4upE4down,STD4upSTD4down,E5upE5down,STD5upSTD5down,E6upE6down,STD6upSTD6down,E7upE7down,STD7upSTD7down]]]>

步骤19：构建聚类标签向量Y＝[y₁，y₂，…，y_j，…y_N]^T，其中y_j＝-1或1，当y_j＝-1时表示 现存虹膜数据库中编号为j的虹膜样本原始图像属于聚类一，当y_j＝1时表示现存虹膜数据库 中编号为j的虹膜样本原始图像属于聚类二；然后将步骤8中所述的虹膜归一化图像特征向 量矩阵X和聚类标签向量Y输入到支持向量机中训练，得到支持向量集合SV和判别函数 f(x)，其中判别函数f(x)中：x表示待检测虹膜样本原始 图像的归一化图像特征向量Feature_A；x_j表示虹膜归一化图像特征向量矩阵X中第j行向量； κ(x_j，x)表示x到x_j的距离；α_j表示x_j的权重系数，且当x_j∈SV时α_j＞0，当时 α_j＝0；b表示支持向量机的冗余度；

步骤20：将待检测虹膜样本原始图像的归一化图像特征向量Feature_A输入到步骤19所 述的判别函数f(x)，计算判别函数f(x)的值，当f(x)＝-1时表示待检测虹膜样本原始图像属 于聚类一；当f(x)＝1时表示待检测虹膜样本原始图像属于聚类二；

通过步骤14至步骤20，就能够得到待检测虹膜样本原始图像在所述虹膜数据库中所属 的一级聚类；若将待检测虹膜样本原始图像在所述虹膜数据库中所属的一级聚类视为新的虹 膜数据库，再采取步骤19至步骤20的操作，就能够得到待检测虹膜样本原始图像在所述虹 膜数据库中所属的二级聚类；以此类推，则能够得到待检测虹膜样本原始图像在所述虹膜数 据库中所属的n级聚类。