CN107609523A - 基于Python的步态周期及三维肢体活动角度算法 - Google Patents
基于Python的步态周期及三维肢体活动角度算法 Download PDFInfo
- Publication number
- CN107609523A CN107609523A CN201710848939.7A CN201710848939A CN107609523A CN 107609523 A CN107609523 A CN 107609523A CN 201710848939 A CN201710848939 A CN 201710848939A CN 107609523 A CN107609523 A CN 107609523A
- Authority
- CN
- China
- Prior art keywords
- joint
- gait cycle
- gait
- angle
- python
- 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
Landscapes
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Abstract
本发明提供了一种基于Python的步态周期及肢体活动角度算法,步骤为:对接三维动作捕捉系统,根据人体重要关节所在的身体部位及关节名称,对原数据进行格式化处理;根据步态周期特征设定跃升阈值,基于髋关节矢状面的角度变化规律确定数据库中所有步态周期;基于数据记录频率及步态周期的起始和结束帧确定步长时间,利用Python在每个步态周期内查找三维关节角的最大、最小值,并计算三维肢体活动角度;根据关节位置生成两个数据矩阵,分别写入每个步态周期的步态和肢体活动指标。本发明可实现大量实时关节角度数据的快速有效处理,适用于各种三维动作捕捉系统,以及人体本身不同负重及服装装备条件,通用性强,灵活性高,扩展性好。
Description
技术领域
本发明涉及一种三维肢体活动角度算法,特别是涉及一种基于Python的步态周期及三维肢体活动角度算法。
背景技术
肢体活动角度(ROM,Range of Motion)是指人体在运动时的肢体活动范围,可用于表征肢体的灵活程度。在三维动作分析中,一般将人体分为3个截面,分别为冠状面、矢状面和水平面。人体在运动时,各个关节在三个截面会分别发生内收/外展、屈曲/伸展、内旋转/外旋转等。
角度计是测量ROM的常用技术手段,但不能实时监测关节在三个截面的活动情况。三维动作捕捉仪的出现实现了对人体关节数据的实时获取,其应用范围非常广泛,包括生物力学、运动学、工效学研究、虚拟现实、电影制作等领域。在生物力学和运动学领域,三维动作捕捉仪主要用于不同负重或身体状态下的行走模式、动作分析及身体平衡的判定,及不同身体截面方向特定关节的运动学特征分析。近年来,三维动作捕捉仪在特种服装装备的工效学测试和性能优化领域逐渐得到了应用。
然而,由于三维动作捕捉仪会实时记录人体各个关节角度在三个维度的变化,监测过程中获得的数据量巨大,因此,需要依靠有效的技术手段实现对关节角度数据的提取和分析。现有技术中还没有这样的有效的技术手段。
发明内容
本发明要解决的技术问题是提供一种三维肢体活动角度算法,与三维动作捕捉仪相结合,用于三维动作捕捉仪获取的大量数据的后处理,计算周期性关节角度,评价肢体灵活性,减少外界因素对人体的影响。
为了解决上述技术问题,本发明的技术方案是提供一种基于Python的步态周期及肢体活动角度算法,其特征在于,步骤为:
步骤1:对接三维动作捕捉系统,根据人体重要关节所在的身体部位及关节名称,对原数据进行格式化处理;
步骤2:根据步态周期特征设定跃升阈值,基于髋关节矢状面的角度变化规律确定数据库中所有步态周期;
步骤3:基于数据记录频率及步态周期的起始和结束帧确定步长时间,利用Python在每个步态周期内查找三维关节角的最大、最小值,并计算三维肢体活动角度;
步骤4:根据关节位置生成两个数据矩阵,分别写入每个步态周期的步态和肢体活动指标。
优选地,所述步骤1包括以下子步骤:
步骤1.1:提取三维动作捕捉仪获取的三维关节角度数据;
步骤1.2:根据关节所处位置,将其划分为左侧、右侧区域;
步骤1.3:纵向以时间顺序、横向以关节名称顺序,形成不同时刻人体重要关节的三维关节角度数据库。
优选地,所述步骤1中,人体重要关节包括肩关节、肘关节、腕关节、髋关节、膝关节、踝关节。
优选地,所述步骤1.3中,关节名称依次排列冠状面、矢状面和水平面角度。
优选地,所述步骤2包括以下子步骤:
步骤2.1:对数据进行步态周期的过滤,删除无效数据;
步骤2.2:以右侧髋关节矢状面的角度变化规律提取所有步态周期。
更优选地,所述步骤2.1中,设定右侧髋关节矢状面角度跃升阈值,当在一个时间区间内角度跃升大于阈值时,说明人体开始行走,删除之前的无效数据。
更优选地,所述步骤2.2中,判断1个步态周期的方法是:当右脚后跟着地,髋关节达到最大屈曲时,作为1个步态周期的起始时刻;当髋关节角度达到最大伸展时,右脚脚尖即将离地;当髋关节角度第二次达到最大屈曲时,右脚脚跟再次着地,第1个步态周期结束,同时第2个周期开始。
优选地,所述步骤3包括以下子步骤:
步骤3.1:根据步态周期的初始帧和结束帧,以及三维动作捕捉仪的数据记录频率,计算每个步长时间;
步骤3.2:在各步态周期内,利用Python查找各个关节在三个维度方向的最大值和最小值,并计算差值;由此,获得所需的步态和肢体活动角度数据。
优选地,所述步骤4具体为:
根据步骤3计算的指标生成步态报告,步态报告写入左、右两个矩阵;输出步态报告中包括各周期的初始帧、结束帧、步长时间在内的步态数据,以及三维关节角度的包括最大值、最小值和差值在内的肢体活动角度数据。
优选地,所述左、右两个矩阵分别表示身体左侧关节和身体右侧关节
本发明提供的方法可实现大量实时关节角度数据的快速有效处理,适用于各种三维动作捕捉系统,以及人体本身的不同负重及服装装备条件,通用性强。可针对不同的人体动作和测试任务,设定不同的数据筛选阈值,灵活性强。可将数据处理结果生成系统性测试报告,有利于进行结果分析。具有一定的扩展性,可在算法基础上增加加速度运算、距离运算等后处理模块,并生成结果图表,从而实现更强的可视化效果。采用本发明提供的步态周期划分方法,可实现对不同外界因素,如负重条件和服装装备等因素对人体肢体活动影响的评价,而该算法的可扩展性又为人体工效学领域的数据分析及符合工效学的产品设计提供了基础。
附图说明
图1为基于Python的步态周期及肢体活动角度算法模块化流程图;
图2为髋关节三维角度示意图;(a)冠状面,其中,1表示内收、关节角度>0°,2表示外展、关节角度<0°;(b)矢状面,其中,3表示屈曲、关节角度>0°,4表示伸展、关节角度<0°;(c)水平面,其中,5表示内旋转、关节角度>0°,6表示外旋转、关节角度<0°;
图3为以右脚为例的一个步态周期示意图;(a)右脚跟着地;(b)右脚掌着地;(c)右脚支撑;(d)左脚掌着地;(e)左脚支撑;(f)右脚跟着地;
图4为右侧髋关节矢状面角度的变化情况示意图。
具体实施方式
下面结合具体实施例,进一步阐述本发明。
在人体正常行走时,人体的运动特征会随步态周期性变化,而该周期一般通过单脚与地面的接触状态判定。在步态周期性变化的同时,人体的各个关节角度均会出现周期性变化。本发明以步态周期为基础,将大量的关节角度数据进行周期性划分,并提取各个周期的关键参数,可有效实现对肢体活动角度的分析。该方法可用于评价外界因素,如服装装备、负重状态对人体活动灵活性的影响,并预测可能对人体造成的肢体损伤。
本发明提供了一种基于Python的步态周期及三维肢体活动角度算法,包括以下步骤:对接三维动作捕捉系统,根据人体重要关节所在的身体部位及关节名称,对原数据进行格式化处理;根据步态周期特征设定跃升阈值,基于髋关节矢状面的角度变化规律确定数据库中所有步态周期;基于数据记录频率及步态周期的起始和结束帧确定步长时间,利用Python在每个步态周期内查找三维关节角的最大最小值,并计算三维肢体活动角度;根据关节位置生成两个数据矩阵,分别写入每个步态周期的步态和肢体活动指标。
图1所示为基于Python的步态周期及肢体活动角度算法模块化流程图,包括原数据处理模块、步态周期过滤模块、步态数据处理模块和步态报告生成模块。具体步骤如下:
步骤1:通过原数据处理模块进行数据标准化处理
原数据处理模块可与三维动作捕捉系统对接,根据人体重要关节所在的身体部位以及关节名称,对原数据进行格式化处理,从而建立不同时刻所有相关部位(肩关节、肘关节、腕关节、髋关节、膝关节、踝关节等)的三维关节角度数据库。
本算法主要用于对三维动作捕捉仪获取的肢体活动角度数据进行后处理,对于不同的三维动作捕捉系统,数据的存储形式不同,因此,首先需要对数据进行格式化处理。人体关节,除头部和躯干部位以外,主要可以分为左右两侧,且一一对应。每个关节处的角度可分为三个维度。图2所示为髋关节三维角度示意图,在冠状面会发生内收或外展、矢状面发生屈曲或伸展,水平面发生内旋转或外旋转。对于所有关节点处的数据,排列方式为纵向以时间顺序排列,横向以关节名称排列,且每个关节处依次排列冠状面、矢状面和水平面角度。
步骤2:对步骤1处理完成的数据矩阵进行步态周期过滤
步态周期过滤模块是对三维肢体角度数据进行后处理的关键。在利用三维动作捕捉仪记录人体行走动作时,由于记录时间的提前或延迟,会出现不满足步态周期形成的数据,因此,需要筛除不能表征标准动作的无效数据。在筛除数据时,需要根据步态周期的特征,设定跃升阈值。标准动作以髋关节矢状面的角度变化为筛选基础,根据角度变化规律,确定处理数据库中的所有步态周期。
图3所示为以右脚为例的一个步态周期。从右脚跟着地开始,至右脚掌着地,然后右脚支撑,至左脚掌着地同时右脚腾空,以右脚跟第二次着地为止。在正常行走时,人体全身尤其是四肢部位的关节角度均会发生周期性变化。在利用三维动作捕捉仪进行数据记录时,由于动作的延迟等因素,会出现不符合实际行走状态的无效数据,因此,需要对数据进行步态周期的过滤。本发明以髋关节矢状面的角度作为判断步态周期的代表性部位。如图4所示为右侧髋关节矢状面角度的变化情况,阶段A代表人体仍在站立,未进行行走的状态,为无效数据。当在一个时间区间内髋关节角度的跃升大于阈值△时,说明人体开始行走,直至角度达到最大时,开始进入步态周期。当右脚后跟着地,髋关节达到最大屈曲(Max1)作为1个步态周期的起始时刻;当髋关节角度达到最大伸展(Min1)时,右脚脚尖即将离地;当髋关节角度为Max2时,右脚脚跟再次着地,第1个步态周期t1结束,同时第2个周期t2开始。阶段B中角度开始增大但未达到最大值时,不属于完整的步态周期,因此,也是要筛除的数据。阶段C和阶段D属于完整的步态周期。在具体操作时,根据步态周期的划分原则,利用Python删除无效数据,并确定所有步态周期。同时,在该步态周期内对所有相关关节的角度数据进行分析。
步骤3:在步骤2确定的所有步态周期中进行参数计算
步态数据处理模块需要对每个步态周期内的数据进行提取和计算。根据三维动作捕捉仪的数据记录频率,以及步态周期的起始和结束帧数,可以确定每个步态周期的步长时间。另外,在每个步态周期内,根据数据变化规律,利用Python查找各个关节在三个维度方向的最大值和最小值,并计算差值,即ROM。由此,获得所需的步态和肢体活动角度数据。
步骤4:以步骤3获得的数据为基础生成步态报告
步态报告生成模块可以输出处理的步态周期和肢体活动角度结果。根据关节所处的身体左侧或右侧,生成两个数据矩阵,分别写入每个步态周期的起始帧、结束帧和步长时间,以及每个关节在三个维度方向的最大值、最小值和ROM。
根据步骤3的计算结果生成步态报告,报告写入左(身体左侧关节)、右(身体右侧关节)两个矩阵。输出报告中包括各周期的初始帧、结束帧、步长时间等步态数据,以及三维关节角度的最大值、最小值和差值(ROM)等肢体活动角度数据。为进一步开展各步态周期之间、以及不同外环境条件下步态之间的分析提供系统数据。
以上所述,仅为本发明的较佳实施例,并非对本发明任何形式上和实质上的限制,应当指出,对于本技术领域的普通技术人员,在不脱离本发明方法的前提下,还将可以做出若干改进和补充,这些改进和补充也应视为本发明的保护范围。凡熟悉本专业的技术人员,在不脱离本发明的精神和范围的情况下,当可利用以上所揭示的技术内容而做出的些许更动、修饰与演变的等同变化,均为本发明的等效实施例;同时,凡依据本发明的实质技术对上述实施例所作的任何等同变化的更动、修饰与演变,均仍属于本发明的技术方案的范围内。
Claims (10)
1.一种基于Python的步态周期及肢体活动角度算法,其特征在于,步骤为:
步骤1:对接三维动作捕捉系统,根据人体重要关节所在的身体部位及关节名称,对原数据进行格式化处理;
步骤2:根据步态周期特征设定跃升阈值,基于髋关节矢状面的角度变化规律确定数据库中所有步态周期;
步骤3:基于数据记录频率及步态周期的起始和结束帧确定步长时间,利用Python在每个步态周期内查找三维关节角的最大、最小值,并计算三维肢体活动角度;
步骤4:根据关节位置生成两个数据矩阵,分别写入每个步态周期的步态和肢体活动指标。
2.如权利要求1所述的一种基于Python的步态周期及肢体活动角度算法,其特征在于:所述步骤1包括以下子步骤:
步骤1.1:提取三维动作捕捉仪获取的三维关节角度数据;
步骤1.2:根据关节所处位置,将其划分为左侧、右侧区域;
步骤1.3:纵向以时间顺序、横向以关节名称顺序,形成不同时刻人体重要关节的三维关节角度数据库。
3.如权利要求1或2所述的一种基于Python的步态周期及肢体活动角度算法,其特征在于:所述步骤1中,人体重要关节包括肩关节、肘关节、腕关节、髋关节、膝关节、踝关节。
4.如权利要求1所述的一种基于Python的步态周期及肢体活动角度算法,其特征在于:所述步骤1.3中,关节名称依次排列冠状面、矢状面和水平面角度。
5.如权利要求1所述的一种基于Python的步态周期及肢体活动角度算法,其特征在于:所述步骤2包括以下子步骤:
步骤2.1:对数据进行步态周期的过滤,删除无效数据;
步骤2.2:以右侧髋关节矢状面的角度变化规律提取所有步态周期。
6.如权利要求5所述的一种基于Python的步态周期及肢体活动角度算法,其特征在于:所述步骤2.1中,设定右侧髋关节矢状面角度跃升阈值,当在一个时间区间内角度跃升大于阈值时,说明人体开始行走,删除之前的无效数据。
7.如权利要求5所述的一种基于Python的步态周期及肢体活动角度算法,其特征在于:所述步骤2.2中,判断1个步态周期的方法是:当右脚后跟着地,髋关节达到最大屈曲时,作为1个步态周期的起始时刻;当髋关节角度达到最大伸展时,右脚脚尖即将离地;当髋关节角度第二次达到最大屈曲时,右脚脚跟再次着地,第1个步态周期结束,同时第2个周期开始。
8.如权利要求1所述的一种基于Python的步态周期及肢体活动角度算法,其特征在于:所述步骤3包括以下子步骤:
步骤3.1:根据步态周期的初始帧和结束帧,以及三维动作捕捉仪的数据记录频率,计算每个步长时间;
步骤3.2:在各步态周期内,利用Python查找各个关节在三个维度方向的最大值和最小值,并计算差值;由此,获得所需的步态和肢体活动角度数据。
9.如权利要求1所述的一种基于Python的步态周期及肢体活动角度算法,其特征在于:所述步骤4具体为:
根据步骤3计算的指标生成步态报告,步态报告写入左、右两个矩阵;输出步态报告中包括各周期的初始帧、结束帧、步长时间在内的步态数据,以及三维关节角度的包括最大值、最小值和差值在内的肢体活动角度数据。
10.如权利要求9所述的一种基于Python的步态周期及肢体活动角度算法,其特征在于:所述左、右两个矩阵分别表示身体左侧关节和身体右侧关节。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710848939.7A CN107609523A (zh) | 2017-09-19 | 2017-09-19 | 基于Python的步态周期及三维肢体活动角度算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710848939.7A CN107609523A (zh) | 2017-09-19 | 2017-09-19 | 基于Python的步态周期及三维肢体活动角度算法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107609523A true CN107609523A (zh) | 2018-01-19 |
Family
ID=61060910
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710848939.7A Pending CN107609523A (zh) | 2017-09-19 | 2017-09-19 | 基于Python的步态周期及三维肢体活动角度算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107609523A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111904427A (zh) * | 2019-05-08 | 2020-11-10 | 株式会社比弗雷斯 | 通过头部加速度传感器的步行因素检测方法及装置 |
CN112435731A (zh) * | 2020-12-16 | 2021-03-02 | 成都翡铭科技有限公司 | 一种判断实时姿势是否满足预设规则的方法 |
CN112998700A (zh) * | 2021-05-26 | 2021-06-22 | 北京欧应信息技术有限公司 | 用于辅助对象运动功能评估的设备、系统和方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103514302A (zh) * | 2013-10-28 | 2014-01-15 | 深圳先进技术研究院 | 人体步态数据库及其建立方法 |
US20160370854A1 (en) * | 2015-06-16 | 2016-12-22 | Wilson Steele | Method and System for Analyzing a Movement of a Person |
CN106933100A (zh) * | 2017-03-19 | 2017-07-07 | 北京工业大学 | 一种基于人体运动捕捉数据的仿人机器人步态规划方法 |
-
2017
- 2017-09-19 CN CN201710848939.7A patent/CN107609523A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103514302A (zh) * | 2013-10-28 | 2014-01-15 | 深圳先进技术研究院 | 人体步态数据库及其建立方法 |
US20160370854A1 (en) * | 2015-06-16 | 2016-12-22 | Wilson Steele | Method and System for Analyzing a Movement of a Person |
CN106933100A (zh) * | 2017-03-19 | 2017-07-07 | 北京工业大学 | 一种基于人体运动捕捉数据的仿人机器人步态规划方法 |
Non-Patent Citations (1)
Title |
---|
田苗 等: "三维动作捕捉仪在服装工效学评价中的应用", 《服装学报》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111904427A (zh) * | 2019-05-08 | 2020-11-10 | 株式会社比弗雷斯 | 通过头部加速度传感器的步行因素检测方法及装置 |
CN112435731A (zh) * | 2020-12-16 | 2021-03-02 | 成都翡铭科技有限公司 | 一种判断实时姿势是否满足预设规则的方法 |
CN112435731B (zh) * | 2020-12-16 | 2024-03-19 | 成都翡铭科技有限公司 | 一种判断实时姿势是否满足预设规则的方法 |
CN112998700A (zh) * | 2021-05-26 | 2021-06-22 | 北京欧应信息技术有限公司 | 用于辅助对象运动功能评估的设备、系统和方法 |
CN112998700B (zh) * | 2021-05-26 | 2021-09-24 | 北京欧应信息技术有限公司 | 用于辅助对象运动功能评估的设备、系统和方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Dorschky et al. | CNN-based estimation of sagittal plane walking and running biomechanics from measured and simulated inertial sensor data | |
CN107609523A (zh) | 基于Python的步态周期及三维肢体活动角度算法 | |
CN105184767B (zh) | 一种运动人体姿态相似性度量方法 | |
CN111383221B (zh) | 一种脊柱侧弯检测模型的生成方法和计算机设备 | |
Bishop et al. | The influence of speed and size on avian terrestrial locomotor biomechanics: predicting locomotion in extinct theropod dinosaurs | |
WO2005002440A1 (en) | Kinematic quantification of gait asymmetry based on bilateral cyclograms | |
Khoshnoudian et al. | A new damage index using FRF data, 2D-PCA method and pattern recognition techniques | |
KR101190840B1 (ko) | 다물체 모델링과 통계 방법론을 이용한 인체 하지부 구속력의 불확실성 추정 방법 | |
Luleci et al. | Generative adversarial networks for data generation in structural health monitoring | |
Chang et al. | Landslide features interpreted by neural network method using a high-resolution satellite image and digital topographic data | |
Obukhov et al. | Human motion capture algorithm for creating digital shadows of the movement process | |
Raspudic et al. | Comparison of the convergence of results on the example of a thin plate | |
Ramey | BIOMECHANICS OF LONG JUMP AND TRIPLE JUMP | |
Qian et al. | Influence of vibration on the grain flow sensor during the harvest and the difference elimination method | |
CN113240044A (zh) | 一种基于多Kinect的人体骨骼数据融合评价方法 | |
Willener et al. | Effect of walking speed on the gait of king penguins: An accelerometric approach | |
Jia | Recognition model of sports athletes’ wrong actions based on computer vision | |
CN110755084A (zh) | 基于主被动、分阶段动作的运动功能评估方法及设备 | |
Arif et al. | An analyze of a backpropagation neural network in the identification of critical land based on ALOS imagery | |
Yao et al. | Research on Trajectory and Modeling of Body Movement | |
Maurice et al. | A subject-specific software solution for the modeling and the visualization of muscles deformations | |
Sheng et al. | Research on Dance Movement Recognition and Assessment through Human Pose Estimation | |
Wang et al. | An intelligent data analysis-based medical management method for lower limb health of football athletes | |
Filippeschi et al. | Online Calibration Procedure for Motion Tracking with Wearable Sensors Using Kalman Filtering | |
van Bijlert et al. | Muscle-controlled physics simulations of the emu (a large running bird) resolve the grounded running paradox |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180119 |