一种结合先验知识的鞋印花纹渐进细化式提取方法与流程

本发明涉及图像处理技术领域，特别是一种结合先验知识的鞋印花纹渐进细化式提取方法。

背景技术：

目前鞋印花纹提取的方式主要有以下两种：

(1)人工手动提取，主要是由经验丰富的专业人员利用图像处理软件，分割出鞋印图像的花纹区域，提取鞋印花纹。

(2)自动提取方法：刘健康运用模糊c-均值聚类算法(fcm)对鞋印图像进行分割，提出了基于二维属性直方图的最大类间方差方法提取鞋印花纹。杨芳将现场鞋印图像按背景进行分类，对不同的背景采用不同的滤波方式进行去噪。利用多尺度top-hat变换的方法来增强鞋印图像的对比度，对于对比度明显的图像运用基于子图的阈值法对图像进行分割，解决光照不均问题。对于对比度不高的图像运用快速模糊c-均值聚类算法，即一种基于分水岭变换的模糊c-均值聚类算法实现花纹的提取。专利号为cn103549960的专利对采集的现场鞋印图像，采用snake模型确定轮廓，然后结合显影图轮廓和显迹图轮廓确定鞋印花纹特征的范围，确定感兴趣区域，用混合高斯模型减影方法提取鞋印花纹。

上述方法均可在一定程度上实现对鞋印花纹的现场提取，但相应的存在的以下问题：(1)人工手动提取的方法工作量大，效率低，并且存在语义的二义性，存在主观性；(2)自动提取方法没有利用鞋印花纹的先验知识以及用户的语义信息，提取的花纹存在花纹丢失或包含较多背景干扰等情况。这些问题将严重影响公安人员侦破案件的效率，还将制约足迹识别技术的应用发展。

技术实现要素：

鉴于已有技术存在的不足，本发明的目的是要提供一种利用花纹特征先验知识并兼顾提取效率的鞋印花纹提取方法。

为了实现上述目的，本发明技术方案如下：

一种结合先验知识的鞋印花纹渐进细化式提取方法，其特征在于，步骤包括：

s1采集原始现场鞋印图像信息，并对鞋印图像进行超像素分割，得到特征一致的图像子区域；

s2对鞋印图像进行粗标记，得到花纹区域f＝{f1,f2,......,fn}以及背景区域b＝{b1,b2,......,bm}，并统计标记花纹区域和背景区域的个数，其中n为花纹区域的个数，m为背景区域的个数；

s3根据用户标记结果对鞋印图像的子区域，利用兼顾方向特征与颜色特征的区域相似性度量方法进行最大相似性区域合并，得到初始的鞋印花纹提取结果p1；

s4对初始的鞋印花纹提取结果p1进行迭代细化修正，得到准确的鞋印花纹提取结果p；

s5对准确的鞋印花纹提取结果p进行光照校正与花纹增强处理；

s6利用otsu算法将处理后的图像花纹提取结果进行二值化，得到最终的鞋印花纹图像二值图。

与现有技术相比，本发明的有益效果：

1)本发明充分考虑了鞋底花纹的特点，利用鞋印花纹的颜色特征和纹线的方向特征的先验知识分别计算区域间的相似性，根据二者的融合得分，并结合用户的交互信息，进行花纹与背景的判定，提高了花纹提取的准确性。

2)将手工标记与自动提取相结合，兼顾了人工标记的语义信息和自动提取的便捷性，使提取的鞋印花纹的更符合主观效果并且降低了人工手动提取的难度。

3)提出了基于knn的区域分类策略，采用渐进细化的交互方式提取鞋印花纹，提高了花纹与背景判定的准确性，减少背景干扰的同时，保留了鞋印花纹细节信息，在提高便捷性的同时，提高了花纹的提取精度。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。

图1为本发明现场鞋印花纹提取流程图；

图2为本发明初始鞋印花纹提取流程图；

图3为本发明初始鞋印花纹提取实施图例；

