首页 > 专利 > 杭州电子科技大学 > 一种基于多特征阈值融合的手指静脉识别方法专利详情

一种基于多特征阈值融合的手指静脉识别方法   0    0

有效专利 查看PDF
专利申请流程有哪些步骤?
专利申请流程图
申请
申请号:指国家知识产权局受理一件专利申请时给予该专利申请的一个标示号码。唯一性原则。
申请日:提出专利申请之日。
2018-01-08
申请公布
申请公布指发明专利申请经初步审查合格后,自申请日(或优先权日)起18个月期满时的公布或根据申请人的请求提前进行的公布。
申请公布号:专利申请过程中,在尚未取得专利授权之前,国家专利局《专利公报》公开专利时的编号。
申请公布日:申请公开的日期,即在专利公报上予以公开的日期。
2018-07-31
授权
授权指对发明专利申请经实质审查没有发现驳回理由,授予发明专利权;或对实用新型或外观设计专利申请经初步审查没有发现驳回理由,授予实用新型专利权或外观设计专利权。
2020-04-07
预估到期
发明专利权的期限为二十年,实用新型专利权期限为十年,外观设计专利权期限为十五年,均自申请日起计算。专利届满后法律终止保护。
2038-01-08
基本信息
有效性 有效专利 专利类型 发明专利
申请号 CN201810015658.8 申请日 2018-01-08
公开/公告号 CN108256456B 公开/公告日 2020-04-07
授权日 2020-04-07 预估到期日 2038-01-08
申请年 2018年 公开/公告年 2020年
缴费截止日
分类号 G06K9/00G06K9/38G06K9/40G06K9/46G06K9/62 主分类号 G06K9/00
是否联合申请 独立申请 文献类型号 B
独权数量 1 从权数量 2
权利要求数量 3 非专利引证数量 0
引用专利数量 4 被引证专利数量 0
非专利引证
引用专利 CN106611168A、CN104239769A、CN105184272A、WO2011052085A1 被引证专利
专利权维持 5 专利申请国编码 CN
专利事件 事务标签 公开、实质审查、授权
申请人信息
申请人 第一申请人
专利权人 杭州电子科技大学 当前专利权人 杭州电子科技大学
发明人 沈雷、蓝师伟、李凡、吕葛梁、杨航 第一发明人 沈雷
地址 浙江省杭州市下沙高教园区2号大街 邮编 310018
申请人数量 1 发明人数量 5
申请人所在省 浙江省 申请人所在市 浙江省杭州市
代理人信息
代理机构
专利代理机构是经省专利管理局审核,国家知识产权局批准设立,可以接受委托人的委托,在委托权限范围内以委托人的名义办理专利申请或其他专利事务的服务机构。
杭州君度专利代理事务所 代理人
专利代理师是代理他人进行专利申请和办理其他专利事务,取得一定资格的人。
朱月芬
摘要
本发明公开了一种基于多特征阈值融合的手指静脉识别方法。本发明先用基于多标准差的二维高斯模板计算手指静脉图像曲率提取图像的背景区域曲率灰度特征,静脉区域曲率灰度特征和曲率细线特征。用相关系数法计算背景区域曲率灰度特征和静脉区域曲率灰度特征的匹配阈值,用MHD算法计算曲率细线特征的匹配阈值,再提出一种多特征阈值融合判决算法,融合以上三种特征的单独识别阈值进行验证判决。本发明提出的基于三种特征阈值融合的算法,由于结合了背景区域曲率灰度特征包含的背景区域可用信息,对模糊的手指静脉图像识别性能下降不明显,因此在低认假率下的拒真率显著低于传统的基于单独细线特征的识别算法,能够合理且高效利用的有效的识别算法。
  • 摘要附图
    一种基于多特征阈值融合的手指静脉识别方法
