CN111144331B - 手肘静脉图像肘正中静脉的识别方法及手肘图像采集装置 - Google Patents

Abstract

本发明公开了手肘静脉图像肘正中静脉的识别方法及手肘图像采集装置，涉及图像识别技术领域，能通过对手肘图像的处理来准确找出肘正中静脉，首先对采集到的手肘图像进行增强，获得高对比度的手肘静脉图像；然后对所述高对比度的手肘静脉图像，进行多尺度自适应滤波获得静脉结构图像；然后对所述静脉结构图像进行阈值分割，在其二值图像中提取静脉特征；最后将所述静脉特征输入到训练好的分类器中完成肘正中静脉的识别。

Description

手肘静脉图像肘正中静脉的识别方法及手肘图像采集装置

技术领域

本发明涉及图像识别技术领域，具体涉及手肘静脉图像肘正中静脉的识别方法及手肘图像采集装置。

背景技术

手肘位置的肘正中静脉，是日常生活中用于静脉采血和注射药物的常用静脉。由于不同的患者皮下脂肪厚度不同，静脉粗壮程度不一，不易找到用于注射的肘正中静脉，容易造成静脉穿刺失败。穿刺失败一方面会给患者带来不良的影响，多次穿刺增加患者的疼痛和焦虑，甚至产生恐惧心理，另一方面会使医护人员产生沮丧，担忧以及缺乏自信等负面情绪。同时，治疗环境拥挤嘈杂，治疗时间受到限制，亦会影响医护人员操作难度。

目前，现有的静脉显像设备和方法，仅仅使医护人员更加容易观察到静脉，没有准确的找到肘正中静脉，这对缺乏经验的医护人员，在吵闹复杂的环境下快速的进行静脉穿刺，依然是一项严峻的挑战。再者，此类设备的售价一般较高，难以普及。

发明内容

本发明是为了解决在手肘静脉穿刺时不易于找到肘正中静脉的问题，提供一种能通过对手肘图像的处理来准确找出肘正中静脉，可靠性高的肘静脉图像肘正中静脉的识别方法及手肘图像采集装置。

以上技术问题是通过下列技术方案解决的：

一种手肘静脉图像肘正中静脉的识别方法，识别方法如下：

对采集到的手肘图像进行增强，获得高对比度的手肘静脉图像；

对所述高对比度的手肘静脉图像，进行多尺度自适应滤波获得静脉结构图像；

对所述静脉结构图像进行阈值分割，在其二值图像中提取静脉特征；

将所述静脉特征输入到训练好的分类器中完成肘正中静脉的识别。

作为优先，所述对采集到的静脉图像进行增强的方法如下：

对采集到的静脉图像进行分块，对每个小块计算直方图，设置裁剪阈值；

对每一灰度级上超过裁剪阈值的直方图进行裁剪，然后均匀的分布到每一灰度级；

对裁剪完成的直方图做直方图均衡化；

将均衡化后的分块图像拼接成完整图像，对左上，左下，右下，右上四个1/4个小块采用直接赋值的方式，对上，左，下，右四条边采用线性插值，图像中剩余部分采用双线性插值方式；所述四条边不包括上述四个小块；

作为优先，所述进行多尺度自适应滤波获得静脉结构图像的方法如下：

根据图像二阶hessian矩阵特征值λ₁，λ₂，|λ₁|≤|λ₂|构建的多尺度自适应滤波函数，如公式(6)所示：

其中

β＝0.5，ω＝50，d由输入图像自适应确定，其形式为公式(7)所示，d_rank为输入图像中像素个数不为0的灰度级，σ_i为当前迭代的尺度因子，滤波器输出0≤v(x，y；σ)≤1；

作为优先，所述提取静脉特征的方法如下：

对所述二值图像中的所有轮廓进行查找，获得轮廓之后，通过求解每个轮廓的最小外包矩形；

然后通过所述最小外包矩形确定每段静脉包括手臂轮廓的宽高及面积，用于表征该段静脉的大小；

然后利用所述最小外包矩形的宽高比来表征这段静脉的弯曲程度；

接着求取每段静脉最小外包矩形的中心距离两侧手臂轮廓的距离；具体方式为通过最小外包矩形的中心向两侧轮廓的最小外包矩形的最外侧边做垂线，分别为左侧距离和右侧距离，用于表征该段静脉的相对位置；