图4为本发明鞋印花纹细化修正流程图；

图5为本发明现场鞋印花纹提取的整体实施图。

具体实施方式

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

本发明设计了一种结合先验知识的鞋印花纹渐进细化式提取方法，下面结合附图以及具体实施例进一步说明本发明的技术方案：

如图1所示为结合先验知识的鞋印花纹渐进细化式提取过程，其步骤包括：

s1采集原始现场鞋印图像信息，并对鞋印图像进行超像素分割，得到特征一致的图像子区域。

s2对鞋印花纹图像进行粗标记，得到花纹区域f＝{f1,f2,......,fn}以及背景区域b＝{b1,b2,......,bm}，并统计标记花纹区域和背景区域的个数，其中n为花纹区域的个数，m为背景区域的个数。

s3对鞋印图像的子区域，根据用户标记，利用兼顾方向特征与颜色特征的区域相似性度量方法进行最大相似性区域合并，得到初始的鞋印花纹提取结果p1。初始鞋印花纹的提取过程见如图2所示，其步骤包括：

s3-1计算步骤s2标记的各图像区域中任意相邻区域之间的颜色相似度s1(r,q)，其中r和q分别表示任意相邻的两区域，s1(r,q)的求取过程如下：

①首先分别在r、g、b三个颜色通道上采用16级量化对图像进行处理，统计量化处理后各区域在16×16×16＝4096维的颜色特征空间内的特征值，作为本发明的较佳实施方式，优选采用颜色特征空间的直方图描述的特征矢量特征，其维度为4096维。

②对于任意相邻区域r和q，利用bhattacharyya系数衡量两区域之间的颜色相似度，即

其中，表示区域r中颜色u出现的概率，表示区域q中颜色u出现的概率。作为本发明的较佳实施方式利用颜色归一化直方图中颜色u的数值表示其出现的概率，和分别表示区域r和区域q的颜色归一化直方图中u的数值。

s3-2计算步骤s2标记的各图像区域中任意相邻区域之间的方向场相似度s2(r,q)，s2(r,q)的求取过程如下：

①对原始鞋印图像利用基于梯度的方法构建方向场，其主要步骤如下：

a.将鞋印图像分成w×w的子图，利用sobel算子计算子图中每一像素点(i,j)沿x轴方向的梯度值和沿y轴方向的梯度值

b.利用下述公式计算子任意子图中心点(i,j)的局部方向值

其中，vx(i,j)为w×w区域内在x方向的平均梯度矢量，vy(i,j)为w×w区域内在y方向的平均梯度矢量θ(i,j)是以像素点(i,j)为中心的子块的局部方向，θ(i,j)∈[0,π]，得到鞋印图像的方向场；

②对于任意相邻区域r和q，利用bhattacharyya系数衡量相邻区域r和q之间的方向场相似度，即

其中，表示区域r中方向v出现的概率，表示区域q中方向v出现的概率。作为本发明的较佳实施方式利用方向场归一化直方图中方向v的数值表示其出现的概率，和分别表示区域r和区域q的方向场归一化直方图中v的数值。

s3-3根据颜色相似度s1(r,q)和方向场相似度s2(r,q)计算相邻区域r和区域q的相似度s(r,q)，

s(r,q)＝λs1(r,q)+(1-λ)s2(r,q)

其中，λ为颜色相似度所占权重，作为本发明的优选，取λ＝0.7。

s3-4统计鞋印花纹图像所有相邻区域的相似度集合s。

s3-5根据所有相邻区域的相似度s，以及步骤s2得到的用户标记区域，进行最大相似性区域合并，得到初始的鞋印花纹提取结果p1。图3给出了初始鞋印花纹提取的一个实施流程。

s4对初始的鞋印花纹提取结果p1进行细化修正，得到准确的鞋印花纹提取结果p，其步骤包括：

