CN104331700B - 基于轨迹能量扩散图的组行为识别方法 - Google Patents

基于轨迹能量扩散图的组行为识别方法 Download PDF

Info

Publication number
CN104331700B
CN104331700B CN201410712945.6A CN201410712945A CN104331700B CN 104331700 B CN104331700 B CN 104331700B CN 201410712945 A CN201410712945 A CN 201410712945A CN 104331700 B CN104331700 B CN 104331700B
Authority
CN
China
Prior art keywords
track
track energy
energy
energy dissipation
fusion
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201410712945.6A
Other languages
English (en)
Other versions
CN104331700A (zh
Inventor
王生生
张同力
温长吉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jilin University
Original Assignee
Jilin University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Jilin University filed Critical Jilin University
Priority to CN201410712945.6A priority Critical patent/CN104331700B/zh
Publication of CN104331700A publication Critical patent/CN104331700A/zh
Application granted granted Critical
Publication of CN104331700B publication Critical patent/CN104331700B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/40Extraction of image or video features
    • G06V10/50Extraction of image or video features by performing operations within image blocks; by using histograms, e.g. histogram of oriented gradients [HoG]; by summing image-intensity values; Projection analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/40Extraction of image or video features
    • G06V10/46Descriptors for shape, contour or point-related descriptors, e.g. scale invariant feature transform [SIFT] or bags of words [BoW]; Salient regional features
    • G06V10/462Salient features, e.g. scale invariant feature transforms [SIFT]
    • G06V10/464Salient features, e.g. scale invariant feature transforms [SIFT] using a plurality of salient features, e.g. bag-of-words [BoW] representations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/40Extraction of image or video features
    • G06V10/50Extraction of image or video features by performing operations within image blocks; by using histograms, e.g. histogram of oriented gradients [HoG]; by summing image-intensity values; Projection analysis
    • G06V10/507Summing image-intensity values; Histogram projection analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V30/00Character recognition; Recognising digital ink; Document-oriented image-based pattern recognition
    • G06V30/10Character recognition
    • G06V30/19Recognition using electronic means
    • G06V30/192Recognition using electronic means using simultaneous comparisons or correlations of the image signals with a plurality of references
    • G06V30/194References adjustable by an adaptive method, e.g. learning
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V2201/00Indexing scheme relating to image or video recognition or understanding
    • G06V2201/07Target detection

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Multimedia (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Databases & Information Systems (AREA)
  • Image Analysis (AREA)

Abstract

本发明公开一种基于轨迹能量扩散图的组行为识别方法,该方法包括以下步骤:一、对原始视频进行多目标检测,多目标跟踪,提取多目标轨迹;二、把多目标轨迹转换为轨迹能量块图;三、利用能量扩散过程由轨迹能量块图生成轨迹能量扩散图;四、由轨迹能量扩散图生成轨迹能量等势线图;五、对轨迹能量等势线图进行特征提取,利用基于典型相关分析的特征融合算法得到融合特征向量;六、对融合特征向量进行分类识别。本发明提出的轨迹能量块图可以很好的表现并保存组行为轨迹的时间信息;轨迹能量扩散图可以平滑因运动随机性产生的噪声;利用基于典型相关分析的特征融合算法实现多特征融合,从而提高识别精度。

Description

基于轨迹能量扩散图的组行为识别方法
技术领域
本发明涉及图像特征提取及模式识别
背景技术
组行为识别主要是对组行为轨迹进行分析,并在此基础上判断运动目标的行为。近年来由于目标检测与跟踪技术取得了显著的进展,为基于轨迹的多目标组行为识别与理解技术提供支持。由于组行为识别在智能视频监控、虚拟现实、视频检索等方面有着广泛的应用前景,因此越来越受到计算机视觉领域研究者的广泛关注。
目前主要的行为识别算法有基于模板匹配的行为识别算法和基于状态空间的行为识别算法。基于模板匹配的算法可以分为帧对帧匹配算法和融合匹配算法。帧对帧匹配算法是指先将待识别视频中图像序列逐帧进行特征提取,然后与已建立的模板的特征序列逐帧进行匹配。融合匹配算法是指整个行为过程融合为一个或者若干个固定数目的模板,通过特征提取建立行为的模板,然后再利用这些模板进行匹配。基于模板匹配的行为识别算法计算虽然相对简单,但是对行为的时间间隔敏感,提取的特征不能有效的包含行为的时间信息,因此在处理复杂的组行为时有很大的局限性。HMM(隐马尔可夫模型)在基于状态空间的识别算法中应用最为广泛,它是一种有效的时变信号的处理算法,隐含了对时间的校正,并提供了学习机制和识别能力,虽然能够保留轨迹的时间信息,却不利于获取大量的训练数据。因此,我们将组行为轨迹转换为轨迹能量块,通过轨迹能量块的灰度变化反映组行为轨迹的时间信息,如图6所示,轨迹能量块的灰度沿着轨迹的方向由白色渐变到黑色。
虽然轨迹能量块可以通过灰度变化反映组行为轨迹的时间信息,但是由于人体运动属于非刚体运动具有随机性,直接使用轨迹能量块图作为特征,识别效果会受运动随机性产生的噪声影响。虽然有些算法利用高斯过程估计或滤波算法来处理这种不确定性,但在处理过程中,并没有考虑轨迹的时间信息问题。因此,我们利用热扩散过程对轨迹能量块图进行扩散处理得到轨迹能量扩散图,如图7所示,由于运动随机性产生的噪声可以被轨迹能量扩散图平滑掉,使得行为识别可获取更高的精度。
综上所述,轨迹能量扩散图不仅可以很好表现并保存组行为轨迹的时间信息,还可以平滑运动随机性产生的噪声,提高识别精度。因此,我们设计了基于轨迹能量扩散图的组行为识别方法。
发明内容:
为解决现有行为识别算法存在的问题,如提取的特征不能很好地保留轨迹的时间信息,运动随机性产生的噪声对识别结果的影响等。本发明提出了一种基于轨迹能量扩散图的组行为识别方法,发明内容主要包括:本识别方法的流程框架,轨迹能量块图,轨迹能量扩散图,轨迹能量等势线图,使用基于典型相关分析的特征融合算法对分别使用AHDH算法和BSM算法提取的能量等势线图特征向量进行特征融合。
一种基于轨迹能量扩散图的组行为识别方法,其特征在于:至少包括以下步骤:
步骤一、对原始视频进行多目标检测,多目标跟踪,提取多目标轨迹;
步骤二、将多目标轨迹转换为轨迹能量块图:把图形窗口分成大小相同且互不重叠的区域块,轨迹经过的区域块定义为轨迹能量块,一组组行为轨迹转换为轨迹能量块图,同时我们为轨迹附加一个衰减系数,轨迹的能量值会随轨迹延伸的方向不断减小,轨迹能量块的灰度也由白色渐变到黑色(如图6所示),通过与灰度条的对比,时间信息可以很好的通过轨迹能量块的灰度信息表现出来;
步骤三、利用能量扩散过程由轨迹能量块图生成轨迹能量扩散图:由于人的运动行为具有随机性,直接使用轨迹能量块图作为特征,识别效果会受运动随机性产生的噪声影响,由能量扩散过程生成的轨迹能量扩散图可以平滑这种噪声,并且保存轨迹的时间信息;
步骤四、由轨迹能量扩散图生成轨迹能量等势线图:
(1)使用追踪法把轨迹能量扩散图转换为轨迹能量等势线图,
(2)使用曲线平滑算法(二次B样条曲线拟合算法等)对轨迹能量等势线做平滑处理;
步骤五、对轨迹能量等势线图进行特征提取:
(1)用形状描述子AHDH算法描述轨迹能量等势线图得到AHDH特征向量,
(2)用形状描述子BSM算法描述轨迹能量等势线图得到BSM特征向量,
(3)用基于典型相关分析的特征融合算法对AHDH特征向量和BSM特征向量进行融合,得到融合特征向量,该算法将两组特征矢量之间的相关性特征作为有效判别信息,既达到了信息融合的目的,又消除了特征之间的信息冗余,为两组特征融合用于分类识别提供了新的途径;
步骤六、对融合特征向量进行分类识别。融合特征向量作为输入,利用支持向量机实现组行为分类识别。
有益效果:
与现有技术相比,采用本发明所述的设计方案,可以达到以下技术效果:
1、提取的特征由于很好的保存了组行为轨迹的时间信息,能够更为准确地反映组行为的时间特性,使得识别精度大幅度提高;
2、利用能量扩散过程生成轨迹能量扩散图,可以平滑由于运动随机性而产生的噪声,减小噪声对识别结果的影响。
3、分别使用AHDH算法和BSM算法对轨迹能量等势线图进行特征提取与描述,可以更准确的获取组行为的轨迹特征,并且完整的保存了组行为轨迹的时间信息,因此可以取得更好的识别效果。
4、使用基于典型相关分析的特征融合算法对特征向量进行融合,保证了组行为轨迹特征信息的多样性,提高识别精度。
附图说明:
图1.方法框架流程图;
图2.图像进行两次迭代分割;
图3.空间区域关键点;
图4.描述物体形状的特征直方图;
图5.相聚行为的原始轨迹图;
图6.相聚行为的能量块图;
图7.相聚行为的轨迹能量扩散图;
图8.相聚行为的轨迹能量等势线图;
图9.分离行为的原始轨迹图;
图10.分离行为的能量块图;
图11.分离行为的轨迹能量扩散图;
图12.分离行为的轨迹能量等势线图。
具体实施方式:
步骤一、对原始视频进行多目标检测,多目标跟踪,提取多目标轨迹
(1)利用帧差法进行多目标检测;
(2)利用粒子滤波算法进行多目标跟踪;
(3)提取多目标轨迹。
步骤二、把多目标轨迹转换为轨迹能量块图
假设在当前组行为中一共有j条轨迹,则某条轨迹经过轨迹能量块i时,该轨迹能量块的能量值Ei为:
其中Q为衰减系数且为某一常数,Q值确定视实际情况而定。E为常数,设置为1,E(f,j)为第j条轨迹初始时的图像帧数,E(l,j)为第j条轨迹离开能量i块时的图像帧数,Judge()用来判断两帧数的差值为几位整数,例如,他们的差值为99,为两位整数。由公式(1)可以看出,沿着轨迹延伸的方向,Ei值越来越小。
步骤三、利用能量扩散过程由轨迹能量块图生成轨迹能量扩散图
得到轨迹能量块的能量值之后,利用能量扩散过程获得轨迹能量扩散图。轨迹能量块i经过能量扩散过程后,能量值为:
其中,N为能量块的总个数,Ek为k处能量块的能量值,P为空间扩散系数,d(i,k)为能量块i和k之间的距离。
步骤四、由轨迹能量扩散图生成轨迹能量等势线图
本文中生成轨迹能量等势线图就是把轨迹能量扩散图中有相同能量值的点用线连接起来。轨迹能量等势线可分为从边界出发到边界结束的轨迹能量等势线和内部封闭的轨迹能量等势线,此处生成轨迹能量等势线我们使用追踪法。追踪法的原理:首先在区域边界或内部网格的边上求得一个等值点,然后由该点出发,判断下一个等值点的坐标,直到下一个等值点落在区域边界上或与起点重合,则完成对该条轨迹能量等势线的追踪。
具体算法描述如下:
(1)依据追踪法的原理,生成轨迹能量等势线图
步骤五、对轨迹能量等势线图进行特征提取。
(1)用形状描述子AHDH算法描述轨迹能量等势线图得到AHDH特征向量
AHDH算法把轨迹能量等势线图看成是黑色像素在白色区域中的分布,在每次迭代中根据像素分布的质心将图像划分为若干个区域(如图2,图像进行两次迭代分割),并计算每阶段子区域与其父区域的密度比例,最终连接所有迭代下的不同子区域的密度值构成整幅图的密度分布特征。
设FV为新的特征向量。经过下列公式:
第L次划分的特征向量:
换成
最后量化可得新的特征向量:
其中,为第L次划分区域产生的子区域的密度和相对密度,的第j个子区域中黑色像素的个数,的第j个子区域的面积。
(2)用形状描述子BSM算法描述轨迹能量等势线图得到BSM特征向量
BSM算法使用一个网格将待识别的目标划分为一系列的空间区域,进而获得这些空间区域的关键点,如图3所示,并计算这些关键点与相邻区域的关系,从而得到一组描述物体形状的特征向量(特征直方图表示,如图4所示),简称BSM算子。具体算法步骤如下:
Step 1:对于给定的图I,获得形状点集合S;
Step 2:将I划分为n×n个子区域,记为R={r1,r2,...,rn×n},且ci是每个区域ri的中心点,i∈[1,...,n2];
Step 3:设N(ri)是区域ri的邻居区域,定义为:
N(ri)={rk|rk∈R,||ck-ci||<2|g|},其中g是元胞的大小;
Step 4:设ri x为包含点x的区域;
Step 5:初始化特征向量v为v(i)=0,
Step 6:For each x∈S,D=0;
For each
di=d(x,ri)=||x-ci||2
End For
更新特征向量v:
End for
Step 7:标准化特征向量v:
(3)用基于典型相关分析的特征融合算法对AHDH特征向量和BSM特征向量融合,得到融合特征向量
首先,抽取同一模式的两组特征矢量,建立描述两组特征矢量之间相关性的判据准则函数;然后,依此准则求取两组典型投影矢量集,通过给定的特征融合策略抽取组合的融合特征向量。
步骤六:将融合特征向量作为输入,利用支持向量机对组行为进行分类识别。
(1)利用训练样本对支持向量机分类器进行训练;
(2)对待分类组行为进行分类识别。

Claims (5)

1.一种基于轨迹能量扩散图的组行为识别方法,其特征在于:至少包括以下步骤:
步骤一、对原始视频进行多目标检测,多目标跟踪,提取多目标轨迹;
步骤二、将多目标轨迹转换为轨迹能量块图:把图形窗口分成大小相同且互不重叠的区域块,轨迹经过的区域块定义为轨迹能量块,一组组行为轨迹转换为轨迹能量块图,同时我们为轨迹附加一个衰减系数,轨迹的能量值会随轨迹延伸的方向不断减小,轨迹能量块的灰度也由白色渐变到黑色,通过与灰度条的对比,时间信息可以很好的通过轨迹能量块的灰度信息表现出来;
步骤三、利用能量扩散过程由轨迹能量块图生成轨迹能量扩散图:由于人的运动行为具有随机性,直接使用轨迹能量块图作为特征,识别效果会受运动随机性产生的噪声影响,由能量扩散过程生成的轨迹能量扩散图可以平滑这种噪声,并且保存轨迹的时间信息;
步骤四、由轨迹能量扩散图生成轨迹能量等势线图:
(1)使用追踪法把轨迹能量扩散图转换为轨迹能量等势线图,
(2)使用二次B样条曲线拟合算法对轨迹能量等势线做平滑处理;
步骤五、对轨迹能量等势线图进行特征提取:
(1)用形状描述子AHDH算法描述轨迹能量等势线图得到AHDH特征向量,
(2)用形状描述子BSM算法描述轨迹能量等势线图得到BSM特征向量,
(3)用基于典型相关分析的特征融合算法对AHDH特征向量和BSM特征向量进行融合,得到融合特征向量,该算法将两组特征矢量之间的相关性特征作为有效判别信息,既达到了信息融合的目的,又消除了特征之间的信息冗余,为两组特征融合用于分类识别提供了新的途径;
步骤六、对融合特征向量进行分类识别;融合特征向量作为输入,利用支持向量机实现组行为分类识别。
2.根据权利要求1所述的一种基于轨迹能量扩散图的组行为识别方法,其特征在于:所述的步骤二中我们把组行为轨迹转换为轨迹能量块图,轨迹能量块通过灰度的变化,可以很好的把时间信息表现出来。
3.根据权利要求1所述的一种基于轨迹能量扩散图的组行为识别方法,其特征在于:所述的步骤三中利用能量扩散过程生成轨迹能量扩散图,用于平滑运动随机性产生的轨迹噪声,并且保存轨迹的时间信息。
4.根据权利要求1所述的一种基于轨迹能量扩散图的组行为识别方法,其特征在于:所述的步骤四中我们提出使用轨迹能量等势线图来获取轨迹信息。
5.根据权利要求1所述的一种基于轨迹能量扩散图的组行为识别方法,其特征在于:所述的步骤五中我们首先独立使用AHDA算法和BSM算法分别描述轨迹能量等势线图的形状特征,然 后使用典型相关分析算法对上述两种形状特征进行特征融合。
CN201410712945.6A 2014-11-28 2014-11-28 基于轨迹能量扩散图的组行为识别方法 Expired - Fee Related CN104331700B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410712945.6A CN104331700B (zh) 2014-11-28 2014-11-28 基于轨迹能量扩散图的组行为识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410712945.6A CN104331700B (zh) 2014-11-28 2014-11-28 基于轨迹能量扩散图的组行为识别方法

Publications (2)

Publication Number Publication Date
CN104331700A CN104331700A (zh) 2015-02-04
CN104331700B true CN104331700B (zh) 2017-08-15

Family

ID=52406422

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410712945.6A Expired - Fee Related CN104331700B (zh) 2014-11-28 2014-11-28 基于轨迹能量扩散图的组行为识别方法

Country Status (1)

Country Link
CN (1) CN104331700B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106875421A (zh) * 2017-01-19 2017-06-20 博康智能信息技术有限公司北京海淀分公司 一种多目标跟踪方法及装置
CN110299063B (zh) * 2018-03-22 2022-01-25 阿里巴巴集团控股有限公司 轨迹数据的可视化展示方法及装置
CN110458782B (zh) * 2019-08-15 2022-05-03 深圳市瑞立视多媒体科技有限公司 三维轨迹平滑处理方法、装置、设备及存储介质
CN111258558B (zh) * 2020-02-13 2023-06-30 苏宁智能终端有限公司 运动轨迹的着色处理方法及装置

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102663452A (zh) * 2012-04-14 2012-09-12 中国人民解放军国防科学技术大学 基于视频分析的可疑行为检测方法
CN103258193A (zh) * 2013-05-21 2013-08-21 西南科技大学 一种基于kod能量特征的群体异常行为识别方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8131011B2 (en) * 2006-09-25 2012-03-06 University Of Southern California Human detection and tracking system

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102663452A (zh) * 2012-04-14 2012-09-12 中国人民解放军国防科学技术大学 基于视频分析的可疑行为检测方法
CN103258193A (zh) * 2013-05-21 2013-08-21 西南科技大学 一种基于kod能量特征的群体异常行为识别方法

Also Published As

Publication number Publication date
CN104331700A (zh) 2015-02-04

Similar Documents

Publication Publication Date Title
CN107657279B (zh) 一种基于少量样本的遥感目标检测方法
CN107679491B (zh) 一种融合多模态数据的3d卷积神经网络手语识别方法
CN107610087B (zh) 一种基于深度学习的舌苔自动分割方法
CN103942577B (zh) 视频监控中基于自建立样本库及混合特征的身份识别方法
CN109376572A (zh) 基于深度学习的交通视频中实时车辆检测与轨迹跟踪方法
CN105740945B (zh) 一种基于视频分析的人群计数方法
CN108665481A (zh) 多层深度特征融合的自适应抗遮挡红外目标跟踪方法
CN110837768B (zh) 一种面向珍稀动物保护的在线检测与识别方法
CN107633226B (zh) 一种人体动作跟踪特征处理方法
CN105160310A (zh) 基于3d卷积神经网络的人体行为识别方法
CN107516316B (zh) 一种在fcn中引入聚焦机制对静态人体图像进行分割的方法
CN103902978B (zh) 人脸检测及识别方法
CN108537751B (zh) 一种基于径向基神经网络的甲状腺超声图像自动分割方法
CN106023257A (zh) 一种基于旋翼无人机平台的目标跟踪方法
CN107767416B (zh) 一种低分辨率图像中行人朝向的识别方法
CN103310194A (zh) 视频中基于头顶像素点梯度方向的行人头肩部检测方法
CN104331700B (zh) 基于轨迹能量扩散图的组行为识别方法
CN107833239B (zh) 一种基于加权模型约束的寻优匹配目标跟踪方法
CN112052772A (zh) 一种人脸遮挡检测算法
CN110991397B (zh) 一种行进方向确定方法及相关设备
CN104966305A (zh) 基于运动矢量划分的前景检测方法
CN106529377A (zh) 一种基于图像的年龄估计方法、装置及系统
Mo et al. Background noise filtering and distribution dividing for crowd counting
CN115527269B (zh) 一种人体姿态图像智能识别方法及系统
CN106529441A (zh) 基于模糊边界分片的深度动作图人体行为识别方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20170815

Termination date: 20171128

CF01 Termination of patent right due to non-payment of annual fee