一种基于改进的主动轮廓算法的指静脉图像分割方法

摘要

本发明属于图像处理领域，具体涉及一种基于改进的主动轮廓算法的指静脉图像分割方法，该方法包括：采用局部二值拟合主动轮廓算法构建指静脉图像分割模型；采用改进的模糊均值聚类算法计算待分割的指静脉图像的初始轮廓，将初始轮廓作为指静脉图像分割模型的初始指静脉轮廓；根据初始指静脉轮廓采用指静脉图像分割模型对待分割的指静脉图像进行分割处理；本发明对每张静脉图像均会进行迭代优化运算，输出最优的手指静脉初始分割结果，然后在该基础上进一步进行手指静脉图像分割，保证分割模型的稳定性及准确率。

说明书

技术领域

本发明属于图像处理领域，具体涉及一种基于改进的主动轮廓算法的指静脉图像分割方法。

背景技术

随着人们对生物识别系统安全性和准确性的要求越来越高，生物特征识别技术得到越来越多的关注。指静脉识别作为众多生物特征识别技术中的一种，由于具有非接触式采集、活体检测、不易伪造、成本较低等优点，成为当前研究的热点。指静脉图像中血管的分割是指静脉识别技术中的一个关键步骤，而拍摄采集到的指静脉图像不仅含有静脉纹路，还含有不规则的噪声、指骨骼和肌肉的不同厚度所产生的阴影。此外，指静脉会随着温度或物理条件的变化而变化，很难精确地提取指静脉血管的细节，因此如何对指静脉图进行准确的分割，直接影响后续识别的精度和准确度。

手指静脉图像在采集过程中，由于手指中的其他组织吸收引起光衰减，捕获的手指静脉图像不仅具有低对比度，还包含由不同厚度的指骨及肌肉产生不规则噪音和阴影；此外，手指静脉特征不仅会随温度、物理条件、外界光照等的变化而暂时变化，还会受到采集过程中平移、旋转、手指压力和采集姿势等的影响。这些影响会导致使用传统算法在特征提取过程中容易出现静脉脉络不连续问题，造成了采集的指静脉图像的准确度低，特征不明显的问题。

发明内容

为解决以上现有技术存在的问题，本发明提出了一种基于改进的主动轮廓算法的指静脉图像分割方法，该方法包括：获取待分割的指静脉图像；采用改进的主动轮廓算法对待分割的指静脉图像进行分割处理，得到分割后的指静脉图像；采用改进的主动轮廓算法对待分割的指静脉图像进行处理的过程包括：

S1：采用局部二值拟合主动轮廓算法构建指静脉图像分割模型；

S2：采用改进的模糊均值聚类算法计算待分割的指静脉图像的初始轮廓，将初始轮廓作为指静脉图像分割模型的初始指静脉轮廓；

S3：根据初始指静脉轮廓采用指静脉图像分割模型对待分割的指静脉图像进行分割处理。

优选的，构建指静脉图像分割模型的过程包括：获取Heaviside函数、水平集函数以及高斯核函数；根据Heaviside函数、水平集函数和高斯核函数计算以x为圆心的圆与零水平集轮廓相交的内部均值f

进一步的，能量泛函模型的表达式为：

ε

其中，ε

优选的，计算待分割的指静脉图像的初始轮廓的过程包括：采用基于果蝇优化算法的模糊C均值聚类算法获取初始轮廓，其步骤包括：

步骤1：初始化参数；初始化参数包括设置初始聚类中心个数n，最大聚类中心个数；

步骤2：采用果蝇优化算法获取全局最优解，最优解为当前聚类中心n条件下的最优聚类中心；

步骤3：记录当前聚类中心数n对应最优值时的味道浓度值Sm；

步骤4：当聚类中心数n大于设置的最大聚类中心个数时，结束果蝇寻优算法；当聚类中心数n小于等于设置的最大聚类中心个数时，设置n＝n+1，返回步骤2继续进行寻优算法；

步骤5：比较所有n的取值，获取最优解时的味道浓度值Sm，即当味道浓度值最大时聚类中心数为最佳的聚类中心数，输出聚类中心对应的位置值；将该该位置值作为初始聚类中心位置；

步骤6：根据初始聚类中心位置采用模糊C均值聚类算法对待分割图像进行分割，输出手指静脉图像的初始分割结果，并将其设置为分割模型的初始轮廓。

进一步的，采用果蝇优化算法获取全局最优解的过程包括：

步骤21：获取果蝇种群位置范围Pr和果蝇单次飞行范围Sr；

步骤22：根据果蝇种群位置范围Pr和果蝇单次飞行范围Sr计算果蝇群体的初始位置x和y；