然后将所述最小外包矩形映射回到还未进行阈值分割时的静脉结构图像，

并计算矩形包围的像素点灰度值之和，用于表征该段静脉易于被观察的程度；

至此每段静脉获得7维特征，将部分特征进行重组，再添加左侧距离与右侧距离之比、右侧距离与左侧距离之比和灰度强度与面积之比，这样使得每段静脉一共获得10维特征。

作为优先，所述分类器包括：随机森林分类器和支持向量机分类器；

对每段静脉添加相应的标签，标注其是否为用于注射的肘正中静脉，将10维特征和标签用于训练分类器；

把采集到的若干幅手肘静脉图像分别独立训练每个分类器后，使每个分类器输出相应类别的概率值，通过软投票的方式将其集成为最终的预测器，完成对肘正中静脉的识别。

一种手肘图像采集装置，包括镜头朝下布置的摄像头，在位于摄像头的正下方设有滤光片，在位于滤光片的正下方设有能放置手肘的肘托，在摄像头的左右两侧分别设有朝下照射的红外光源，在每个红外光源的正下方都分别设有散光片。

作为优先，两个散光片所处位置都高于镜头所处位置，并且两个散光片位于同一高度处。

作为优先，位于摄像头左右两侧方的红外光源位于同一高度处。

作为优先，在位于摄像头的镜头正下方的肘托上设有放肘槽。

本发明能够达到如下效果：

本发明能通过对手肘图像的处理来准确找出肘正中静脉，可靠性高。

附图说明

图1所示为本发明的总体流程图。

图2为手肘图像采集装置的原理示意图，

具体实施方式

下面结合附图与实施例对本发明作进一步的说明。

实施例1，肘静脉图像肘正中静脉的识别方法，参见图1-2所示，包括手肘图像采集模块、图像预处理模块、肘静脉特征提取模块、肘正中静脉预测模块。图像预处理模块包括手肘图像增强、静脉结果提取、阈值分割和静脉分段。

完成手肘图像采集之后，进行手肘图像处理和识别，其识别方法如下：

对采集到的手肘图像进行灰度化增强，获得单通道高对比度的手肘静脉图像；

再对所述高对比度的手肘静脉图像，进行多尺度自适应滤波获得静脉结构图像；

对所述静脉结构图像进行阈值分割，在其二值图像中提取静脉特征；

对静脉结构图像的阈值分割方式为：先选取固定尺寸的局部滑动窗口，通过求取窗口加权和的方式得到分割的阈值；对静脉结构图像分割完成后得到二值图像。

对获得的二值图像，做一次形态学开运算，使得每段静脉之间的连接点断开，并且保留手臂两侧轮廓。对二值图像中的所有轮廓进行查找，获得所有轮廓后，接着求解每个轮廓的最小外包矩形。

将所述静脉特征输入到训练好的分类器中完成肘正中静脉的识别。

所述对采集到的静脉图像进行增强的方法如下：

对采集到的静脉图像进行分块，对每个小块计算直方图，设置裁剪阈值；

对每一灰度级上超过裁剪阈值的直方图进行裁剪，然后均匀的分布到每一灰度级；

对裁剪完成的直方图做直方图均衡化；

对裁剪完成的直方图做限制对比度的自适应直方图的均衡化处理；增强静脉区域与图像背景之间的对比度，并采用直接赋值，线性插值，双线性插值结合的方式加快运行速度。在手肘静脉图像中，静脉直径存在变化，不适合采用在单一尺度下的增强方法。多尺度空间理论可用于表示同一空间不同尺度的信号特征。在该理论基础上，将Hessian矩阵的差分运算与高斯函数结合；求解Hessian矩阵的二阶微分运算转换为原始图像数据与高斯函数二阶导数的卷积。卷积完成后能得到原图像中每个像素点所对应的二阶Hessian矩阵，求解Hessian矩阵的特征值和特征向量，并将特征值按其绝对值大小进行排序，获得特征值与几何形状之间的一一对应关系，当该点所对应为原图像中的管状结构即静脉区域时，较小特征值的大小接近于0，较大特征值的绝对值远大于较小特征值的绝对值。

对应为块状区域时即阴影区域时，较小特征值与较大特征值的绝对值都为较大的值；对应为噪声区域时，较小特征值与较大特征值的绝对值都为较小值。

