CN103268500B - 一种对行走状态变化具有鲁棒性的步态识别方法 - Google Patents
一种对行走状态变化具有鲁棒性的步态识别方法 Download PDFInfo
- Publication number
- CN103268500B CN103268500B CN201310206820.1A CN201310206820A CN103268500B CN 103268500 B CN103268500 B CN 103268500B CN 201310206820 A CN201310206820 A CN 201310206820A CN 103268500 B CN103268500 B CN 103268500B
- Authority
- CN
- China
- Prior art keywords
- circletimes
- prime
- sample
- gait
- matrix
- 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
Links
- 230000005021 gait Effects 0.000 title claims abstract description 91
- 238000000034 method Methods 0.000 title claims abstract description 29
- 230000008859 change Effects 0.000 title abstract description 18
- 239000011159 matrix material Substances 0.000 claims abstract description 57
- 238000012549 training Methods 0.000 claims abstract description 24
- 238000012360 testing method Methods 0.000 claims abstract description 22
- 230000009466 transformation Effects 0.000 claims description 17
- 230000006870 function Effects 0.000 claims description 14
- 239000013598 vector Substances 0.000 claims description 14
- 230000008878 coupling Effects 0.000 claims description 13
- 238000010168 coupling process Methods 0.000 claims description 13
- 238000005859 coupling reaction Methods 0.000 claims description 13
- 238000000354 decomposition reaction Methods 0.000 claims description 11
- 238000004422 calculation algorithm Methods 0.000 claims description 10
- 238000005259 measurement Methods 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000013507 mapping Methods 0.000 claims description 4
- 238000005457 optimization Methods 0.000 claims description 4
- 230000003595 spectral effect Effects 0.000 claims description 4
- 238000012545 processing Methods 0.000 claims description 3
- 230000002427 irreversible effect Effects 0.000 claims description 2
- 230000007704 transition Effects 0.000 claims description 2
- 238000013519 translation Methods 0.000 claims description 2
- 238000004088 simulation Methods 0.000 claims 1
- 230000017105 transposition Effects 0.000 claims 1
- 238000010801 machine learning Methods 0.000 abstract description 7
- 230000001149 cognitive effect Effects 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 230000000717 retained effect Effects 0.000 description 3
- 210000003462 vein Anatomy 0.000 description 3
- 210000002318 cardia Anatomy 0.000 description 2
- 230000001186 cumulative effect Effects 0.000 description 2
- 230000007123 defense Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- FKOQWAUFKGFWLH-UHFFFAOYSA-M 3,6-bis[2-(1-methylpyridin-1-ium-4-yl)ethenyl]-9h-carbazole;diiodide Chemical compound [I-].[I-].C1=C[N+](C)=CC=C1C=CC1=CC=C(NC=2C3=CC(C=CC=4C=C[N+](C)=CC=4)=CC=2)C3=C1 FKOQWAUFKGFWLH-UHFFFAOYSA-M 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000003542 behavioural effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 210000003128 head Anatomy 0.000 description 1
- 235000008216 herbs Nutrition 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 238000003909 pattern recognition Methods 0.000 description 1
- 210000001525 retina Anatomy 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Landscapes
- Image Analysis (AREA)
Abstract
本发明公开了一种对行走状态变化具有鲁棒性的步态识别方法,建立不同行走状态的步态特征矩阵之间的距离度量表达式,通过机器学习手段,在训练阶段将同一个体的不同行走状态的样本构成联系并耦合到同一像空间中,训练完成后可得到不同行走状态样本各自的投影矩阵;在识别阶段,当测试样本与注册集样本库的行走状态不一致时,在训练得到的不同行走状态的投影轴上分别投影,采用最近邻分类器进行分类。本发明不需要从一种行走状态到另外一种行走状态进行预测估计,可直接采用机器学习解决行走状态变化的步态识别问题。
Description
技术领域
本发明属于模式识别、机器学习领域,具体地说是一种对行走状态变化具有鲁棒性的步态识别方法。
背景技术
“9·11”的悲剧事件引起了全世界对增强国防、恐怖袭击安全以及恐怖袭击之后的自动防护能力的格外重视。生物特征识别技术已被成功地应用在身份确认、门禁系统中,还有可能被应用在机场以及其他的安全敏感场所的恐怖分子的确认上。美国国防高级研究项目署在2000年资助的HID(humanidentificationatadistance)研究项目,就是开发并提高当前远距离下大规模身份识别系统的性能,以具有高可靠性、鲁棒性的身份识别能力。2003年,世界民航组织公布了护照中加入生物特征信息便于进入其他国家边境时确认个体身份的应用规划,这一规划已经在美国、欧盟、南非、日本、澳大利亚、韩国等国家投入运行。生物特征身份证在不久的将来会深入到人们的社会生活中,例如,好多款笔记本都需要人脸识别、指纹识别认证,自动提款机也可以通过生物特征来确认身份。现有的生物特征包括生理特征和行为特征。生理特征包括人脸、指纹、静脉(手指静脉、手背静脉)、虹膜、人耳、手形、掌纹、视网膜、唇形、DNA、体味等;行为特征包括步态、语音、击键、签名笔迹等,这些特征都具有可采集性、普遍性、唯一性和稳定性。步态识别是远距离下唯一可用的生物特征识别技术,它具有低分辨率、受环境影响小、易采集等优点,而且是在个体没有觉察的情况下采集到身份信息的,个体行走姿势也不易模仿和伪装,这些都使得步态识别越来越受到研究人员的关注。
鞋子、背包、衣着变化、受伤、年龄、行走速度、拍摄视角变化等联合因素是严重影响步态识别的,Bouchrika等人[1]分析了这些因素单独影响识别性能的程度。贲晛烨等人[2]提出一种线性插值方法对速度不等的步态识别具有鲁棒性。Enokida等人[3]提出3种预测模型用于将鞋的影响从时间和鞋的共同影响中分离出来,使用正常鞋步态来估计拖鞋步态的预测模型,结果使得带预测的一段时间以后的拖鞋步态比不带预测模型的具有更好的识别性能。Lee等人[4]通过双线性模型解决步速变化的步态识别。彭彰等人[5]利用脚间距的计算方法和动态身体分割方法,用曲线拟合方法得到场景的转换参数,将头的高度、躯干的高度、两腿的长度(取最大值为hl)、人的图像高度HI和一个周期内的脚间距之和这5个参数估算出各部分的实际参数,实现了与视角无关的步态识别算法。解决视角变化的步态识别的方法还有三维步态识别[6]、全景摄像机方法[7]。解决衣着变化的步态识别问题可以采用模型法,如:腿部的钟摆模型[8]、骨架模型[9]、恢复静态身体参数模型[10]、椭圆模型[11]、5棒模型[12]和三维模型[13]。模型法受外界干扰小、特征短、能够描述身体各部位的变化情况,在建模准确时,识别效果好,但是,建模复杂程度大,匹配过程复杂。本专利从机器学习角度来解决行走状态变化(衣着变化、背包)的步态识别问题,大大降低了计算复杂度。
与本发明相关的公开报道包括:
[1]ImedBouchrika,MarkSNixion.Exploratoryfactoranalysisofgaitrecognition.Proc.ofthe8thIEEEInternationalConferenceonAutomaticFaceandGestureRecognition.2008:1-6P
[2]王科俊,贲晛烨.基于线性插值的步态识别算法,华中科技大学学报(自然科学版),2010,38(2):41-44页。
[3]EnokidaShuichi,ShimomotoRyo,WadaTomohito,etal.APredictiveModelforGaitRecognition.Proc.of2006BiometricsSymposium:SpecialSessiononResearchattheBiometricConsortiumConference.2006:1-6P。
[4]ChanSuLee,AhmedElgammal.Gaitstyleandgaitcontent:bilinearmodelsforgaitrecognitionusinggaitre-sampling.Proc.ofthesixthIEEEInternationalConferenceonAutomaticFaceandGestureRecognition(FGR’04).2004:147-152P。
[5]彭彰,吴晓娟,杨军.基于肢体长度参数的多视角步态识别算法.自动化学报.2007,33(2):210-213页。
[6]GuoyingZhao,GuoyiLiu,HuaLi,etal.3Dgaitrecognitionusingmultiplecameras.Proc.ofthe7thInternationalConferenceonAutomaticFaceandGestureRecognition(FGR'06).2006:529-534P。
[7]SugiuraK,MakiharaY,YagiY.Gaitrecognitionbasedonmulti-viewobservationusingomnidirectionalcamera.InternationalComferenceonComputerVision(ACCV).2007:452-461P
[8]CunadoD,NixonMS,CarterJN.UsingGaitasaBiometricviaPhaseweightedMagnitudePectual.LectureNotesinComputerScienceProc.ofAVBPA97.1997:95-102P。
[9]NashJM,CarterJN,NixonMS.ExtractionofMovingArticulatedObjectsbyEvidenceGathering.Proc.oftheNinthBritishMachineVisionConference(BMVC’98).1998:609-618P。
[10]BobickAF,JohnsonAY.Gaitrecognitionusingstatic,activityspecificparameters.Proc.of2001IEEEComputerSocietyConferenceonComputerVisionandPatternRecognition.2001:423-430P。
[11]LeeL,GrimsonWEL.Gaitanalysisforrecognitionandclassification.Proc.ofIEEEInt.Conf.AutomaticFaceandGestureRecognition.2002:148-155P。
[12]ZhangRong,ChristianV,DimitrisM.Humangaitrecognition.Proc.of2004IEEEComputersocietyConferenceonComputerVisionandPatternRecognitionWorkshops(CVPRW'04).2004:1-8P。
[13]UrtasunR,FuaP.3DTrackingforGaitCharacterizationandRecognition.Proc.oftheSixthIEEEInternationalConferenceonAutomaticFaceandGestureRecognition.2004:17-22P。
发明内容
本发明的目的在于提供一种对行走状态变化具有鲁棒性的步态识别方法,是通过机器学习手段,在训练阶段将同一个体的不同行走状态的样本构成联系并耦合到同一像空间中,训练完成后可得到不同行走状态样本各自的投影矩阵;在识别阶段,当测试样本与注册集样本库的行走状态不一致时,在训练得到的不同行走状态的投影轴上分别投影,采用最近邻分类器进行分类。该方法对行走状态变化具有强鲁棒性,具有很好的识别性能。
为了实现上述目的,本发明采用如下技术方案。
本发明是通过如下技术方案实现的:
一种对行走状态变化具有鲁棒性的步态识别方法。主要包括:建立不同行走状态的步态特征矩阵之间的距离度量表达式、训练阶段和识别阶段。训练时,在训练样本集合中建立相似度矩阵,建立目标优化函数,然后对该目标函数解耦得到不同行走状态的步态的各自的变换矩阵,不同行走状态的步态分别在这两个变换矩阵上投影得到矩阵空间下不同行走状态的步态样本特征集合,然后将所有样本特征向量化处理,再采用向量空间下的耦合度量学习得到向量空间下的不同行走状态的步态的各自的变换矩阵,最后得到注册集样本的最终的特征集合。测试样本与注册集样本的行走状态不一致时,利用训练阶段得到的矩阵空间下的变换矩阵和向量空间下的变换矩阵投影变换两次,最终采用最近分类器判定该步态样本所属的类别。
本发明不需要从一种行走状态到另外一种行走状态进行预测估计,可直接采用机器学习解决行走状态变化的步态识别问题。
对于两个不同行走状态的步态特征矩阵集合 和分别表示维数为Dxm×Dxn和Dym×Dyn的空间,其中的步态样本和的距离定义为κC:R表示实数空间。常规的实现步骤是:①通过映射函数fx和fy将Xi和Yj映射到同一个像空间中:fx: 表示维数为Dc×Dr的空间;②在像空间中再进行传统的矩阵距离度量κC:因此
其中,||||C是定义计算κC的距离度量关系,T表示转置。C是定义一个矩阵以代替距离度量关系||||C,C是一个半正定阵。
令C=WcWc T, 则
再令Px=WxmWc,Qx=Wxn,Py=WymWc,Qy=Wyn,有
求取κC(Xi,Yj)中的Px,Qx,Py和Qy即可得到两个不同行走状态的步态样本之间的距离,并且需要利用监督信息(边信息、局部保留)来指导学习得到这四个投影矩阵。
边信息相似约束条件下的不同行走状态的步态样本的匹配学习依靠图模型:若和相似,那么样本i和j的相似性关系其中,为相似关系集合。设和都包含M个步态样本,S表示相似矩阵,其中的每一个元素
局部保留相似约束条件下的不同行走状态的步态样本的匹配学习与边信息相似约束条件下的不同之处在于式(4)中的S矩阵为
或者,
其中,xi是xj分别是Xi和Xj的向量化版本,而且还满足:xi是xj的k-近邻,或者xj是xi的k-近邻,否则,Sij=0。Sij为相似度函数S的元素。式(5)是余弦相似度,式(6)是高斯相似度。t是尺度因子。
此时的不同行走状态的步态样本矩阵的匹配学习准则应该满足相似对应关系点对在像空间中仍然具有相似性:
其中,||||F表示F范数。
初始化 其中,为维数为Dr×Dr的单位阵,和表示维数为下角标所示的全零矩阵。中的样本X1,X2,…,XM和Y1,Y2,…,YM分别向方向为Qx、Qy的空间上投影,令投影后得到的样本集为集合X1=[X1Qx0,X2Qx0,…,XMQx0],其中每一个样本的维数为Dxm×Dr,集合投影后得到的样本集为集合Y1=[Y1Qy0,Y2Qy0,…,YMQy0],其中每一个样本的维数为Dym×Dr,那么,此时的学习准则转换为
其中:代表克罗尼克积。Fh(S)和Fv(S)都是对角阵,其中的每一元素都是S的列或者行的和,
再令: 此时的学习准则函数为
为了使式(9)具有唯一解,这里还需要加入尺度不变约束:(I表示单位阵)和平移不变约束:PTZ1e=0,其中,e为全1、大小为2MDr×1的列向量。求此优化问题,采用谱分解求取EP=λFP,其中 λ为特征值,将求得的特征值按照从小到大排列,取第2到第Dc+1个特征值对应的特征向量按列的顺序拼接排列得到P,这里没有使用第1个特征值对应的特征向量是通过抛弃零特征值对应的特征向量实现中心化。P的维数为(Dxm+Dym)×Dc,那么Px的维数为Dxm×Dc,Py的维数为Dym×Dc,通常是不可逆的,可通过规整化来消除奇异性,其中,τ取接近于0的正数,可取τ=10-6。
和中的数据分别向Px和Py投影之后,只相当于集合和在一个方向上构造了数据间的耦合关系,这里还需要在另外一个方向也构造出数据间的耦合关系。那么,集合 中的样本分别再向方向为Px、Py的空间上投影,得到的数据和分别记作X′1,X′2,…,X′M和Y′1,Y′2,…,Y′M,集合中的样本依次拼接得到X2=[X′1,X′2,…,X′M],其中每一个样本的维数为Dc×Dr,集合中的样本也同样拼接得到Y2=[Y′1,Y′2,…,Y′M],其中每一个样本的维数为Dc×Dr,从新排列集合X2和Y2,有X′2=[(X′1)T,(X′2)T,…,(X′M)T],Y′2=[(Y′1)T,(Y′2)T,…,(Y′M)T]。那么此时的学习准则转换为
再令: 此时的学习准则函数为
Q的求解方式同前谱分解,其维数为Dr×2Dr,Qx、Qy的维数分别为Dr×Dr和Dr×Dr。当然,为了确保特征的紧凑性与有效性,Qx、Qy的维数也可保留到Dr×d′和Dr×d′(d′≤Dr),d′代表列的维数。
在求得投影矩阵Px,Qx,Py和Qy之后,很容易将不同行走状态的步态特征矩阵集合的元素进行耦合距离度量,可完全采用批处理方式,如图2所示。它的主体思想就是将的样本集按列拼接在一起形成一个新的矩阵,Qx初始化为Qx0,将的样本集投影到Qx0上,获得样本集的维数大小为Dxm×(Dr·M),将其向方向为Px的空间上投影,得到新矩阵,其维数大小为Dc×(Dr·M),再将此矩阵从新排列,排列的原则为变换后的样本集中的每个样本取转置,并按照列来拼接,那么,此时样本集构成新矩阵的维数为Dr×(Dc·M),将它向方向为Qx的空间投影得到最终样本集构成的矩阵大小为Dr×(Dc·M),此时的每一个样本的维数变为Dc×Dr,的样本集也进行类似的操作,只不过,初始化投影矩阵为Qy0,两个投影矩阵分别是Py和Qy,最终变换得到的样本集中每一个样本的大小亦为Dc×Dr,这样两个样本集就耦合到同一空间中,而且数据的行、列方向都分别进行了耦合。
该算法的特征提取过程如下:
(1)Qx初始化为 Qy初始化为 两个不同行走状态的步态特征矩阵集合中的样本分别向方向为Qx0、Qy0的空间上投影,集合投影后得到的样本集为集合X1=[X1Qx0,X2Qx0,…,XMQx0],集合投影后得到的样本集为集合Y1=[Y1Qy0,Y2Qy0,…,YMQy0]。
(2)求式(9)的广义特征分解,Px对应于P的第1行到第Dxm行,Py对应于P的第Dxm+1行到Dxm+Dym行。
(3)集合中的样本分别向方向为Px、Py的空间上投影,得到的数据 和 分别记作X′1,X′2,…,X′M和Y′1,Y′2,…,Y′M,集合中的样本依次拼接得到X2=[X′1,X′2,…,X′M],集合中的样本也同样拼接得到从新排列集合X2和Y2,有X′2=[(X′1)T,(X′2)T,…,(X′M)T],Y′2=[(Y′1)T,(Y′2)T,…,(Y′M)T]。
(4)求式(11)的广义特征分解,Qx和Qy分别对应于Q的第1行到第Dr行和第Dr+1行到第2·Dr行。
(5)样本集X′1,X′2,…,X′M和Y′1,Y′2,…,Y′M分别向方向为Qx、Qy的空间上投影,则投影后的集合和分别为X3=[X′1Qx,X′2Qx,…,X′MQx],Y3=[Y′1Qy,Y′2Qy,…,Y′MQy]。
(6)将X3和Y3中的各个样本都向量化并按照列拼接成矩阵a和b,然后再采用向量空间下的耦合度量学习算法得到变换矩阵P1和P2,X3样本集投影得到ya=P1 Ta。
该算法的分类识别过程如下:
假设注册样本集的特征为yak(k=1,2,…,M),对于测试样本B',首先经初始化的Qy0矩阵以及训练得到的Py和Qy矩阵投影得到特征YB′
YB′向量化得到b',它的特征y′b为
如果
其中,yaj是注册样本集的特征,j=1,2,…,M,则B'属于yak所在的类别。
本发明的有益效果:本发明不需要从一种行走状态到另外一种行走状态进行预测估计,可直接采用机器学习解决行走状态变化的步态识别问题。本发明可有效地解决不同分辨率图像的直接匹配、同一个体但不同模态图像的直接匹配、同一个体但由不同设备采集的图像的直接匹配问题。
附图说明
图1是本发明识别方法的流程图。
图2是不同行走状态的步态特征矩阵集合间耦合的示意图。
图3步态能量图,(a)正常行走状态,(b)背包,(c)外套变化的。
图4使用边信息时的特征值分布情况。
图5使用边信息时的注册样本集和测试样本集的特征。
图6高斯形式下随近邻数k和尺度因子t变化的训练识别率。
图7使用高斯相似度时的特征值分布情况。
图8使用高斯相似度时的注册样本集和测试样本集的特征。
图9使用余弦相似度的注册样本集和测试样本集的特征。
具体实施方式
下面结合附图与实施例对本发明作进一步说明。
实验的数据库选择CASIA(B)步态数据库,该库含有3种行走状态,分别为正常步态(记为’nm’),带有背包的步态(记为‘bg’)以及外套变化的步态(记为‘cl’)。以侧面视角下的步态识别为例来说明本专利所提供的方法的实验效果,对于每段步态视频图像都采用步态能量图(gaitenergyimage,GEI)的形式来表达步态特征,且GEI的大小均为64×64像素。如图3所示。
四组对行走状态变化具有鲁棒性的步态识别实验:
(1)注册集由每人第一个正常步态组成;训练集由每人第一个正常步态和第一个带有背包的步态组成;测试集由每人第二个带有背包的步态组成。该数据集记为‘nm&bg’;
(2)注册集由每人第一个正常步态组成;训练集由每人第一个正常步态和第一个外套变化的步态组成;测试集由每人第二个带有外套变化的步态组成。该数据集记为‘nm&cl’;
(3)注册集由每人第一个背包步态组成;训练集由每人第一个背包步态和第一个外套变化的步态组成;测试集由每人第二个外套变化的步态组成。该数据集记为‘bg&cl’;
(4)注册集由每人第一个外套变化的步态组成;训练集由每人第一个外套变化的步态和第一个带有背包的步态组成;测试集由每人第二个带有背包的步态组成。该数据集记为‘cl&bg’;
首先测试边信息相似约束条件下的对行走状态变化具有鲁棒性的步态识别方法的性能,在数据集‘nm&bg’中,初始化然后通过式(9)和(10)的广义特征分解得到投影矩阵Px、Py、Qx和Qy,并通过查看特征值分布情况来确定投影矩阵的维数。具体来讲,式(9)的广义特征分解得到的特征值分布如图4中上面的子图,可以看出特征值由小到大排序,直到第30个特征值以后才迅速上升,但是特征值分布的数量级为100,因为累计的特征值的和越大,我们用来建立样本集之间联系的差别程度就越大,但是保留的特征值个数较少时,累计的特征值的和虽然小,保留的特征又不足以用来分类。均衡两者,这里选择从小到大排列后的第2个到第11个特征值对应的特征向量来组成变换矩阵Px和Py。图4中的下面的子图对应的是通过式(11)的广义特征分解的特征值分布情况,数量级是10-24,十分接近于零,所以Qx和Qy的维数可以保留稍微多一些,这里可以选择从小到大排列后的第2个到第31个特征值对应的特征向量来组成Qx和Qy。这样,两个样本集就在行方向和列方向都进行了耦合,且样本的维数为300维。然后再采用向量空间下的耦合度量学习算法,特征仅保留到1维时,识别率就可以达到100%。
如图5所示的从上到下的子图分别是测试bg集合的124个个体的一维特征,注册nm集合的124个个体的一维特征以及两者之差,其横轴是不同的样本个体,可以看出对于同一个体,即使注册集是正常步态样本,测试集是带有背包的步态样本,采用本专利所提供的方法,两者的一维特征数值大小十分相近。
下面测试局部保留相似约束条件下的对行走状态变化具有鲁棒性的步态识别方法的识别性能,当相似度矩阵选择为式(6)时,将nm与bg集进行耦合,并测试识别的效果,定义训练识别率为测试用于训练的bg识别率,而注册样本集依然为nm。首先,依据式(8)和(9)将nm和bg在一个方向上耦合,且相似度矩阵选择为S3=min{S1,S2},其中S1和S2分别为通过nm集合、bg集合计算出来的相似度矩阵,当k=5,t=300时最佳的训练识别率为0.5645。如图6所示为高斯形式下随近邻数k和尺度因子t变化的训练识别率。然后,再依据式(10)和(11)将nm和bg集在另一个方向上耦合,此时的训练识别率为0.9839。图7的上、下的两个子图分别为通过这两个式子计算得到的特征值分布情况,它与图4的分布走势相类似,数量级亦为100和10-24,这里的Px和Py同样也取从小到大排列后的第2个到第11个特征值对应的特征向量;Qx和Qy也取选择从小到大排列后的第2个到第31个特征值对应的特征向量。这样,两个样本集就在行方向和列方向都进行了耦合,且样本的维数为300维。然后再采用向量空间下的耦合度量学习算法,同样也是特征保留到1维时,识别率就可以达到100%。
如图8所示的从上到下的子图分别是测试bg集合的124个个体的一维特征,注册nm集合的124个个体的一维特征以及两者之差,其横轴是不同的样本个体,可以看出:经特征提取后两者的一维特征数值大小十分相近。
若相似度矩阵为余弦相似度矩阵时,通过实验确定k=1,通过式(9)和(11)的广义特征分解得到的投影变换矩阵也将nm和bg数据集的样本维数降维到300维,然后再采用向量空间下的耦合度量学习算法,同样也是将特征保留到1维时,识别率就可以达到100%。
如图9所示的从上到下的子图分别是测试bg集合的124个个体的一维特征,注册nm集合的124个个体的一维特征以及两者之差,其横轴是不同的样本个体,可以看出:经特征提取后两者的一维特征数值大小十分相近。
将‘nm&bg’的实验结果同‘nm&cl’、‘bg&cl’和‘cl&bg’的都总结在表1中,无论是基于边信息的抑或余弦相似度和高斯相似度下的对行走状态变化具有鲁棒性的步态识别方法都在矩阵耦合阶段将特征维数保留到10×30维,最终用于识别的特征维数仅为1维时,四个数据集下的识别率大小均能达到100%。
表1本专利所提供方法的识别性能
上述虽然结合附图对本发明的具体实施方式进行了描述,但并非对本发明保护范围的限制,所属领域技术人员应该明白,在本发明的技术方案的基础上,本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。
Claims (3)
1.一种对行走状态变化具有鲁棒性的步态识别方法,其特征是,主要包括:建立不同行走状态的步态特征矩阵之间的距离度量表达式、训练阶段和识别阶段;训练时,在训练样本集合中建立相似度矩阵,建立目标优化函数,然后对该目标函数解耦得到不同行走状态步态的各自变换矩阵,不同行走状态的步态分别在这两个变换矩阵上投影得到矩阵空间下不同行走状态的步态样本特征集合,然后将所有样本特征向量化处理,再采用向量空间下的耦合度量学习得到向量空间下的不同行走状态步态的各自变换矩阵,最后得到注册集样本的最终特征集合;测试样本与注册集样本的行走状态不一致时,利用训练阶段得到的矩阵空间下的变换矩阵和向量空间下的变换矩阵投影变换两次,最终采用最近邻分类器判定该步态样本所属的类别;所述不同行走状态的步态特征矩阵之间的距离度量表达式的建立方法为:对于两个不同行走状态的步态特征矩阵集合和分别表示维数为Dxm×Dxn和Dym×Dyn的空间,其中的步态样本和的距离定义为κC:R表示实数空间;①通过映射函数fx和fy将Xi和Yj映射到同一个像空间中: 表示维数为Dc×Dr的空间;②在像空间中再进行传统的矩阵距离度量κC:因此
其中,||||C是定义计算κC的距离度量关系,T表示转置;C是定义一个矩阵以代替距离度量关系||||C,C是一个半正定阵;
令 则
再令Px=WxmWc,Qx=Wxn,Py=WymWc,Qy=Wyn,有
求取κC(Xi,Yj)中的Px,Qx,Py和Qy即得到两个不同行走状态的步态样本之间的距离,并且需要利用监督信息来指导学习得到这四个投影矩阵;
其中:Px,Qx,Py和Qy均为投影矩阵,Wc为C的分解矩阵,Wxm、Wxn、Wym、Wyn为投影变换矩阵。
2.如权利要求1所述的对行走状态变化具有鲁棒性的步态识别方法,其特征是,所述训练阶段的具体步骤为:
边信息相似约束条件下的不同行走状态的步态样本的匹配学习依靠图模型:若和相似,那么样本i和j的相似性关系其中,为相似关系集合;设和都包含M个步态样本,S表示相似矩阵,其中的每一个元素
局部保留相似约束条件下的不同行走状态的步态样本的匹配学习与边信息相似约束条件下的不同之处在于式(4)中的S矩阵为
或者,
其中,xi和xj分别是Xi和Xj的向量化版本,而且还满足:xi是xj的k-近邻,或者xj是xi的k-近邻,否则,Sij=0;Sij为相似度函数S的元素;式(5)是余弦相似度,式(6)是高斯相似度;t是尺度因子;
此时的不同行走状态的步态样本矩阵的匹配学习准则应该满足相似对应关系点对在像空间中仍然具有相似性:
其中,||||F表示F范数;
初始化 其中,为维数为Dr×Dr的单位阵,和表示维数为下角标所示的全零矩阵;中的样本X1,X2,…,XM和Y1,Y2,…,YM分别向方向为Qx、Qy的空间上投影,令投影后得到的样本集为集合X1=[X1Qx0,X2Qx0,…,XMQx0],其中每一个样本的维数为Dxm×Dr,集合投影后得到的样本集为集合Y1=[Y1Qy0,Y2Qy0,…,YMQy0],其中每一个样本的维数为Dym×Dr,那么,此时的学习准则转换为
其中:代表克罗尼克积;Fh(S)和Fv(S)都是对角阵,其中的每一元素都是S的列或者行的和,
再令: 此时的学习准则函数为
为了使式(9)具有唯一解,需要加入尺度不变约束:和平移不变约束:PTZ1e=0,其中,I表示单位阵,e为全1、大小为2MDr×1的列向量;求此优化问题,采用谱分解求取EP=λFP,其中λ为特征值,将求得的特征值按照从小到大排列,取第2到第Dc+1个特征值对应的特征向量按列的顺序拼接排列得到P;P的维数为(Dxm+Dym)×Dc,那么Px的维数为Dxm×Dc,Py的维数为Dym×Dc,通常是不可逆的,通过规整化来消除奇异性,其中,τ取接近于0的正数,可取τ=10-6;
集合中的样本分别再向方向为Px、Py的空间上投影,得到的数据和分别记作X′1,X′2,…,X′M和Y′1,Y′2,…,Y′M,集合中的样本依次拼接得到X2=[X′1,X′2,…,X′M],其中每一个样本的维数为Dc×Dr,集合中的样本也同样拼接得到Y2=[Y′1,Y′2,…,Y′M],其中每一个样本的维数为Dc×Dr,从新排列集合X2和Y2,有X′2=[(X′1)T,(X′2)T,…,(X′M)T],Y′2=[(Y′1)T,(Y′2)T,…,(Y′M)T];那么此时的学习准则转换为
再令: 此时的学习准则函数为
Q的求解方式同前谱分解,其维数为Dr×2Dr,Qx、Qy的维数分别为Dr×Dr和Dr×Dr;
样本集X′1,X′2,…,X′M和Y′1,Y′2,…,Y′M分别向方向为Qx、Qy的空间上投影,则投影后的集合和分别为X3=[X′1Qx,X′2Qx,…,X′MQx],Y3=[Y′1Qy,Y′2Qy,…,Y′MQy];
将X3和Y3中的各个样本都向量化并按照列拼接成矩阵a和b,然后再采用向量空间下的耦合度量学习算法得到变换矩阵P1和P2,X3样本集投影得到
3.如权利要求1所述的对行走状态变化具有鲁棒性的步态识别方法,其特征是,所述识别阶段的具体步骤为:
假设注册样本集的特征为yak(k=1,2,…,M),对于测试样本B',首先经初始化的Qy0矩阵以及训练得到的Py和Qy矩阵投影得到特征YB′
YB′向量化得到b',它的特征y′b为
如果
其中,yaj是注册样本集的特征,j=1,2,…,M,则B'属于yak所在的类别。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310206820.1A CN103268500B (zh) | 2013-05-29 | 2013-05-29 | 一种对行走状态变化具有鲁棒性的步态识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310206820.1A CN103268500B (zh) | 2013-05-29 | 2013-05-29 | 一种对行走状态变化具有鲁棒性的步态识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103268500A CN103268500A (zh) | 2013-08-28 |
CN103268500B true CN103268500B (zh) | 2016-04-06 |
Family
ID=49012127
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310206820.1A Expired - Fee Related CN103268500B (zh) | 2013-05-29 | 2013-05-29 | 一种对行走状态变化具有鲁棒性的步态识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103268500B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107103296A (zh) * | 2017-04-20 | 2017-08-29 | 四川大学 | 基于完全矫正Boosting和子空间学习算法的步态识别方法 |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103745205B (zh) * | 2014-01-22 | 2017-01-11 | 山东大学 | 一种基于多线性均值成分分析的步态识别方法 |
CN104217211B (zh) * | 2014-08-07 | 2017-05-24 | 哈尔滨工程大学 | 基于最优判别耦合投影的多视角步态识别方法 |
CN104463099B (zh) * | 2014-11-05 | 2017-10-31 | 哈尔滨工程大学 | 一种基于图的半监督耦合度量的多角度步态识别方法 |
CN106096528B (zh) * | 2016-06-03 | 2019-10-01 | 山东大学 | 一种基于二维耦合边距Fisher分析的跨视角步态识别方法 |
CN106203255B (zh) * | 2016-06-24 | 2019-04-23 | 华中科技大学 | 一种基于时间对齐的行人重识别方法及系统 |
CN108921114B (zh) * | 2018-07-10 | 2021-09-28 | 福州大学 | 基于多非线性局部保持投影的两视角步态识别方法 |
CN109241870B (zh) * | 2018-08-16 | 2022-02-22 | 中国矿业大学(北京) | 基于步态识别的煤矿井下人员身份识别方法 |
CN109871817B (zh) * | 2019-02-27 | 2021-08-27 | 河北工业大学 | 基于足底压力和自适应集成学习的行走状态识别方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102982323A (zh) * | 2012-12-19 | 2013-03-20 | 重庆信科设计有限公司 | 快速步态识别方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010118233A2 (en) * | 2009-04-08 | 2010-10-14 | University Of Southern California | Cadence analysis of temporal gait patterns for seismic discrimination |
-
2013
- 2013-05-29 CN CN201310206820.1A patent/CN103268500B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102982323A (zh) * | 2012-12-19 | 2013-03-20 | 重庆信科设计有限公司 | 快速步态识别方法 |
Non-Patent Citations (5)
Title |
---|
Improved gait recognition by multiple-projections normalization;Ekinci M等;《Machine Vision and Application,2010》;20100228;第143-161页 * |
基于Trace变换的步态识别算方法;贲晛烨等;《吉林大学学报(工学版)》;20120115;第156-160页 * |
基于步态能量图像和2维主成分分析的步态识别方法;王科俊;《中国图象图形学报》;20091215;第2503-2509页 * |
基于线性插值的步态识别算法;王科俊等;《华中科技大学学报(自然科学版)》;20100215;第41-44页 * |
采用Radon变化和二维主成分分析的步态识别算法;王科俊等;《智能系统学报》;20100615;第266-271页 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107103296A (zh) * | 2017-04-20 | 2017-08-29 | 四川大学 | 基于完全矫正Boosting和子空间学习算法的步态识别方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103268500A (zh) | 2013-08-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103268500B (zh) | 一种对行走状态变化具有鲁棒性的步态识别方法 | |
US9117105B2 (en) | 3D face recognition method based on intermediate frequency information in geometric image | |
Hu et al. | Gait-based gender classification using mixed conditional random field | |
Luo et al. | Robust arbitrary view gait recognition based on parametric 3D human body reconstruction and virtual posture synthesis | |
CN103745205B (zh) | 一种基于多线性均值成分分析的步态识别方法 | |
CN105138980A (zh) | 基于身份证件信息和人脸识别的身份验证方法及系统 | |
Jalal et al. | Daily human activity recognition using depth silhouettes and transformation for smart home | |
Kusakunniran et al. | Automatic gait recognition using weighted binary pattern on video | |
Shu et al. | Age progression: Current technologies and applications | |
Ouarda et al. | MLP Neural Network for face recognition based on Gabor Features and Dimensionality Reduction techniques | |
CN106709508A (zh) | 一种利用特征信息的加权典型相关分析方法 | |
Li et al. | Adapting geometric attributes for expression-invariant 3D face recognition | |
CN104217211B (zh) | 基于最优判别耦合投影的多视角步态识别方法 | |
CN102592150A (zh) | 基于模糊理论决策的双向二维主成分分析的步态识别方法 | |
Du et al. | Face recognition under pose variations | |
Mahalingam et al. | Face verification of age separated images under the influence of internal and external factors | |
Fotiadou et al. | Activity-based methods for person recognition in motion capture sequences | |
Kare et al. | Using bidimensional regression to assess face similarity | |
Vasconcelos et al. | Methods to automatically build point distribution models for objects like hand palms and faces represented in images | |
Mráček et al. | 3D face recognition on low-cost depth sensors | |
Wang et al. | Robust face recognition from 2D and 3D images using structural Hausdorff distance | |
Nayak et al. | Efficient face recognition with compensation for aging variations | |
Kumar et al. | Palmprint Recognition in Eigen-space | |
Deng et al. | Expression-robust 3D face recognition based on feature-level fusion and feature-region fusion | |
Sagonas et al. | Back to the future: A fully automatic method for robust age progression |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160406 Termination date: 20170529 |
|
CF01 | Termination of patent right due to non-payment of annual fee |