专利名称：：一种基于核聚类的虹膜分类方法
技术领域：
：本发明属于图像处理
技术领域：
，主要涉及生物特征鉴别中的虹膜身份识别技术。
背景技术：
：目前，基于虹膜特征的身份识别技术己经渗透到日常生活的每一个方面，虹膜数据库的样本量也越来越大。例如英国政府发起的IRIS计划要求至少一百万个经常在英国出入境的旅行者通过视频头采集虹膜进行身份识别，而不需要随身附带其他的身份鉴别证明；同时，英国政府正试图将生物信息加入到公民身份证中，如果成功采集到英国大约六千万公民的虹膜图像，那虹膜数据库的样本将达到海量。而我国人口是英国的约20倍，出入境人数更是达到了3亿多。所以，在海量虹膜数据库的情况下对于特定虹膜样本的识别过程中，需要在虹膜数据库中快速而准确的找到其分类检索，这对于特定虹膜样本的识别具有十分必要意义的。在海量虹膜数据库中建立适当标准的分类检索不仅有利于提高虹膜识别准确率，更加有利于提高虹膜识别的效率。目前虹膜识别
技术领域：
的研究主要集中在虹膜的定位与匹配识别方面，尚没有针对海量虹膜数据库进行分类索引的研究。
发明内容本发明的任务是提供一种基于核聚类的虹膜分类方法，该方法首先针对现有虹膜数据库(尤其是海量数据库)建立一种便于分类索引的聚类，然后针对特定的待识别虹膜样本进行分类。本发明有利于提高虹膜识别准确率和识别效率。本发明详细技术方案如下一种基于核聚类的虹膜分类方法，如图1所示，包括以下两个过程一、针对现存虹膜数据库建立聚类的过程，具体包括以下步骤步骤1:设现存虹膜数据库中虹膜样本原始图像的数量为iV，对现存虹膜数据库的所有虹膜样本原始图像进行编号；步骤2:提取虹膜数据库中一张经编号的虹膜样本原始图像/;步骤3:将步骤2中的虹膜样本原始图像/进行去噪、定位和归一化处理，得到其大小为64x256像素的虹膜归一化图像A;步骤4:将步骤3所得的虹膜归一化图像A进行光照补偿和直方图拉伸后，等分成大小为32x256像素的上下两层图像上层图像A_up和下层图像A—down;步骤5:对上层图像A一up采用Daubechies-4小波进行三层分解，删除三层小波分解过程中得到的三个高频通道下的系数，保留其余7个通道并计算每个通道的平均能量五"P和标准差mr;然后将各个通道上的平均能量和标准差S7Z)up按照Z￡3,Z//3，/tt3，丄//2,/ffi2,丄//1，/ffi1的顺序排列，得到上层图像A—up的分类特征向量Feature—A—up:Feat鹏—A—up-[^,節r，五2"P,5T《,《，Sraf，《，，,^,即5"P，E7，Sn^,^,，7"P];步骤6:对下层图像A—down采用与步骤5相同的处理方法，得到下层图像A—down的分类特征向量Feature—A—down:Feature—A—down=步骤7:将上层图像A一up的分类特征向量Feature—A—up与下层图像A—down的分类特征向量Feature一B—down—一对应相除，得到标准虹膜原始图像J内外半圆在小波系数统计量上的比值，作为虹膜归一化图像的特征向量Feature—A=67Zf1S7Zf!57Zf《!57Zf!6TZf^57Zf一步骤8:反复执行步骤2至步骤7，对虹膜数据库中的所有编号的虹膜样本原始图像，计算其虹膜归一化图像A的特征向量，得到iV个虹膜归一化图像的特征向量，并组成一个iVx14的虹膜归一化图像特征向量矩阵X。步骤9:定义一个7Vx2的聚类矩阵B，聚类矩阵B的行号与虹膜数据库中虹膜样本原始图像的编号一一对应，其中每一行中的两个元素中一个为O，另一个为l，"l"表示对应行号的虹膜样本原始图像属于"1"所在列所对应的聚类。7步骤10:定义对角矩阵D-I"0o1/f;Bi.，其对角元素依次是聚类矩阵B的第一、二列所有元素之和的倒数；计算质心矩阵C-X'BD，其中X'为虹膜归一化图像特征向量矩阵X的转置矩阵，质心矩阵C的第it列C;，A:=l，2，表示第A:类聚类的质心。步骤lh计算虹膜数据库中所有虹膜样本原始图像的归一化图像的特征向量Feature—A与两聚类质心C;的距离Ar^(力，C;);其中-("表示归一化图像的特征向量Feature—A;x表示虹膜样本原始图像的编号，即聚类矩阵B的行号;然后比较K(冉x),C,)与A:W(x),q)的大小，把编号为x的虹膜样本原始图像分到距离聚类的质心最小的那个聚类，即对应修改聚类矩阵B的x行中"1"的位置；得到更新后的聚类矩阵B'。步骤12:比较聚类矩阵B'与聚类矩阵B，若二者对应元素完全相同，则执行步骤13;若二者对应元素不完全相同，则令B-B'并返回步骤IO。步骤13:聚类矩阵B'的行号与虹膜数据库中虹膜样本原始图像的编号一一对应，其中每一行中的元素"1"所在列的列号代表对应行号的虹膜样本原始图像所属的聚类，这样就把所述虹膜数据库的虹膜样本原始图像分成了两个聚类。通过步骤1至步骤12，就能够得到所述虹膜数据库的第一级分类的两个聚类；若将第一级分类的两个聚类分别视为两个新的虹膜数据库，再分别采取步骤1至步骤13的操作，就能够得到所述虹膜数据库的第二级分类的四个聚类；以此类推，则能够得到所述虹膜数据库的第"级分类的2"个聚类。二、待检测虹膜样本原始图像的归类过程，具体包括以下步骤步骤14:对于待检测虹膜样本原始图像，进行去噪、定位和归一化处理，得到其大小为64x256像素的虹膜归一化图像A;步骤15:将步骤14所得的虹膜归一化图像A进行光照补偿和直方图拉伸后，等分成大小为32x256像素的上下两层图像上层图像A—up和下层图像A一down;步骤16:对上层图像A—up采用Daubechies-4小波进行三层分解，删除三层小波分解过程中得到的三个高频通道下的系数，保留其余7个通道并计算每个通道的平均能量五"p和标然后将各个通道上的平均能量五"p和标准差S7P"p按照丄￡3,￡//3,/ffi3，￡//2,/ft2,Z//1，/ffi1的顺序排列，得到上层图像A—up的分类特征向量Feature—A—up:Feature—A—up=[《，即严，￡2"p,即7,《,S7P3"P，五7,S7Df,￡7,S7Z)，，《，S7D6"P,￡，,節7];步骤17:对下层图像A—down采用与步骤16相同的处理方法，得到下层图像A—down的分类特征向量Feature—A—down:Feature—A—down=步骤18:将上层图像A—up的分类特征向量Feature—A—up与下层图像A—down的分类特征向量Feature—B—down—一对应相除，得到标准虹膜原始图像/内外半圆在小波系数统计量上的比值，作为虹膜归一化图像的特征向量Feature—A=.157Zf!57Zf!$71^!>S7If!S7Zf《5TIf!步骤19:构建聚类标签向量Yz[x，;;2,…,力,…;v]r，其中力=-1或1，当R--l时表示现存虹膜数据库中编号为J'的虹膜样本原始图像属于聚类一，当^=1时表示现存虹膜数据库中编号为J'的虹膜样本原始图像属于聚类二；然后将步骤8中所述的虹膜归一化图像特征向量矩阵X和聚类标签向量Y输入到支持向量机中训练，得到支持向量集合^和判别函数/(;c)，其中/(x)-sgnZaj,，:c)+6卜判别函数/(x)中x表示待检测虹膜样本原始图像的归一化图像特征向量Feature—A;.表示虹膜归一化图像特征向量矩阵X中第7'行向量；K(jc";c)表示jc到^.的距离；^表示、的权重系数，且当^.eSr时>0，当、.gSr时=0;力表示支持向量机的冗余度。步骤20:将待检测虹膜样本原始图像的归一化图像特征向量;c输入到步骤19所述的判别函数/(x)，计算判别函数/(x)的值，当/(;c)=-l时表示待检测虹膜样本原始图像属于聚类一；当/00=1时表示待检测虹膜样本原始图像属于聚类二。通过步骤14至步骤20，就能够得到待检测虹膜样本原始图像在所述虹膜数据库中所属的一级聚类；若将待检测虹膜样本原始图像在所述虹膜数据库中所属的一级聚类视为新的虹膜数据库，再采取步骤19至步骤20的操作，就能够得到待检测虹膜样本原始图像在所述虹；以此类推，则能够得到待检测虹膜样本原始图像在所述虹膜数据库中所属的w级聚类。需要说明的是步骤11和步骤19中计算所述距离时，具体计算方法可以采用范数距离或核空间距离的计算方法。若采用核空间距离的计算方法，具体方法可以采用高斯核空间距离计算方法。但不管采用什么计算方法，步骤11和步骤19应当采用同样的距离计算方法。本发明的有益效果是本发明创新性地提出了一种基于核聚类的虹膜分类的方法，将海量的虹膜数据库中的虹膜样本进行无监督地聚类并对待检测的虹膜样本原始图像进行分类，找到其所属的聚类。本发明利用虹膜纹理能量和标准差在方向、频率、位置上的不同，使用Daubechies-4小波对虹膜纹理提取并构建特征向量；在无监督聚类算法中运用了核方法，得到比经典聚类方法(如C均值聚类)更接近数据特性、性能更为良好的分类结果；利用聚类中核算法的多样性，机动地、灵活地选取适合表达虹膜纹理的核算法，而不是经典聚类算法中确定一种算法，就只能得到一种聚类结果；最后采用支持向量机对待测试虹膜样本进行分类，使分类成为凸最优化问题，更容易得到全局最优解。本发明能够快速而准确的找到待检测的虹膜样本原始图像所属的聚类，从而减少虹膜识别时待检测的虹膜样本原始图像在海量虹膜库中匹配的时间，有利于提高虹膜识别准确率和识别效率。图1:本发明提供的一种基于核聚类的虹膜分类方法的流程示意图。图2:归一化示意图。图3:小波分解示意图。步骤4中所提及的7个通道分别为"3,￡//3，/tt3，￡//2,itt2,i^1，H丄1。具体实施方式采用本发明的方法，使用Matlab语言编写虹膜识别软件，然后用中科院自动化所采集的CASIA2.0虹膜数据库为源数据。从中提取出40个人的虹膜图片，分别对其进行高斯噪声(选取25组不同的参数)的攻击，乘积性噪声(选取25组不同的参数)的攻击和不同角度(角度范围0度到10度)的旋转，一共2400张图片进行分类。正确分类结果为98.3%,各个噪声的正确分类率为96.6%，97.8%和96.3%。10同时从CASIA2.0库中选取40个人,每人9张拍摄于不同时期的虹膜图片进行分类训练，成功将40个人的虹膜分成大小合适的3类，并得到支持向量机的各个参数。另外选取同样40个人，每人3张的虹膜图片作为测试库。经过提取特征向量后，输入支持向量机进行分类测试，得到平均92%的分类正确率。权利要求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个通道并计算每个通道的平均能量Eup和标准差STDup；然后将各个通道上的平均能量Eup和标准差STDup按照LL3，LH3，HL3，LH2，HL2，LH1，HL1的顺序排列，得到上层图像A_up的分类特征向量Feature_A_up<mathsid="math0001"num="0001"><math><![CDATA[<mrow><mi>Feature</mi><mo>_</mo><mi>A</mi><mo>_</mo><mi>up</mi><mo>=</mo><mo>[</mo><msubsup><mi>E</mi><mn>1</mn><mi>up</mi></msubsup><mo>,</mo><msubsup><mi>STD</mi><mn>1</mn><mi>up</mi></msubsup><mo>,</mo><msubsup><mi>E</mi><mn>2</mn><mi>up</mi></msubsup><mo>,</mo><msubsup><mi>STD</mi><mn>2</mn><mi>up</mi></msubsup><mo>,</mo><msubsup><mi>E</mi><mn>3</mn><mi>up</mi></msubsup><mo>,</mo><msubsup><mi>STD</mi><mn>3</mn><mi>up</mi></msubsup><mo>,</mo><msubsup><mi>E</mi><mn>4</mn><mi>up</mi></msubsup><mo>,</mo><msubsup><mi>STD</mi><mn>4</mn><mi>up</mi></msubsup><mo>,</mo><msubsup><mi>E</mi><mn>5</mn><mi>up</mi></msubsup><mo>,</mo><msubsup><mi>STD</mi><mn>5</mn><mi>up</mi></msubsup><mo>,</mo><msubsup><mi>E</mi><mn>6</mn><mi>up</mi></msubsup><mo>,</mo><msubsup><mi>STD</mi><mn>6</mn><mi>up</mi></msubsup><mo>,</mo><msubsup><mi>E</mi><mn>7</mn><mi>up</mi></msubsup><mo>,</mo><msubsup><mi>STD</mi><mn>7</mn><mi>up</mi></msubsup><mo>]</mo><mo>;</mo></mrow>]]></math></maths>步骤6对下层图像A_down采用与步骤5相同的处理方法，得到下层图像A_down的分类特征向量Feature_A_down<mathsid="math0002"num="0002"><math><![CDATA[<mrow><mi>Feature</mi><mo>_</mo><mi>A</mi><mo>_</mo><mi>down</mi><mo>=</mo></mrow>]]></math></maths><mathsid="math0003"num="0003"><math><![CDATA[<mrow><mo>[</mo><msubsup><mi>E</mi><mn>1</mn><mi>down</mi></msubsup><mo>,</mo><msubsup><mi>STD</mi><mn>1</mn><mi>down</mi></msubsup><mo>,</mo><msubsup><mi>E</mi><mn>2</mn><mi>down</mi></msubsup><mo>,</mo><msubsup><mi>STD</mi><mn>2</mn><mi>down</mi></msubsup><mo>,</mo><msubsup><mi>E</mi><mn>3</mn><mi>down</mi></msubsup><mo>,</mo><msubsup><mi>STD</mi><mn>3</mn><mi>down</mi></msubsup><mo>,</mo><msubsup><mi>E</mi><mn>4</mn><mi>down</mi></msubsup><mo>,</mo><msubsup><mi>STD</mi><mn>4</mn><mi>down</mi></msubsup><mo>,</mo><msubsup><mi>E</mi><mn>5</mn><mi>down</mi></msubsup><mo>,</mo><msubsup><mi>STD</mi><mn>5</mn><mi>down</mi></msubsup><mo>,</mo><msubsup><mi>E</mi><mn>6</mn><mi>down</mi></msubsup><mo>,</mo><msubsup><mi>STD</mi><mn>6</mn><mi>down</mi></msubsup><mo>,</mo><msubsup><mi>E</mi><mn>7</mn><mi>down</mi></msubsup><mo>,</mo><msubsup><mi>STD</mi><mn>7</mn><mi>down</mi></msubsup><mo>]</mo><mo>;</mo></mrow>]]></math></maths>步骤7将上层图像A_up的分类特征向量Feature_A_up与下层图像A_down的分类特征向量Feature_B_down一一对应相除，得到标准虹膜原始图像I内外半圆在小波系数统计量上的比值，作为虹膜归一化图像的特征向量<mathsid="math0004"num="0004"><math><![CDATA[<mrow><mi>Feature</mi><mo>_</mo><mi>A</mi><mo>=</mo></mrow>]]></math></maths><mathsid="math0005"num="0005"><math><![CDATA[<mrow><mo>[</mo><mfrac><msubsup><mi>E</mi><mn>1</mn><mi>up</mi></msubsup><msubsup><mi>E</mi><mn>1</mn><mi>down</mi></msubsup></mfrac><mo>,</mo><mfrac><msubsup><mi>STD</mi><mn>1</mn><mi>up</mi></msubsup><msubsup><mi>STD</mi><mn>1</mn><mi>down</mi></msubsup></mfrac><mo>,</mo><mfrac><msubsup><mi>E</mi><mn>2</mn><mi>up</mi></msubsup><msubsup><mi>E</mi><mn>2</mn><mi>down</mi></msubsup></mfrac><mo>,</mo><mfrac><msubsup><mi>STD</mi><mn>2</mn><mi>up</mi></msubsup><msubsup><mi>STD</mi><mn>2</mn><mi>down</mi></msubsup></mfrac><mo>,</mo><mfrac><msubsup><mi>E</mi><mn>3</mn><mi>up</mi></msubsup><msubsup><mi>E</mi><mn>3</mn><mi>down</mi></msubsup></mfrac><mo>,</mo><mfrac><msubsup><mi>STD</mi><mn>3</mn><mi>up</mi></msubsup><msubsup><mi>STD</mi><mn>3</mn><mi>down</mi></msubsup></mfrac><mo>,</mo><mfrac><msubsup><mi>E</mi><mn>4</mn><mi>up</mi></msubsup><msubsup><mi>E</mi><mn>4</mn><mi>down</mi></msubsup></mfrac><mo>,</mo><mfrac><msubsup><mi>STD</mi><mn>4</mn><mi>up</mi></msubsup><msubsup><mi>STD</mi><mn>4</mn><mi>down</mi></msubsup></mfrac><mo>,</mo><mfrac><msubsup><mi>E</mi><mn>5</mn><mi>up</mi></msubsup><msubsup><mi>E</mi><mn>5</mn><mi>down</mi></msubsup></mfrac><mo>,</mo><mfrac><msubsup><mi>STD</mi><mn>5</mn><mi>up</mi></msubsup><msubsup><mi>STD</mi><mn>5</mn><mi>down</mi></msubsup></mfrac><mo>,</mo><mfrac><msubsup><mi>E</mi><mn>6</mn><mi>up</mi></msubsup><msubsup><mi>E</mi><mn>6</mn><mi>down</mi></msubsup></mfrac><mo>,</mo><mfrac><msubsup><mi>STD</mi><mn>6</mn><mi>up</mi></msubsup><msubsup><mi>STD</mi><mn>6</mn><mi>down</mi></msubsup></mfrac><mo>,</mo><mfrac><msubsup><mi>E</mi><mn>7</mn><mi>up</mi></msubsup><msubsup><mi>E</mi><mn>7</mn><mi>down</mi></msubsup></mfrac><mo>,</mo><mfrac><msubsup><mi>STD</mi><mn>7</mn><mi>up</mi></msubsup><msubsup><mi>STD</mi><mn>7</mn><mi>down</mi></msubsup></mfrac><mo>]</mo></mrow>]]></math></maths>步骤8反复执行步骤2至步骤7，对虹膜数据库中的所有编号的虹膜样本原始图像，计算其虹膜归一化图像A的特征向量，得到N个虹膜归一化图像的特征向量，并组成一个N×14的虹膜归一化图像特征向量矩阵X；步骤9定义一个N×2的聚类矩阵B，聚类矩阵B的行号与虹膜数据库中虹膜样本原始图像的编号一一对应，其中每一行中的两个元素中一个为0，另一个为1，“1”表示对应行号的虹膜样本原始图像属于“1”所在列所对应的聚类；步骤10定义对角矩阵<mathsid="math0006"num="0006"><math><![CDATA[<mrow><mi>D</mi><mo>=</mo><mfencedopen='['close=']'><mtable><mtr><mtd><mn>1</mn><mo>/</mo><munderover><mi>Σ</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>B</mi><mrow><mi>i</mi><mo>,</mo><mn>1</mn></mrow></msub></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>1</mn><mo>/</mo><munderover><mi>Σ</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>B</mi><mrow><mi>i</mi><mo>,</mo><mn>2</mn></mrow></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math>id="icf0006"file="A2009100593580003C1.tif"wi="43"he="19"top="69"left="71"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>其对角元素依次是聚类矩阵B的第一、二列所有元素之和的倒数；计算质心矩阵C＝X′BD，其中X′为虹膜归一化图像特征向量矩阵X的转置矩阵，质心矩阵C的第k列Ck，k＝1，2，表示第k类聚类的质心；步骤11计算虹膜数据库中所有虹膜样本原始图像的归一化图像的特征向量Feature_A与两聚类质心Ck的距离κ(φ(x)，Ck)；其中φ(x)表示归一化图像的特征向量Feature_A；x表示虹膜样本原始图像的编号，即聚类矩阵B的行号；然后比较κ(φ(x)，C1)与κ(φ(x)，C2)的大小，把编号为x的虹膜样本原始图像分到距离聚类的质心最小的那个聚类，即对应修改聚类矩阵B的x行中“1”的位置；得到更新后的聚类矩阵B′；步骤12比较聚类矩阵B′与聚类矩阵B，若二者对应元素完全相同，则执行步骤13；若二者对应元素不完全相同，则令B＝B′并返回步骤10；步骤13聚类矩阵B′的行号与虹膜数据库中虹膜样本原始图像的编号一一对应，其中每一行中的元素“1”所在列的列号代表对应行号的虹膜样本原始图像所属的聚类，这样就把所述虹膜数据库的虹膜样本原始图像分成了两个聚类；通过步骤1至步骤12，就能够得到所述虹膜数据库的第一级分类的两个聚类；若将第一级分类的两个聚类分别视为两个新的虹膜数据库，再分别采取步骤1至步骤13的操作，就能够得到所述虹膜数据库的第二级分类的四个聚类；以此类推，则能够得到所述虹膜数据库的第n级分类的2n个聚类；二、待检测虹膜样本原始图像的归类过程，具体包括以下步骤步骤14对于待检测虹膜样本原始图像，进行去噪、定位和归一化处理，得到其大小为64×256像素的虹膜归一化图像A；步骤15将步骤14所得的虹膜归一化图像A进行光照补偿和直方图拉伸后，等分成大小为32×256像素的上下两层图像上层图像A_up和下层图像A_down；步骤16对上层图像A_up采用Daubechies-4小波进行三层分解，删除三层小波分解过程中得到的三个高频通道下的系数，保留其余7个通道并计算每个通道的平均能量Eup和标准差STDup；然后将各个通道上的平均能量Eup和标准差STDup按照LL3，LH3，HL3，LH2，HL2，LH1，HL1的顺序排列，得到上层图像A_up的分类特征向量Feature_A_up<mathsid="math0007"num="0007"><math><![CDATA[<mrow><mi>Feature</mi><mo>_</mo><mi>A</mi><mo>_</mo><mi>up</mi><mo>=</mo><mo>[</mo><msubsup><mi>E</mi><mn>1</mn><mi>up</mi></msubsup><mo>,</mo><msubsup><mi>STD</mi><mn>1</mn><mi>up</mi></msubsup><mo>,</mo><msubsup><mi>E</mi><mn>2</mn><mi>up</mi></msubsup><mo>,</mo><msubsup><mi>STD</mi><mn>2</mn><mi>up</mi></msubsup><mo>,</mo><msubsup><mi>E</mi><mn>3</mn><mi>up</mi></msubsup><mo>,</mo><msubsup><mi>STD</mi><mn>3</mn><mi>up</mi></msubsup><mo>,</mo><msubsup><mi>E</mi><mn>4</mn><mi>up</mi></msubsup><mo>,</mo><msubsup><mi>STD</mi><mn>4</mn><mi>up</mi></msubsup><mo>,</mo><msubsup><mi>E</mi><mn>5</mn><mi>up</mi></msubsup><mo>,</mo><msubsup><mi>STD</mi><mn>5</mn><mi>up</mi></msubsup><mo>,</mo><msubsup><mi>E</mi><mn>6</mn><mi>up</mi></msubsup><mo>,</mo><msubsup><mi>STD</mi><mn>6</mn><mi>up</mi></msubsup><mo>,</mo><msubsup><mi>E</mi><mn>7</mn><mi>up</mi></msubsup><mo>,</mo><msubsup><mi>STD</mi><mn>7</mn><mi>up</mi></msubsup><mo>]</mo><mo>;</mo></mrow>]]></math></maths>步骤17对下层图像A_down采用与步骤16相同的处理方法，得到下层图像A_down的分类特征向量Feature_A_down<mathsid="math0008"num="0008"><math><![CDATA[<mrow><mi>Feature</mi><mo>_</mo><mi>A</mi><mo>_</mo><mi>down</mi><mo>=</mo></mrow>]]></math></maths><mathsid="math0009"num="0009"><math><![CDATA[<mrow><mo>[</mo><msubsup><mi>E</mi><mn>1</mn><mi>down</mi></msubsup><mo>,</mo><msubsup><mi>STD</mi><mn>1</mn><mi>down</mi></msubsup><mo>,</mo><msubsup><mi>E</mi><mn>2</mn><mi>down</mi></msubsup><mo>,</mo><msubsup><mi>STD</mi><mn>2</mn><mi>down</mi></msubsup><mo>,</mo><msubsup><mi>E</mi><mn>3</mn><mi>down</mi></msubsup><mo>,</mo><msubsup><mi>STD</mi><mn>3</mn><mi>down</mi></msubsup><mo>,</mo><msubsup><mi>E</mi><mn>4</mn><mi>down</mi></msubsup><mo>,</mo><msubsup><mi>STD</mi><mn>4</mn><mi>down</mi></msubsup><mo>,</mo><msubsup><mi>E</mi><mn>5</mn><mi>down</mi></msubsup><mo>,</mo><msubsup><mi>STD</mi><mn>5</mn><mi>down</mi></msubsup><mo>,</mo><msubsup><mi>E</mi><mn>6</mn><mi>down</mi></msubsup><mo>,</mo><msubsup><mi>STD</mi><mn>6</mn><mi>down</mi></msubsup><mo>,</mo><msubsup><mi>E</mi><mn>7</mn><mi>down</mi></msubsup><mo>,</mo><msubsup><mi>STD</mi><mn>7</mn><mi>down</mi></msubsup><mo>]</mo><mo>;</mo></mrow>]]></math></maths>步骤18将上层图像A_up的分类特征向量Feature_A_up与下层图像A_down的分类特征向量Feature_B_down一一对应相除，得到标准虹膜原始图像I内外半圆在小波系数统计量上的比值，作为虹膜归一化图像的特征向量<mathsid="math0010"num="0010"><math><![CDATA[<mrow><mi>Feature</mi><mo>_</mo><mi>A</mi><mo>=</mo></mrow>]]></math></maths><mathsid="math0011"num="0011"><math><![CDATA[<mrow><mo>[</mo><mfrac><msubsup><mi>E</mi><mn>1</mn><mi>up</mi></msubsup><msubsup><mi>E</mi><mn>1</mn><mi>down</mi></msubsup></mfrac><mo>,</mo><mfrac><msubsup><mi>STD</mi><mn>1</mn><mi>up</mi></msubsup><msubsup><mi>STD</mi><mn>1</mn><mi>down</mi></msubsup></mfrac><mo>,</mo><mfrac><msubsup><mi>E</mi><mn>2</mn><mi>up</mi></msubsup><msubsup><mi>E</mi><mn>2</mn><mi>down</mi></msubsup></mfrac><mo>,</mo><mfrac><msubsup><mi>STD</mi><mn>2</mn><mi>up</mi></msubsup><msubsup><mi>STD</mi><mn>2</mn><mi>down</mi></msubsup></mfrac><mo>,</mo><mfrac><msubsup><mi>E</mi><mn>3</mn><mi>up</mi></msubsup><msubsup><mi>E</mi><mn>3</mn><mi>down</mi></msubsup></mfrac><mo>,</mo><mfrac><msubsup><mi>STD</mi><mn>3</mn><mi>up</mi></msubsup><msubsup><mi>STD</mi><mn>3</mn><mi>down</mi></msubsup></mfrac><mo>,</mo><mfrac><msubsup><mi>E</mi><mn>4</mn><mi>up</mi></msubsup><msubsup><mi>E</mi><mn>4</mn><mi>down</mi></msubsup></mfrac><mo>,</mo><mfrac><msubsup><mi>STD</mi><mn>4</mn><mi>up</mi></msubsup><msubsup><mi>STD</mi><mn>4</mn><mi>down</mi></msubsup></mfrac><mo>,</mo><mfrac><msubsup><mi>E</mi><mn>5</mn><mi>up</mi></msubsup><msubsup><mi>E</mi><mn>5</mn><mi>down</mi></msubsup></mfrac><mo>,</mo><mfrac><msubsup><mi>STD</mi><mn>5</mn><mi>up</mi></msubsup><msubsup><mi>STD</mi><mn>5</mn><mi>down</mi></msubsup></mfrac><mo>,</mo><mfrac><msubsup><mi>E</mi><mn>6</mn><mi>up</mi></msubsup><msubsup><mi>E</mi><mn>6</mn><mi>down</mi></msubsup></mfrac><mo>,</mo><mfrac><msubsup><mi>STD</mi><mn>6</mn><mi>up</mi></msubsup><msubsup><mi>STD</mi><mn>6</mn><mi>down</mi></msubsup></mfrac><mo>,</mo><mfrac><msubsup><mi>E</mi><mn>7</mn><mi>up</mi></msubsup><msubsup><mi>E</mi><mn>7</mn><mi>down</mi></msubsup></mfrac><mo>,</mo><mfrac><msubsup><mi>STD</mi><mn>7</mn><mi>up</mi></msubsup><msubsup><mi>STD</mi><mn>7</mn><mi>down</mi></msubsup></mfrac><mo>]</mo></mrow>]]></math></maths>步骤19构建聚类标签向量Y＝[y1，y2，…，yj，…yN]T，其中yj＝-1或1，当yj＝-1时表示现存虹膜数据库中编号为j的虹膜样本原始图像属于聚类一，当yj＝1时表示现存虹膜数据库中编号为j的虹膜样本原始图像属于聚类二；然后将步骤8中所述的虹膜归一化图像特征向量矩阵X和聚类标签向量Y输入到支持向量机中训练，得到支持向量集合sv和判别函数f(x)，其中<mathsid="math0012"num="0012"><math><![CDATA[<mrow><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><mi>sgn</mi><mo>{</mo><munder><mi>Σ</mi><mrow><msub><mi>x</mi><mi>j</mi></msub><mo>∈</mo><mi>SV</mi></mrow></munder><msub><mi>α</mi><mi>j</mi></msub><msub><mi>y</mi><mi>j</mi></msub><mi>κ</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>j</mi></msub><mo>,</mo><mi>x</mi><mo>)</mo></mrow><mo>+</mo><mi>b</mi><mo>}</mo><mo>,</mo></mrow>]]></math>id="icf0012"file="A2009100593580004C6.tif"wi="60"he="13"top="218"left="40"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>判别函数f(x)中x表示待检测虹膜样本原始图像的归一化图像特征向量Feature_A；xj表示虹膜归一化图像特征向量矩阵X中第j行向量；κ(xj，x)表示x到xj的距离；αj表示xj的权重系数，且当xj∈SV时αj＞0，当<mathsid="math0013"num="0013"><math><![CDATA[<mrow><msub><mi>x</mi><mi>j</mi></msub><mo>∉</mo><mi>SV</mi></mrow>]]></math>id="icf0013"file="A2009100593580004C7.tif"wi="13"he="4"top="247"left="170"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>时αj＝0；b表示支持向量机的冗余度；步骤20将待检测虹膜样本原始图像的归一化图像特征向量x输入到步骤19所述的判别函数f(x)，计算判别函数f(x)的值，当f(x)＝-1时表示待检测虹膜样本原始图像属于聚类一；当f(x)＝1时表示待检测虹膜样本原始图像属于聚类二；通过步骤14至步骤20，就能够得到待检测虹膜样本原始图像在所述虹膜数据库中所属的一级聚类；若将待检测虹膜样本原始图像在所述虹膜数据库中所属的一级聚类视为新的虹膜数据库，再采取步骤19至步骤20的操作，就能够得到待检测虹膜样本原始图像在所述虹膜数据库中所属的二级聚类；以此类推，则能够得到待检测虹膜样本原始图像在所述虹膜数据库中所属的n级聚类。全文摘要一种基于核聚类的虹膜分类方法，属于图像处理
技术领域：
，涉及虹膜身份识别技术。针对虹膜数据库建立聚类时，首先将虹膜归一化图像等分成上下两层，分别对上下层做三层小波分解，选取其中7个通道下的小波系数计算能量和标准差；然后将上下两层的小波能量和标准差相除得到虹膜库样本的特征向量；将特征向量进行基于核方法的聚类，得到聚类结果。对待检测虹膜进行分类时，先将虹膜库所有特征向量输入到支持向量机中进行训练，得到支持向量和判别函数；最后将待检测虹膜样品的相应的特征向量输入判别函数中检测，得到分类结果。本发明能快速而准确的找到待检测虹膜样本所属的聚类，有利于提高虹膜识别准确率和识别效率。文档编号G06K9/00GK101556646SQ200910059358公开日2009年10月14日申请日期2009年5月20日优先权日2009年5月20日发明者梅解,韬郑申请人:电子科技大学