一种基于运动学数据确定步态站立中期的方法与流程

本发明涉及一种确定步态站立中期的方法。特别是涉及一种基于运动学数据确定步态站立中期的方法。
背景技术：
：步态是人类步行的行为特征，正常步态具有明显的周期性。通常使用七个主要事件即首次触地、对侧足趾离地、足跟抬起、对侧首次触地、足趾离地、双足邻近和胫骨垂直将步态周期划分为承重反应期、站立中期、站立末期、摆动前期、摆动初期、摆动中期和摆动末期共七个阶段。人类步行常简化为倒立摆模型，即假定人体质心(com)相对人体位置不变，其绕支撑侧踝关节做单摆运动。其中摆长使用com至外踝尖的距离估计。将站立中期定义为在前进方向即com瞬时速度vcom方向上com经过外踝尖的时刻。该时刻对内外向干扰极为敏感，可导致跌倒或额外迈步，故将其作为关键事件。确定站立中期，有利于下一步进行内外稳定性分析、步态特征提取以及生物力学计算等。传统的站立中期使用com轨迹来进行估计，即com达到其最高点时认为达到站立中期，并以此区别跑步。但实际操作中发现其误差较大，且不符合其原始定义。技术实现要素：本发明所要解决的技术问题是，提供一种能够较为准确的站立中期的基于运动学数据确定步态站立中期的方法。本发明所采用的技术方案是：一种基于运动学数据确定步态站立中期的方法，包括如下步骤：1)提取人体在步行中骨盆和足上的6个部位的三维位置坐标；2)估计人体在步行中质心的三维位置坐标；3)计算人体在步行中质心的速度；4)确定步态站立中期。步骤1)所述的6个部位是：右侧髂前上棘、左侧髂前上棘、右侧髂后上棘、左侧髂后上棘、右侧外踝尖和左侧外踝尖。步骤2)所述的估计人体在步行中质心的三维位置坐标是由下式获得：其中，midpsis是左侧髂后上棘lpsis与右侧髂后上棘rpsis的中点，lpsis是人体在步行中左侧髂后上棘的三维位置坐标，rpsis是人体在步行中右侧髂后上棘的三维位置坐标，com是人体在步行中质心的三维位置坐标，lasis是人体在步行中左侧髂前上棘的三维位置坐标，rasis是人体在步行中右侧髂前上棘的三维位置坐标。步骤3)所述的计算人体在步行中质心的速度是采用下式计算：其中，vcom(1)是第1帧质心的瞬时速度；vcom(i)是第i帧质心的瞬时速度；vcom(end)是最后一帧质心的瞬时速度；δt为采样间隔时间，即采样频率的倒数；com(1)是第1帧质心的三维位置坐标；com(2)是第2帧质心的三维位置坐标；com(i)是第i帧质心的三维位置坐标；com(end)是最后一帧质心的三维位置坐标。步骤4)所述确定步态站立中期是采用下式：min(|coslmidstanceangle|)对应的lmidstanceangle→90°式中，lmidstanceangle是左侧站立期内摆长地面投影与质心瞬时速度地面投影的夹角；lpendulumlengthprojection是左侧站立期内摆长地面投影；vcomprojection是质心瞬时速度地面投影。以左侧为例，连结com与左侧外踝尖，即左侧摆长，左侧摆长地面投影lpendulumlengthprojection与质心瞬时速度地面投影vcomprojection的内积除以二者模的乘积得到lmidstanceangle的余弦值，当该余弦值的绝对值取最小值的时候，有lmidstanceangle最接近90°，该时刻即站立中期。本发明的一种基于运动学数据确定步态站立中期的方法，利用合理的简化方法，通过一个夹角来逼近站立中期时刻，严格的数学运算，较为准确的确定步态站立中期这一关键事件。本发明原理可靠，计算便捷，实施简单，精度较高，从而为进一步的内外稳定性分析、步态特征提取以及生物力学计算等奠定基础。本技术较传统方法更符合定义并更加准确。附图说明图1是本发明一种基于运动学数据确定步态站立中期的方法的流程图；图2是站立中期的确定示意图。具体实施方式下面结合实施例和附图对本发明的一种基于运动学数据确定步态站立中期的方法做出详细说明。本发明的一种基于运动学数据确定步态站立中期的方法为：首先提取采集到的右侧髂前上棘(rasis)、左侧髂前上棘(lasis)、右侧髂后上棘(rpsis)、左侧髂后上棘(lpsis)、右侧外踝尖(rlateralmalleolus)和左侧外踝尖(llateralmalleolus)共6个标记点的三维运动捕捉数据，其次根据髂嵴处的四个标记点估计人体质心(com)的三维坐标，然后计算人体质心瞬时速度vcom及其地面投影vcomprojection，最后连结com与外踝尖，确定其地面投影与vcomprojection正交之时即为站立中期。但由于采样频率限制，其正交时刻可能难以捕捉，故计算该侧站立期内摆长地面投影与质心瞬时速度地面投影的夹角midstanceangle，取其最接近90°的时刻作为站立中期。如图1所示，本发明的一种基于运动学数据确定步态站立中期的方法，包括如下步骤：1)提取人体在步行中骨盆和足上的6个部位的三维位置坐标；所述的6个部位如表1所示为：右侧髂前上棘、左侧髂前上棘、右侧髂后上棘、左侧髂后上棘、右侧外踝尖和左侧外踝尖。表1.标记点位置说明标记点名称位置标记点名称位置rasis右侧髂前上棘lasis左侧髂前上棘rpsis右侧髂后上棘lpsis左侧髂后上棘rlateralmalleolus右侧外踝尖llateralmalleolus左侧外踝尖此处默认各标记点三维位置坐标(x，y，z)，其中，x轴为前后方向，y轴为左右方向，z轴正方向为竖直向上，坐标系服从右手定则。本发明实施例是使用运动捕捉系统如vicon(viconmotionsystemltd,oxford,uk)或motion(cortex,motionanalysis,santarosa,ca)采集人类步态的运动学数据。然后将采集的三维运动捕捉数据导入数据处理软件如matlab(themathworksinc.,natick,ma,usa)或labview(nationalinstruments,austintxusa)。选取合适的滤波器进行零相位滤波，以移除噪声，避免相移。然后从中提取出所需的lasis、rasis、lpsis、rpsis、llateralmalleolus、rlateralmalleolus共6个标记点。2)估计人体在步行中质心的三维位置坐标；质心坐标由右侧髂前上棘(rasis)、左侧髂前上棘(lasis)、右侧髂后上棘(rpsis)、左侧髂后上棘(lpsis)四个标记点的坐标估计得到，即先计算左侧髂后上棘lpsis与右侧髂后上棘rpsis的中点midpsis，再计算其与左侧髂前上棘lasis和右侧髂前上棘rasis构成的三角形的重心，即人体质心com，所述的估计人体在步行中质心的三维位置坐标是由下式获得：其中，midpsis是左侧髂后上棘lpsis与右侧髂后上棘rpsis的中点，lpsis是人体在步行中左侧髂后上棘的三维位置坐标，rpsis是人体在步行中右侧髂后上棘的三维位置坐标，com是人体在步行中质心的三维位置坐标，lasis是人体在步行中左侧髂前上棘的三维位置坐标，rasis是人体在步行中右侧髂前上棘的三维位置坐标。3)计算人体在步行中质心的速度；对人体质心坐标com进行中心差分，计算从第1帧到最后一帧(end)质心的瞬时速度vcom，亦即步行速度。所述的计算人体在步行中质心的速度是采用下式计算：其中，vcom(1)是第1帧质心的瞬时速度；vcom(i)是第i帧质心的瞬时速度；vcom(end)是最后一帧质心的瞬时速度；△t为采样间隔时间，即采样频率的倒数；com(1)是第1帧质心的三维位置坐标；com(2)是第2帧质心的三维位置坐标；com(i)是第i帧质心的三维位置坐标；com(end)是最后一帧质心的三维位置坐标。4)确定步态站立中期。连结com与外踝尖，即倒立摆模型的摆长，当摆长地面投影pendulumlengthprojection与com瞬时速度地面投影vcomprojection正交之时即为站立中期。但由于采样频率限制，其正交时刻可能难以捕捉，故计算该侧站立期内摆长地面投影与com瞬时速度地面投影的夹角midstanceangle，取其最接近90°的时刻作为站立中期。以左侧为例，如图2与公式(3)所示。其中二维向量lpendulumlengthprojection与vcomprojection的内积除以二者模的乘积得到lmidstanceangle的余弦值，当该余弦值的绝对值取最小值的时候，有lmidstanceangle最接近90°，该时刻即站立中期。所述确定步态站立中期是采用下式：min(|coslmidstanceangle|)对应的lmidstanceangle→90°式中，lmidstanceangle是左侧站立期内摆长地面投影与质心瞬时速度地面投影的夹角；lpendulumlengthprojection是左侧站立期内摆长地面投影；vcomprojection是质心瞬时速度地面投影。本发明的一种基于运动学数据确定步态站立中期的方法，站立中期作为步态的一个关键事件，本技术发明将之确定后可用于步态特征的提取，然后进行后续的步态分析。站立中期时刻对内外向干扰极为敏感，可导致跌倒或额外迈步，故使用本技术发明确定站立中期后，可计算该时刻的内外稳定性，作为步态稳定性评估的一个重要的量化指标。对于一些需要计算步态中的生物力学等场景，如假肢、矫形器、内固定、关节假体等，使用本技术发明确定站立中期后，可重点对该关键事件时刻进行力学分析，从而实现不同个体间的对比分析。当前第1页12