利用对应关系，构建多尺度自适应滤波器用于提取静脉结构，抑制非静脉区域，获得静脉结构图像。对获得的二值图像，做一次形态学开运算，使得每段静脉之间的连接点断开，并且保留手臂两侧轮廓。对二值图像中的所有轮廓进行查找，获得轮廓之后，通过求解每个轮廓的最小外包矩形。

将均衡化后的分块图像拼接成完整图像，对左上，左下，右下，右上四个1/4个小块采用直接赋值的方式，对上，左，下，右四条边采用线性插值，图像中剩余部分采用双线性插值方式；所述四条边不包括上述四个小块。

所述进行多尺度自适应滤波获得静脉结构图像的方法如下：

根据图像二阶hessian矩阵特征值λ₁，λ₂，|λ₁|≤|λ₂|构建的多尺度自适应滤波函数，如公式(6)所示：

其中

β＝0.5，ω＝50，d由输入图像自适应确定，其形式为公式(7)所示，d_rank为输入图像中像素个数不为0的灰度级，σ_i为当前迭代的尺度因子，滤波器输出0≤v(x，y；σ)≤1；

所述提取静脉特征的方法如下：

对所述二值图像中的所有轮廓进行查找，获得轮廓之后，通过求解每个轮廓的最小外包矩形；

然后通过所述最小外包矩形确定每段静脉包括手臂轮廓的宽高及面积，用于表征该段静脉的大小；

然后利用所述最小外包矩形的宽高比来表征这段静脉的弯曲程度；

接着求取每段静脉最小外包矩形的中心距离两侧手臂轮廓的距离；具体方式为通过最小外包矩形的中心向两侧轮廓的最小外包矩形的最外侧边做垂线，分别为左侧距离和右侧距离，用于表征该段静脉的相对位置；

然后将所述最小外包矩形映射回到还未进行阈值分割时的静脉结构图像，并计算矩形包围的像素点灰度值之和，用于表征该段静脉易于被观察的程度；

至此每段静脉获得7维特征，将部分特征进行重组，再添加左侧距离与右侧距离之比、右侧距离与左侧距离之比和灰度强度与面积之比，这样使得每段静脉一共获得10维特征。

所述分类器包括：随机森林分类器和支持向量机分类器；

对每段静脉添加相应的标签，标注其是否为用于注射的肘正中静脉，将10维特征和标签用于训练分类器；

把采集到的若干幅手肘静脉图像分别独立训练每个分类器后，使每个分类器输出相应类别的概率值，通过软投票的方式将其集成为最终的预测器，完成对肘正中静脉的识别。

一种手肘图像采集装置，包括镜头朝下布置的摄像头，在位于摄像头的正下方设有滤光片，在位于滤光片的正下方设有能放置手肘的肘托，在摄像头的左右两侧分别设有朝下照射的红外光源，在每个红外光源的正下方都分别设有散光片。两个散光片所处位置都高于镜头所处位置，并且两个散光片位于同一高度处。位于摄像头左右两侧方的红外光源位于同一高度处。在位于摄像头的镜头正下方的肘托上设有放肘槽。

采用200万像素的CCD摄像头；摄像头前装的滤光片其中心波长为850±10nm，峰值透射率＞86％；摄像头两侧的红外光源为850nm的LED红外光源；在红外光源前装有散光片。摄像头正下方为肘托，采集图像时将手肘置于肘托的放肘槽内，采集到的图像大小为640*480。

本实施例1能通过对手肘图像的处理来准确找出肘正中静脉，可靠性高。

实施例2，本发明为一种肘静脉图像肘正中静脉的识别方法及手肘图像采集装置，目的是能够通过图像识别方法快速准确的找到用于手肘处静脉注射的肘正中静脉。附图1所示为本发明的总体流程图，其包括手肘静脉图像采集模块，图像预处理模块，肘静脉特征提取模块，肘正中静脉预测模块。

附图2为图像采集装置的原理图，包括肘托，200万像素CCD摄像头，850nm红外LED灯，中心波长为850±10hm，峰值透射率＞86％的红外滤光片，散光片。采集图像时，将手肘置于摄像头正下方的肘托上，采集到的图像大小为640*480。

图像预处理模块包括静脉图像增强，静脉结果提取，阈值分割，静脉分段。

