CN101789075A - 基于特征值归一化双向加权的手指静脉识别方法 - Google Patents

Abstract

本发明提供的是基于特征值归一化双向加权的手指静脉识别方法。(1)通过图像采集装置进行手指静脉图像的采集；(2)对采集的手指静脉图像进行预处理，所述预处理包括：彩色图像进行灰度化、手指区域提取、采用组合滤波器分别消除椒盐噪声和高斯噪声、采用局部动态阈值算法分割图像并二值化、然后采用面积消除法去噪、按照手指轮廓标记提取手指静脉脉络图像，最后将图像的大小标准化为统一的图像；(3)通过在行列两个方向上都加权的二维主成分分析算法提取手指静脉特征；(4)通过最近邻分类器进行匹配与识别。本发明用于手指静脉身份识别系统。明显地提高手指静脉的识别速度，识别率稳定、而且高。

Description

基于特征值归一化双向加权的手指静脉识别方法

技术领域

本发明涉及一种生物特征身份识别技术，是一种手指静脉特征身份认证技术。

背景技术

生物特征识别技术(Biometric Identification Technology)是指利用人体生物特征进行身份认证的一种技术。与传统的方法不同在于，生物特征识别方法依据的是我们人体本身所拥有的个体特性。生物特征分为基于身体特征和基于行为特点两类。身体特征包括：指纹、掌型、眼睛(视网膜和虹膜视网膜)、人体气味、脸形、皮肤毛孔、手的血管纹理和DNA等；行为特点包括：签名、语音、行走的步态、击键的力度等。

目前较为成熟的生物特征识别系统绝大多数都是基于单模态的，在实际应用中还存在一些问题，影响其使用的可靠性，对这些问题的解决，如仅仅从算法和采集设备上来进行改进，成本很高或者根本无法解决。

在众多生物特征识别技术中，由于指纹识别具有很强的普遍性、唯一性、稳定性、易用性等特点，受到广泛的关注、研究和应用。指纹识别技术是目前较为理想的身份确认技术，实用程度很高，但在实际应用中发现存在一些瓶颈问题：

(1)由于手指有污物、过湿、过干，或指纹仪不清洁，或采集指纹的照相机聚焦不准等原因，导致采集的指纹图像不够清晰，引起图像质量问题，直接影响指纹的识别精度和结果；

(2)根据NIST(美国国家标准技术研究所)的报告指出，由于手指受伤(疤痕、磨损)，或手指爆皮大约有2％的人，不能提供好质量的待验和注册图像，因此这些人就不能通过指纹来进行识别。

(3)指纹采集过程中，手指按压的扭转和拉伸，按压力度等因素会使指纹产生形变、旋转，导致特征的描述缺乏恒定性，会使注册的特征值的描述和验证时的特征值描述不同，这样成功比对的可能性就会降低，这类情况在现实中是比较常见的。因此设计验证系统时，必须将这些因素考虑在内，这样系统的可靠性会变差和复杂度会增加，会增加拒识率。

(4)虽然很难窃取生物特征，但是这种可能性是存在的。目前已经出现了利用硅树脂制造的克隆指纹。

指纹识别方式还仅停留在观察人的生物特征的“表象”上，安全系数较低。近年，出现了一种新的生物特征识别技术——静脉识别。现在主要的研究有手背静脉识别、手掌静脉识别以及手指静脉识别。同指纹一样，静脉也具有很强的普遍性和唯一性，而且具有指纹无法比拟的优势：

(1)在获取静脉图像时，可以利用透射光或者反射光两种方式进行获取。这里我们采用透视光获取手指静脉图像，因为该法穿透手指获取内部静脉图像特征，这样可以有效避免反射光获取图像时，因皮肤表面的皱纹、褶皱、粗糙、干裂或潮湿等造成的获取精确图像特征的障碍。

(2)对于指纹识别装置而言，只要被识别对象的纹路符合识别规则，那么便被识别装置判断为正确，识别过程即告完成。因此被识别对象的可仿造程度很高，可以是贴有正确纹路的人造手指或真人手指，也可以是伪造好纹路的各种无生命载体。静脉识别装置的关键在于“活体识别”，也就是说被识别对象必须是活着的人，才能达到“可识别”的第一步。很难伪造或是手术改变。