法律状态
序号 法律状态公告日 法律状态 法律状态信息
1 2020-04-07 授权
2 2018-07-31 实质审查的生效 IPC(主分类): G06K 9/00 专利申请号: 201810015658.8 申请日: 2018.01.08
3 2018-07-06 公开
权利要求
权利要求书是申请文件最核心的部分,是申请人向国家申请保护他的发明创造及划定保护范围的文件。
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为标准差取σ1的二维高斯模板提取的曲率结果,kσ2为标准差取σ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)中Ka,Kb为两个由曲率灰度特征矩Cov意为计算两个曲率灰度特征Ka,Kb的协方差,Var意为计算两个曲率灰度特征Ka,Kb的方差;记背景曲率灰度特征识别阈值为T0,静脉曲率灰度特征识别阈值为T1;
步骤7、曲率细线特征单独识别阈值计算
曲率细线特征阈值采用修正的Hausdorff距离(modified Hausdorff distance,MHD)算法计算;对于A与B两个曲率细线特征,MHD计算公式如式(10)(11):
T2(A,B)=max(d(A,B),d(B,A))         (11)
式(10)中d为A到B的最小平均距离,式(11)中T2为A与B的MHD阈值;先计算A中某点Am到B中所有点Bn的距离,其中1≤n≤N共N个点;然后求出最小值,遍历A中所有M个点并求平均值;再从B到A重复一遍,取两者较大值作为A与B间的MHD阈值;
步骤8、三种特征阈值T0,T1和T2融合,获取融合后的最终阈值T;
步骤9、匹配识别:
阈值越大表示图像相似度越高,识别验证表达式如式(13):
式(13)中Ts为根据图像库计算得到的识别阈值,不同的认假率对应有不同的识别阈值,T为两幅手指静脉图像计算出的匹配的最终阈值,根据式(13)能够判断两幅手指静脉图像是否来源于同一个手指。

2.根据权利要求1所述的一种基于多特征阈值融合的手指静脉识别方法,其特征在于步骤8所述的融合的实现,具体如下:
由于曲率灰度特征和曲率细线特征是不同维度的特征,其阈值计算方式不同,因此三种阈值融合顺序是先进行两区域曲率灰度特征的阈值融合,得到二者最佳的权重α及新阈值,再与曲率细线特征阈值融合得到权重β及最终阈值T,表达式如式(12):
T=(1-β)(αT0+(1-α)T1)+βT2         (12)
式(12)中T0,T1,T2分别为背景曲率灰度特征,静脉曲率灰度特征和曲率细线特征的阈值,α为阈值T0,T1的融合权重,β为阈值T2与T0,T1融合阈值的融合权重。

3.根据权利要求1所述的一种基于多特征阈值融合的手指静脉识别方法,其特征在于步骤曲率细线特征提取,具体实现如下:
二值化处理采用局部动态阈值分割算法,利用窗口内的均值和方差得到阈值,对曲率灰度图K进行二值化,采用连通区域滤波去噪,再用查表细化法得到曲率细线特征。
说明书

技术领域

[0001] 本发明属于生物特征识别及信息安全技术领域,特别涉及到手指静脉识别技术领域。

背景技术

[0002] 手指静脉识别技术属于生物特征识别技术,它利用手指内部静脉拓扑结构进行人体身份验证。因其识别速度快、性能好以及特征安全不易伪造而成为国内外诸多学者的研究热点。手指静脉识别系统主要包括手指静脉图像采集,图像预处理,特征提取和特征匹配识别。考虑到目前采集设备的工艺水平,存在着一些人的手指在采集过程中受到设备不稳定或者采集环境的影响,导致采集到的手指静脉图像存在静脉纹路模糊,静脉区域不稳定等问题。因此从手指静脉图像中提取出稳定的可用特征以及如何利用特征进行识别算法设计是手指静脉识别技术的关键环节。
[0003] 管凤旭等人运用主成分分析改进算法提取特征,识别率较高,但其针对静脉图像特点改进的算法复杂度高,提取时间长。王科俊等人先确定静脉方向,再用该方向的模板提取细线特征。提取效率高,且抑制了其它方向模板提取的噪声。但其设计的模板对低质量图像提取效果不佳。林剑等人使用高斯模板二阶导数卷积图像得到图像的Hessian矩阵,再根据Hessian矩阵的迹和特征值得到细线特征,并通过改变高斯模板方差来针对不同粗细的静脉获得最佳的提取效果。高斯模板比方向模板对低质量图像的适应性更好,提取结果稳定。孙晓琳根据图像中垂直静脉纹路方向呈谷形分布的特点,计算多方向上离散曲率,根据曲率值提取静脉中心点。但其曲率计算采用离散算法,提取的特征存在较多噪声。Naoto Miura等人根据曲率值对静脉中心进行跟踪寻找,准确提取静脉中心点,但其跟踪算法复杂度较高。由于曲率综合了一阶导数和二阶导数信息,对图像增强效果比二阶导数好。目前的手指静脉识别算法中,大多数使用单个特征进行识别,且从静脉拓扑结构来设计算法时只提取细线特征或者端点和交叉点特征进行识别,而手指静脉纹路模糊的图像在特征提取环节如果只提取细线特征,会使特征包含的信息少且不稳定,当数据库中静脉纹路模糊图像累计到一定数量会造成不同类手指静脉图像之间相似度提高,使图像模糊的用户识别困难或者出现误识。

