手指静脉特征提取分析与对比
梁雪慧5'b赵菲6程云泽6张瑞杰6
(天津理工大学a.天津市复杂系统控制理论与应用重点实验室;b.电气电子工程学院)
摘要针对手指静脉特征提取，首先通过对比度受限自适应直方图均衡化(CLAHE)的方法，对图像进行预处理，增强对比度，再对采用Gabor滤波器、重复线性跟踪和最大曲率法3种特征提取方法提取的静脉图像进行了对比，从结果看出Gabor滤波器效果较好%
关键词特征提取图像增强手指静脉生物识别
中图分类号TP391.41文献标识码A文章编号1000-3932(2020)03-0256-04
手指静脉识别在现代生物识别领域中具有重要的地位［1］，相比其他生物识别(如人脸、指纹、虹膜及语音等)，手指静脉识别具有不易俘获、识别数据小及安全性高等特点。
静脉识别首先要从近红外图像中获得较为完整的静脉结构，然后提取静脉完成身份识别。因此，从近红外图像中获得较为完整的静脉线结构成为至关重要的工作，静脉提取算法的优劣直接关系到匹配识别的效果，甚至对整个静脉识别系统的性能产生重大影响。因此，笔者深入研究了静脉提取的方法，以提高识别系统的性能。
现有的静脉特征提取方法中一类是基于机器学习的方法，另一类是通过研究图像的静脉纹路。基于机器学习的方法其主要思想是从原始灰度图像中直接提取特征，已有的研究有：通过减少维数，获得静脉图像主要成分分量特征的主成分分析法(PCA)［2］；利用HOL特征提取获得静脉的线性响应和方向的基于2DHOL特征与2D_2FPCA的方法［3］；通过重新定义类间离散度矩阵的双边二维线性判别分析(MB2D-LDA)方法［4］；利用流形学习的基本思想进行降维的正交邻域保持投影(ONPP&方法［5］；还有一些利用稀疏表达等方法来对静脉特征进行提取的方法［6］(研究图像的静脉纹路这类方法是将图像中静脉纹路分割出来进行特征提取，有Gabor滤波器⑺、重复线性追踪［:］及最大曲率［9］等方法。笔者主要是针对第2类的静脉特征提取方法进行研究分析。1图像预处理
1.1图像归一化
通过采集设备获取的原始手指静脉图像中由于包含了大量的背景信息，例如平移、旋转、缩放和不均匀照明，这些信息对于图像的分析完全没有意义，且会增加算法的处理时间，所以首先应提取出采集图像的感兴趣区域(ROI)，再进行尺度和灰度的归一化，提供更加准确的特征区域给后续的操作，从而提高静脉识别系统的性能。
1.2图像增强
成像设备通常采用近红外(NIR)透射照明，采集过程中多种光散射占主导地位，因此手指静脉图像的质量
始终很差，因为散射会大幅降低静脉和非静脉区域之间的对比度，直接对静脉纹理进行提取，会导致一部分静脉信息丢失，也可能会误识背景信息，产生伪静脉，这对之后的匹配和识别造成巨大的影响，因此在特征提取前要对静脉图像进行增强处理。
对比度受限自适应直方图均衡化［10］(Contrast Limited Adaptive Histogram Equalization，CLAHE)方法，最初是为医学成像而开发的，并已被证明
基金项目：天津市科技项目(18YFZCNC01120)。
作者简介：梁雪慧(1970-)，副教授，从事智能控制理论与应用的研究。
通讯作者：赵菲(1994-)，硕士研究生，从事图像处理的研究，*****************0
可以成功地增强低对比度图像，例如射野验证片(Portal Film,PF)。它实际上是改进的插值自适应直方图均衡方法，其基本思想是对比度限幅，即限制局部直方图的高度，用以限制过分放大的噪声和过度增强的局部对比度。算法步骤如下：
a.将原始图像分成互不重叠且大小相等的M x N个子区域，并计算每一个子区域的像素个数Num%和直方图&%(')。
b.求子区域的裁剪极限！，氏二[1+"()-1))・Num%/I%+其中，"为截断系数,"![0,1);)为裁剪斜率;'%为第%个子区域的灰度级数。
C.按照特定的步长对&%(')进行循环搜索，如果某一个灰度级的像素个数大于其0%，则将它裁剪到0%以内，并将剩余的像素个数分配给像素个数小于0%的灰度级，遍历所有灰度级$
%.对每个子区域对比度进行限制后，将其直方图进行均衡化处理$
e.利用双线性插值求每个像素点的灰度值$将每个子区域中心位置像素点的灰度值作为参考点$
从图1可以看出，处理后的图像暗区域内的模糊效果得到很好的增强，使得图像更加清晰，对比度更高$
图1高低质量图像增强后效果
2静脉特征提取
2.1Gabor滤波器
Gabor滤波器已被广泛用于模式分析，因为它可以有效地捕获图像中具有不同频率和方向的有用特征$Gabor滤波器的实部称为偶数对称Gabor滤波器，通常用于从生物特征图像中提取特征：
2・22
+(，,-；#,$,%，&,')=exp(---—-—)cos(2^-+%)(1)
2&2入
其中,,$(,cos$+-si;$,-$二一,si;$+-cos$,$表示法线相对于Gabor函数的平行条纹的方向,#是正弦因子的波长,%是相位偏移，&是高斯包络的标准偏差,'是空间长宽比$
假设/(,,-)是输入图像,F k(,,-)表示方向索引处对应的Gabor滤波图像之一$F k(,,-)可以这样获得：
F k(,,-)(
G k(，,-)*0(，,-)(2)
其中,231,2,3，…,3；*表示卷积运算符。
2.2重复线性跟踪
重复线性跟踪的思想是根据水平和垂直方向上的预定义概率按照选定的方向跟踪图像中的静脉，并随机选择起始种子，整个过程重复执行一定次数。
该方法基于在各个位置开始的线性跟踪。首先识别出局部暗线，并通过沿线逐像素移动来执行线性跟踪。当无法检测到黑线时，新的跟踪操作将从另一个位置开始。通过重复执行这样的局部线性跟踪操作，可以跟踪图像中的所有暗线+最终，线的位置重叠并且统计获得手指静脉的图案。
首先确定线性跟踪的起点和移动方向。
线性跟踪的起点是(，s，-s)，即当前跟踪点(，c，-c&的初始值为(，s，-S)，这是从手指轮廓内像素点的集合4•中选择的一对均匀随机数+此后，确定移动方向属性6lr、6<d$6lr和几是防止跟踪点跟随曲率过大的路径参数$久和分别确定如下：
屛(1，0)凡(2X1⑶
丨(一1,0),otherwise
((0,1),4W7(2)<1
6ud3(4)
(0，一1),otherwise
其中4,3)是0〜3之间的随机数$
其次检测暗线的方向和移动跟踪点：
a.轨迹位置表匚的初始化。跟踪点移动到的位置存储在轨迹位置表中。
b.确定当前跟踪点可以移动到的像素组化, N e=T e/$N:(,c，-c)，其中(光c，-c)是(光c，-c)的相邻像素的集合$
C.检测当前跟踪点附近的暗线方向，为了确 定当前跟踪点(，c,-c)应该移动到的像素，需要计
算线性评估函数，该函数反映了当前跟踪点周围的横截面轮廓中的谷底深度。
d.在轨迹位置表匚中记录轨迹并移动跟踪点，当前跟踪点(亿，儿)被添加到轨迹位置表中。之后，如果J为正，则将(##,$#)更新为J最大的(衝，$%)。
最后重复执行步骤b-d&
步骤d中的线性评估函数为：
0尸max{((#c+!cos!%-----sin!%,$c+rsin!%+—cos!%)+F(#c+rcos!%+—sin!%,$c+rsin!%------cos!%)-2F(#c+rcos!%,$c+rsin!%)}(5)
其中,!是轮廓的宽度，"是(#c,$J与横截面之间的距离,!%是线段(#,$〉(#c+l,$c)和(#,$〉(#%,$%)之间的角度&
2.3最大曲率法
最大曲率法主要通过搜索静脉图像横截面轮廓的曲率局部最大位置来检测中心线&
首先，提取具有各种宽度和亮度的静脉中心线，检查手指静脉图像的横截面轮廓&因为静脉比周围区域
暗，所以静脉的横截面轮廓看起来像是凹痕&这些凹曲线具有较大的曲率，因此可以通过计算横截面轮廓中的局部最大曲率来获得静脉的中心位置&
假设'为手指静脉图像,F(#,$)为对应像素点(#,$)的灰度值，将F(#,$)在位置)处的剖面图定义为*())，则曲率可以表示为：
，一吟))压2(6)
2[l+(d存())/d))2]|
其次，为了连接静脉中心点并消除噪声，检查中心点像素右侧的两个相邻像素和左侧的两个相邻像素&如果(#,$)与两侧的像素值都较大，则水平绘制一条线；如果(#，刃的值较小且两侧的像素值较大，在(#,$)处会产生间隙，因此应增加(#,$)的值以连接线路；当(#,$)具有较大的值而两侧的像素具有较小的值时，噪声点位于(#,$)，因此应减小(#，刃的值以消除噪声&以相同的方式对4个方向均进行该操作&
最后，使用阈值对静脉图案G(#,$)进行二值化&值小于阈值的像素被标记为背景，值大于或等于阈值的像素被标记为静脉区域&假设G(#,$)中值的直方图在形式上是两相的，因此确定阈值以使G(#,$)中值的组之间的离散最大&
3实验对比
将采用Gabor滤波器、重复线性跟踪和最大曲率法0种特征方法提取的静脉图像进行对比，结果如图2所示&
图2不同特征提取方法的效果
由图2可以看出,3种特征方法表现都较为优良，其中重复线性跟踪法受预先设置的参数和迭代次数影响较大，在实验中迭代的次数为200000次，仍有部分细小的静脉特征被遗漏，且运算效率低，不适用于实时静脉识别系统中。最大曲率法的静脉结构得到了增强，但是背景噪声也被突出显示了，这可能会从生成的静脉结构中产生虚假匹配&Gabor滤波器提取出来的静脉较为完整，提取的静脉结构非常适合原始图像的血管拓扑结构，并且很好地抑制了伴随的噪声&
4结束语
特征提取是静脉识别最重要的一环，提取出的特征好坏直接影响到识别效果&笔者通过对3种不同特征提取方法提取的图像进行对比，结果 表明使用Gabor滤波器进行特征提取，提取出的静脉线具有更高的完整性，但没有对运算效率进行研究，今后可以将提高特征提取算法的效率作为一个新的研究方向&
参考文献
[1]Jain A K.Technology:Biometric Recognition[J].Nature,
2007,449(7158)：38-40.
[2]Wu J D,Liu C T.Finger-vein Pattern Identification
Using Principal Component Analysis and the Neural
Network Technique[J].Expert Systems with Applic­
ations:An International Journal,2011,38(5)：5423〜
5427.
[3]张丽萍，李卫军，宁欣，等.一种基于2DHOL特征与
(2D)〜2FPCA结合的手指静脉识别方法[J].计算机
辅助设计与图形学学报,2018,30(2):254-261. [4]王贺，邓茂云，姜守坤，等.改进的双边二维线性判别
分析的手背静脉识别[J].吉林大学学报(信息科学
版),2017,35(1):32-36.
[5]Liu Z,Yin Y,Wang H,et al.Finger Vein Recognition
with Manifold Learning[J].Journal of Network and
Computer Applications,2010,33(3):275-282.
[6]Xin Y,Liu Z,Zhang H,et al.Finger Vein Verification
System Based on Sparse Representation[J].Applied
Optics,2012,51(25):6252-6258.
[7]Mohsin A,Juboori A,Bu W,et al.Palm Vein Verification
Using Gabor Filter[J].International Journal of Computer Science,2013,10(1):678-683.[8]Miura N,Nagasaka A,Miyatake T.Feature Extraction
of Finger Vein Patterns Based on Iterative Line
Tracking and Its Application to Personal Identification
[J].Systems&Computers in Japan,2004,35(7)：61-
71.
[9]Miura N,Nagasaka A,Miyatake T.Extraction of Finger­
Vein Patterns Using Maximum Curvature Points in
Image Profiles[J].IEICE Transactions on Information
and Systems,2005,E90D(8):347-350.
[10] Cromartie R,Pizer S M.Structure-sensitive Adaptive
Contrast Enhancement Methods and Their Evaluation
[J].Image and Vision Computing,1993,11(8):460-
467.
(收稿日期:2019-11-13，修回日期:2020-02-28)
Study on Finger Vein Feature Extraction
LIANG Xue-hui a'b,ZHAO Fei b,CHENG Yun-ze b,ZHANG Rui-jie b
(a.Tianjin Key Laboratory of Control Theory and Applications in Complicated Systems；
b.School of Electrical and Electronic Engineering,Tianjin University o,Technology)
Abstract In view of the finger vein feature extraction,the method of contrast-limited adaptive histogram equalization(CLAHE)was adopted to preprocess images and to enhance their contrast,and then three feature extraction methods of Gabor filter,repeated linear tracking and maximum curvature were compared to show that the Gabor filter works well.
Key words feature extraction,image enhancement (上接第204页)
裂解炉自动投料专家控制系统，保证了非平稳操作的安全性和操作一致性，为实现工业裂解炉的全工况自动化提供了一条新的有效途径。
参考文献
[1]华丰，方舟，邱彤.乙烯裂解炉反应与传热耦合的智能
混合建模与模拟[J].化工学报,2018,69(3):923-930.finger vein,biometric recognition
[2]刘俊杰，郭莹，张利军.基于Aspen Plus用户模型的
裂解炉对流段流程模拟[J].天津化工,2009,23(6)：
直方图均衡化方法25-30.
[3]张伦，张悦.基于广义预测控制的混杂系统控制[J].
计算机仿真,2017,34(8)：138-141.
(收稿日期:2019-11-22，修回日期:2020-03-23)
Expert Control System Simulation for Automatic Feeding
Process of Industrial Cracking Furnace
JIANG Ming-jing1,LU Xiang-dong1,ZHAO Zhong2,CHEN De-kun2
(1.Sinopec Zhenhai Refining&Chemical Company；
2.College of Information Science and Technology,Beijing University of Chemical Technology)
Abstract Considering the feeding operation of BA109furnace at olefin section in Sinopec Zhenhai Refining&Chemical Company and existing problems and requirements there,an expert control system was developed and simulated for the auto-feeding process of industrial cracking furnace,including the detailed description of the system)s structure,control strategy and control scheme；meanwhile,aiming at BA109 furnace*s requirements for the auto-feeding process,an expert control software was simulated.
Key words industrial cracking furnace,auto-feeding,expert control,DCS,COT

版权声明：本站内容均来自互联网，仅供演示用，请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198，我们将在24小时内删除。