(3)静脉是身体内部的血管特征，很难伪造或者是手术改变，是非接触性的信息采集，手指无须与仪器接触，不会造成采集界面的污染，没有手指接触设备时的不卫生、及特征可能被复制所带来的安全问题，并避免了被当作审查对象的心理不适。

很明显，手指静脉识别克服了指纹识别等的诸多缺点，有着广阔的应用前景，因此手指静脉识别是近几年来生物认证技术开辟的新领域。

手指静脉识别研究最先始于日立公司，日立公司Information & Telecom研究小组针对手指静脉认证算法的作了较多的研究。2000年，日立公司工程师等首次提出手指静脉用于个人身份的鉴定的方法。2004年，日立公司发表了NaotoMiura等几位工程师关于提取手指静脉特征的成果。2005年，日立公司成立了全球指静脉识别技术推广中心，开发内容包括ATM、电脑、楼宇系统、汽车安全等等领域。据报道，在日本已有上万家企业采用了手指静脉识别技术的安全保密系统，东京各大银行采用了ATM及柜台身份认证、保险箱及电脑系统的登录；目前，该技术已经登陆我国。

在国内关于手指静脉的研究和资助的项目相对比较少。中国民航大学天津市智能信号与图像处理重点实验室杨金锋博士获得了国家自然科学基金资助，这是目前国内手指静脉识别技术研究领域第一项获得国家自然科学基金资助的研究项目。哈尔滨工程大学模式识别与智能系统实验室的王科俊教授对手指静脉识别进行了研究，制作了采集设备和提出了相应的手指静脉提取和匹配方法。吉林大学张忠波等在指静脉识别算法上也进行了研究。另外还有其他一些学者进行研究，文献对手指静脉采集装置进行了设计。不仅手指静脉图像预处理作了一些深入而细致的研究。也分别采用不用方法提取手指静脉特征，均取得令人满意的效果。

发明内容

本发明的目的在于提供一种能提高手指静脉的识别速度，识别率稳定、而且高的基于特征值归一化双向加权的手指静脉识别方法。

本发明的目的是这样实现的：

(1)通过图像采集装置进行手指静脉图像的采集；

(2)对采集的手指静脉图像进行预处理，所述预处理包括：彩色图像进行灰度化、手指区域提取、采用组合滤波器分别消除椒盐噪声和高斯噪声、采用局部动态阈值算法分割图像并二值化、然后采用面积消除法去噪、按照手指轮廓标记提取手指静脉脉络图像，最后将图像的大小标准化为统一的图像；

(3)通过在行列两个方向上都加权的二维主成分分析算法提取手指静脉特征；

(4)通过最近邻分类器进行匹配与识别。

所述通过在行列两个方向上都加权的二维主成分分析算法提取手指静脉特征的方法为：首先采用二维主成分分析方法，分别计算训练样本A在行方向上和列方向上的特征值和对应的特征向量组；其次将两组特征值分别按从大到小排序，并计算特征值累积贡献率，分别获取行和列方向上的最优维数d和t；然后按照最优维数d和t，分别选择其特征值所对应的特征向量组，构成行和列方向上的最佳投影矩阵X和B^T；按照行和列方向上的特征值λ和λ′的大小，将X和B^T分别按照行和列两个方向上的进行特征值归一化加权处理，得到投影加权矩阵X_W和B^T _W；最后将训练样本和测试样本分别在行和列方向上的加权矩阵X_W和B^T _W进行投影，获得训练样本和测试样本的图像特征矩阵C_W＝B^T _WAX_W。

行方向上特征向量组X的加权方法为：

对于行方向通过在行列两个方向上都加权的二维主成分分析算法提取手指静脉特征的方法中前d个最大的特征值λ₁≥λ₂≥...≥λ_d，采用累积特征值贡献率方式，获取行方向上满足贡献率为percet％的特征值λ_k作为其他所有特征值归1化分母，其行方向上累积特征值贡献率：

$\frac{\underset{k=1}{\overset{d}{\mathrm{\Σ}}}{\mathrm{\λ}}_k}{\underset{i=1}{\overset{M}{\mathrm{\Σ}}}{\mathrm{\λ}}_i}\≥\mathrm{percet}\%$

