一种外骨骼设备穿戴误差校正方法与流程

1.本技术属于运动检测技术领域，具体涉及一种外骨骼设备穿戴误差校正方法。


背景技术：

2.近年来，随着科学技术的逐步发展和国家医疗卫生水平的显著提高，穿戴式外骨骼设备在健康养老、医疗康复和运动健身等领域得到了飞速的发展。外骨骼设备在正常工作时需要准确地获得穿戴者的关节运动情况，作为外骨骼控制器的参考信息或者是用于分析穿戴者的运动、康复等情况。
3.在所有穿戴式运动测量方案中，惯性测量单元(以下简称imu)以其低成本、便携性好、动态响应快以及不受使用场景限制等优点成为主要采用的传感器。但是使用imu进行关节运动角度监测时存在以下几点问题：(1)如专利申请公布号为cn110522458a和cn108852360a的专利所述，几乎所有使用imu的产品都直接或者间接地要求传感器必须固定于标准平面内(一般平行或者垂直于人体矢状面)。然而人体不是一种规则的几何体，尤其是涉及到髋、膝等多关节的运动测量，在穿戴设备时无法保证imu的轴向能够完全沿着人体运动方向；(2)imu以欧拉角作为姿态测量的输出，然而只使用三个参数的欧拉角表示法并非是一种全局的、非奇异的刚体运动描述方法，况且由于欧拉角存在“万向节锁”和复杂的周期性等问题，使用欧拉角来表示关节的运动角度并非那么直观，尤其是在(1)中所述问题存在的情况下，使用欧拉角进行人体运动姿态描述会变得非常不直观。
4.基于imu使用存在的问题，现有技术中也提出了一些相应解决方案，例如专利申请公布号为cn111166346a的专利文献中描述了一种仅使用角速度传感器进行人体运动检测的装置，并且根据人体下肢运动的特点，可以在不对穿戴有任何要求的前提下测量出膝关节的运动角度。但是该专利所述方法需要先离线优化，通过对人体运动数据进行分析得到运动轴，之后再进行实时角度解算，复杂度高且解算费时。又如专利申请公布号为cn110646014a的专利文献描述了一种基于人体关节位置捕捉设备进行imu穿戴校准的方法，该方法可以直接在线使用，但是需要借助除外骨骼以外的人体关节位置捕捉设备对imu的穿戴偏差进行校准，增加了设备复杂度，且需要进行设备间的数据转换计算。


技术实现要素：