在具体实施时，将采集到的图像进行灰度化，使其成为单通道灰度图像。再对图像进行限制对比度自适应直方图均衡化，具体的方式为将图像分割为10*10的小块，对每个小块求解灰度直方图，并且设置裁剪值，当某一灰度级的像素个数超过裁剪值，将其均匀的分布到每一灰度级，将灰度直方图重构，接重构后的直方图对每一个小块做直方图均衡。上述过程可由公式(1)，(2)，(3)完成。

p_r(r_k)＝n_k，k＝0，1，2，...，L-1 公式(1)

其中n_k是灰度为r_k的像素个数，L是灰度级的数量，在本实施例中为256。α_clip为裁剪系数，M，N为每个小块的高和宽。T(r_k)为直方图均衡化的转换函数。完成上述过程后，需要将10*10个小块图像重新合成。具体合成方式为，对左上，左下，右下，右上四个1/4个小块采用直接赋值的方式，对上，左，下，右四条边(不包括上述四个小块)采用线性插值，图像中剩余部分采用双线性插值方式。

完成对比度限制自适应直方图均衡后的图像，进行多尺度自适应滤波操作。首先求解图像的二阶hessian矩阵，求解方式可用公式(4)表示。

其中I为原图像，I_xy为原图像与二元高斯函数的二阶偏导数卷积后的结果即原图像的二阶偏导数，同理可求解I_xx，I_yy，则原图像上(x，y)点所对应的二阶hessian矩阵为公式(5)所示。

滤波器的局部特性分析窗口的半宽设置为3σ，通过对σ取不同的值，设置不同的窗口半宽和不同的二元高斯偏导数达到多尺度的目的，筛选宽度不一的静脉，在本实施例中σ＝1，3，5，7，9，一共迭代计算五次。对得到的hessian矩阵求解特征值和特征向量，并将特征值按其绝对值大小进行排序|λ₁|≤|λ₂|。可以获得与之对应的几何形状之间的关系，当该点所对应为原图像中的管状结构即静脉区域时，|λ₁|接近于0，|λ₂|远大于|λ₁|。对应为块状区域时即阴影区域，|λ₁|与|λ₂|均为较大的值。对应为噪声区域时，|λ₁|与|λ₂|都为较小值。应用此关系并根据近红外手肘静脉图像的特点可以构建多尺度自适应滤波器，其具体形式为公式(6)所示

其中

在本实施例中β＝0.5，ω＝50，d由输入图像自适应确定，其形式为公式(7)所示，d_rank为输入图像中像素个数不为0的灰度级，σ_i为当前迭代的尺度因子。滤波器输出0≤v(x，y；σ)≤1。

接着，经过5轮迭代之后，每个像素点对应最大的滤波器输出作为最终的输出如公式(8)所示：

V(x，y)＝max(v(x，y；σ_i))，i＝1，2，3，4，5 公式(8)

进一步地，对完成多尺度自适应滤波后的图像即静脉结构图像进行自适应阈值分割，在本实施例中自适应方法为领域加权和，领域大小为25*25，权重为一个高斯窗口，获得二值图像。

进一步地，对获得的静脉二值图像，做一次形态学开运算，使得每段静脉之间的连接点断开，并且保留手臂两侧轮廓，在本实施例中形态学运算所选取的结构元素为椭圆形，长短轴为5和3。

所述特征提取模块包括对各段静脉几何特征，位置特征及其组合特征的提取，具体实施方式为：

对所述二值图像中的所有轮廓进行查找，获得轮廓之后，通过求解每个轮廓的最小外包矩形。

进一步地，通过所述最小外包矩形确定每段静脉包括手臂轮廓的宽高及面积，用于表征该段静脉的大小。

进一步地，利用所述最小外包矩形的宽高比来表征这段静脉的弯曲程度。

进一步地，求取每段静脉最小外包矩形的中心距离两侧手臂轮廓的距离，具体方式为通过最小外包矩形的中心向两侧轮廓的最小外包矩形的最外侧边做垂线，分别为左侧距离和右侧距离，用于表征该段静脉的相对位置。

进一步地，将所述最小外包矩形映射回所述静脉结构图像即未进行阈值分割的图像，计算矩形包围的像素点灰度值之和，用于表征该段静脉易于被观察的程度。

