一种步态分析方法与流程

本发明涉及步态分析技术领域，尤其涉及一种步态分析方法。

背景技术：

步态是人体结构与功能、运动调节系统、行为及心理活动在行走时的外在表现，但是，其中某个系统或某些方面的功能障碍，都可引起步态异常。足内旋本身是一种正常的动作，它是脚部为了防止脚内侧足弓受伤害而自然产生的一系列缓冲动作。但是长期的内旋不足或者内旋过度都有可能引起跑步损伤。脚内旋会破坏正常负重路径，导致腿向内过度翻转。脚过度内旋是一种不良形态，可能导致不适和损伤。外旋则与内旋相反，指脚向外翻卷。脚过度内旋和外旋引起的常见疾病包括：足弓疼痛、脚跟疼痛、扁平足、膝盖疼痛、踝关节扭伤、肌腱炎、关节疼痛、腰背痛、胫夹疼痛和/或应力性骨折。因此，有必要提供一种方便，快捷，精准的步态分析方法，尤其是分析足内旋和外旋的方法。

技术实现要素：

本发明的目的在于克服上述不足，提供一种步态分析方法。

为实现上述目的，本发明的技术解决方案是：一种步态分析方法，对信号采集系统采集传输的数据进行分析，其具体分析过程如下：

1)定位完全着地的时间点:设置初始阈值为0.5g，方波信号为0，将原始信号放大3倍，消除干扰，对信号数据进行低通滤波处理，利用y轴加速度生成方波信号，若检测到方波是为下降边沿，开始检测x轴角速度极大值，判断极大值是否大于500dps，若否，则以该极大值的时间点为着地时间点，若是，则继续寻找以角速度为零的时间点为完全着地时间点；

2)对信号进行步态检测：检测信号的极大值，从极大值开始寻找与极大值相差大于阈值的实时信号数据，实时信号数据与极大值的差值大于阈值时将方波信号置为-1，在方波信号为-1时开始检测y轴加速度的最小值，并检测完全着地时间点，在找到完全着地时间点后，方波信号置为0，以极大值和最小值的差值的0.8倍作为下一步的阈值；

3)单步信号划分：检测生成的方波信号的下降沿，在实时信号大于最小值0.2g时(减少过程误差)开始检测角速度的极大值判断极大值的大小，小于零时以该极值点作为一步结束的时间，大于零时，从极值点位置开始寻找角速度为零的时间点，作为一步的结束；

4)计算角度的大小：定位着地后脚掌翻转的时间，对角速度进行积分，求解出角度θ，θ＝∫(w*t),式中w为角速度，t为采样时间；

5)判断着地方式：判断该角度的正负关系以及绝对值大小，若θ绝对值小于α，则当前步态为全掌着地；若θ绝对值大于α，且θ方向为正，则当前步态为前脚掌着地；若θ绝对值大于α，且θ方向为负，则当前步态为后脚掌着地；

6)判断足内旋和足外旋：若角度θ大于β，则当前步态为足外旋；若角度θ小于或等于0°，则当前步态为足内旋；当角度θ属于0°到β之间，则为正常内旋；

7)判断走跑方式：定位x轴离地时间和着地时间，判断触地腾空比是否小于1，若小于1，则为跑步状态，反之则为走路状态，其中触地腾空比为单步着地时间/腾空时间。

优选的，所述采集系统为六轴加速度传感器，至少采集在空间x、y、z三个方向上加速度数据，所述x、y、z方向为设定的方向。

优选的，所述的滤波方式采用二阶巴特沃斯低通滤波器进行滤波。

优选的，所述步态检测还包括极大值的更新：在检测到极大值后极大值的值会锁定，在之后若检测到10个以上的极大值时该极大值的值重新确定。或若后面的极大值大于该极大值，将极大值替换为较大的极大值。

本发明由于采用了以上技术方案，具有显著的技术效果：本方法通过对用户在行走或运动过程中基于足部的各种参数进行实时采集和算法处理，能够实现实时检测与分析用户步态，并且分析过程方便快捷，准确度高；在鞋子上添加六轴加速度传感器配合移动终端，整个足内旋和足外旋的判断方法更加全面、准确，成本低廉，有利于运动及早发现进行矫正。

附图说明

图1为本发明步态分析方法的流程示意图；

图2为本发明定位完全着地的时间点流程示意图；

图3为本发明的步态检测信号示意图；

图4为本发明的单步划分信号示意图；

图5为本发明x轴角速度示意图；

图6为本发明判断判断着地方式示意图；

图7为本发明判断判断足内旋和足外旋示意图。

图中：1、y轴机速度信号；2、生成的方波信号；3、单步划分信号；4、y轴加速度；5、方波信号下降沿；6、直线一；7、直线二。

具体实施方式