5.本技术的目的在于提供一种外骨骼设备穿戴误差校正方法，在不借助于外部设备的情况下，准确测量穿戴者的关节角度，为外骨骼设备或者康复监测设备提供准确的运动信息反馈。
6.为实现上述目的，本技术所采取的技术方案为：
7.一种外骨骼设备穿戴误差校正方法，所述外骨骼设备安装有imu，所述外骨骼设备穿戴误差校正方法，包括：
8.步骤1、在外骨骼设备的穿戴者直立情况下进行imu自校准；
9.步骤2、在外骨骼设备的穿戴者平躺情况下进行下肢平躺校正，包括：
10.步骤2.1、检测穿戴者是否满足下肢平躺校正条件，所述下肢平躺校正条件为：穿戴者下肢脚尖朝上平躺于水平面上并保持静止状态至少2秒；
11.步骤2.2、若检测到穿戴者不满足下肢平躺校正条件，则返回步骤2.1继续检测；否则执行步骤2.3；
12.步骤2.3、获取imu测量得到的加速度值a
x
和ay；
13.步骤2.4、根据加速度值a
x
和ay计算imu安装平面与人体矢状面的夹角α；
14.步骤2.5、根据步骤2.4得到的夹角α计算并得到角速度向量ωv，计算公式如下：
[0015][0016]
式中，δt表示imu的采样时间，单位为秒，角速度向量ωv的单位为弧度/秒；
[0017]
步骤3、在外骨骼设备的穿戴者正常运动情况下，将imu测量得到的欧拉角转换为四元数形式，并利用记录的角速度向量ωv对四元数进行校正，根据校正后的四元数换算得到补偿后的欧拉角作为最终的关节角度，完成外骨骼设备穿戴的误差校正。
[0018]
以下还提供了若干可选方式，但并不作为对上述总体方案的额外限定，仅仅是进一步的增补或优选，在没有技术或逻辑矛盾的前提下，各可选方式可单独针对上述总体方案进行组合，还可以是多个可选方式之间进行组合。
[0019]
作为优选，所述根据加速度值a
x
和ay计算imu安装平面与人体矢状面的夹角α，包括：
[0020]
若外骨骼设备穿戴在穿戴者下肢外侧面，则计算得到的夹角α为：
[0021][0022]
若外骨骼设备穿戴在穿戴者下肢正面，则计算得到的夹角α为：
[0023][0024]
式中，夹角α为弧度制。
[0025]
作为优选，所述将imu测量得到的欧拉角转换为四元数形式，并利用记录的角速度向量ωv对四元数进行校正，根据校正后的四元数换算得到补偿后的欧拉角作为最终的关节角度，包括：
[0026]
步骤3.1、将imu测量得到的欧拉角转换为四元数形式如下：
[0027][0028]
式中，q为四元数，ψ、θ、φ为欧拉角；
[0029]
步骤3.2、将四元数q以及角速度向量ωv代入下式：
[0030][0031]
式中，为角速度向量ωv的四元素形式，为校正后的四元数；
[0032]
步骤3.3、根据校正后的四元数换算得到补偿后的欧拉角。
[0033]
本技术提供的外骨骼设备穿戴误差校正方法，先获取imu传感器安装平面相对于人体矢状面倾斜的度数，然后在数值层面对imu传上的角度数据进行偏差消除。在校正偏差的过程中不需要借助于任何外部测量仪器，仅仅需要穿戴者进行站立校准和下肢平躺校正两步流程，就可以准确测量出穿戴者髋关节的伸展/弯曲、外展/内收以及膝关节的伸展/弯曲等运动的关节角度。
附图说明
[0034]
图1为本技术的外骨骼设备穿戴误差校正方法的流程图；
[0035]
图2为本技术外骨骼设备穿戴的俯视图。
具体实施方式
[0036]
下面将结合本技术实施例中的附图，对本技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是本技术一部分实施例，而不是全部的实施例。基于本技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本技术保护的范围。
[0037]
除非另有定义，本文所使用的所有的技术和科学术语与属于本技术的技术领域的技术人员通常理解的含义相同。本文中在本技术的说明书中所使用的术语只是为了描述具体的实施例的目的，不是在于限制本技术。
[0038]
目前外骨骼设备在穿戴过程中存在无法保证imu的轴向能够完全沿着人体运动方向，也不能保证imu的安装平面完全贴合于人体解剖面的问题。在设备穿戴出现上述偏差时，imu给出的欧拉角并不能直观反映出髋、膝关节沿某一轴向的运动角度。然而无论是助力控制领域还是康复监测领域，普遍都是对人体下肢沿特定轴向的运动角度感兴趣，需要获取关于特定轴向的精确角度，这导致目前外骨骼设备在应用上存在较大问题。
[0039]
针对上述问题，其中一个实施例中，提供一种外骨骼设备穿戴误差校正方法，在imu设备穿戴出现偏差的情况下，在不借助于外部设备的条件下对imu的穿戴偏差进行校正，从而准确测量穿戴者髋、膝关节沿特定运动轴向的运动角度，为外骨骼助力设备或者康复监测设备提供准确的运动信息反馈。
[0040]
如图1所示，本实施例中的外骨骼设备穿戴误差校正方法，包括：
[0041]
步骤1、在外骨骼设备的穿戴者直立情况下进行imu自校准。
[0042]
imu自校准是指目前市场上的惯性测量单元(简称imu)自带的校准程序，主要是为了消除陀螺仪和加速度计的偏差。根据不同型号的imu启动对应的校准程序即可，这里不对imu自带的校准程序进行详述。
[0043]
并且在穿戴者直立情况下进行imu自校准中对直立情况的判断可以是主动触发或者是被动触发方式，主动触发例如对穿戴者的状态进行检测，满足直立状态时主动触发进
行imu自校准，被动触发例如在穿戴者切换至直立状态后收到触发imu进行自校准。
[0044]
步骤2、在外骨骼设备的穿戴者平躺情况下进行下肢平躺校正。
[0045]
在穿戴者平躺情况下进行imu校正，是为了得到imu安装平面相对于人体矢状面(将身体分为左右两部分的纵切面)倾斜的角度，以作为后期对imu获取的关节角度数据的补偿基础。让imu以其所安装的平面与人体矢状面平行的方式来解算欧拉角，达到偏差矫正的目的。
[0046]
根据外骨骼设备的穿戴情况，本实施例将外骨骼设备穿戴于人体下肢外侧面且平行于人体矢状面的状态定义为正确佩戴的状态，即此时经过imu自校准后即可得到有效的运动角度数据。其中下肢的外侧面理解为下肢的外侧，即两腿相对的一侧作为下肢的内侧，两腿相背的一侧作为下肢的外侧。
[0047]
本实施例中记下肢坐标系为b系，xb轴方向定义为垂直于人体冠状面(将身体分为前后两部分的纵切面)且指向人体后方的方向，定义yb轴为位于人体冠状面内与xb轴垂直且指向人体右侧的方向，zb轴根据右手定则(指三维坐标系定义所使用的右手定则)确定。记imu坐标系为s系，zs轴与zb平行且指向同一个方向，xs与xb位于人体横断面内且夹角为α，ys轴根据右手定则确定。且取imu安装平面为平行于xs轴且垂直于ys轴。
[0048]
如图2所示，本实施例以穿戴于人体左侧下肢且由左侧下肢的外侧面起围绕zb轴的负方向偏斜(偏斜后同样视为穿戴于下肢外侧面)的外骨骼设备为例说明本技术的误差校正方法。
[0049]
具体的，本实施例中下肢平躺校正具体包括以下步骤：
[0050]
步骤2.1、检测穿戴者是否满足下肢平躺校正条件。为了保证校正效果，本实施例中的下肢平躺校正条件设为：穿戴者下肢脚尖朝上平躺于水平面上并保持静止状态至少2秒。
[0051]
若外骨骼设备佩戴于人体左侧下肢的外侧面，则检测是否满足下肢平躺校正条件时，可以基于外骨骼设备的陀螺仪和imu进行检测，若陀螺仪数据未发生变化且膝关节的imu和髋关节的imu同时满足采集的加速度向量与[-1,0,0]这一向量平行，保持至少2秒后即认为满足下肢平躺校正条件。
[0052]
但考虑到实际应用过程中外骨骼设备存在佩戴倾斜的情况，因此为了较为准确的判断是否下肢平躺校正条件，可设定为直接接收用户输入的满足信号，当接收到满足信号即判断为满足下肢平躺校正条件，否则不满足下肢平躺校正条件。
[0053]
步骤2.2、若检测到穿戴者不满足下肢平躺校正条件，则返回步骤2.1继续检测；否则执行步骤2.3。由于下肢平躺校正为必须进行的校准步骤，因此若穿戴者姿势不满足将一直重新执行检测，直到满足条件为止。
[0054]
步骤2.3、获取imu测量得到的加速度值a
x
和ay。在穿戴者静止位于地面上时，其中的加速度值a
x
为当地重力加速度沿xs轴的分量，加速度值ay为当地重力加速度沿ys轴的分量。
[0055]
容易理解的是，当穿戴者处于运动状态或者相对地面不静止时，施加在imu上的加速度则不一定是当地重力加速度，则加速度值a
x
为实际施加在imu上的加速度沿xs轴的分量，加速度值ay为实际施加在imu上的加速度沿ys轴的分量。
[0056]
步骤2.4、根据加速度值a
x
和ay计算imu安装平面与人体矢状面的夹角α(弧度制)：
[0057][0058]
步骤2.5、根据步骤2.4得到的夹角α计算角速度向量ωv并记录。角速度向量ωv的计算公式如下：
[0059][0060]
式中，δt表示imu的采样时间，单位为秒，角速度向量ωv的单位为弧度/秒。
[0061]
步骤3、在外骨骼设备的穿戴者正常运动情况下，将imu测量得到的欧拉角转换为四元数形式，并利用记录的角速度向量ωv对四元数进行校正，根据校正后的四元数换算得到补偿后的欧拉角作为最终的关节角度，完成外骨骼设备穿戴的误差校正。
[0062]
为避免欧拉角解算引起的奇异性问题，在角度补偿阶段利用四元数对imu解算得到的角度进行校正。因此本实施例中对欧拉角进行补偿的过程如下：
[0063]
定义四元数(quaternions)定义如下：
[0064][0065]
式中，q0表示四元数的标量部分，q表示四元数的矢量部分，q＝[q
0 q
1 q
2 q3]
t
。四元数表示刚体旋转时，其动力学方程有解析解，离散形式如下：
[0066][0067]
其中，为四元数乘法运算的符号，q
ω
＝[0，ω
t
]
t
为imu旋转角速度ω的四元素(four-component quantity)表达形式，q(t)表示t时刻的四元数，q(t+δt)表示t+δt时刻的四元数，δt表示imu的采样时间。这里四元数的指数运算可以根据旋转矢量和四元数之间的映射公式进行简便运算，令：
[0068][0069]
则有如下映射成立：
[0070][0071]
为了对设备的穿戴偏差进行校正，我们必须对原始的四元数(经过imu上传的欧拉角转换得到的四元数)进行校正，本实施例根据四元数的计算原理提出四元数的校正公式如下：
[0072][0073]
式中，q表示原始的四元数，表示校正后的四元数，为角速度
向量ωv的四元素形式，只需要将公式(2)中计算的角速度向量ωv带入到的计算公式中即可得到校正后四元数。符号ωv中的下标表示虚拟含义，用于将imu旋转角速度与角速度向量区分表示。
[0074]
因此在穿戴者进行髋、膝关节运动或者康复治疗时，对imu测量的欧拉角进行校正的具体步骤为：
[0075]
步骤3.1、将imu测量得到的欧拉角转换为四元数形式如下：
[0076][0077]
式中，q为四元数，ψ、θ、ψ为欧拉角，ψ为航向角，θ为俯仰角，φ为翻滚角。
[0078]
步骤3.2、将四元数q以及角速度向量ωv代入式(7)中得到校正后的四元数。
[0079]
步骤3.3、根据校正后的四元数换算得到补偿后的欧拉角。其中将四元数换算成欧拉角可以根据以下式子进行换算：
[0080][0081]
经过矫正以后的欧拉角将能够准确、直接地反映出髋、膝关节沿特定轴向的运动角度。
[0082]
需要说明的是，本实施例以穿戴于人体左侧下肢且由左侧下肢的外侧面起围绕zb轴的负方向偏斜的外骨骼设备为例进行说明，但这并不限定本技术的误差校正方法仅能应用于该穿戴方式。
[0083]
本技术的误差校正方法适用于四肢的任一穿戴方式，以上述穿戴于人体左侧下肢且由左侧下肢的外侧面起围绕zb轴的负方向偏斜的外骨骼设备的误差校正方式作为基础校正方式，其中穿戴于人体左侧下肢且由左侧下肢的外侧面起围绕zb轴的正方向偏斜时的误差校正方式与基础校正方式相同；穿戴于人体右侧下肢且由右侧下肢的外侧面起围绕zb轴的正/负方向偏斜时的误差校正方式与基础校正方式相同。
[0084]
在外骨骼设备的穿戴中还有一种较为常规的穿戴方式为穿戴在穿戴者下肢的正面，因此可以将外骨骼设备穿戴于人体下肢正面且平行于人体冠状面的状态同样定义为正确佩戴的状态，即此时经过imu自校准后即可得到有效的运动角度数据。其中下肢的正面理解为下肢位于穿戴者面部的一面。当外骨骼设备穿戴在人体下肢正面时同样可能存在佩戴倾斜的情况，即存在由下肢的正面起围绕zb轴的正/负方向偏斜的情况，当外骨骼设备为穿戴于人体下肢正面的正确佩戴状态或佩戴倾斜的状态下时，其误差校正方法相对于基础校正方式而言仅为计算夹角α的公式有区别，其他步骤或公式均相同，且该情况下夹角α的公式为：
[0085]
至此本技术提出的误差校正方法适用于外骨骼设备穿戴在人体下肢的任一穿戴方式中。其中外骨骼设备佩戴在下肢的外侧面还是正面可借助外部传感器测量判断，也可
以直接由穿戴者输入。
[0086]
对于外骨骼设备穿戴在人体上肢的情况与外骨骼设备穿戴在人体下肢的情况相同，本实施例中不再进行重复说明。
[0087]
另外，由于外骨骼设备中的imu可采用任意型号，而不同型号之间的imu可能其坐标系的坐标轴不一样，但其坐标轴均可以相对于imu安装平面进行分析。例如本实施例中所取的坐标轴中xs轴为平行于imu安装平面，ys轴为垂直于imu安装平面。则当imu的坐标轴变换时，取坐标轴中与imu安装平面平行的轴等价于本技术中的xs轴，即该轴上的加速度分量作为本技术中的a
x
进行计算，取坐标轴中与imu安装平面垂直的轴等价于本技术中的ys轴，即该轴上的加速度分量作为本技术中的ay进行计算即可，因此本技术的误差校正方法适用于具有任意imu型号的外骨骼设备。
[0088]
本实施例显著降低了外骨骼设备的穿戴难度，现有技术要求在穿戴阶段需要严格保持imu安装平面位于人体解剖面内，或者要求增加额外的校准设备，或者需要离线的复杂运算，增加设备的使用复杂度。而本实施例所提供方法仅仅要求穿戴者按照正常的穿戴方式，将外骨骼设备绑定于人体下肢，在穿戴过程中不需要进行任何校正操作。而在穿戴完成后，仅仅需要执行简单的校正程序就可以实现精确的髋、膝关节运动角度测量。尤其适合于现在的使用imu进行运动角度测量的设备使用。
[0089]
以上所述实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。
[0090]
以上所述实施例仅表达了本技术的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本技术构思的前提下，还可以做出若干变形和改进，这些都属于本技术的保护范围。因此，本技术专利的保护范围应以所附权利要求为准。