步骤23：群体的果蝇根据食物气味进行随机搜索，得到新的位置坐标(x

步骤24：计算个体位置与原点的距离D

步骤25：根据味道浓度判定值s

步骤26：在群体中选择味道浓度值最优的果蝇个体，记录浓度值bestSm和位置bestId，并求极大值max(Sm)，该极大值为味道浓度值最优的位置；

步骤27：不断重复步骤23～步骤26，判断迭代后的的气味浓度是否比迭代更新前的浓度值更优，若是，则执行步骤28，否则继续进行迭代操作；

步骤28：其他果蝇个体根据获取的最优浓度值位置进行位置更新。

优选的，采用模糊C均值聚类算法对待分割图像进行分割的过程包括：

步骤61：根据初始聚类中心位置获取模糊C均值聚类算法目标函数；

步骤62：确定隶属度矩阵U＝{u

其中，n表示聚类中心数，u

步骤63：采用拉格朗日乘数法获取目标函数的极小值函数；

步骤64：对极小值函数的隶属度u

步骤65：根据约束条件对目标函数的极小值函数进行化简，将化简后的表达式带入到隶属度偏导函数中，得到隶属度公式；

步骤66：对目标函数的极小值函数的c

步骤67：重复计算u

进一步的，目标函数表达式为：

其中，

优选的，隶属度公式为：

其中，n表示聚类中心数，l表示样本对聚类中心的加权值，x

优选的，聚类中心表达式为：

其中，

优选的，采用指静脉图像分割模型对待分割的指静脉图像进行分割处理的过程包括：

步骤1：初始化参数α

步骤2：获取初始指静脉轮廓，根据初始指静脉轮廓初始化指静脉图像分割模型中的水平集；

步骤3：采用梯度下降算法求解能量泛函模型对于水平集函数φ的极小化；即计算与零水平集轮廓相交的内部均值f

步骤4：采用变分法表示能量泛函模型的梯度下降流，对梯度下降流公式进行正则化，得到正则化后的梯度流方程；

步骤5：根据正则化后的梯度流方程对水平集函数进行更新；

步骤6：判断正则化后的梯度流方程是否收敛，若不收敛，则返回步骤3继续迭代，若收敛，则输出更新后的水平集函数。

本发明的有益效果：

本发明不仅可以解决手指静脉图像分割过程出现静脉脉络不连续的问题，而且整个模型不需要大量样本训练来提升静脉图像的分割精度，针对每张静脉图像，该模型都会进行迭代优化运算输出最优的手指静脉初始分割结果，然后在该基础上进一步进行手指静脉图像分割，保证分割模型的稳定性及准确率。

附图说明

图1为本发明的主动轮廓手指静脉图像分割算法流程图；

图2为本发明的主动轮廓手指静脉图像分割算法中初始轮廓获取流程图；

图3为本发明的主动轮廓手指静脉图像分割算法中果蝇优化算法流程图；

图4为本发明的主动轮廓手指静脉图像分割算法中模糊C均值聚类算法流程图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

一种基于改进的主动轮廓算法的指静脉图像分割方法，如图1所示，该方法包括：获取待分割的指静脉图像；采用改进的主动轮廓算法对待分割的指静脉图像进行分割处理，得到分割后的指静脉图像；采用改进的主动轮廓算法对待分割的指静脉图像进行处理的过程包括：

S1：采用局部二值拟合主动轮廓算法构建指静脉图像分割模型；

S2：采用改进的模糊均值聚类算法计算待分割的指静脉图像的初始轮廓，将初始轮廓作为指静脉图像分割模型的初始指静脉轮廓；

S3：根据初始指静脉轮廓采用指静脉图像分割模型对待分割的指静脉图像进行分割处理。

构建指静脉图像分割模型的过程包括：获取Heaviside函数、高斯核函数以及水平集函数；根据Heaviside函数、高斯核函数以及水平集函数计算图像的内外局部均值；根据内外局部均值构建能量泛函模型。当水平集函数φ处于零水平集时，可将图像域分成两个互不相交的内外轮廓区域，内部轮廓表示为H

H

其中，φ表示水平集函数。

K

其中，σ控制核函数的尺度大小，n表示聚类中心数。

根据Heaviside函数、高斯核函数以及水平集函数计算图像的内外局部均值的公式为：

以x为圆心的圆与零水平集轮廓相交的内部均值的表达式为：

以x为圆心的圆与零水平集轮廓相交的外部均值的表达式为：

其中，K

能量泛函模型由核函数定义，其公式为：

ε

+α

其中，ε

如图2所示，计算待分割的指静脉图像的初始轮廓的过程包括：采用基于果蝇优化算法的模糊C均值聚类算法获取初始轮廓，其步骤包括：

步骤1：初始化各项参数，设置初始聚类中心个数n，最大聚类中心个数；

步骤2：依照果蝇优化算法步骤，利用果蝇优化算法步骤寻找全局最优解，找到的最优结果就是当前聚类中心n条件下的最优聚类中心；

步骤3：记录当前聚类中心数n对应最优值时的味道浓度值Sm；

步骤4：当聚类中心数n大于设置的最大聚类中心个数时，结束果蝇寻优算法；当聚类中心数n小于等于设置的最大聚类中心个数时，设置n＝n+1，返回步骤2继续进行寻优算法；

步骤5：通过比较所有n的取值情况，找到最优解时的味道浓度值Sm，值最大时的聚类中心数就判定为最佳的聚类中心数，而此时输出的聚类中心对应的位置值，就作为接下来执行模糊C均值聚类算法的初始聚类中心位置；

步骤6：根据果蝇优化算法输聚类中心；

步骤7：执行模糊C均值聚类算法，输出手指静脉图像的初始分割结果，并将其设置为分割模型的初始轮廓。

如图3所示，果蝇优化算法的实现过程包括：用Pr表示果蝇种群位置范围、用Sr表示果蝇单次飞行范围，用(x,y)表示各个果蝇个体的位置信息，设置初始条件，其步骤为：

步骤1：随机获取果蝇群体的初始位置x和y；

x＝rand(Pr)

y＝rand(Pr)

其中，rand()表示随机生成函数。

步骤2：群体里的每个果蝇会按食物味道随机搜寻，得到新的位置坐标(x

x

y

步骤3：计算个体位置与原点的距离D

步骤4：将s

Sm

其中，Fiction表示判定函数。

步骤5：在群体中选择味道浓度值最优的果蝇个体，记录浓度值bestSm和位置bestId，并求极大值max(Sm)；

[bestSm,bestId]＝max(Sm)

步骤6：不断重复步骤2～步骤5，判断迭代后的气味浓度是否比迭代更新前的浓度值更优，若是，则执行步骤7，否则继续进行迭代操作；

步骤7：获得最优浓度值位置后，其它果蝇个体则利用视觉根据位置信息行动，公式如下：

Smbest＝bestSm

x＝x(bestId)

y＝y(bestId)

如图4所示，采用模糊C均值聚类算法对待分割图像进行分割的过程包括：将一幅手指静脉图像看作所有像素点的集合，若图中所有像素点数为n，将总的像素点划分为c类，那么总的聚类中心就有c个，将其表示为c

步骤1：根据初始聚类中心位置获取模糊C均值聚类算法目标函数；目标函数的表达式为：

其中，

步骤2：确定隶属度矩阵U＝{u

步骤3：采用拉格朗日乘数法获取目标函数的极小值函数；其表达式为：

其中，λ拉格朗日乘子，F表示目标函数极小值。

步骤4：对极小值函数的隶属度u

步骤5：将约束条件

步骤6：将步骤5所得式代入式步骤4得隶属度公式，如下式所示：

步骤7：对F关于c

步骤8：反复计算u

采用指静脉图像分割模型对待分割的指静脉图像进行分割处理的过程包括：

步骤1：初始化参数α

步骤2：获取初始指静脉轮廓，根据初始指静脉轮廓初始化指静脉图像分割模型中的水平集；

步骤3：采用梯度下降算法求解能量泛函模型对于水平集函数φ的极小化；即计算与零水平集轮廓相交的内部均值f

步骤4：采用变分法表示能量泛函模型的梯度下降流，对梯度下降流公式进行正则化，得到正则化后的梯度流方程；

步骤5：根据正则化后的梯度流方程对水平集函数进行更新；

步骤6：判断正则化后的梯度流方程是否收敛，若不收敛，则返回步骤3继续迭代，若收敛，则输出更新后的水平集函数。

具体的，采用梯度下降算法，得出能量泛函模型ε

其中，

利用变分法进一步表示能量泛函的梯度下降流，其表达式为：

其中，e

e

e

其中，K

δ

其中，τ表示有效宽度。

为加强水平集函数φ演化的稳定性，在利用变分法表示的能量泛函梯度下降流方程中引入正则化，则引入正则化后梯度流方程如下式所示，并依据下述公式持续演化，直至收敛，则目标边缘由对应的曲线或轮廓定位，实现手指静脉图像分割。

其中，α

以上所举实施例，对本发明的目的、技术方案和优点进行了进一步的详细说明，所应理解的是，以上所举实施例仅为本发明的优选实施方式而已，并不用以限制本发明，凡在本发明的精神和原则之内对本发明所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。