以下结合附图和具体实施例来进一步说明本发明。

如图1-6所示，一种步态分析方法，对信号采集系统采集传输的数据进行分析，所述采集系统为六轴加速度传感器，至少采集在空间x、y、z三个方向上加速度数据，所述x、y、z方向为设定的方向，其具体分析过程如下：

1)参照图2，定位完全着地的时间点：设置初始阈值为0.5g，方波信号为0，将原始信号放大3倍，消除干扰，对信号数据进行低通滤波处理，利用y轴加速度生成方波信号，若检测到方波是为下降边沿，开始检测x轴角速度极大值，判断极大值是否大于500dps，若否，则以该极大值的时间点为着地时间点，若是，则继续寻找以角速度为零的时间点为完全着地时间点；

还可以采用通带滤波处理方式，优选的，本发明使用二阶巴特沃斯低通滤波器进行滤波。滤波公式如下：

data_fil＝

(0.0201*data1+0.0402*data2+0.0201*data3+1.5610f*data_fil1-0.6414*data_fil2)

其中datafil为某一时刻滤波后的信号，datafil1、datafil2分别为上一时刻和上上一时刻滤波后的信号，data1、data2、data3分别为某一时刻的原始数据、上一时刻的原始数据和上上一时刻的原始数据。

2)参照图3，对信号进行步态检测：本发明分别用加速度计的y轴加速度的特性进行步态进行检测和单步信号划分(运动过程中y轴加速度信号和陀螺仪x轴角速度的主要特征不会因使用者的不同而出现较大变化，具有良好的普适性)。本发明使用动态阈值的方法生成方波信号进行步态检测。相对于静态阈值，动态阈值检测的准确性更高。检测信号的极大值，从极大值开始寻找与极大值相差大于阈值的实时信号数据，实时信号数据与极大值的差值大于阈值时将方波信号置为-1，在方波信号为-1时开始检测y轴加速度的最小值，并检测完全着地时间点，在找到完全着地时间点后，方波信号置为0，以极大值和最小值的差值的0.8倍作为下一步的阈值；在检测到极大值后极大值的值会锁定，在之后若检测到10个以上的极大值时该极大值的值重新确定。或若后面的极大值大于该极大值，将极大值替换为较大的极大值；

3)参照图4，进行单步信号划分：检测生成的方波信号的下降沿，在实时信号大于最小值0.2g时(减少过程误差)开始检测角速度的极大值判断极大值的大小，小于零时以该极值点作为一步结束的时间，大于零时，从极值点位置开始寻找角速度为零的时间点，作为一步的结束，一步的结束即为下一步的开始；

4)计算角度的大小：定位着地后脚掌翻转的时间，对角速度进行积分，求解出角度θ，θ＝∫(w*t),式中w为角速度，t为采样时间(以0.005s为例)；

参照图5，角度的计算原理：脚掌着地瞬间由于鞋底与地面的缓冲，陀螺仪x轴角速度会产生一个微小的突变(图中直线一6的位置)，以该时间点为脚掌翻转的起点，直线二7的位置为脚掌完全着地的时间点，该时间点由着地时后y轴加速度趋近于零的极值点来确定，积分该过程陀螺仪x轴角速度得到角度变化量。

5)参照图6所示，判断着地方式：判断该角度的正负关系以及绝对值大小，若绝对值小于10°，则当前步态为全掌着地；若绝对值大于10°，且方向为正，则当前步态为前脚掌着地；若绝对值大于10°，且方向为负，则当前步态为后脚掌着地；

6)参照图7所示，判断足内旋和足外旋转：若角度大于10°，则当前步态为足外旋；若角度小于或等于0°，则当前步态为足内旋；反之，则为正常内旋；正常情况下着地时脚掌向内翻转的角度在0°～10°时为正常内旋，当内旋角度大于10°时为足外旋，小于或等于0°时为足内旋。

7)判断走跑方式：定位x轴离地时间和着地时间，判断触地腾空比是否小于1，若小于1，则为跑步状态，反之则为走路状态，其中触地腾空比为单步着地时间/腾空时间。

当使用者在穿着过程中，鞋子通过六轴传感器采集到信号数据传输到移动终端，移动终端对信号数据进行处理，分析出步态信息，反馈到鞋子的控制芯片上，鞋子的控制芯片根据不同的步态信息，控制气囊的变化，当该气压值大于预设的最高气压值时，气囊放气，直至当前气囊的气压值与预设的最高气压值接近或相同；当该气压值小于预设的最高气压值时，气囊充气，直至当前气囊的气压值与预设的最高气压值接近或相同。

以上所述的，仅为本发明的较佳实施例而已，不能限定本发明的范围，凡是依本发明申请专利范围所作的均等变化与装饰，皆应仍属于本发明涵盖的范围内。