[发明专利]一种基于多特征阈值融合的手指静脉识别方法

权利要求书

1.一种基于多特征阈值融合的手指静脉识别方法，其特征在于包括以下步骤：

步骤1、构造多个标准差的二维高斯模板，高斯模板表达式如下：

式(1)中，x和y取值由二维高斯函数标准差σ与高斯模板中标准差因子t共同决定，二维高斯模板正方形窗口边长为p＝2×t×σ+1，因此x和y取值均为[-p,p]，σ取值可变，并且决定着窗口的大小，不同窗口大小对应提取不同粗细的静脉；

步骤2、先求出二维高斯模板的一阶偏导数h′_x(x,y)、h′_y(x,y)，二阶偏导数h″_xx(x,y)、h″_yy(x,y)，混合偏导数h″_xy(x,y),并用这5个偏导数模板分别卷积图像矩阵得到图像矩阵的5个偏导数H′_x,H′_y,H′_xx,H′_yy,H′_xy，再用方向导数公式求四个角度的一阶、二阶方向导数；坐标轴x轴方向为0°，y轴为90°；

式(2)(3)中取0°,45°,90°,135°，将方向导数代入曲率公式(4)，得到静脉图像在各方向上对应的曲率值k′，并求四个方向最大曲率值保存为静脉图像的最终曲率值，如式(5):

式(5)中k为最终提取的曲率特征，为四个方向上最大曲率值；

步骤3、改变标准差σ，选取多个标准差的提取结果进行曲率特征融合；融合公式如式(6)所示：

式(6)中k_σ1为标准差取σ₁的二维高斯模板提取的曲率结果，k_σ2为标准差取σ₂的二维高斯模板提取的曲率结果，K为两个标准差提取结果的融合结果；

步骤4、曲率灰度特征分解，按照式(7)(8)将曲率特征按照正负区域分别归一化得到背景曲率灰度特征K0′和静脉曲率灰度特征K1′；

式(7)(8)中，K是曲率灰度图像，K_min和K_max分别是将曲率灰度图像矩阵值按照从小到大和从大到小排序的两个向量；取K_min的前10个值的平均值作为负区域归一化的上限，取K_max的前10个值的平均值作为正区域归一化的上限，减轻图像边缘的曝光突变对曲率结果的影响；

步骤5、曲率细线特征提取，曲率灰度图像K经过二值化，滤波，细化得到细线特征；

步骤6、曲率灰度特征单独识别阈值计算

背景曲率灰度特征和静脉曲率灰度特征同属于曲率灰度特征，本文计算矩阵的相关系数来表示曲率灰度特征的相似度；相关系数计算公式如式(9):

式(9)中K_a,K_b为两个由曲率灰度特征矩Cov意为计算两个曲率灰度特征K_a,K_b的协方差，Var意为计算两个曲率灰度特征K_a,K_b的方差；记背景曲率灰度特征识别阈值为T₀，静脉曲率灰度特征识别阈值为T₁；

步骤7、曲率细线特征单独识别阈值计算

曲率细线特征阈值采用修正的Hausdorff距离(modified Hausdorff distance,MHD)算法计算；对于A与B两个曲率细线特征，MHD计算公式如式(10)(11):

T₂(A,B)＝max(d(A,B),d(B,A)) (11)

式(10)中d为A到B的最小平均距离，式(11)中T₂为A与B的MHD阈值；先计算A中某点A_m到B中所有点B_n的距离，其中1≤n≤N共N个点；然后求出最小值，遍历A中所有M个点并求平均值；再从B到A重复一遍，取两者较大值作为A与B间的MHD阈值；

步骤8、三种特征阈值T₀,T₁和T₂融合，获取融合后的最终阈值T；

步骤9、匹配识别：

阈值越大表示图像相似度越高，识别验证表达式如式(13):

式(13)中T_s为根据图像库计算得到的识别阈值，不同的认假率对应有不同的识别阈值，T为两幅手指静脉图像计算出的匹配的最终阈值，根据式(13)能够判断两幅手指静脉图像是否来源于同一个手指。