发明内容

[0004] 对于现有的算法如基于方向滤波提取细线特征的算法和基于Hessian矩阵提取细线特征的算法对模糊图像提取特征不稳定导致的在低认假率情况下拒真率低的问题,本发明提供一种基于多特征阈值融合的手指静脉识别算法。
[0005] 本发明所采用的技术方案为:
[0006] 一种基于多特征阈值融合的手指静脉识别算法,包括以下步骤:
[0007] S1、构造多个标准差的二维高斯模板,高斯模板表达式如下:
[0008]
[0009] 式(1)中x和y取值由二维高斯函数标准差σ与高斯模板中标准差因子t共同决定,二维高斯模板正方形窗口边长为p=2×t×σ+1。因此x和y取值均为[-p,p]。σ取值可变,并且决定着窗口的大小,不同窗口大小对应提取不同粗细的静脉。
[0010] S2、先求出二维高斯模板的一阶偏导数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,再用方向导数公式(2)、(3)求四个角度的一阶、二阶方向导数。坐标轴x轴方向为0°,y轴为90°。
[0011]
[0012]
[0013] 式(2)(3)中 取0°,45°,90°,135°。将方向导数代入曲率公式(4),得到静脉图像在各方向上对应的曲率值。并求四个方向最大曲率值保存为静脉图像的最终曲率值,如式(5):
[0014]
[0015]
[0016] 式(5)中k为最终提取的曲率特征,为四个方向上最大曲率值。
[0017] S3、改变标准差σ,选取合适的多个标准差的提取结果进行曲率特征融合。实验表明选取两个标准差提取结果进行融合效果最好。融合公式如式(6)所示:
[0018]
[0019] 式(6)中kσ1为标准差取σ1的二维高斯模板提取的曲率结果,kσ2为标准差取σ2的二维高斯模板提取的曲率结果,K为两个标准差提取结果的融合结果。
[0020] S4、曲率灰度特征分解,按照式(7)(8)将曲率特征按照正负区域分别归一化得到背景曲率灰度特征K0和静脉曲率灰度特征K1。
[0021]
[0022]
[0023] 式(7)(8)中,K是曲率灰度图像,K_min和K_max分别是将曲率灰度图像矩阵值按照从小到大和从大到小排序的两个向量。取K_min的前10个值的平均值作为负区域归一化的上限,取K_max的前10个值的平均值作为正区域归一化的上限,减轻图像边缘的曝光突变对曲率结果的影响。背景曲率灰度特征和静脉曲率灰度特征是本发明提取的两种用于融合识别的特征。
[0024] S5、曲率细线特征提取,曲率细线特征是静脉图像中静脉区域的最精简的特征。曲率灰度图像K经过二值化,滤波,细化得到细线特征。本文二值化处理采用局部动态阈值分割算法,利用窗口内的均值和方差得到阈值,对曲率灰度图K进行二值化,采用连通区域滤波去噪,再用查表细化法得到曲率细线特征。
[0025] S6、曲率灰度特征单独识别阈值计算
[0026] 背景曲率灰度特征和静脉曲率灰度特征同属于曲率灰度特征,本文计算矩阵的相关系数来表示曲率灰度特征的相似度。相关系数计算公式如式(9):
[0027]
[0028] 式(9)中Ka,Kb为两个由曲率灰度特征矩Cov意为计算两个曲率灰度特征Ka,Kb的协方差,Var意为计算两个曲率灰度特征Ka,Kb的方差。记背景曲率灰度特征识别阈值为T0,静脉曲率灰度特征识别阈值为T1。
[0029] S7、曲率细线特征单独识别阈值计算
[0030] 曲率细线特征阈值采用修正的Hausdorff距离(modified Hausdorff distance,MHD)算法计算。对于A与B两个曲率细线特征,MHD计算公式如式(10)、(11):
[0031]
[0032] T2(A,B)=max(d(A,B),d(B,A))   (11)
[0033] 式(10)中d为A到B的最小平均距离,式(11)中T2为A与B的MHD阈值。先计算A中某点Am到B中所有点Bn的距离,其中1≤n≤N共N个点。然后求出最小值,遍历A中所有M个点并求平均值。再从B到A重复一遍,取两者较大值作为A与B间的MHD阈值。
[0034] S8、三种特征阈值融合识别算法
[0035] 由于曲率灰度特征和曲率细线特征是不同维度的特征,其阈值计算方式不同,因此三种阈值融合顺序是先进行两区域曲率灰度特征的阈值融合,得到二者最佳的权重α及新阈值,再与曲率细线特征阈值融合得到权重β及最终阈值T,表达式如式(12):
[0036] T=(1-β)(αT0+(1-α)T1)+βT2   (12)
[0037] 式(12)中T0,T1,T2分别为背景曲率灰度特征,静脉曲率灰度特征和曲率细线特征的归一化百分制阈值,α为阈值T0,T1的融合权重,β为阈值T2与T0,T1融合阈值的融合权重。
[0038] S9、匹配识别。本文中定义阈值越高表示图像相似度越高,识别验证表达式如式(13):
[0039]
[0040] 式(13)中Ts为根据图像库计算得到的识别阈值,不同的认假率对应有不同的识别阈值,T为两幅手指静脉图像计算出的匹配阈值,根据式(13)可以判断两幅手指静脉图像是否来源于同一个手指。
[0041] 本发明有益效果如下:
[0042] 本发明提出一种基于多特征阈值融合的手指静脉识别算法,首先利用多标准差的二维高斯模板提取手指静脉图像的静脉区域曲率灰度特征,背景区域曲率灰度特征及曲率细线特征并且采用阈值融合算法确定各特征阈值的最佳融合权重融合三种特征进行识别。实验表明,本发明提出的阈值融合算法综合了手指静脉图像不同维度的特征所包含的可用信息,特别是融合了背景区域的可用信息后,拒真率比传统的基于单独细线特征的识别算法在低认假率下有较大提升,说明本发明提出的多特征阈值融合算法是一种有效的手指静脉识别算法。
[0043] 本发明提取的手指静脉图像的三种特征从不同维度上刻画了手指静脉图像中的不同程度可用信息,而所提阈值融合算法融合了三种特征的可用信息进行验证判决,识别性能明显高于传统的基于单独细线特征识别的算法。特别地,对于手指静脉纹路模糊图像,传统算法提取的细线特征包含的可用信息少且不稳定,在低认假率下识别性能明显降低,而本发明提出的基于三种特征阈值融合的算法,由于结合了背景区域曲率灰度特征包含的背景区域可用信息,对模糊的手指静脉图像识别性能下降不明显,因此在低认假率下的拒真率显著低于传统的基于单独细线特征的识别算法。因此,基于多特征阈值融合的手指静脉识别算法是一种对手指静脉图像信息合理且高效利用的有效的识别算法。