至此每段静脉获得7维的特征，将部分特征进行重组，添加左侧距离与右侧距离之比，右侧距离与左侧距离之比，灰度强度与面积之比，一共10维的特征。

肘正中静脉识别模块包括模型的选择，训练，集成方式。其具体实施方式为：

分类器1为随机森林，在本实施列2中，随机森林分类器集成了500个决策树。

分类器2为支持向量机，在本实施例中，支持向量机的核函数为高斯RBF核函数。

对每段静脉添加相应的标签，标注其是否为用于注射的肘正中静脉。将特征和标签用于训练分类器。本实施例中一共选取150幅手肘静脉图像，其中包括75幅左手，75幅右手，分别独立训练每个分类器后，使每个分类器输出相应类别的概率值，通过软投票的方式将其集成为最终的预测器，完成对肘正中静脉的识别，其中随机森林的识别准确率为86.4％，支持向量机的准确率为85.2％，将其集成后的分类器准确率为95.1％，能够准确识别肘正中静脉。

本实施例2能通过对手肘图像的处理来准确找出肘正中静脉，可靠性高。

Claims (3)

1.一种手肘静脉图像肘正中静脉的识别方法，其特征在于，识别方法如下：

对采集到的手肘图像进行增强，获得高对比度的手肘静脉图像；

对所述高对比度的手肘静脉图像，进行多尺度自适应滤波获得静脉结构图像；

对所述静脉结构图像进行阈值分割，在其二值图像中提取静脉特征；

将所述静脉特征输入到训练好的分类器中完成肘正中静脉的识别；

所述对采集到的静脉图像进行增强的方法如下：

对采集到的静脉图像进行分块，对每个小块计算直方图，设置裁剪阈值；

对每一灰度级上超过裁剪阈值的直方图进行裁剪，然后均匀的分布到每一灰度级；

对裁剪完成的直方图做直方图均衡化；

将均衡化后的分块图像拼接成完整图像，对左上，左下，右下，右上四个1/4个小块采用直接赋值的方式，对上，左，下，右四条边采用线性插值，图像中剩余部分采用双线性插值方式；四条边不包括上述四个小块；

所述提取静脉特征的方法如下：

对所述二值图像中的所有轮廓进行查找，获得轮廓之后，通过求解每个轮廓的最小外包矩形；

然后通过所述最小外包矩形确定每段静脉包括手臂轮廓的宽高及面积，用于表征该段静脉的大小；

然后利用所述最小外包矩形的宽高比来表征这段静脉的弯曲程度；

接着求取每段静脉最小外包矩形的中心距离两侧手臂轮廓的距离；具体方式为通过最小外包矩形的中心向两侧轮廓的最小外包矩形的最外侧边做垂线，分别为左侧距离和右侧距离，用于表征该段静脉的相对位置；

然后将所述最小外包矩形映射回到还未进行阈值分割时的静脉结构图像，

并计算矩形包围的像素点灰度值之和，用于表征该段静脉易于被观察的程度；

至此每段静脉获得7维特征，将部分特征进行重组，再添加左侧距离与右侧距离之比、右侧距离与左侧距离之比和灰度强度与面积之比，这样使得每段静脉一共获得1Q维特征。

2.根据权利要求1所述的手肘静脉图像肘正中静脉的识别方法，其特征在于，

所述进行多尺度自适应滤波获得静脉结构图像的方法如下：

根据图像二阶hessian矩阵特征值λ₁，λ₂，|λ₁|≤|λ₂|构建的多尺度自适应滤波函数，如公式(6)所示：

其中

β＝0.5，ω＝50，d由输入图像自适应确定，其形式为公式(7)所示，d_rank为输入图像中像素个数不为0的灰度级，σ_i为当前迭代的尺度因子，滤波器输出0≤v(x，y；σ)≤1；

3.根据权利要求2所述的手肘静脉图像肘正中静脉的识别方法，其特征在于，所述分类器包括：随机森林分类器和支持向量机分类器；

对每段静脉添加相应的标签，标注其是否为用于注射的肘正中静脉，将10维特征和标签用于训练分类器；

把采集到的若干幅手肘静脉图像分别独立训练每个分类器后，使每个分类器输出相应类别的概率值，通过软投票的方式将其集成为最终的预测器，完成对肘正中静脉的识别。

Images