一种基于隐马尔科夫模型的步态去噪方法与流程

本发明属于计算机视觉和模式识别领域，具体地说涉及一种基于隐马尔科夫(hmm)模型的步态去噪方法。

背景技术：

在公共场合中，人流量比较大，大量的人流同时就会产生难以估计的行为，如果能快速准确的识别公共场合中的行人，并理解其产生的行为，对维护公共秩序和安全将会有极大的帮助，因此基于生物特征的识别研究成为现在公共场合安防的主流技术。

作为生物特征识别技术的重要组成部分，步态识别正在被广泛应用。步态识别主要的优点如下：1)步态识别是一种可以实现远距离身份认证的生物识别技术，在操作上具有很大的灵活性和便捷性；2)步态识别可以在不影响被观测者的情况下，完成特征的提取和识别过程，是一种真正的非侵扰式的生物识别技术；3)步态识别可以在被观测者毫不知情的情况下，完成整个身份识别过程，具有很强的隐蔽性；4)步态识别对视频分辨率要求不高，即便在分辨率相对较低的情况下，也可以完成身份识别的任务。因此，需要一种有效的步态识别算法。

在现实场景下,步态噪声，如被观察者携带的背包、挎包、大衣等，会对步态识别的准确性造成影响，有时甚至是严重影响。

技术实现要素：

发明目的：为了克服现有技术中存在的不足，本发明提供了一种有效提供步态识别效率和准确率的基于隐马尔科夫(hmm)模型的步态去噪方法。

发明内容：为解决上述技术问题，本发明提出一种基于隐马尔科夫模型的步态去噪方法，包含以下步骤：

步骤1：对采集到的步态轮廓高度进行归一化处理；

步骤2：对采集到的步态序列长度进行归一化处理；

步骤3：根据步骤2得到的归一化后的步态序列建立步态集合概率分布序列；

步骤4：采用隐马尔可夫模型计算获得每个步态集合概率的阈值，根据得到的每个步态集合概率的阈值获得步态的先验估计的结果；

步骤5：将步态集合和对应的步态集合概率分布相乘，再做二值化处理，得到步态集合的似然估计的结果；

步骤6：结合步骤4中获得先验估计的结果和步骤5中获得的步态集合的似然估计的结果对步态进行后验估计；

步骤7：在步态后验估计的基础上进行步态回归。

进一步，步骤2中所述对步态序列长度进行归一化处理，将长度为l的步态序列聚类成长度为n的步态序列。

进一步，步骤3中所述步态集合概率分布序列建立方法为，通过累加的方式得到每个步态集合的累加和，然后对相对应的图像做归一化操作得到每个步态集合的概率分布。

进一步，根据公式获得第j个步态集合的累加和其中，d为步态集合的长度，sj,d表示第j个步态集合sj的第d个步态轮廓；根据公式获得第j个步态集合的概率分布pj；其中，sj,d(x,y)表示第j个步态集合sj的第d个步态轮廓图片中的像素点坐标，x表示步态轮廓图片归一化后的像素高度，y表示步态轮廓图片归一化后的像素宽度。

进一步，步骤7中所述在步态后验估计的基础上进行步态回归，对步态后验估计做逆归一化操作。即将归一化之后的步态还原到原始的步态尺寸，这样原始步态的静态属性得以保留。

工作原理：本发明建立步态轮廓概率分布，并将其看作步态的先验概率，使用隐马尔科夫模型的解码算法，找出每个步态图像所对应的步态先验概率，进而得到步态的后验估计并结合原始步态的静态属性的恢复，有效的去除了背包、大衣、肢体等缺失引入的步态噪声。

有益效果：与现有技术相比，本发明主要是在步态去噪步骤中使用静态属性结合动态属性的方式提高步态去噪性能；显著的提高了步态识别的准确率，减少了携带物、大衣、裙子等因素对步态识别效果的影响。在所述建立步态概率分布序列中，目的是对一个二值化的步态图像，得到其所有像素值为1的像素属于该步态的概率，然后取一个阈值，将步态图像中低概率的像素去掉，从而实现步态去噪，有助于显著提高步态识别的准确率。

附图说明

图1为本发明的流程图；

图2为观察步态和隐状态的间距与混淆概率之间的关系图；

图3位本发明提供的方法对步态识别结果的影响。

具体实施方式

下面结合附图对本发明的技术方案作进一步解释。

如图1所示，本实施例公开的基于hmm模型的步态去噪方法，包括以下步骤：

步骤1：步态轮廓高度归一化：为了消除步态轮廓高度引起的差异，按照步态轮廓高度，对所有的步态序列中的步态轮廓做了归一化处理。将所有的步态轮廓的高度都缩放到120像素，步态的宽度做等比例缩放，并将步态轮廓图片归一化成为120×80px的图片。本实施例中共有k个步态序列。

步骤2：步态序列长度归一化：在步态序列完成高度归一化之后，对于一个连续的步态序列，需要将其抽象成为n个离散的步态集合，并用这n个离散的步态集合来表示一个完整的步态周期。

使用knn算法将长度为l(l>n)的步态序列聚类成为长度为n的步态序列。首先将一个步态序列等距划分为n个步态集合,然后根据步态轮廓的距离在相邻步态集合之间做聚类操作，保证任意一个步态集合中至少有一个步态轮廓。

步骤3：建立步态集合概率分布序列：在对所有的步态序列完成高度归一化和长度归一化之后，对于每一个步态序列都可以得到n个步态集合，每个步态集合中有d个步态轮廓，其中，对于第j个归一化的步态集合为sj，通过累加的方式求得该步态集合的累加和然后对其做图像归一化操作，即得到第j个步态集合的概率分布pj。其中x表示步态轮廓图片归一化后的像素高度，y表示步态轮廓图片归一化后的像素宽度：

其中d为步态集合的长度，sj,d表示第j个步态集合sj的第d个步态轮廓。

步骤4：初始化隐马尔可夫模型：为了得到步态集合所对应的步态集合概率的阈值，因此使用了已有的隐马尔科夫模型viterbi算法。

hmm表示为：λ＝(n,m,π,a,b)，参数具体意义如下：

1)n是隐状态的数目，本实施例中，n为归一化后的每个步态序列长度。

2)m是观察状态的数目，即输入的一个步态周期的长度，不同的步态周期其对应的m不一定相同。

3)π是隐状态的初始向量。

4)a为隐状态的状态转移矩阵。

5)b为隐状态的混淆矩阵。

隐状态的状态转移矩阵a和混淆矩阵b则需要基于现有的步态数据集估计，下面是对矩阵a和b的估计过程。

步骤5：计算混淆矩阵：隐马尔科夫模型的混淆矩阵b中的元素bj,k是指在第t次迭代时给定hmm中的一个隐状态qt时，得到第t次迭代时的观察状态ot的概率，所以：

采用混合高斯密度函数,初始估计混淆矩阵b：

其中的是ot到第j个隐状态集合均值的距离。如图2所示，μj＝0.5时上式的函数曲线图，观察步态和隐状态之间的间距越小，其对应的混淆概率b就越大，即观察状态ot的隐状态qt为pj的概率就越大。

其中μj是一个和步态隐状态相关的参数，其定义为：

其中|sj|为隐状态集合sj的势。

步骤6：计算转移矩阵：对输入的步态序列做循环边界处理，即对于输入的步态序列o＝{s1,s2,…,sn}，在其序列的最末端添加第一个步态集合，得到扩展后的步态序列eo＝{s1,s2,…,sn,s1}，并基于扩展步态序列eo对转移矩阵a做初始化，得到初始化转移矩阵a中的元素a⁰(s，v)：

其中，a⁰(s，v)表示在初始状态下第s个隐状态集合转移为第v个隐状态集合的概率，|st(s,v)|为隐状态转移集合st(s,v)的势，|ss|表示第s个隐状态集合的势。转移矩阵的行数为n，列数为n+1。

st(s,v)＝{ss∈eot&&sv∈eot+1}；

其中，eot表示第t次迭代时扩展后的步态序列，eot+1表示第t+1次迭代时扩展后的步态序列，ss表示第s个隐状态集合，sv表示第v个隐状态集合。

因此在扩展步态序列下得到的结果表示为：

其中，k是步态序列的总数，|sv|表示第v个隐状态集合的势，转移矩阵a的初始值如下，由于对步态序列的边界做了扩展，所以转移矩阵a的初始值中最后一行第一个元素是非零的。

为了对转移矩阵a进行迭代估计，需要定义前向局部变量和后向局部变量

因此可得出前向局部变量计算公式为:

后向局部变量为：

其中，上标k指的是所有步态序列中第k个步态序列，t是迭代次数。pj是第j个隐状态，即第j个步态集合概率分布。表示在第k个步态序列中第t次迭代时的观察状态；表示在第k个步态序列中第m次迭代时的观察状态。

因此可得出后向局部变量计算公式为：

其中是第t次迭代时转移矩阵a的元素，当转移矩阵a不再变化时，停止迭代。

步骤7：对步态进行后验估计：相当于对步态的重构，将二值化图像看做是似然估计，即将与人相连的标记为1，背景为0；将步骤1-6算得步态轮廓概率分布看做是步态的先验估计，然后将两者相乘，就可以对存在缺失的步态做填补。如果不进行后验估计，那么手套被作为噪声去除，得到的步态就是没有手。主要是通过步骤4到步骤6中完成了对hmm模型λ＝(n,m,π,a,b)的参数估计，任意一个步态轮廓都有一个与之对应的步态概率分布，然后采用hmm解码算法，得到步态序列中每个步态轮廓所对应的步态概率分布，看作是步态轮廓的先验概率。将步态轮廓和对应的步态概率分布相乘，再做二值化处理，即得到了去噪之后的步态轮廓，看作是步态轮廓的似然估计。将先验概率和似然估计相乘，得到的就是步态的后验估计，相对于直接提取到的似然估计而言，步态后验估计可以去除背包、大衣等人为噪声，更好的描述步态的原始特征。

步骤8：在步态后验估计的基础上进行步态回归：基于步骤7得到经过后验估计后的步态，步态后验估计在去除了步态中的大部分噪声的同时，由于步态是经过了尺寸上的归一化操作，所以后验估计之后的步态所得到的是步态的动态属性，步态的静态属性被归一化处理掉。

为了同时将步态的动态属性和静态属性用于步态识别之中，在步态后验估计的基础上，进行步态回归。步态回归的主要思想是对步态后验估计做逆归一化操作，即将归一化之后的步态还原到原始的步态尺寸，这样原始步态的静态属性得以保留。具体算法如下：

步骤801：输入步骤7中得到的步态后验估计图片，步态轮廓对应的原始高度h；

步骤802：计算高度缩放比例α，这里归一化的高度h＝120；

α＝h/h；

步骤803：计算等比例缩放之前的宽度w，这里归一化的宽度为w：

步骤804：计算步态质心c，并按照xc将步态对齐在逆归一化之后的尺寸为h×w的图片中，这里q是逆归一化之后步态轮廓图片中前景像素点的数目，p(xw,yh)是第i个像素点及其坐标；xw表示逆归一化之后步态轮廓图片中前景像素点的横坐标，yh表示逆归一化之后步态轮廓图片中前景像素点的纵坐标。

通过步态质心c就可以还原原始步态的静态属性。

如图3所示，采用本发明提供的方法步态去噪结果得到显著提升。

以上所述仅是本发明的优选实施方式，应当指出：对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。