本发明涉及指静脉识别领域，特别地，涉及一种指静脉特征比对方法、装置、存储介质和处理器。
背景技术：
指静脉识别是依据血液中的血红蛋白会吸收近红外线的特性而发展起来的一种新的生物识别技术。
当前生物识别技术中：指纹技术虽然很成熟，但指纹特征很容易被模仿，人脸也存在易被模仿的问题，所以他们在安全上的表现缺少依靠并有待改进，而且人脸识别技术还存在大姿态、大表情变化下结果不可靠的风险；虹膜识别技术表现稳定，但其接触式的操作方法不被广泛人群接受，不适合大规模推广；声纹识别虽然有采集、应用方面的优势，但首先个人声纹特征会随年龄、情绪、身体状况而变化，其次所采集声纹信号质量与采集设备密切相关，甚者，背景环境等噪声信号会对声纹信号形成严重的干扰；上述生物识别技术，虽然各有其优点，但使用上的局限亦不可小觑。
指静脉识别则易于采集，受环境因素影响较小，更关键的它是一种天然的活体特征，安全性上有保障，所以，近年来静脉识别技术在生物识别应用领域获得越来越多的关注。
指静脉识别技术的成熟依赖两方面的进步：采集设备和比对算法；采集设备需具备获取到清晰静脉信息的能力，比对算法则要能合理区分比对目标，给出正确的目标识别结果。显而易见，有效的静脉描述特征能为比对算法的准确性提供可靠的支撑。而如何获得有效的稳定的静脉特征则是目前指静脉识别技术的研究关键。
技术实现要素：
本发明提供了一种指静脉特征比对方法、装置、存储介质和处理器，以获得更稳定的描述特征，为指静脉识别提供保障。
本发明采用的技术方案如下：
本发明一方面提供了一种指静脉特征比对方法，包括以下步骤：
对输入的指静脉图像进行预处理；
采用Gabor函数对预处理后的图像进行处理获得Gabor矢量特征；
对获得的Gabor矢量特征进行降维处理；
对降维后的Gabor矢量特征进行二值化Gabor梯度特征提取：通过计算不同频率参数下对应方向处Gabor系数的梯度，并依据梯度值的正负符号将梯度值二值化为“0”、“1”值，然后将对应的值置于由区块坐标、频率、方向参数索引的分量位置对应的全新矢量中，生成二进制表达的全新矢量特征；
对二进制表达的全新矢量特征进行特征比对并输出比对结果。
进一步地，对输入的指静脉图像进行预处理的步骤包括：
对输入的指静脉图像进行滤波处理和手指边缘定位；
依据手指边缘位置信息对输入图像进行区域裁剪以获得有效静脉纹理区域的子图；
对裁剪后的子图做增强获得曲率增强图。
进一步地，对预处理后的图像采用Gabor滤波器进行处理获得Gabor矢量特征的步骤包括：
将曲率增强图归一化到标准模型尺寸大小；
定义Gabor函数的频率参数和方向参数，将该组参数代入Gabor函数公式，获得对应参数下的Gabor卷积滤波系数；
定义特征提取区块在水平方向和垂直方向的步进参数；
在步进参数引导下遍历曲率增强图中各区块，并在各区块内采用预先获得的Gabor卷积滤波系数进行卷积计算对应频率、方向参数下的Gabor特征；
将各区块计算获得的Gabor特征串接起来生成Gabor矢量特征。
优选地，频率参数设置为7；方向参数设置为16。
进一步地，对获得的Gabor矢量特征进行降维处理的步骤包括：
对曲率增强图中各区块按均匀分布方式进行选择，滤除部分区块特征；
对剩余区块特征按照设定方向参数进行下采样，构成新的降维后的Gabor矢量特征。
进一步地，对降维后的Gabor矢量特征进行二值化Gabor梯度特征提取的步骤包括：
计算同一频率参数下相邻两个系数的梯度，并二值化；
计算不同频率参数下对应近邻两个系数和的复合梯度值，并二值化；
通过依序计算各频率参数中单个频率下相邻系数的二值化Gabor梯度值，再组合各频率参数中近邻频率参数下相邻系数和的复合梯度二进制值，获得由区块坐标、频率参数、方向参数索引的二进制表达的全新矢量特征。
对二进制表达的全新矢量特征进行特征比对并输出比对结果的步骤包括：
对两幅识别目标图像进行二值化后的全新矢量特征中对应位置中“0”、“1”出现的相同次数进行计数，在归一化操作下转变为两幅识别目标图像的相似度值；
判断并输出比对结果，将相似度值转换为实数表达的0到1间的得分值，形成最终的比对结果得分，再将得分值与预设阈值比较，大于预设阈值则认为是同一目标，否则为异类目标。
根据本发明的另一方面，还提供了一种指静脉特征比对系统，包括：
预处理模块，用于对输入的指静脉图像进行预处理；
Gabor特征提取模块，用于对预处理后的图像采用Gabor函数进行处理获得Gabor矢量特征；
Gabor特征降维模块，用于对获得的Gabor矢量特征进行降维处理；
二值化Gabor梯度特征提取模块，用于对降维后的Gabor矢量特征进行二值化Gabor梯度特征提取：通过计算不同频率参数下对应方向处Gabor系数的梯度，并依据梯度值的正负符号将梯度值二值化为“0”、“1”值，然后将对应的值置于由区块坐标、频率、方向参数索引的分量位置对应的全新矢量中，生成二进制表达的全新矢量特征；
比对输出模块，用于对二进制表达的全新矢量特征进行特征比对并输出比对结果。
本发明还提供了一种存储介质，存储介质包括存储的程序，程序运行时控制存储介质所在的设备执行上述的指静脉特征比对方法。
本发明还提供了一种处理器，处理器用于运行程序，处理器运行时执行上述的指静脉特征比对方法。
本发明通过二值化Gabor梯度特征提取，对初期特征矢量进行再次提取，在低阶特征基础上做高阶抽象，计算开销相对增加较少，无需大规模调整系统流程，且能获得更稳定的表述特征，进一步提高指静脉识别算法的性能。
除了上面所描述的目的、特征和优点之外，本发明还有其它的目的、特征和优点。下面将参照附图，对本发明作进一步详细的说明。
附图说明
构成本申请的一部分的附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：
图1是本发明优选实施例的指静脉特征比对方法的流程图；
图2是同频率下相邻系数梯度二值化的示意图；
图3是不同频率下相邻系数复合梯度二值化的示意图；
图4是本发明优选实施例的指静脉特征比对系统的结构框图。
具体实施方式
需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本发明。
参照图1，本发明的优选实施例提供了一种指静脉特征比对方法，包括以下步骤：
步骤S100，对输入的指静脉图像进行预处理；
步骤S200，采用Gabor函数对预处理后的图像进行处理获得Gabor矢量特征；
步骤S300，对获得的Gabor矢量特征进行降维处理；
步骤S400，对降维后的Gabor矢量特征进行二值化Gabor梯度特征提取：通过计算不同频率参数下对应方向处Gabor系数的梯度，并依据梯度值的正负符号将梯度值二值化为“0”、“1”值，然后将对应的值置于由区块坐标、频率、方向参数索引的分量位置对应的全新矢量中，生成二进制表达的全新矢量特征；
步骤S500，对二进制表达的全新矢量特征进行特征比对并输出比对结果。
Gabor函数是一个线性滤波器，在空间域中，一个二维Gabor滤波器是一个由正弦平面波调制的高斯核函数。Gabor滤波器的频率和方向表达同人类视觉系统类似，它可以很好地近似单细胞的感受野函数(光强刺激下的传递函数)。
静脉识别系统比对操作中，即使是同一目标手指，因为时间、个人动作习惯、采集环境变化，设备参数调整等原因必定导致所采集图出现差异，但无论如何，该目标手指的静脉纹理肯定是大致确定的。本发明首先优选了Gabor特征做为指静脉低阶特征，然后在其基础上做高阶特征的提取，新特征平滑了上述原因导致的对应特征处的这种差异，更容易捕捉到静脉纹理的固有特性，由此获得一种更稳定的描述特征，使指静脉比对结果鲁棒性更强。
下面介绍本发明优选实施例的指静脉比对方法的具体处理过程。
对输入的指静脉图像进行预处理的步骤S100具体包括：首先对输入的指静脉图像进行滤波处理和手指边缘定位。依据手指边缘位置信息对输入图像进行区域裁剪以获得有效静脉纹理区域的子图。接着对裁剪后的子图做增强获得曲率增强图。
对预处理后的图像采用Gabor滤波器进行处理获得Gabor矢量特征的步骤S200包括：
将曲率增强图归一化到标准模型尺寸大小。
定义Gabor函数的频率参数和方向参数，将该组参数代入Gabor函数公式，获得对应参数下的Gabor卷积滤波系数。本优选实施例中，考虑到本发明最终采用的是基于Gabor之上的高阶特征，本发明将方向参数做了更精细的设置，将人们常用的8方向扩展为16方向。此外，本优选实施例中，频率参数设置为7。
定义特征提取区块在水平方向和垂直方向的步进参数。
在步进参数引导下遍历曲率增强图中各区块，并在各区块内采用预先获得的Gabor卷积滤波系数(滤波核)进行卷积计算对应频率、方向参数下的Gabor特征；
将各区块计算获得的Gabor特征串接起来生成Gabor矢量特征。
上述步骤中，定义了7个频率和16个方向Gabor参数，这样计算下来，原始Gabor系数组成的矢量特征的维度是相当可观的。本发明考虑到图像近邻像素的灰度缓慢过度分布，它们在近邻参数下(频率近邻或方向近邻)的特征一定会包含相当多的冗余信息，我们完全可以对它们做下采样，一定的稀疏密度下的采样数据足可以保全原目标的特性。
具体地，对获得的Gabor矢量特征进行降维处理的步骤S300包括：
首先对曲率增强图中各区块按均匀分布方式进行选择，滤除部分区块特征。
然后对剩余区块特征按照设定方向参数进行下采样，滤除偶数方向系数，保留奇数方向系数，经过这样的截取操作构成新的降维后的Gabor矢量特征。
为了增强特征表达的稳定性，本发明中对降维后的Gabor矢量特征进行二值化Gabor梯度特征提取的步骤S400包括：
计算同一频率参数下相邻两个系数的梯度，并二值化。参照图2，具体计算过程如下：频率参数1下，计算方向参数1与方向参数2对应Gabor系数的梯度，计算方向参数2与方向参数3对应Gabor系数的梯度，以此类推，计算方向参数16系数反转回来与方向参数1系数的梯度，遍历16个方向获得16个梯度值，然后根据各值正负符号分布，正梯度值用“1”代替，负梯度值用“0”代替，将16个梯度值转换为16位二进制值。
计算不同频率参数下对应近邻两个系数和的复合梯度值，并二值化。参照图3，具体计算过程如下：先求得频率参数1下，方向参数1对应系数与邻近方向参数2系数的和，再先求频率参数2下，方向参数1对应系数与邻近方向参数2系数的和，再计算前面两个和值相减获得复合梯度值，以此类推，遍历16个方向获得8个复合梯度值，同样低，依据梯度值的符号将梯度值转换为“0”、“1”表达的8位二进制值。
通过依序计算7个频率参数中单个频率下相邻系数的二值化Gabor梯度值，再组合7个频率参数中近邻频率参数下相邻系数和的复合梯度二进制值，获得由区块坐标、频率参数、方向参数索引的二进制表达的全新矢量特征。
对二进制表达的全新矢量特征进行特征比对并输出比对结果的步骤S500包括：
对两幅识别目标图像进行二值化后的全新矢量特征中对应位置中“0”、“1”出现的相同次数进行计数，在归一化操作下转变为两幅识别目标图像的相似度值；
判断并输出比对结果，将相似度值转换为实数表达的0到1间的得分值，形成最终的比对结果得分，再将得分值与预设阈值比较，大于预设阈值则认为是同一目标，否则为异类目标。
参照图4，本发明还提供了一种指静脉特征比对装置，包括：
预处理模块100，用于对输入的指静脉图像进行预处理；具体的预处理过程参照上文中指静脉特征比对方法中的描述。
Gabor特征提取模块200，用于对预处理后的图像采用Gabor函数进行处理获得Gabor矢量特征；具体的特征提取过程参照上文中指静脉特征比对方法中的描述。
Gabor特征降维模块300，用于对获得的Gabor矢量特征进行降维处理；具体的降维过程参照上文中指静脉特征比对方法中的描述。
二值化Gabor梯度特征提取模块400，用于对降维后的Gabor矢量特征进行二值化Gabor梯度特征提取：通过计算不同频率参数下对应方向处Gabor系数的梯度，并依据梯度值的正负符号将梯度值二值化为“0”、“1”值，然后将对应的值置于由区块坐标、频率、方向参数索引的分量位置对应的全新矢量中，生成二进制表达的全新矢量特征；具体的计算过程参照上文中指静脉特征比对方法中的描述。
比对输出模块500，用于对二进制表达的全新矢量特征进行特征比对并输出比对结果。具体的比对过程参照上文中指静脉特征比对方法中的描述。
本发明实施例还提供了一种存储介质，存储介质包括存储的程序，程序运行时控制存储介质所在的设备执行上述的指静脉特征比对方法。
本发明实施例还提供了一种处理器，处理器用于运行程序，处理器运行时执行上述的指静脉特征比对方法。
本发明先利用函数变换获得静脉纹理的分块Gabor特征，再依序将分块Gabor特征组合成对应指静脉图像的识别特征矢量，与常见直接采用Gabor矢量特征进行目标比对分类不同，本发明对Gabor矢量特征进行再次加工，对不同参数区间的Gabor值，将频率和尺度(方向)对应下的分量值求取梯度，并依据梯度值正负结果转换为“0”、“1”二值化数据(如：A分量大于B分量，则置结果为“1”，否则置“0”)，且将该“0/1”值放置在与该分量索引位置对应的全新矢量中，最后获得的是二进制表达的全新特征矢量；因为全新特征矢量记录的不是单纯的Gabor特征，而是Gabor特征分量间的梯度信息，所以，本发明将指静脉成像中由于设备、环境、采集操作不规范等不可控因素所带来的干扰大大压缩了，呈现的特征更能捕捉目标的真实特性，也必然大大提升特征表达的稳定性。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。