附图说明

[0044] 图1为一幅手指静脉原始灰度图像;
[0045] 图2为0°方向的模板提取出的曲率灰度图像;
[0046] 图3为45°方向的模板提取出的曲率灰度图像;
[0047] 图4为90°方向的模板提取出的曲率灰度图像;
[0048] 图5为135°方向的模板提取出的曲率灰度图像;
[0049] 图6为标准差1.5的模板提取出的曲率灰度图像;
[0050] 图7为标准差2.5的模板提取出的曲率灰度图像;
[0051] 图8为标准差1.5和2.5的模板提取的曲率灰度图像融合图像;
[0052] 图9为融合曲率灰度图像分解出的静脉区域曲率灰度特征图像;
[0053] 图10为融合曲率灰度图像分解出的背景区域曲率灰度特征图像;
[0054] 图11为融合曲率灰度图像提取出的曲率细线特征图像;
[0055] 图12为基于本发明提取的三种特征单独识别算法ROC曲线和其他论文中的基于方向滤波细线特征识别算法和基于Hessian矩阵细线特征识别算法ROC曲线对比图;
[0056] 图13为静脉曲率灰度特征和背景曲率灰度特征融合识别算法的拒真率随阈值融合权重α的变化曲线;
[0057] 图14为基于静脉曲率灰度特征融合背景曲率灰度特征识别算法的ROC曲线和基于单独静脉曲率灰度特征识别算法的ROC曲线对比图;
[0058] 图15为三种特征阈值融合识别算法的拒真率随阈值融合权重β的变化曲线;
[0059] 图16为三种特征阈值融合的识别算法ROC曲线和基于单独静脉曲率灰度特征识别算法的ROC曲线对比图;具体实施例
[0060] 下面结合附图对本发明的具体实施例作进一步的说明。
[0061] 本实施例的基于多特征阈值融合的手指静脉识别方法,包括以下步骤:
[0062] S1、式(1)中σ取1.5,t取4,模板窗口边长即为p=13,构造成一个13×13的正方形二维高斯模板。
[0063] S2、分别求出二维高斯函数x方向一阶导数h′x(x,y)、y方向的一阶导数h′y(x,y)、x方向二阶偏导数h″xx(x,y)、y方向二阶偏导数h″yy(x,y)和二阶混合偏导数h″xy(x,y)。
[0064] S3、将模板窗口坐标代入步骤S2中的各导数表达式得到模板的5个偏导数。
[0065] S4、用模板的5个偏导数卷积手指静脉原始图像,原始图像如图1所示,卷积过程中,原始图像矩阵扩充采用边界值复制法进行扩充。卷积得到的结果为图像矩阵的5个偏导数。
[0066] S5、将步骤S4得到的图像矩阵的5个偏导数代入方向导数公式(2)(3),得到图像矩阵的4个方向导数。
[0067] S6、将步骤S5得到的4个方向导数逐个代入曲率计算公式即式(4)得到图像矩阵在四个方向上的曲率特征图,如图2-5所示。
[0068] S7、将步骤S6得到的四个方向上的曲率代入式(5),求得图像矩阵的四个方向上最大曲率图,如图6所示。图6的曲率灰度图像是二维高斯模板方差为1.5时的提取结果,由图6可见,方差小时提取的曲率灰度图像中静脉细节更多,但是粗静脉轮廓并不能完整提取出来。
[0069] S8、式(1)中σ取2.5,方差因子t取4,模板窗口边长即为p=21,构造成一个21×21的正方形二维高斯模板,并重复步骤S2-S7,得到如图7所示的方差为2.5的曲率灰度图像提取结果。由图7可见,与图6相比,大方差的提取结果中粗静脉轮廓信息更加完整,但是细节静脉纹路变得不清晰。
[0070] S9、将S8与S7得到的标准差为1.5和2.5的曲率灰度特征图按照式(6)进行融合,融合结果如图8所示。图8中所提取的曲率灰度大的区域对应于手指静脉图像的静脉区域,可见融合了两个标准差的结果与单个标准差提取结果相比较,即包含了粗静脉的完整轮廓,又包含了较多的细节静脉纹路,取得了较好的融合效果。
[0071] S10、将融合的曲率灰度特征按照式(7)(8)进行分解,由式(8)得到静脉区域曲率灰度特征如图9所示,由式(7)得到背景区域曲率特征如图10所示。图9和图10中白色区域即为曲率值较大的区域,经过式(7)的归一化处理,背景区域曲率灰度特征的曲率值范围也为0-255,归一化处理是为了方便后续步骤的阈值计算和阈值融合。
[0072] S11、将融合的曲率灰度特征经过局部动态阈值分割二值化,连通区域滤波,查表细化得到曲率细线特征如图11所示。由图11可见,曲率细线特征是反映一个手指静脉图像静脉拓扑结构的最简单的特征,只包含了静脉区域的单像素主干信息。
[0073] S12、阈值计算和匹配验证。将两幅手指静脉图像均执行完步骤S1-S11得到两幅手指静脉图像的三对特征即静脉曲率灰度特征,背景曲率灰度特征和曲率细线特征。静脉曲率灰度特征和背景曲率灰度特征按照式(9)进行阈值计算得到阈值T1和T0,曲率细线特征按照式(10)(11)进行阈值计算得到阈值T2,再按照阈值融合公式即式(12)得到融合阈值T,并通过比较当前比对阈值与由图像库计算出的某个认假率下的识别阈值Ts来判断两幅手指静脉图像是否来自同一根手指。
[0074] 以下是基于图像数据库对本发明所提算法的实验数据分析,确定阈值权重α与β。并将本发明所提算法的性能与其他手指静脉识别算法的性能进行比较的实验结果。
[0075] 本次算法仿真图像是本实验室研发的静脉采集设备采集的500类指静脉图像,每类图像有3幅图,图像原始尺寸为500×220,为了节约图像处理时间,先将图像采用双线性插值归一化为160×64。采集时被采集人员正常使用设备,一根手指使用三次,设备自动保存三幅图像。采集过程时间跨度长,采集环境变化复杂,被采集人员数量多,符合目前大部分的实际应用场景,因此该图像库所得结论具有实际应用意义。算法程序在MATLAB R2014a上编写和运行,仿真电脑内存4G,主频率3.4GHZ,操作系统为64位Windows7。认假率(False Acceptance Rate,FAR)为不同手指分别采集的指静脉图像,在进行1:1匹配时被判定为同一手指所占的比例;拒真率(False Rejection Rate,FRR)为同一手指分别采集的指静脉图像,在进行1:1匹配时被判定为不同手指所占的比例。本次仿真实验1:1计算的不同手指类间图像匹配次数为1122750次,同一手指类内图像匹配次数为1500次。
[0076] 通过比较不同特征的FAR-FRR曲线即ROC(receiver operator characteristic curve,ROC)曲线来进行识别性能的比较,ROC曲线靠近坐标轴,性能越好。本发明对所提出的基于多特征阈值融合的静脉识别算法做了仿真,同时为了比较,对基于方向滤波的细线特征识别算法和基于Hessian矩阵的细线特征识别算法做了仿真。
[0077] 图12是本文所提取三种特征单独识别算法和其他的识别算法ROC曲线。由图12可知,本发明提出的基于静脉曲率灰度特征识别算法和基于曲率细线特征识别算法比基于方向滤波细线特征识别算法和基于Hessian矩阵细线特征识别算法具有更好的性能。这是由于用二维高斯模板计算曲率对图像有更好的增强作用。由于静脉曲率灰度特征包含的有用信息比单像素曲率细线特征多,所以静脉曲率灰度特征性能略高于曲率细线特征。由于背景曲率灰度特征只包含了少量的静脉图像可用信息,所以背景曲率灰度特征性能低于静脉曲率灰度特征和曲率细线特征,在融合识别算法中起辅助作用。
[0078] 按照式(12),先进行两区域曲率灰度特征的阈值融合,图13是在FAR为10-6时两区域曲率灰度特征阈值融合后识别曲线。由图13可知α在20%-30%内FRR提升幅度最大。结合其它FAR时的性能情况,权重α选取为23%,在FAR为10-6时FRR从24.19%降低至19.73%,可见融合识别性能提高了,在低认假率下两区域曲率灰度特征阈值融合后拒真率明显降低。图14是两区域曲率灰度特征阈值融合识别算法的ROC曲线,由图14可见,经过23%的背景曲率灰度特征阈值加权后,两区域曲率灰度特征融合识别算法的ROC曲线更贴近坐标轴,说明融合后识别性能提升了。融合后的新阈值称为两曲率灰度特征融合阈值。
[0079] 再将两曲率灰度特征融合阈值与曲率细线特征阈值进行融合,确定权重β。图15是FAR为10-6时的最终融合FRR随权重β的变化曲线,由图15可见,β在40%左右FRR最低,结合其他FAR时的性能情况,选择β为40%。在FAR为10-6时三种特征阈值融合的FRR从19.73%降低到了8.13%。图16是三种特征阈值融合识别算法的ROC曲线与静脉曲率灰度特征单独识别算法的ROC曲线对比。由图16可知阈值融合后识别性能提升显著,特别是在认假率为10-6时,拒真率从单静脉曲率灰度特征识别算法的24.19%提升到基于多特征阈值融合识别算法的8.13%,说明低认假率下识别性能提升明显。最终的融合阈值为T=T0×0.138+T1×0.462+T2×0.4,各阈值融合权重符合理论分析。这是因为静脉曲率灰度特征和曲率细线特征包含了静脉图像的主要信息,在阈值融合识别算法中起到主要作用,所以阈值T1,T2在T中所占比例高达46.2%和40%;而背景曲率灰度特征一般只有静脉图像的少量信息,只有在静脉区域模糊的低质量图像中才会起主要识别作用,而在基于整个图像库计算性能的融合识别时只起辅助用,所以阈值T0在T中所占比例只有13.8%。本发明实验所采用的图像库是模拟了采集设备实际使用环境采集而来的,因此该权重对基于本实验室研发的指静脉采集设备的手指静脉识别系统具有适用性。
[0080] 上面结合附图对本发明的实施例作了详细说明,但是本发明并不限于上述实施例,在本领域普通技术人员所具备的知识范围内,还可以在不脱离本发明宗旨的前提下作出得各种变化,也应视为本发明的保护范围。
版权所有:盲专网 ©2023 zlpt.xyz  蜀ICP备2023003576号