s4-1对初始的鞋印花纹提取结果p1不满意的部分进行再次标记，不满意的鞋印花纹图像区域为r，r＝r^f∪r^b，其中r^f＝{r^f1,r^f2,......,r^fk}为本应属于花纹错被判定为背景的区域，r^b＝{r^b1,r^b2,......,r^bl}为本应属于背景错被判定为花纹的区域，k表示用户标记的本应属于花纹错被判定为背景的区域个数，l表示本应属于背景错被判定为花纹的区域个数；

s4-2计算区域ri与区域q的相似度s(ri,q)＝{s1,s2,......,sn+m}，其中表示与任意一个不满意的鞋印花纹图像标记区域rm的相邻区域集合，q＝{f1,f2,......,fn,b1,b2,......bm}表示在步骤s2中标记的所有区域；

s4-3根据计算得到的s(ri,q)，利用knn分类算法判断区域ri的区域种类，即判断区域ri属于花纹还是背景区域，具体地是对由步骤s4-2得到的相似度集合s(ri,q)中的相似度进行降序排序，统计前k个相似度对应区域q中的区域标记种类，区域ri自动标记为图像种类数目较多的区域种类，其中k为奇数，作为本发明的较佳实施方式，优选k＝7。

s4-4如果区域ri与区域rm的种类相同，则将区域ri合并至区域rm，即rm＝rm∪ri，否则不做处理；

s4-5重复步骤s4-2～s4-4，直至遍历所有标记区域rm的邻接区域ri；

s4-6重复步骤s4-2～s4-5，直至遍历步骤s4-1中标记的不满意的鞋印花纹图像区域集合r的所有区域rm，得到区域合并结果t＝{t^f,t^b}，其中t^f表示由标记区域r^f得到的合并结果，t^b表示由标记区域r^b得到的合并结果。

s4-7将步骤s4-6得到的区域合并结果t与初始的鞋印花纹提取结果p1按照下述规则进行合并：

p1＝p1∪t^f

以得到细化修正的鞋印花纹提取结果p1。

s4-8观察步骤s4-7得到的细化修正的鞋印花纹提取结果p1，如果对结果仍然不满意，重复步骤s4-1～s4-7，否则准确的鞋印花纹提取结果p＝p1。

图4给出了准确鞋印花纹提取过程的实施流程。

s5对准确的鞋印花纹提取结果p进行光照校正与花纹增强处理，其具体步骤包括：

s5-1对准确的鞋印花纹提取结果p进行灰度化处理，得到灰度图像i；

s5-2灰度图像i进行光照校正，光照校正的步骤包括：

①将灰度图像i分割为大小相同的若干图像区域，计算任意区域内的像素灰度均值μ和标准差σ，并以max(min,μ-0.8σ)作为该区域的近似背景灰度，其中min为该区域像素的最小灰度值；

②利用所有区域的近似背景灰度构成近似背景灰度矩阵；

③利用插值法将上述近似背景灰度矩阵扩展为大小与灰度图像i相同的背景灰度矩阵；

④从灰度图像i中减去步骤③得到的背景灰度矩阵，以得到光照均匀图像i1；

s5-3利用对数变换s＝clog(1+r)对图像i1进行图像增强以增强图像中较暗部分的细节，其中c为常数，s为增强后的图像像素值，r为图像i1像素值。

s6利用otsu算法将处理后的图像花纹提取结果进行二值化，得到最终的鞋印花纹图像二值图。

图5给出了本发明鞋印花纹现场提取的实施流程。

本发明充分考虑了鞋底花纹的特点，利用鞋印花纹的颜色特征和纹线的方向特征的先验知识分别计算区域间的相似性，并提出了基于knn的区域分类策略，采用渐进细化的交互方式提取鞋印花纹，提高了花纹与背景判定的准确性，减少背景干扰的同时，保留了鞋印花纹细节信息，在提高便捷性的同时，提高了花纹的提取精度。同时将手工标记与自动提取相结合，兼顾了人工标记的语义信息和自动提取的便捷性，使提取的鞋印花纹的更符合主观效果并且降低了人工手动提取的难度。

最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。