CN112115923A - 一种基于直接特征提取的多通道时序步态分析算法 - Google Patents
一种基于直接特征提取的多通道时序步态分析算法 Download PDFInfo
- Publication number
- CN112115923A CN112115923A CN202011083624.6A CN202011083624A CN112115923A CN 112115923 A CN112115923 A CN 112115923A CN 202011083624 A CN202011083624 A CN 202011083624A CN 112115923 A CN112115923 A CN 112115923A
- Authority
- CN
- China
- Prior art keywords
- gait
- lambda
- observing
- alpha
- phase
- 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.)
- Pending
Links
- 230000005021 gait Effects 0.000 title claims abstract description 162
- 238000004458 analytical method Methods 0.000 title claims abstract description 21
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 21
- 238000000605 extraction Methods 0.000 title claims abstract description 18
- 238000012549 training Methods 0.000 claims abstract description 30
- 238000004364 calculation method Methods 0.000 claims abstract description 9
- 210000002683 foot Anatomy 0.000 claims description 27
- 238000007637 random forest analysis Methods 0.000 claims description 16
- 230000033001 locomotion Effects 0.000 claims description 15
- 238000000034 method Methods 0.000 claims description 13
- 210000003414 extremity Anatomy 0.000 claims description 10
- 210000003127 knee Anatomy 0.000 claims description 8
- 230000004044 response Effects 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 6
- 210000003423 ankle Anatomy 0.000 claims description 4
- 239000012535 impurity Substances 0.000 claims description 3
- 210000002414 leg Anatomy 0.000 claims description 3
- 210000001699 lower leg Anatomy 0.000 claims description 3
- 238000002360 preparation method Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 210000001694 thigh bone Anatomy 0.000 claims description 3
- 210000002303 tibia Anatomy 0.000 claims description 3
- 210000000689 upper leg Anatomy 0.000 claims description 3
- 238000005457 optimization Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 4
- 201000010099 disease Diseases 0.000 description 4
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 4
- 238000012935 Averaging Methods 0.000 description 2
- 208000018737 Parkinson disease Diseases 0.000 description 2
- QBWCMBCROVPCKQ-UHFFFAOYSA-N chlorous acid Chemical compound OCl=O QBWCMBCROVPCKQ-UHFFFAOYSA-N 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 206010000117 Abnormal behaviour Diseases 0.000 description 1
- 238000012952 Resampling Methods 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 239000011324 bead Substances 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000013399 early diagnosis Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 210000001624 hip Anatomy 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000000474 nursing effect Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000008447 perception Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 210000002832 shoulder Anatomy 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V40/00—Recognition of biometric, human-related or animal-related patterns in image or video data
- G06V40/20—Movements or behaviour, e.g. gesture recognition
- G06V40/23—Recognition of whole body movements, e.g. for sport training
- G06V40/25—Recognition of walking or running movements, e.g. gait recognition
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
- G06F18/232—Non-hierarchical techniques
- G06F18/2321—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
- G06F18/23213—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions with fixed number of clusters, e.g. K-means clustering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
- G06F18/2413—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on distances to training or reference patterns
- G06F18/24147—Distances to closest patterns, e.g. nearest neighbour classification
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/243—Classification techniques relating to the number of classes
- G06F18/24323—Tree-organised classifiers
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N20/00—Machine learning
- G06N20/20—Ensemble learning
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/40—Extraction of image or video features
- G06V10/44—Local feature extraction by analysis of parts of the pattern, e.g. by detecting edges, contours, loops, corners, strokes or intersections; Connectivity analysis, e.g. of connected components
- G06V10/443—Local feature extraction by analysis of parts of the pattern, e.g. by detecting edges, contours, loops, corners, strokes or intersections; Connectivity analysis, e.g. of connected components by matching or filtering
- G06V10/449—Biologically inspired filters, e.g. difference of Gaussians [DoG] or Gabor filters
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Software Systems (AREA)
- Multimedia (AREA)
- Psychiatry (AREA)
- Molecular Biology (AREA)
- Biomedical Technology (AREA)
- Probability & Statistics with Applications (AREA)
- Social Psychology (AREA)
- Human Computer Interaction (AREA)
- Medical Informatics (AREA)
- Biodiversity & Conservation Biology (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Abstract
本发明公开了一种基于直接特征提取的多通道时序步态分析算法,涉及步态分析方法领域,该方法包括以下步骤:拍摄若干受试者的步态视频信息,将所有视频信息分割成单个分离的步态周期,调整为同一长度;计算标准化后的步态参数Sz的距离,使用k‑nn算法和参数距离方程对已经标记好的训练数据做聚类;对步态相位进行分类并重建步态相位,选择当前时间序列中相关系数最大的点作为步态事件,并以此重新分割步态,再次迭代上面的操作:以上一次迭代计算得到的步态事件分割而成的步态分类,计算相邻步态间相关系数,找出全局最大相关系数的点作为步态事件优化上一次迭代的步态分类,直到分类收敛。本发明能够有效提高提取效率和计算效率。
Description
技术领域
本发明涉及步态分析方法领域,具体涉及一种基于直接特征提取的多通道时序步态分析算法。
背景技术
随着人口老龄化的日趋严重,老年人常见的疾病如帕金森综合征患者的数量急剧增加,该疾病严重时会导致病人完全丧失运动能力,给患者和家人带来极大的经济和护理压力,对帕金森综合征早期诊断和治疗能够有效缓解病情的恶化。
步态作为人体固有的生理特征,可在智能假肢、智能监控、临床医学、康复治疗、运动分析等众多领域发挥其作用。通过步态检测与识别,可以对运动人体的行为进行分析,进而实现对特殊人群的异常行为或状态进行跟踪检测,判断患者患病程度或康复程度。
在进行步态分析时,需要采集大量的视屏数据,所采集的数据具有数据量大、范围大且个体差异大的特点,在采集数据后需要进行详细的分析才能得到较准确的数据,难以进行实时测评,且由于个体差异性大,步态分析难以到统一的识别。
早期的步态分析装置主要基于光学的方法,通常在人体上粘贴大量的主动光源(LED)或反射红外的小球作为标记,采用多相机的人体动作捕捉***,如VICON(维康动作捕捉***)、Optitrack(一种全身动作捕捉)等。该方法对步态的测量准确,但成本高且需要专业的测量场地;还有一些***使用惯性传感器和测力板配合的方法,该方法需要准确的传感器进行校准和同步,对传感器要求较高,成本较高。
目前出现一类基于视觉方式的方法,该方法通过对受试者的视频进行轮廓或骨架运动进行提取,再使用移动平台进行计算得到相应的结果,这种方法对不仅对设备提取运动要求较高,且对移动平台的计算能力要求较高,随着科技的发展,现有的便携式深度相机能够有效估计人体姿态,进而极大的减少了图像处理的计算,降低了对移动平台计算能力的要求,使得基于视觉方式的步态分析装置得到有效发展,但是,现有的装置提取效率和计算效率均较低。
发明内容
针对现有技术中存在的缺陷,本发明的目的在于提供一种基于直接特征提取的多通道时序步态分析算法,能够有效提高提取效率和计算效率。
为达到以上目的,本发明采取的技术方案是:
一种基于直接特征提取的多通道时序步态分析算法,包括以下步骤:
S1、拍摄若干受试者的步态视频信息,将所有视频信息分割成单个分离的步态周期,对分割后的每个步态周期进行插值和重新采样,将不同运动速度/采集帧率的步态周期调整为同一长度;
S2、定义标准化后的步态参数为Sz,Sz的距离方程为:
所述公式1中,Pi,x,Pi,y均为1到K的整数,K取1到9中的整数,Pi,x,Pi,y为相对应的步态相位,为在第i个周期内第x帧对应的标准步态参数Sz值,为在第j个周期内第y帧对应的标准步态参数Sz值,φ为标签距离;
标签距离φ(a,b)定义为:
其中,bool(x)为二值函数,如果x为非0则其输出为1;
S3、使用k-nn(k最近邻)的算法和参数距离方程对已经标记好的训练数据做聚类,对任意一个聚类的距离方程为:
其中,M为聚类的大小;
S4、在x帧周围选取任意一个小窗(u,v)的时间序列,其中u,v均∈[-0.5,0.5),根据公式4计算该时间序列的特征值:
其中,L为标准化后的步态周期长度,u,v均为{-0.5,-0.49,-0.48,……,0.49}中的值;
为避免分类器过拟合训练数据,通过基尼杂质Q(Cn)来度量本实施例中时间序列的质量,
u,v均为{-0.5,-0.49,-0.48,……,0.49}中的值,选取Cn值最高的B个(u,v)时间序列作为训练数据;
采用增强随机森林作为训练的模型,在训练过程中我们将目标帧x的标签与这些时间序列映射起来,Cn值最高的B个(u,v)时间序列作为随机森林训练的输入,而(u,v)所在帧x的步态分类标签作为训练的期待输出;
S5、步态相位重建
增强随机森林的输出为当前帧x分类的概率向量ρ={ρ1,……,ρ9},其中每一个概率向量的值表示对应的步态相位的概率,为了对任一时间片段的特征值fs,……,fe进行步态输出,通过公式6来检测当前时间序列中的步态相位:
公式6中,a、b为两个相邻步态相位,ηf(a、b)为当前时间片段在相邻步态相位a、b之间的相关系数;
选择当前时间序列中相关系数最大的点作为步态事件,并以此重新分割步态,再次迭代上面的操作:以上一次迭代计算得到的步态事件分割而成的步态分类,计算相邻步态间相关系数,找出全局最大相关系数的点作为步态事件优化上一次迭代的步态分类,直到分类收敛。
进一步的,步骤S1中所述将不同运动速度/采集帧率的步态周期调整为同一长度具体包括以下步骤:使用2阶贝塞尔曲线回归出连续的运动曲线,并均匀地在每个步态周期样本中采样N个采样点。
进一步的,步骤S4和S5之间还包括以下步骤:对输入序列进行嵌入编码,将训练数据编码后在输入到随机森林,编码方式按公式5进行计算:
γ(f)=(sin(20πf),cos(20πf),...,sin(2L-1πf),cos(2L-1πf)) 公式5
其中,f=fx(u,v),L为4。
进一步的,每个所述步态周期包括9种步态相位:第一加载响应相位P1、第二加载响应相位P2、中间步态相位P3、第一终结步态相位P4、第二终结步态相位P5、摇摆准备相位P6、摇摆开始相位P7、摇摆中间相位P8和摇摆结束相位P9。
进一步的,每一帧的步态参数为z,z∈(λ1,……,λ12),z通过表格1中定义的公式计算得到,定义为在第j个周期内第x帧对应的步态参数z值,其中,λ1至λ12的定义如下:足间距为:λ1=d2/d1,用于观测脚的切换状态;膝盖角率为:λ2=a1/π,用于观测肢体支撑状态;大腿平面角率为:λ3=a2/π,用于观测股骨摇摆状态;右侧脚趾抬高率为:λ4=d3/d1,用于观测右侧脚趾接触状态;右侧脚跟抬高率为:λ5=d4/d1,用于观测右侧脚跟接触状态;左侧脚趾抬高率为:λ6=d5/d1,用于观测左侧脚趾接触状态;左侧脚跟抬高率为:λ7=d6/d1,用于观测左侧脚跟接触状态;腿平面脚率为:λ8=a3/π,用于观测脚的摇摆状态;脚踝角率为:λ9=a4/π,用于观测肢体支撑状态;小腿平面角率为:λ10=a5/π,用于观测胫骨摇摆状态;右脚角率为:λ11=a6/π,用于观测右脚支撑状态;左脚角率为:λ12=a7/π,用于观测左脚支撑状态;
d1为臀到地面的距离,d2为左右脚之间的距离,d3为右侧脚趾到地面的距离,d4为右侧脚跟到地面的距离,d5为左侧脚趾到地面的距离,d6为左侧脚跟到地面的距离;a1至a7为关节或肢体之间的夹角。
与现有技术相比,本发明的优点在于:
(1)本发明中基于直接特征提取的多通道时序步态分析算法,先将所有周期均调整为相同的帧率和速度,使其标准化,在每个步态周期样本中采样N个采样点,然后使用k-nn算法和参数距离方程对已经标记好的训练数据做聚类,将聚类进行平均,由此可以得到聚类平均值,使用随机森林作为训练模型对步态相位的时间序列进行分类,最后根据训练结果重建步态相位,本发明使用k-nn算法和随机森林模型进行计算,计算量较小,效率较高,对平台计算能力要求较低,能够适用于绝大多数移动平台,用于各类感知方式,动作捕捉、视频提取或传感器等***,使用范围较广。
附图说明
图1为本发明实施例中步态周期的结构示意图;
图2为发明实施例中的姿态参数示意图。
具体实施方式
以下结合附图对本发明的实施例作进一步详细说明。
参见图1所示,本发明实施例提供一种基于直接特征提取的多通道时序步态分析算法,该算法基于以下情况:
虽然不同个体的步态存在较大的差异,但是人类标准的步行姿态周期(步态周期)基本分为8到9种步态,本发明实施例中基于各步态中分布的均匀性,选用9步态分割法。
参见图1所示,将步态分为以下9种相位:加载响应1(P1)、加载响应2(P2)、中间步态(P3)、终结步态1(P4)、终结步态2(P5)、摇摆准备(P6)、摇摆开始(P7)、摇摆中间(P8)和摇摆结束(P9),通过分析步态中关节的轨迹可知,每个步态均会以某个关节如脚跟或脚趾接触地面或离开地面为当前步态的起始标志,本实施例中将起始标志称为步态事件。
参见图1中人体骨架示意图,自上而下的四个点分别为:肩、臀、膝盖、脚踝;以脚跟接触地面的骨架图为例,下方左侧为左脚趾和脚跟,右侧为右脚趾和脚跟,不同相位对应的膝盖角比例不相同。
参见图2所示,d1为臀到地面的距离,d2为左右脚之间的距离,d3为右侧脚趾到地面的距离,d4为右侧脚跟到地面的距离,d5为左侧脚趾到地面的距离,d6为左侧脚跟到地面的距离;a1至a7为关节或肢体之间的夹角。
定义足间距为:λ1=d2/d1,用于观测脚的切换状态;
膝盖角率为:λ2=a1/π,用于观测肢体支撑状态;
大腿平面角率为:λ3=a2/π,用于观测股骨摇摆状态;
右侧脚趾抬高率为:λ4=d3/d1,用于观测右侧脚趾接触状态;
右侧脚跟抬高率为:λ5=d4/d1,用于观测右侧脚跟接触状态;
左侧脚趾抬高率为:λ6=d5/d1,用于观测左侧脚趾接触状态;
左侧脚跟抬高率为:λ7=d6/d1,用于观测左侧脚跟接触状态;
腿平面脚率为:λ8=a3/π,用于观测脚的摇摆状态;
脚踝角率为:λ9=a4/π,用于观测肢体支撑状态;
小腿平面角率为:λ10=a5/π,用于观测胫骨摇摆状态;
右脚角率为:λ11=a6/π,用于观测右脚支撑状态;
左脚角率为:λ12=a7/π,用于观测左脚支撑状态。
P1至P9组成了一个完整的步态周期,参见图1可知每一个步态周期对应不同的轨迹曲线,图1中以膝盖角比例(即膝盖角率λ2)为例,示意出步态参数随着步态周期的推进而变化的曲线,同一步态周期内不同时间(帧数)所对应的膝盖角比例不相同,且构成一个完整连续的曲线,其他步态参数随步态周期也呈现出连续完整的曲线。
基于此,本发明以步态参数来对不同时间的步态特征进行表征,每一帧的步态参数为z,z∈(λ1,……,λ12),z通过λ1至λ12相应的公式计算得到,定义为在第j个周期(j∈(P1,……,P9)内第x帧对应的步态参数z值。
由于步态参数收到年龄、健康程度、身体情况等因素的影响,不同个体的步态参数差异较大,本发明通过以下方式提取步态模式:
S1、步态周期标准化
拍摄若干受试者的步态视频信息,将所有视频信息分割成单个分离的步态周期(包括P1到P9的周期),由于不同受试者的速度不同,且同一受试者步同时间的速度也不相同,需要对分割后的每个步态周期进行插值和重新采样,将不同运动速度或采集帧率的步态数据调整为同一长度(即每个周期都调整为相同的帧率和速度)。
具体的:使用2阶贝塞尔曲线回归出连续的运动曲线,并均匀地在每个步态周期样本中采样N个采样点,这样所有步态周期就被统一在相同的数据尺度下,步态数据被调整后,每一帧的步态参数z也随之变化,以Sz表示标准化后的步态参数,其计算方式与z相同。
S2、定义标准化后的参数距离方程
对标准化后的步态参数Sz进行距离定义
公式1中,Pi,x,Pi,y∈(1,……,K),为相对应的步态相位的标签(Pi,x,Pi,y即P1,……,P9,K取1到9中的整数),为在第i个周期内第x帧对应的标准后的步态参数Sz值,为在第j个周期内第y帧对应的标准化后的步态参数Sz值,φ为标签距离。
标签距离φ(a,b)定义为:
其中,bool(x)为二值函数,如果x为非0则其输出为1。
S3、对已经标准化的周期数据做聚类
使用k-nn(k最近邻)的算法和参数距离方程对已经标记好的训练数据做聚类,对任意一个聚类的距离方程为:
其中,M为聚类的大小,M的具体数值根据需要进行选择,为大于等于1的整数。
S4、对步态相位的时间序列进行分类
由于步态是一个时序概念,对单独的一帧进行分类往往会导致输出的观测结果有非常大的扰动或噪声,所以本发明实施例使用时间序列进行分类。
在x帧周围选取任意一个小窗(u,v)的时间序列,其中u,v均∈[-0.5,0.5),根据公式4定义该时间序列的特征值为:
其中,L为标准化后的步态周期长度,本实施例中,L=100,u,v均为{-0.5,-0.49,-0.48,……,0.49}中的值,通过公式4能够消除时间序列长度对特征的影响。
为避免分类器过拟合训练数据,通过基尼杂质Q(Cn)来度量本实施例中时间序列的质量,
u,v均为{-0.5,-0.49,-0.48,……,0.49}中的值
本实施例选取Q(Cn)最高的B个(u,v)时间序列作为训练数据,B的具体值根据实际需要确定。
采用增强随机森林作为训练的模型,在训练过程中将目标帧x的标签与这些时间序列映射起来,时间序列(u,v)作为随机森林训练的输入,而(u,v)所在帧x的步态分类标签作为训练的期待输出。
此外,为了增加随机森林的泛化能力,对输入序列做了嵌入编码,将训练数据编码后在输入到随机森林,编码方式按公式5进行计算:
γ(f)=(sin(20πf),cos(20πf),...,sin(2L-1πf),cos(2L-1πf)) 公式5
其中,f=fx(u,v),L为4
随机森林的计算量小,适用于绝大多数移动平台,可以快速的训练与回归,为算法的便携性和可移植性提供了保障。
S5、步态相位重建
增强随机森林的输出为当前帧x分类的概率向量ρ={ρ1,……,ρ9},其中每一个概率向量的值表示对应的步态相位的概率,为了对任一时间片段的特征值fs,……,fe进行步态输出,通过公式6来检测当前时间序列中的步态相位:
公式6中,a、b为两个相邻步态相位,ηf(a、b)为当前时间片段在相邻步态相位a、b之间的相关系数。
选择当前时间序列中相关系数最大的点作为步态事件,并以此重新分割步态,再次迭代上面的操作:以上一次迭代计算得到的步态事件分割而成的步态分类,计算相邻步态间相关系数,找出全局最大相关系数的点作为步态事件优化上一次迭代的步态分类,直到分类收敛。
本发明不仅局限于上述最佳实施方式,任何人在本发明的启示下都可得出其他各种形式的产品,但不论在其形状或结构上作任何变化,凡是具有与本发明相同或相近似的技术方案,均在其保护范围之内。
Claims (5)
1.一种基于直接特征提取的多通道时序步态分析算法,其特征在于,包括以下步骤:
S1、拍摄若干受试者的步态视频信息,将所有视频信息分割成单个分离的步态周期,对分割后的每个步态周期进行插值和重新采样,将不同运动速度/采集帧率的步态周期调整为同一长度;
S2、定义标准化后的步态参数为Sz,Sz的距离方程为:
所述公式1中,Pi,x,Pi,y均为1到K的整数,K取1到9中的整数,Pi,x,Pi,y为相对应的步态相位,为在第i个周期内第x帧对应的标准步态参数Sz值,为在第j个周期内第y帧对应的标准步态参数Sz值,φ为标签距离;
标签距离φ(a,b)定义为:
其中,bool(x)为二值函数,如果x为非0则其输出为1;
S3、使用k-nn(k最近邻)的算法和参数距离方程对已经标记好的训练数据做聚类,对任意一个聚类的距离方程为:
其中,M为聚类的大小;
S4、在x帧周围选取任意一个小窗(u,v)的时间序列,其中u,v均∈[-0.5,0.5),根据公式4计算该时间序列的特征值:
其中,L为标准化后的步态周期长度,u,v均为{-0.5,-0.49,-0.48,……,0.49}中的值;
为避免分类器过拟合训练数据,通过基尼杂质Q(Cn)来度量本实施例中时间序列的质量,
u,v均为{-0.5,-0.49,-0.48,……,0.49}中的值,选取Cn值最高的B个(u,v)时间序列作为训练数据;
采用增强随机森林作为训练的模型,在训练过程中我们将目标帧x的标签与这些时间序列映射起来,Cn值最高的B个(u,v)时间序列作为随机森林训练的输入,而(u,v)所在帧x的步态分类标签作为训练的期待输出;
S5、步态相位重建
增强随机森林的输出为当前帧x分类的概率向量ρ={ρ1,……,ρ9},其中每一个概率向量的值表示对应的步态相位的概率,为了对任一时间片段的特征值fs,……,fe进行步态输出,通过公式6来检测当前时间序列中的步态相位:
公式6中,a、b为两个相邻步态相位,ηf(a、b)为当前时间片段在相邻步态相位a、b之间的相关系数;
选择当前时间序列中相关系数最大的点作为步态事件,并以此重新分割步态,再次迭代上面的操作:以上一次迭代计算得到的步态事件分割而成的步态分类,计算相邻步态间相关系数,找出全局最大相关系数的点作为步态事件优化上一次迭代的步态分类,直到分类收敛。
2.如权利要求1所述的一种基于直接特征提取的多通道时序步态分析算法,其特征在于:步骤S1中所述将不同运动速度/采集帧率的步态周期调整为同一长度具体包括以下步骤:使用2阶贝塞尔曲线回归出连续的运动曲线,并均匀地在每个步态周期样本中采样N个采样点。
3.如权利要求1所述的一种基于直接特征提取的多通道时序步态分析算法,其特征在于:步骤S4和S5之间还包括以下步骤:对输入序列进行嵌入编码,将训练数据编码后在输入到随机森林,编码方式按公式5进行计算:
γ(f)=(sin(20πf),cos(20πf),...,sin(2L-1πf),cos(2L-1πf)) 公式5
其中,f=fx(u,v),L为4。
4.如权利要求1所述的一种基于直接特征提取的多通道时序步态分析算法,其特征在于:每个所述步态周期包括9种步态相位:第一加载响应相位P1、第二加载响应相位P2、中间步态相位P3、第一终结步态相位P4、第二终结步态相位P5、摇摆准备相位P6、摇摆开始相位P7、摇摆中间相位P8和摇摆结束相位P9。
5.如权利要求1所述的一种基于直接特征提取的多通道时序步态分析算法,其特征在于:
每一帧的步态参数为z,z∈(λ1,……,λ12),z通过表格1中定义的公式计算得到,定义(x)为在第j个周期内第x帧对应的步态参数z值,其中,λ1至λ12的定义如下:足间距为:λ1=d2/d1,用于观测脚的切换状态;膝盖角率为:λ2=a1/π,用于观测肢体支撑状态;大腿平面角率为:λ3=a2/π,用于观测股骨摇摆状态;右侧脚趾抬高率为:λ4=d3/d1,用于观测右侧脚趾接触状态;右侧脚跟抬高率为:λ5=d4/d1,用于观测右侧脚跟接触状态;左侧脚趾抬高率为:λ6=d5/d1,用于观测左侧脚趾接触状态;左侧脚跟抬高率为:λ7=d6/d1,用于观测左侧脚跟接触状态;腿平面脚率为:λ8=a3/π,用于观测脚的摇摆状态;脚踝角率为:λ9=a4/π,用于观测肢体支撑状态;小腿平面角率为:λ10=a5/π,用于观测胫骨摇摆状态;右脚角率为:λ11=a6/π,用于观测右脚支撑状态;左脚角率为:λ12=a7/π,用于观测左脚支撑状态;
d1为臀到地面的距离,d2为左右脚之间的距离,d3为右侧脚趾到地面的距离,d4为右侧脚跟到地面的距离,d5为左侧脚趾到地面的距离,d6为左侧脚跟到地面的距离;a1至a7为关节或肢体之间的夹角。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011083624.6A CN112115923A (zh) | 2020-10-12 | 2020-10-12 | 一种基于直接特征提取的多通道时序步态分析算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011083624.6A CN112115923A (zh) | 2020-10-12 | 2020-10-12 | 一种基于直接特征提取的多通道时序步态分析算法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112115923A true CN112115923A (zh) | 2020-12-22 |
Family
ID=73798082
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011083624.6A Pending CN112115923A (zh) | 2020-10-12 | 2020-10-12 | 一种基于直接特征提取的多通道时序步态分析算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112115923A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112800908A (zh) * | 2021-01-19 | 2021-05-14 | 中国科学院心理研究所 | 一种基于视频中个体步态分析的焦虑感知模型的建立方法 |
CN113504778A (zh) * | 2021-07-26 | 2021-10-15 | 广东工业大学 | 基于融合概率模型的足式机器人控制方法、***及设备 |
CN114663913A (zh) * | 2022-02-28 | 2022-06-24 | 电子科技大学 | 一种基于Kinect的人体步态参数提取方法 |
CN115211846A (zh) * | 2022-03-14 | 2022-10-21 | 北京中科睿医信息科技有限公司 | 一种基于深度传感器的重复座椅起立测试定量评估方法和装置 |
-
2020
- 2020-10-12 CN CN202011083624.6A patent/CN112115923A/zh active Pending
Non-Patent Citations (1)
Title |
---|
MINXIANG YE 等: "GAIT PHASE CLASSIFICATION FOR IN-HOME GAIT ASSESSMENT", 《IEEE》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112800908A (zh) * | 2021-01-19 | 2021-05-14 | 中国科学院心理研究所 | 一种基于视频中个体步态分析的焦虑感知模型的建立方法 |
CN112800908B (zh) * | 2021-01-19 | 2024-03-26 | 中国科学院心理研究所 | 一种基于视频中个体步态分析的焦虑感知模型的建立方法 |
CN113504778A (zh) * | 2021-07-26 | 2021-10-15 | 广东工业大学 | 基于融合概率模型的足式机器人控制方法、***及设备 |
CN113504778B (zh) * | 2021-07-26 | 2023-09-19 | 广东工业大学 | 基于融合概率模型的足式机器人控制方法、***及设备 |
CN114663913A (zh) * | 2022-02-28 | 2022-06-24 | 电子科技大学 | 一种基于Kinect的人体步态参数提取方法 |
CN114663913B (zh) * | 2022-02-28 | 2023-10-31 | 电子科技大学 | 一种基于Kinect的人体步态参数提取方法 |
CN115211846A (zh) * | 2022-03-14 | 2022-10-21 | 北京中科睿医信息科技有限公司 | 一种基于深度传感器的重复座椅起立测试定量评估方法和装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP3468450B1 (en) | Method and system for analyzing human gait | |
CN112115923A (zh) | 一种基于直接特征提取的多通道时序步态分析算法 | |
Gu et al. | Markerless gait analysis based on a single RGB camera | |
Rast et al. | Systematic review on the application of wearable inertial sensors to quantify everyday life motor activity in people with mobility impairments | |
Leightley et al. | Automated analysis and quantification of human mobility using a depth sensor | |
CN109815858B (zh) | 一种日常环境中的目标用户步态识别***及方法 | |
Jalloul et al. | Activity recognition using complex network analysis | |
Hossain et al. | Deepbbwae-net: A cnn-rnn based deep superlearner for estimating lower extremity sagittal plane joint kinematics using shoe-mounted imu sensors in daily living | |
Auvinet et al. | Lower limb movement asymmetry measurement with a depth camera | |
WO2014153665A1 (en) | System and method for monitoring a subject | |
US20230355135A1 (en) | Intelligent gait analyzing apparatus | |
CN114052718A (zh) | 步态评估***及步态评估方法 | |
Wei et al. | Using sensors and deep learning to enable on-demand balance evaluation for effective physical therapy | |
Perez et al. | A smartphone-based system for clinical gait assessment | |
CN117883074A (zh) | 一种基于人体姿态视频的帕金森病步态量化分析方法 | |
Ali et al. | Detection and analysis of transitional activity in manifold space | |
Li et al. | Quantitative analysis of abnormal and normal gait based on inertial sensors | |
CN116999057A (zh) | 基于穿戴式传感器的偏瘫步态识别及偏瘫步态评估方法 | |
KR102194313B1 (ko) | 검출정보별 신경망 분석을 실시하여 개인을 식별하는 장치 및 방법 | |
Natarajan et al. | Reha@ home-a vision based markerless gait analysis system for rehabilitation at home | |
Dong et al. | An “optical flow” method based on pressure sensors data for quantification of Parkinson's disease characteristics | |
Low et al. | Lower extremity kinematics walking speed classification using long short-term memory neural frameworks | |
KR20210011097A (ko) | 검출정보별 신경망 분석을 실시하여 보행 종류를 분류하는 장치 및 방법 | |
Khokhlova et al. | Kinematic covariance based abnormal gait detection | |
Sant'Anna et al. | Developing a motion language: Gait analysis from accelerometer sensor systems |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20201222 |
|
RJ01 | Rejection of invention patent application after publication |