特征值归1化公式为：λ^* _i＝λ_i/λ_k，i＝1，2，…，d

其特征向量组X的加权策略如下：

$\left\{\begin{array}{c}{X}_W=X\&times;\mathrm{diag}({{\mathrm{\&lambda;}}^{*}}_1,{{\mathrm{\&lambda;}}^{*}}_2,...,{{\mathrm{\&lambda;}}^{*}}_d)\\ {{\mathrm{\&lambda;}}^{*}}_i=\left\{\begin{array}{cc}{{\mathrm{\&lambda;}}^{*}}_i& {{\mathrm{\&lambda;}}^{*}}_i<1\\ {{\mathrm{\&lambda;}}^{*}}_i^{\mathrm{\&omega;}}& {{\mathrm{\&lambda;}}^{*}}_i\&GreaterEqual;1\end{array}\right.\\ \begin{array}{cc}i=\mathrm{1,2},...,d,& 0<\mathrm{\&omega;}<1\end{array}\end{array}\right.$

其中diag(λ^* ₁，λ^* ₂，...，λ^* _d)为经过归1化的新特征值λ^* ₁，λ^* ₂，...，λ^* _d组成的对角阵，ω为加权因子；

大于1的特征值经加权因子ω加权后，减小特征值，对于等于和小于1的特征值则没有变化。

本发明提出的的手指静脉识别方法，包括手指静脉图像获取、图像预处理、获取尺寸标准化的手指静脉脉络图像、特征提取、身份识别等步骤。

本发明的主要贡献和特点在于：(1)手指静脉脉络是隐藏在皮肤下面，只能通过红外图像采集装置才能获取手指静脉图像。(2)对采集的手指静脉图像进行预处理，包括：彩色图像进行灰度化、手指区域提取、采用组合滤波器分别消除椒盐噪声和高斯噪声、采用局部动态阈值算法分割图像并二值化、然后采用面积消除法去噪、按照手指轮廓标记提取手指静脉图像，最后将图像的大小标准化为统一的图像。通过预处理，可以有效的去除各种噪声，以获取准确的手指静脉脉络图像。(3)对尺寸标准化的脉络图像采取(2D)²PCA方法，并结合累积特征值贡献率的分析，可以极大的降低特征矩阵的维数，进而提高识别速度；而且在保留了手指静脉结构特征的同时，从而使得最终识别精度得到显著提高。(4)基于行列两个方向上的特征值的特性，对手指静脉图像特征矩阵进行行列两个方向上加权的方法，可以滤掉更多的冗余信息，从而获得更加稳定的识别率，从而明显提高特征矩阵的维数对手指静脉识别的鲁棒性。

附图说明

图1手指静脉识别算法的流程图；

图2手指静脉图像预处理，其中图2(a)为原始图像、图2(b)为滤波、图2(c)为分割、图2(d)为去噪、图2(e)为静脉提取与归一化；

图3在2DPCA和(2D)²PCA方法下累积特征贡献率与识别率关系图；

图4在(2D)²PCA和(OW2D)²PCA方法下累积特征贡献率与识别率关系图。

具体实施方式

下面结合附图举例对本发明做更详细地描述：

1.手指静脉红外图像的采集

选用的手指静脉采集装置的基本原理是利用近红外线照射手指，并由图像传感器感应手指透射过来的光线。其中的关键在于流到静脉红血球中的血红蛋白因照射会失脱氧份，而这个还原的血红蛋白对波长760nm附近的近红外线会有吸收，本实施方式中采用的是850nm的红外光源，导致静脉部分的透射较少，在成像设备上就会产生静脉图案。所以，选用的手指静脉采集器就是利用透射的近红外线的强弱来凸现静脉血管的。

2.手指静脉图像的预处理

为了提取手指静脉脉络，首先将原始手指静脉图像进行灰度变换(如图2(a))；通过判别连通区域面积的大小确定手指轮廓并标记，去除漏光干扰区域；采用组合滤波器分别消除椒盐噪声和高斯噪声(如图2(b))；采用局部动态阈值算法分割图像并二值化(如图2(c))，然后采用面积消除法去噪(如图2(d))；按照手指轮廓标记提取手指静脉图像，最后将图像的大小归一化为60×150的图像(如图2(e))。

3.行列相结合二维线性判别分析

3.1 2DPCA

设模式类别有w个，每类有样本d个，训练样本总数为M＝w×d，每个样本大小为m×n的图像A_k(k＝1，2，…，M)。

x是一个n维的列向量，考虑将图像样本矩阵A向方向为x的空间上投影，y就是A到方向为x的空间上的投影，y为一个m维的列向量，如下式表示：

y＝Ax                          (1)

实际上，最重要的工作是选择好x的方向，x方向的不同，将使图像样本投影后的可分离程度不同，从而直接影响到图像识别的效果。因此，最主要的任务就是寻找最优的投影向量x，使得所有样本投影到x后，其投影样本的总体散布矩阵达到最大。投影样本的总体散布矩阵可以用投影特征向量的协方差矩阵的迹来表述。从这个观点来看，可以采用下面的准则：

J(x)＝trace(S_x)                (2)

其中，S_x表示训练样本A到方向为x空间上的投影后得到特征向量y的协方差，trace(S_x)表示S_x的迹。协方差矩阵S_x表示如下：

S_x＝E[(y-Ey)(y-Ey)^T]＝E{[Ax-E(Ax)][Ax-E(Ax)]^T}         (3)

＝E[x^T(A-EA)^T(A-EA)x]

由于矩阵的迹是一个数的性质，可有

trace(S_x)＝x^TE[(A-EA)^T(A-EA)]x                         (4)

已知图像矩阵的协方差矩阵为：

$G=E\left[{(A-\mathrm{EA})}^T(A-\mathrm{EA})\right]=\frac{1}{M}\underset{k=1}{\overset{M}{\mathrm{\&Sigma;}}}{(A_k-\stackrel{\&OverBar;}{A})}^T(A_k-\stackrel{\&OverBar;}{A})---\left(5\right)$

其中

表示所有训练样本的平均图像。

因此，(2)式中准则可表示成：

J(x)＝trace(x^TGx)                                      (6)

最大化该准则J(x)的物理意义是寻找到一个最优投影轴x，这就意味着投影样本的总体散布矩阵在图像矩阵集合均投影到x上后被最大化了。这样，最优投影轴x也就是对应于G的最大特征值的特征向量。一般情况下，只有一个最优投影轴是远远不够的，所以需要选取一组投影向量x₁，x₂，...，x_d构成二维最佳投影矩阵，使得J(x)取极大值。当然，这组投影轴要满足相互正交的约束条件，同时也要最大化准则J(x)，也即：

$\left\{\begin{array}{c}\{x_1,x_2,...,x_d\}=\arg \max J\left(x\right)\\ x_i^T{x}_j=0,i\&NotEqual;j,i,j=\mathrm{1,2},...,d\end{array}\right.---\left(7\right)$

也就是说，使J(x)取极大值的前d(d＜n)个最大的特征值所对应的相互正交特征向量组就是最佳投影矩阵，令X＝[x₁，x₂，…，x_d]，则有：

Y＝[y₁，y₂，…，y_d]＝[Ax₁，Ax₂，…，Ax_d]AX                      (8)

Y就是尺寸为m×n维图像矩阵A在最佳投影矩阵X投影后得到的一个尺寸为m×d维投影特征矩阵。

3.2列方向2DPCA

上述2DPCA实际是图像矩阵A寻找行方向上的最佳投影矩阵X，同理在列方向上也可以寻找一个最佳投影矩阵。

假设A′是一个尺寸为m×d的图像，x′^T是一个m维的行向量，将A′向方向为x′^T的空间上投影，得到一个d维行向量y′，如下式表示：

y′＝x′^TA′                                (9)

其图像矩阵A′的协方差矩阵估计值G′：

$G^{\&prime;}=\frac{1}{M}\underset{k=1}{\overset{M}{\mathrm{\&Sigma;}}}(A_k-\stackrel{\&OverBar;}{A}){(A_k-\stackrel{\&OverBar;}{A})}^T---\left(10\right)$

最大化准则为：

J(x′^T)＝trace(x′^TG′x′)                  (11)

同样，需要选取一组向量构成最佳投影矩阵，使得J(x′^T)取极大值，这组最佳投影轴为：

$\left\{\begin{array}{c}{{\{x^{\&prime;}}_1,{x^{\&prime;}}_2,...,{x^{\&prime;}}_t\}}^T=\arg \max J\left(x^{\&prime;T}\right)\\ {x^{\&prime;}}_i^T{x^{\&prime;}}_j=0,i\&NotEqual;j,i,j=\mathrm{1,2},...,t\end{array}\right.---\left(12\right)$

也就是说，使J(x′^T)取极大值的前t(t＜m)个最大的特征值所对应的特征向量组就是最佳投影矩阵，令B＝[x′₁，x′₂，…，x′_t]，则有：

Y′＝[y′₁，y′₂，…，y′_t]＝[x′₁ ^TA′，x′₂ ^TA′，…，x′_t ^TA′]＝B^TA′         (13)

Y′就是尺寸为m×d维图像矩阵A′在最佳投影矩阵B^T投影后得到的一个尺寸为t×d维投影特征矩阵。

3.3行列方向相结合的2DPCA

行方向上的2DPCA和列方向上的2DPCA分别只在图像的行方向和列方向上发挥作用。行2DPCA能从一组训练图像中生成一个反应图像行之间信息的最优矩阵X，相似地，列2DPCA也能生成一个反应图像列之间信息的最优矩阵B^T，然后将尺寸为m×n的矩阵A同时投影到X和B^T上，产生一个大小为t×d的矩阵C：

C＝B^TAX                           (14)

矩阵C称为特征矩阵，当把每幅图像A_k投影到X和B^T后，可以得到图像的特征矩阵C_k(k＝1，2，...，M)。

4.特征值归一的双向加权(2D)²PCA

尽管(2D)²PCA算法虽然减少了图像特征矩阵的维数，但(2D)²PCA仍然平等对待了每一维特征。然而不同的特征值所对应的特征向量对识别的贡献是不同的，大的特征值所对应的特征向量对识别贡献较大，因此，本文针对行列两个方向上特征投影矩阵分别进行加权，以突出不同特征值所对应的特征向量对识别的贡献，提出了行列两个方向上的特征值归1化并加权(2D)²PCA算法((OW2D)²PCA)。

对于行方向2DPCA算法中式(7)中前d个最大的特征值λ₁≥λ₂≥...≥λ_d，经实验检验，特征值λ_i均大于1，为实现本文加权策略，特采用累积特征值贡献率方式，获取行方向上满足贡献率为percet％的特征值λ_k作为其他所有特征值归1化分母，其行方向上累积特征值贡献率：

$\frac{\underset{k=1}{\overset{d}{\mathrm{\&Sigma;}}}{\mathrm{\&lambda;}}_k}{\underset{i=1}{\overset{M}{\mathrm{\&Sigma;}}}{\mathrm{\&lambda;}}_i}\&GreaterEqual;\mathrm{percet}\%---\left(15\right)$

特征值归1化公式为：

λ^* _i＝λ_i/λ_k，i＝1，2，…，d                     (16)

其特征向量组X的加权策略如下：

$\left\{\begin{array}{c}{X}_W=X\&times;\mathrm{diag}({{\mathrm{\&lambda;}}^{*}}_1,{{\mathrm{\&lambda;}}^{*}}_2,...,{{\mathrm{\&lambda;}}^{*}}_d)\\ {{\mathrm{\&lambda;}}^{*}}_i=\left\{\begin{array}{cc}{{\mathrm{\&lambda;}}^{*}}_i& {{\mathrm{\&lambda;}}^{*}}_i<1\\ {{\mathrm{\&lambda;}}^{*}}_i^{\mathrm{\&omega;}}& {{\mathrm{\&lambda;}}^{*}}_i\&GreaterEqual;1\end{array}\right.\\ \begin{array}{cc}i=\mathrm{1,2},...,d,& 0<\mathrm{\&omega;}<1\end{array}\end{array}\right.---\left(17\right)$

其中diag(λ^* ₁，λ^* ₂，...，λ^* _d)为经过归1化的新特征值λ^* ₁，λ^* ₂，...，λ^* _d组成的对角阵，ω为加权因子。

当λ≥1时，加权因子ω可以适当地减小权值对特征向量的放大作用。λ＜1时由于权值将缩小特征向量的作用，因此在λ＜1的条件下，不对特征向量进行加权。

同理针对列方向2DPCA算法中式(12)中前t个最大的特征值进行归1化处理，得到λ′₁≥λ′₂≥...≥λ′_t，其特征向量组B^T的加权策略如下：

$\left\{\begin{array}{c}{B^T}_W={(B\&times;\mathrm{diag}({{\mathrm{\&lambda;}}^{\&prime;}}_1,{{\mathrm{\&lambda;}}^{\&prime;}}_2,...,{{\mathrm{\&lambda;}}^{\&prime;}}_t))}^T\\ {{\mathrm{\&lambda;}}^{\&prime;}}_j=\left\{\begin{array}{cc}{{\mathrm{\&lambda;}}^{\&prime;}}_j& {{\mathrm{\&lambda;}}^{\&prime;}}_j<1\\ {{\mathrm{\&lambda;}}^{\&prime;}}_j^{\mathrm{\&omega;}}& {{\mathrm{\&lambda;}}^{\&prime;}}_j\&GreaterEqual;1\end{array}\right.\\ \begin{array}{cc}j=\mathrm{1,2},...,t,& 0<\mathrm{\&omega;}<1\end{array}\end{array}\right.---\left(18\right)$

于是，式(14)可以重新表示为：

C_W＝(B×diag(λ′_j))^T×A×(X×diag(λ^* _i))＝B^T _WAX_W          (19)

其中i＝1，2，…，d；j＝1，2，…，t

C_W就是图像矩阵A经过行和列两个方向上特征值归1化并加权的(2D)²PCA后的图像特征矩阵。

5.通过最近邻分类器进行识别

最近邻分类器：

$d(Y_{\mathrm{ij}},Y)=\underset{c=1}{\overset{t}{\mathrm{\&Sigma;}}}\underset{r=1}{\overset{s}{\mathrm{\&Sigma;}}}{\left|\right|{\left(A_{\mathrm{ij}}\right)}_{\mathrm{rc}}{P}_{\mathrm{rc}}{Q}_{\mathrm{rc}}-A_{\mathrm{rc}}{P}_{\mathrm{rc}}{Q}_{\mathrm{rc}}\left|\right|}_F=\underset{c=1}{\overset{t}{\mathrm{\&Sigma;}}}\underset{r=1}{\overset{s}{\mathrm{\&Sigma;}}}\sqrt{\mathrm{trace}{W}^{T}W}---\left(17\right)$

其中，W＝(A_ij)_rcP_rcQ_rc-A_rcP_rcQ_rc

如果那么Y属于C类.

6.实验

为了验证本方法的识别性能，通过自行研制的手指静脉红外图像采集装置，以手指(考虑使用者的方便，主要是采集每个人的食指和中指)为单位，一共采集132个手指的静脉图像，每个手指采集5次，一共采集了132×5＝660幅静脉图像，构成手指静脉数据库。

选取每类中4幅图像作为训练样本，剩余1幅作为测试样本，表1为PCA、2DPCA、(2D)²PCA几种算法，在1:n识别模式下的执行效率，从中可以最重要信息就是PCA的训练时间太长，竟然都超过6分钟，要远远高于其他2种算法的训练时间。这主要是由于PCA算法将二维图像矩阵展开为一维向量，而导致的计算中要处理的协方差矩阵的维数远高于二维矩阵的协方差矩阵，造成特征提取的时间过长。另外二维图像数据运算2DPCA和(2D)²PCA方法中，而采用(2D)²PCA的训练时间略低，因此采用(2D)²PCA方法可以提高识别速度。

表1.不同算法的识别时间

算法	维数	总训练时间(s)	总测试时间(s)	识别率(％)
					PCA	向量维数4000	414.329	16.015	92.42
2DPCA	矩阵维数40×100	13.156	4.062	92.42
					(2D)2PCA	矩阵维数40×100	12.469	3.812	92.42

其次，将经过图像预处理并标准化手指静脉图像，按照一定比例降低分辨率，分别处理成图像尺寸为80×200、60×150、50×125、40×100、20×50、10×25，在累积特征贡献率为90％条件下，研究图像分辨率对特征提取与识别的影响。选取每类中4幅图像作为训练样本，剩余1幅作为测试样本，训练时间和测试时间均定义为一个样本进行数据处理的平均时间，得到1:n的识别模式下的识别率如表2所示。

表2.不同分辨率图像的识别率

通过从表2可知，手指静脉图像一旦确定，随着图像尺寸逐渐变小，其相应的图像质量也随之逐渐变得越差，但训练出来的图像特征矩阵的表征性并不一定随之下降。因而，当图像大小分别为80×200、60×150、50×125和40×100的手指静脉图像识别率相差不多，但从识别率、识别时间和图像质量等几方面因素进行综合考虑，选择图像大小为60×150作为后续研究对象。

为了不受某个特定的累积特征贡献率对识别率的影响，基于2DPCA和(2D)2PCA两种算法，针对上述两种尺寸的图像进行对比研究。选择累积特征贡献率在30％～100％变化，选取每类中4幅图像作为训练样本，剩余1幅作为测试样本，得到1:n的识别模式下的识别率如图3所示，平均识别率如表3所示。

表3.不同算法的识别率

算法	最优识别率	平均识别率％
			2DPCA	93.18	92.13
(2D)2PCA	95.45	93.46

从图3和表3中可以看出，(2D)²PCA的平均识别率略高于2DPCA，尤其当累积特征贡献率在50％～70％之间，(2D)²PCA算法的识别率比2DPCA平均高出2个百分点，但是当累积特征贡献率超过70％之后，(2D)²PCA算法的识别率不稳定。

接下来，为研究加权对(OW2D)²PCA的影响，针对(OW2D)²PCA(归一行列两方向加权)算法，在没有加权值以及权值为0.05～1之间变化进行对比研究。选择累积特征贡献率选择80％，每类中4幅图像作为训练样本，剩余1幅作为测试样本，得到1:n的识别模式下的识别率如表4所示。

表4权重与识别率关系

权重	0	0.05	0.1	0.15	0.2	0.25	0.3
								识别率	0.932	0.955	0.947	0.947	0.947	0.947	0.947
权重	0.35	0.4	0.45	0.5	0.55	0.6	0.65
								识别率	0.939	0.939	0.939	0.939	0.939	0.939	0.939
权重	0.7	0.75	0.8	0.85	0.9	0.95	1
								识别率	0.939	0.939	0.939	0.939	0.939	0.939	0.939

发现当加权值在0.05时，识别率会落在波峰处，比没有加权值的同类型算法高出2个多百分比。

为了进一步了解加权值对(2D)²PCA的影响，选取权值ω＝0.05时，针对60×150尺寸图像，研究累积特征贡献率在30％～100％之间变化时对识别率的影响。选取每类中4幅图像作为训练样本，剩余1幅作为测试样本，得到1:n的识别模式下的识别率如图4所示，平均识别率如表4所示。

表5.权重值取0.05时2种算法的识别率

算法	最优识别率	平均识别率％
			(2D)2PCA	95.45	93.46
(OW2D)2PCA	95.45	94.34

从图4和表4可以看出，当累积特征贡献率较高(约超过70％)时，(OW2D)²PCA明显取得更加稳定，而且识别率高于(2D)²PCA。这说明(OW2D)²PCA对(2D)²PCA提取特征向量中的冗余信息有很强的抑制作用，(OW2D)²PCA整体平均识别率高于(2D)²PCA。

PCA是将图像矩阵转化为高维向量，导致计算量非常大，而2DPCA和(2D)²PCA直接基于图像矩阵进行PCA，从而越过将图像矩阵转化为向量的过程。这样，在取得效率的同时，避免了高维数带来的运算困难。其中2DPCA仅在行方向上进行特征提取的问题，特征提取不全面。(2D)²PCA以从行列两个方向上降低图像特征的矩阵维数，使得识别速度进一步加快。然而(2D)²PCA提取特征向量中存在大量冗余信息，对识别率影响很大，如果不能很好的选择降维比例，识别率不仅不高，而且不稳定。

因此本发明采用行列两个方向上的特征值归一化加权的二维主成分分析(OW2D)²PCA的处理方法，不仅可以很容易确定(2D)²PCA的降维比例，从而提高识别速度，而且使得识别率高、而且稳定。

Claims (3)

1.一种基于特征值归一化双向加权的手指静脉识别方法，其特征是：

(1)通过图像采集装置进行手指静脉图像的采集；

(2)对采集的手指静脉图像进行预处理，所述预处理包括：彩色图像进行灰度化、手指区域提取、采用组合滤波器分别消除椒盐噪声和高斯噪声、采用局部动态阈值算法分割图像并二值化、然后采用面积消除法去噪、按照手指轮廓标记提取手指静脉脉络图像，最后将图像的大小标准化为统一的图像；

(3)通过在行列两个方向上都加权的二维主成分分析算法提取手指静脉特征；

(4)通过最近邻分类器进行匹配与识别。

2.根据权利要求1所述的基于特征值归一化双向加权的手指静脉识别方法，其特征是所述通过在行列两个方向上都加权的二维主成分分析算法提取手指静脉特征的方法为：首先采用二维主成分分析方法，分别计算训练样本A在行方向上和列方向上的特征值和对应的特征向量组；其次将两组特征值分别按从大到小排序，并计算特征值累积贡献率，分别获取行和列方向上的最优维数d和t；然后按照最优维数d和t，分别选择其特征值所对应的特征向量组，构成行和列方向上的最佳投影矩阵X和B^T；按照行和列方向上的特征值λ和λ′的大小，将X和B^T分别按照行和列两个方向上的进行特征值归一化加权处理，得到投影加权矩阵X_W和B^T _W；最后将训练样本和测试样本分别在行和列方向上的加权矩阵X_W和B^T _W进行投影，获得训练样本和测试样本的图像特征矩阵C_W＝B^T _WAX_W。

3.根据权利要求2所述的基于特征值归一化双向加权的手指静脉识别方法，其特征是行方向上特征向量组X的加权方法为：对于行方向通过在行列两个方向上都加权的二维主成分分析算法提取手指静脉特征的方法中前d个最大的特征值λ₁≥λ₂≥...≥λ_d，采用累积特征值贡献率方式，获取行方向上满足贡献率为percet％的特征值λ_k作为其他所有特征值归1化分母，其行方向上累积特征值贡献率：

$\frac{\underset{k=1}{\overset{d}{\mathrm{\&Sigma;}}}{\mathrm{\&lambda;}}_k}{\underset{i=1}{\overset{M}{\mathrm{\&Sigma;}}}{\mathrm{\&lambda;}}_i}\&GreaterEqual;\mathrm{percet}\%$

特征值归1化公式为：λ^* _i＝λ_i/λ_k，i＝1，2，…，d

其特征向量组X的加权策略如下：

$\left\{\begin{array}{c}{X}_W=X\&times;\mathrm{diag}({{\mathrm{\&lambda;}}^{*}}_1,{{\mathrm{\&lambda;}}^{*}}_2,...,{{\mathrm{\&lambda;}}^{*}}_d)\\ {{\mathrm{\&lambda;}}^{*}}_i=\left\{\begin{array}{cc}{{\mathrm{\&lambda;}}^{*}}_i& {{\mathrm{\&lambda;}}^{*}}_i<1\\ {{\mathrm{\&lambda;}}^{*}}_i^{\mathrm{\&omega;}}& {{\mathrm{\&lambda;}}^{*}}_i\&GreaterEqual;1\end{array}\right.\\ i=\mathrm{1,2},...d,0<\mathrm{\&omega;}<1\end{array}\right.$

其中diag(λ^* ₁，λ^* ₂，...，λ^* _d)为经过归1化的新特征值λ^* ₁，λ^* ₂，..，λ^* _d组成的对角阵，ω为加权因子；

大于1的特征值经加权因子ω加权后，减小特征值，对于等于和小于1的特征值则没有变化。

Images