CN106778535B - 一种基于小波包分解的虹膜特征提取与匹配方法 - Google Patents
一种基于小波包分解的虹膜特征提取与匹配方法 Download PDFInfo
- Publication number
- CN106778535B CN106778535B CN201611065197.2A CN201611065197A CN106778535B CN 106778535 B CN106778535 B CN 106778535B CN 201611065197 A CN201611065197 A CN 201611065197A CN 106778535 B CN106778535 B CN 106778535B
- Authority
- CN
- China
- Prior art keywords
- iris
- image
- feature
- sub
- wavelet packet
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V40/00—Recognition of biometric, human-related or animal-related patterns in image or video data
- G06V40/10—Human or animal bodies, e.g. vehicle occupants or pedestrians; Body parts, e.g. hands
- G06V40/18—Eye characteristics, e.g. of the iris
- G06V40/193—Preprocessing; Feature extraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V40/00—Recognition of biometric, human-related or animal-related patterns in image or video data
- G06V40/10—Human or animal bodies, e.g. vehicle occupants or pedestrians; Body parts, e.g. hands
- G06V40/18—Eye characteristics, e.g. of the iris
- G06V40/197—Matching; Classification
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Ophthalmology & Optometry (AREA)
- Human Computer Interaction (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Multimedia (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Collating Specific Patterns (AREA)
Abstract
本发明公开了一种虹膜特征提取与匹配方法,该方法包括:将规范化虹膜图像进行分块和去噪处理,并将去噪后的所有子块进行虹膜特征提取,形成特征模板;分别对各子块图像进行一层小波包分解,提取虹膜的径向纹理特征和弧向纹理特征,形成环向特征向量;分别对各子块图像进行多层小波包分解,提取虹膜的总体能量分布特征和局部能量分布特征,形成能量特征向量;将子块的环向特征向量和能量特征向量组成子块图像的虹膜特征向量,将各个子块图像的虹膜特征向量组成虹膜图像的特征模板;对两幅虹膜图像的特征模板进行匹配识别,并基于相似度和汉明距离逐级进行判别分类。本方案在保证识别精度的基础上能够有效提高虹膜特征提取与匹配的实时性和适应性。
Description
技术领域
本发明涉及虹膜图像处理与识别领域,特别是涉及一种用于虹膜生物识别***的虹膜特征提取方法。
背景技术
虹膜识别具有准确性、稳定性、安全性和实用性等显著优势,是目前已知最先进的生物识别技术,拥有广阔的市场前景。在实际的虹膜识别***中,虹膜特征提取与匹配是虹膜识别的核心环节,只有提取到尽可能多的稳定而独立的特征点信息,并采取合理的分类判别策略,同时兼顾算法的适应性,降低对虹膜图像采集与预处理的要求,才能有效提升虹膜识别应用***的性能,兼顾识别的精度和速度。
目前在实际应用的自动虹膜识别***中,最常用的虹膜特征提取与匹配方法包括:Daugman提出的基于二维Gabor纹理滤波相位编码和Hamming距离判别的方法;Wildes提出的基于Gaussian-Laplace金字塔多分辨分解编码和Fisher线性判据的方法;谭铁牛提出的基于圆形对称空间滤波编码和最近邻分类决策的方法。这些方法虽然能提取到尽可能多的稳定的虹膜特征,但计算量大、冗余特征信息多,并且对虹膜图像成像质量以及预处理要求较高,算法复杂,适应性差,不利于满足虹膜识别***的实时性要求。
发明内容
本发明要解决的技术问题是提供一种基于小波包分解的虹膜特征提取与匹配方法,在保证识别精度的基础上解决虹膜特征提取与匹配实时性和适应性差的问题。
为解决上述技术问题,本发明采用下述技术方案:
一种基于小波包分解的虹膜特征提取与匹配方法,该方法的步骤包括:
S1、将规范化虹膜图像进行分块和去噪处理,并将去噪后的所有子块进行虹膜特征提取,形成特征模板;
S2、分别对各子块图像进行一层小波包分解,提取虹膜的径向纹理特征和弧向纹理特征,形成环向特征向量;
S3、分别对各子块图像进行多层小波包分解,提取虹膜的总体能量分布特征和局部能量分布特征,形成能量特征向量;
S4、将子块的环向特征向量和能量特征向量组成子块图像的虹膜特征向量,将各个子块图像的虹膜特征向量组成虹膜图像的特征模板;
S5、对两幅虹膜图像的特征模板进行匹配识别,并基于相似度和汉明距离逐级进行判别分类。
优选地,所述步骤S1包括:
S11、将获取的规范化虹膜图像进行分块处理;
S12、根据噪声掩膜图像剔除含噪声的图像块,并将不含噪声的图像块组成虹膜图像子块序列;
S13、分别对各子块进行虹膜特征提取,形成特征模板。
优选地,所述步骤S5包括:
S51、根据加权环向特征向量的相似度进行初筛分类,如果属于异类虹膜,则匹配失败,进入下一对虹膜特征模板的匹配,如果属于同类虹膜,则初筛成功执行步骤S52;
S52、根据加权能量特征向量的汉明距离进行判别分类。
本发明的有益效果如下:
本发明所述技术方案可以快速提取到稳定的虹膜特征,同时,通过获取不含噪声的虹膜区域,就能有效地提取到稳定的虹膜特征,包含径向和弧向纹理特征,以及总体和局部能量特征,增强了特征提取的适应性。同时利用初筛和精分相结合的特征匹配分类策略,可以快速判别异类虹膜,并充分考虑不同特征的相似度占比,有效地提高了虹膜识别的平均判别时间和准确性,从而兼顾了识别的速度和精度。
附图说明
下面结合附图对本发明的具体实施方式作进一步详细的说明;
图1示出本方案所述虹膜特征提取与匹配方法的示意图。
具体实施方式
为了更清楚地说明本发明,下面结合优选实施例和附图对本发明做进一步的说明。附图中相似的部件以相同的附图标记进行表示。本领域技术人员应当理解,下面所具体描述的内容是说明性的而非限制性的,不应以此限制本发明的保护范围。
本发明的目的在于提供一种基于小波包分解的虹膜特征提取与匹配方法,兼顾虹膜特征提取与匹配的速度和精度,提升虹膜识别的实时性和适应性。
以分辨率为480×640的标准虹膜测试图像为例,假设规范化虹膜图像的分辨率为64(行)×512(列)。一种基于小波包分解的虹膜特征提取与匹配方法的具体步骤是:
第一步,将规范化虹膜图像进行分块,并根据噪声掩膜图像剔除含噪声的图像块,将不含噪声的图像块组成虹膜图像子块序列,然后分别对各子块进行虹膜特征提取。
虹膜特征提取的输入是规范化的虹膜图像I(x,y)和噪声掩膜图像M(x,y),其中,0≤x≤m-1,0≤y≤n-1,m=64、n=512分别为虹膜图像径向和弧向的采样点数,M中值为1的点对应I中的噪声点,值为0的点对应I中的非噪声点。为了便于并行处理,提升虹膜特征提取算法的效率,对I进行均匀分块。每个子块Bk(i,j)的大小为32×64,其中,0≤i≤31,0≤j≤63,0≤k≤15,共分割成16个子块。将不含噪声的子块按列方向组成虹膜图像子块序列BS={Bks},ks对应子块索引向量IB中非1元素的索引值,给定IB(ks)=0,Bks表示非噪声子块,非噪声子块的个数num=16-∑IB。
然后,对虹膜图像子块序列BS中的每个子块Bks分别进行特征提取,最后形成特征模板FT。
第二步,分别对各子块图像进行一层小波包分解,提取虹膜的径向纹理特征和弧向纹理特征,形成环向特征向量。
一层小波包分解可以将原始图像分解成概貌子图、水平细节子图、垂直细节子图和斜向细节子图等四个子图,其中概貌子图和三方向细节子图分别反映了原始图像的总体特征分布和局部特征分布。分别对各子块图像Bks进行一层小波包分解,Bks的小波包垂直细节子图BHks反映了虹膜径向纹理特征的分布,Bks的小波包水平细节子图BVks反映了虹膜弧向纹理特征的分布,提取子块图像Bks的径向纹理特征和弧向纹理特征。Bks的径向纹理特征向量BFRks为:
式中,0≤i≤511。Bks的弧向纹理特征向量BFAks为:
式中,0≤i≤511。定义Bks的环向特征向量BFLks=[BFRks,BFAks]。
第三步,分别对各子块图像进行多层小波包分解,提取虹膜的总体能量分布特征和局部能量分布特征,形成能量特征向量。
多层小波包分解逐层在上一层的概貌子图和三方向细节子图上分别进行小波包分解,得到更精细的特征分布集合。为了提取虹膜的总体特征并进一步补充局部特征,并且考虑到算法的实现效率,分别对各子块图像Bks进行3层小波包分解,并在第3层小波包概貌子图BON ks上提取虹膜的总体能量分布特征,在第3层小波包斜向细节子图BDN ks上提取虹膜的局部能量分布特征。小波包子图的能量分布利用非零奇异值表征,而能量主要集中在前M个非零奇异值,因此采用奇异值截断变换获取能量系数。Bks的总体能量特征向量BFOks为:
式中,0≤i≤MO-1,svdcM表示MO-奇异值截断变换,λ表示非零奇异值。Bks的局部能量特征向量BFDks为:
式中,0≤i≤MD-1,svdcM表示MD-奇异值截断变换,λ表示非零奇异值。定义Bks的能量特征向量BFEks=[BFOks,BFDks]。奇异值截断个数MO和MD由下式确定:
第四步,将子块的环向特征向量和能量特征向量组成子块图像的虹膜特征向量,将各个子块图像的虹膜特征向量组成虹膜图像的特征模板。
定义Bks的虹膜特征向量FTks为:
FTks=[BFLks,BFEks]=[BFRks,BFAks,BFOks,BFDks] (6)
由此可得,虹膜图像的特征模板FT={FTks}0≤ks≤num-1。
第五步,在对两幅虹膜图像的特征模板进行匹配时,首先根据加权环向特征向量的相似度进行初筛分类,如果属于异类虹膜,则匹配失败,进入下一对虹膜特征模板的匹配;如果属于同类虹膜,则初筛成功,再根据加权能量特征向量的汉明距离进行判别分类。
为了提升虹膜识别的速度和适应性,在进行特征匹配分类时采用初筛和精分相结合的判别策略。在提取到虹膜的环向特征向量FL=[FR,FA]={BFLks}=[{BFRks},{BFAks}]后,根据加权环向特征向量的相似度SIMw进行初筛分类:
式中,SIMR表示径向特征向量的相似度,SIMA表示弧向特征向量的相似度,α=0.63是相似度权值,其中,FL1和FL2仅由子块索引向量IB1和IB2相同坐标均为0的元素对应的索引子块的环向特征向量组成。设初筛判别阈值为ta=0.65,如果SIMw<ta,则匹配失败,判定两幅虹膜图像来自于不同个体,继续输入下一幅虹膜图像进行识别,减少虹膜识别的耗时;如果SIMw≥ta,则初筛匹配成功,继续进行精分判别,增强虹膜识别的准确性。
在提取到虹膜的能量特征向量FE=[FO,FD]={BFEks}=[{BFOks},{BFDks}]后,根据加权能量特征向量的汉明距离HDw进行精确判别分类:
式中,HDO表示总体能量特征向量的汉明距离,HDD表示局部能量特征向量的汉明距离,β=0.77是汉明距离权值,其中,FE1和FE2仅由子块索引向量IB1和IB2相同坐标均为0的元素对应的索引子块的能量特征向量组成,每个子块的能量特征向量元素由前min{MO1,MO2}个总体能量特征向量的元素和前min{MD1,MD2}个局部能量特征向量的元素组成。设精分判别阈值为td=0.28,如果HDw>td,则匹配失败,判定两幅虹膜图像来自于不同个体;如果HDw≤td,则匹配成功,判定两幅虹膜图像来自于同一个体。
根据以上步骤可以从不含噪声的虹膜区域快速有效地提取到稳定的虹膜特征,包含径向和弧向纹理特征,以及总体和局部能量特征;同时利用初筛和精分相结合的特征匹配分类策略,可以快速判别异类虹膜,并充分考虑不同特征的相似度占比,有效地提高了虹膜识别的平均判别时间和准确性,从而提升了虹膜识别的实时性和适应性。
显然,本发明的上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定,对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动,这里无法对所有的实施方式予以穷举,凡是属于本发明的技术方案所引伸出的显而易见的变化或变动仍处于本发明的保护范围之列。
Claims (3)
1.一种基于小波包分解的虹膜特征提取与匹配方法,其特征在于,该方法的步骤包括:
S1、将规范化虹膜图像进行分块和去噪处理,并将去噪后的所有子块进行虹膜特征提取,形成特征模板;
S2、分别对各子块图像进行一层小波包分解,提取虹膜的径向纹理特征和弧向纹理特征,形成环向特征向量;
S3、分别对各子块图像进行多层小波包分解,提取虹膜的总体能量分布特征和局部能量分布特征,形成能量特征向量;
S4、将子块的环向特征向量和能量特征向量组成子块图像的虹膜特征向量,将各个子块图像的虹膜特征向量组成虹膜图像的特征模板;
S5、对两幅虹膜图像的特征模板进行匹配识别,并基于相似度和汉明距离逐级进行判别分类;
其中,在所述S3中,多层小波包分解逐层在上一层的概貌子图和三方向细节子图上分别进行小波包分解,得到特征分布集合,分别对各子块图像进行三层小波包分解,并在第三层小波包概貌子图上提取虹膜的总体能量分布特征,在第三层小波包斜向细节子图上提取虹膜的局部能量分布特征。
2.根据权利要求1所述的虹膜特征提取与匹配方法,其特征在于,所述步骤S1包括:
S11、将获取的规范化虹膜图像进行分块处理;
S12、根据噪声掩膜图像剔除含噪声的图像块,并将不含噪声的图像块组成虹膜图像子块序列;
S13、分别对各子块进行虹膜特征提取,形成特征模板。
3.根据权利要求1所述的虹膜特征提取与匹配方法,其特征在于,所述步骤S5包括:
S51、根据加权环向特征向量的相似度进行初筛分类,如果属于异类虹膜,则匹配失败,则进入下一对虹膜特征模板的匹配,如果属于同类虹膜,则初筛成功执行步骤S52;
S52,根据加权能量特征向量的汉明距离进行判别分类。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611065197.2A CN106778535B (zh) | 2016-11-28 | 2016-11-28 | 一种基于小波包分解的虹膜特征提取与匹配方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611065197.2A CN106778535B (zh) | 2016-11-28 | 2016-11-28 | 一种基于小波包分解的虹膜特征提取与匹配方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106778535A CN106778535A (zh) | 2017-05-31 |
CN106778535B true CN106778535B (zh) | 2020-06-02 |
Family
ID=58904492
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611065197.2A Active CN106778535B (zh) | 2016-11-28 | 2016-11-28 | 一种基于小波包分解的虹膜特征提取与匹配方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106778535B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109117880A (zh) * | 2018-08-06 | 2019-01-01 | 广东工业大学 | 基于小波包分解选择系数加权重构的瓷砖图像分类算法 |
CN112270271A (zh) * | 2020-10-31 | 2021-01-26 | 重庆商务职业学院 | 一种基于小波包分解的虹膜识别方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1426760A (zh) * | 2001-12-18 | 2003-07-02 | 中国科学院自动化研究所 | 基于活体虹膜的身份识别方法 |
CN101344913A (zh) * | 2007-07-10 | 2009-01-14 | 电子科技大学中山学院 | 一种通过提取虹膜纹理特征进行身份识别的方法 |
CN102201063A (zh) * | 2011-06-13 | 2011-09-28 | 中国科学院自动化研究所 | 基于纹理图像局部特征的形变虹膜匹配方法 |
US9412022B2 (en) * | 2012-09-06 | 2016-08-09 | Leonard Flom | Iris identification system and method |
CN105975960A (zh) * | 2016-06-16 | 2016-09-28 | 湖北润宏科技有限公司 | 基于纹理方向能量特征的虹膜识别方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100453943B1 (ko) * | 2001-12-03 | 2004-10-20 | 주식회사 세넥스테크놀로지 | 개인 식별을 위한 홍채 영상의 처리 및 인식방법과 시스템 |
CN1209732C (zh) * | 2002-10-09 | 2005-07-06 | 中国科学院自动化研究所 | 基于小波分析和过零点描述的虹膜识别方法 |
CN101540000B (zh) * | 2008-03-20 | 2011-07-27 | 中国科学院自动化研究所 | 基于纹理基元统计特性分析的虹膜分类方法 |
CN102708362A (zh) * | 2012-05-14 | 2012-10-03 | 南京光华科技发展有限公司 | 一种虹膜识别法 |
GB201209390D0 (en) * | 2012-05-28 | 2012-07-11 | Optos Plc | Improvements in or relating to image processing |
-
2016
- 2016-11-28 CN CN201611065197.2A patent/CN106778535B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1426760A (zh) * | 2001-12-18 | 2003-07-02 | 中国科学院自动化研究所 | 基于活体虹膜的身份识别方法 |
CN101344913A (zh) * | 2007-07-10 | 2009-01-14 | 电子科技大学中山学院 | 一种通过提取虹膜纹理特征进行身份识别的方法 |
CN102201063A (zh) * | 2011-06-13 | 2011-09-28 | 中国科学院自动化研究所 | 基于纹理图像局部特征的形变虹膜匹配方法 |
US9412022B2 (en) * | 2012-09-06 | 2016-08-09 | Leonard Flom | Iris identification system and method |
CN105975960A (zh) * | 2016-06-16 | 2016-09-28 | 湖北润宏科技有限公司 | 基于纹理方向能量特征的虹膜识别方法 |
Non-Patent Citations (3)
Title |
---|
Texture feature extraction using 2D Gabor Filters;Rosniza Roslan et al;《2012 International Symposium on Computer Applications and Industrial Electronics》;20130321;全文 * |
基于小波包分析的虹膜特征提取方法;王琪等;《计算机工程与应用》;20060614;第3节 * |
基于小波包分解的虹膜图像识别方法研究;刘磊等;《科技通报》;20120830;第28卷(第8期);第2-3节 * |
Also Published As
Publication number | Publication date |
---|---|
CN106778535A (zh) | 2017-05-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Casanova et al. | Plant leaf identification using Gabor wavelets | |
CN109543548A (zh) | 一种人脸识别方法、装置及存储介质 | |
Cope et al. | Plant texture classification using gabor co-occurrences | |
Doegar et al. | Cnn based image forgery detection using pre-trained alexnet model | |
CN105825183B (zh) | 基于部分遮挡图像的人脸表情识别方法 | |
CN109409384A (zh) | 基于细粒度图像的图像识别方法、装置、介质及设备 | |
CN103366180A (zh) | 一种基于自动特征学习的细胞图像分割方法 | |
CN103853724B (zh) | 多媒体数据分类方法及装置 | |
CN101996328B (zh) | 一种木材识别方法 | |
CN107135664A (zh) | 一种人脸识别的方法以及人脸识别装置 | |
CN106203356B (zh) | 一种基于卷积网络特征提取的人脸识别方法 | |
Alsmadi et al. | Fish recognition based on robust features extraction from color texture measurements using back-propagation classifier | |
CN109389585A (zh) | 一种基于全卷积神经网络的脑组织提取方法 | |
Mohamed et al. | An improved LBP algorithm for avatar face recognition | |
CN112950477A (zh) | 一种基于双路径处理的高分辨率显著性目标检测方法 | |
CN112861605A (zh) | 基于时空混合特征的多人步态识别方法 | |
CN106778535B (zh) | 一种基于小波包分解的虹膜特征提取与匹配方法 | |
Yin et al. | Dynamic difference learning with spatio-temporal correlation for deepfake video detection | |
Abaza et al. | Ear recognition: a complete system | |
Yuan et al. | A review of recent advances in ear recognition | |
CN108805183A (zh) | 一种融合局部聚合描述符和局部线性编码的图像分类方法 | |
Lawgali | An evaluation of methods for Arabic character recognition | |
Saranya et al. | An approach towards ear feature extraction for human identification | |
Li | Novel face recognition algorithm based on adaptive 3D local binary pattern features and improved Singular Value Decomposition method | |
CN113688715A (zh) | 面部表情识别方法及*** |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |