本发明专利技术公开了一种融合目标空间分布和时序分布特征子空间的视觉跟踪算法，主要包含步骤：首先将目标区域特征从三个角度进行展开，分别进行子空间学习；其次运用一种更新策略来在线地更新三个模式下目标的均值与子空间。最后通过张量重构将三个子空间模式统一到一个似然函数中用于候选观测图像的评价。本发明专利技术实现了有效的目标跟踪方法，是一种通用的方法。实验结果表明，相对于其他经典子空间跟踪算法，该发明专利技术更加有效和鲁棒，具有很好的应用前景。
【技术实现步骤摘要】
本专利技术涉及计算机视觉领域，特别是涉及一种。
技术介绍
目标跟踪是计算机视觉领域中非常重要的研究问题，因为它是高层视觉问题如运动分析和行为识别等研究的基础。总的来说，目前的目标跟踪算法主要两个关键性的问题:(I)表观模型；(2)跟踪框架。一般来说，表观模型就是如何对目标物体进行有效的表达，并且进行实时的更新。 目标模板，是一种最直接的目标建模方法，不过该模型缺乏判别性和鲁棒性。尽管目标区域的颜色直方图对于目标尺度、旋转以及非刚性形变较为鲁棒，但是由于其忽略了目标表观的颜色空间分布信息，存在一定的缺陷。虽然基于核密度估计的表观模型很好的解决了这一缺陷，不过换来的代价是计算与存储复杂度的增加。另外，基于条件随机场的表观模型通过马尔可夫随机场来建模邻近像素之间的内在关系，但是其训练代价非常巨大。近年来，基于子空间学习的表观模型由于其子空间不变假设更为合理，因而被广泛地应用于视觉跟踪领域。但是该模型在训练时需要足够多的样本，在实际运用中很难达到实时性的要求。基于此，Levy和Lindenbaum提出了序列KL(Sequential Karhunen-Loeve)变换算法用于增量地学习图像的特征基。Lim等扩展了序列KL变换算法， 同时对目标图像的均值和特征基进行增量更新，并将该算法首次应用于目标的视觉跟踪。 而后，鲁棒估计策略，Yang的基于数据驱动的加强自适应方法，Liao的基于鲁棒卡尔曼滤波的跟踪方法以及Gai和Stevenson基于动态模型的方法，虽然在某些特定的场景中获得了较好的跟踪性能，但是有一定的不足之处上述所有基于子空间的跟踪算法首先要将图像展成一维向量，目标表观的空间分布信息几乎完全丢失，从而使得模型对目标表观的全局性变化以及噪音非常敏感。针对这一缺点，Hu等引入张量思想，在一定程度上起到了效用。不过由于其在使用R-SVD更新过程只保留了前R个较大特征值所对应的特征向量带来了误差，并且随着跟踪进行时，误差会逐步累计，导致模型漂移。虽然一种基于动态张量分析的模型避免了上述误差，得到了更加精确的结果，但是由于小样本问题使得计算得到的协方差矩阵无法描述样本的分布情况，从而导致子空间的计算退化。
技术实现思路
为了克服传统的基于子空间学习的表观模型只关注于目标区域特征之间的时间关联性，而忽略了目标区域的空间分布信息。本专利技术提供一种，相对于其它的经典子空间跟踪算法，该跟踪方法更加有效和鲁棒。为了实现上述目的，本专利技术采用了以下的技术方案一种，包括以下步骤(I)运用张量理论将目标图像观测的训练样本表示为一个三阶张量，并将此张量从图像的宽度、图像的高度和图像时间轴这三个不同的方向展开成三个分解矩阵；其中，前两个分解矩阵的列样本子空间对于目标的空间分布信息；第三个分解矩阵的行样本子空间对应于目标的时间分布信息；(2)对于前两个分解矩阵，计算各分解矩阵的样本协方差矩阵，并通过对协方差矩阵做特征分解得到目标水平和垂直方向的模式子空间；对于第三个分解矩阵，采用奇异值分解算法获得目标在时间轴上的模式子空间；(3)对于给定的目标候选区域，计算其在目标表观空间上的重构误差，并由此定义目标候选区域的观测似然函数；(4)根据定义的观测似然函数，运用粒子滤波器计算目标状态的后验概率，采用最大后验估计策略从一系列候选观测值中挑选后验概率最大的作为当前时刻的目标状态；(5)利用跟踪的结果，对对步骤(2)中获得的目标模式子空间进行更新，并在更新过程中加入遗忘因子；其中，遗忘因子是一种介于O和I之间的权重，在子空间模型中越大的遗忘因子意味着模型越关注近期的数据。进一步说，所述的步骤(I)具体包括以下子步骤首先，将目标图像观测的训练样本用一个三阶张量表示；其次，将获得的张量从宽度、高度、时间这三个不同的方向展开成三个分解矩阵；更进一步说，所述的步骤(2)具体包括以下子步骤首先，对于前两个分解矩阵，计算各分解矩阵的样本协方差矩阵，通过对协方差矩阵做特征分解得到各个模式的子空间；其次，对于第三个分解矩阵，采用奇异值分解算法获得目标在时间轴上的模式子空间。更进一步说，所述的步骤(3)具体包括以下子步骤首先，求出重构误差；其次，由重构误差定义似然函数。更进一步说，所述的步骤(4)具体包括以下子步骤首先，假定目标在连续两帧图像之间的运动是仿射运动，目标的状态由两帧之间的目标仿射运动参数表征，并假定目标状态动态转移概率模型为高斯模型；其次，利用似然函数计算目标状态的后验概率；最后，采用最大后验估计策略从从一系列候选观测值中挑选后验概率最大的作为当前时刻的目标状态。更进一步说，所述的步骤(5)具体包括以下子步骤对于前两个分解矩阵，首先，对步骤⑵中获得的样本矩阵的协方差矩阵进行更新；其次，将遗忘因子加入更新方程中，得到更符合实际情况的协方差矩阵。对于第三个分解矩阵，首先，将新获得的数据矩阵进行QR分解；其次，对融合新旧数据的数据矩阵进行奇异值分解。5本专利技术的有益效果是I、本专利技术所提出的跟踪算法是一种通用的方法，对任意跟踪目标类型都适用；2、本专利技术所采用的融合目标空间分布和时序分布特征子空间的视觉跟踪算法，充分挖掘了目标区域的垂直和水平方向上的空间分布信息以及目标区域特征之间的时间关联性，使得算法相较于经典算法降低了对目标表观的全局性变化以及噪音的敏感程度；3、本专利技术提出一种更新策略，对于水平和垂直方向的两个空间模式，通过更新协方差矩阵在线地更新目标的空间模式分布，有效地提高了更新效率；4、本专利技术提出一种更新策略，对于时间分布模式，采用增量奇异值分解方法在线地更新目标的时序模式分布，不但解决了由于样本个数远小于样本维数而产生的子空间退化问题，同时也一定程度上缓解了计算中引入的模型漂移问题。附图说明下面结合附图和实施例对本专利技术进一步说明。图I是本专利技术的跟踪系统的整体框架；图2是本专利技术的图像观测分解示意图；图3是本专利技术的三阶张量展开示意图。具体实施例方式下面通过实施例对本专利技术进行具体的描述，只用于对本专利技术进行进一步说明，不能理解为对本专利技术保护范围的限定，该领域的技术人员可根据上述专利技术的内容对本专利技术作出一些非本质的改进和调整。如图I所示，图I为本专利技术整体框架。本专利技术为一种基于多角度子空间融合的目标跟踪方法，本专利技术的方法具体运行的硬件和编程语言并不限制，用任何语言编写都可以完成，为此其它工作模式不再赘述。本专利技术的实施例采用一台具有3. 2G赫兹中央处理器和IG字节内存的奔腾4计算机并用Matlab语言编制了序列粒子群优化跟踪框架的工作程序，实现了本专利技术的方法，本专利技术的包括以下步骤计算基于张量子空间的展开、协方差矩阵的更新、R-SVD算法即增量的奇异值分解算法、粒子滤波器Condensation以及最大后验估计(MAP)等模块，具体步骤如下所述(I)将目标图像观测的训练样本表示为一个三阶张量Oe分"2"3来表示，其中I1, 12和I3分别表示为图像观测的宽度、高度，以及在时间轴上的长度。根据张量理论，该三阶图像观测张量可以按I1, 12和13方向分别展开成三个分解矩阵Αω，Α⑵，A(3)与图2进行对应。对于Αω和A⑵的列样本协方差矩阵，Α(3)的行样本协方差矩阵分布计算如下=Σ ~ "I )(^(1) " "I Yi乙2 — 2、— "2 為2、- "2本文档来自技高网...
【技术保护点】
【技术特征摘要】
【专利技术属性】
技术研发人员：张笑钦，
申请(专利权)人：温州大学，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。