一种踝关节三维运动空间的量化算法
张云超;肖金壮;王洪瑞
【摘 要】目的 结合踝关节在各方向上运动之间的约束关系,对其三维转动角度进行积分,得到三维运动空间的大小,用于更加准确全面地量化其灵活度和康复程度.方法 选取25例大学生志愿者,其中无踝关节损伤的健康者20例,近期有过踝关节损伤的患者5例.使用运动捕捉系统实时记录被测者脚上标记点的三维坐标.利用三角函数将这些坐标转化为转动角度并进行处理,通过使用多项式拟合,得到踝关节在各个方向上最大转动角度随跖屈和背屈角度变化的方程.对拟合得到的方程进行积分,得到踝关节三维运动空间的大小.结果 健康人群踝关节的三维运动空间为41.256,个体间差距小于10%.患者患侧三维运动空间的大小为33.163,由于康复程度的不同,个体间差距小于15%.两组之间存在显著性差异(t=8.804,P<0.01).在同一跖屈和背屈角度下,内翻的最大转动角度比外翻大,内旋的最大转动角度比外旋大.结论 通过此算法能更准确地量化踝关节三维运动空间的大小,实现对踝关节灵活度和康复程度更全面的评估.%Objective To obtain the three-dimensional (3D) motion space of ankle by an integral for the 3D rotation angle with mutual relationships of rotation angles in different directions, which is us
ed to evaluate the flexibility and recovery of ankle more accurately and ful-ly. Methods Twenty-five graduate students were recruited to participate in this study, one group consisted of 20 healthy individuals without ankle injury, and the other was compod of 5 patients with ankle injury. A motion capturing system was ud to simultaneously capture the 3D coordinates of the marked points on the foot. Next, the 3D coordinates were converted into rotation angles through trigonometric func-tions procesd. The maximum rotation angles of adduction, abduction, varus, and eversion change with the rotation angles of plantar flex-ion and dorsal flexure were measured by using polynomial fitting. At last, the size of the ankle's 3D motion space was obtained by integrat-ing the fitting functions. Results The ankle's 3D motion space of the healthy people was 41.256, and the variation among them was less than 10%. The ankle's 3D motion space of the patient side was 33.163, and the variation among them was less than 15%becau of the different degrees of rehabilitation. There was significant difference between two groups (t=8.804, P<0.01). With the same rotation angles of plantar flexion or dorsal flexure, the maximum rotation angle of varus was larger than that of eversion, and the maximum rotat
ion angle of adduc-tion was larger than that of abduction. Conclusion The algorithm can be ud to quantify the 3D motion space of ankle more accurately, and to achieve a more comprehensive evaluation of ankle's flexibility and rehabilitation.
【期刊名称】《中国康复理论与实践》儿子
【年(卷),期】2017(023)010
【总页数】7页(P1209-1215)
【关键词】踝关节;康复评定;三维运动空间;多项式拟合;积分;运动捕捉
【作 者】张云超;肖金壮;王洪瑞
【作者单位】河北大学电子信息工程学院,河北保定市071000;河北大学电子信息工程学院,河北保定市071000;河北大学电子信息工程学院,河北保定市071000
【正文语种】中 文
工作内容【中图分类】R322.7
踝关节的受损越来越常见，它的恢复周期较长，恢复效果常欠理想，影响人类日常行为[1-3]。随着科技的发展，踝关节康复机器人的种类越来越多，康复训练的方式也多种多样[4-7]。评价其康复程度的方法一般有两种：一种是根据医生的经验进行判断，这种评价方法依赖于医生的个人经验，主观性强且缺乏统一的标准；另一种是通过测量关节活动度进行判断。关节活动度指踝关节在各方向上运动时一维的转动范围。虽然它考虑到用准确数值表达足踝的灵活度，但并没有考虑各方向上运动的综合表达。因此，这两种评估踝关节康复程度的方法都存在不足。
笃笃笃卖糖粥足踝运动是一个立体空间内的三维运动，而且它在各方向上的运动都相互影响[8-9]。为了找到一个数值来描述踝关节在三维空间内运动空间的大小，既需要考虑踝关节在各个方向上的转动范围又要考虑踝关节在各个方向上转动之间的影响[10-11]。本研究的目的是通过一种算法量化踝关节三维运动空间的大小，并用来评价踝关节的灵活性和康复程度[12]，为建立更加合理的踝关节康复标准提供准确的数据和方法，也为踝关节康复训练提供反馈指导；同时为踝关节康复训练仪器以及此方面假肢机器人提供更加完善的数据，使其运动轨迹或运动空间更加接近人体，以提高设备的舒适性[13-15]。
悬崖歌词
人体踝关节是由胫骨的下关节面、内踝关节面和腓骨的外踝关节面共同形成的叉状关节窝，以及距骨滑车的关节头构成[16-17]。足可绕此关节做屈伸运动，足向下为跖屈，向上为背屈。由于关节窝和关节头之间存在缝隙，所以当足跖屈或背屈时，足还可以做内旋、外旋、内翻和外翻运动[18]。踝关节三维运动空间的大小应为三维转动角度的一种积分。由于内旋、外旋、内翻和外翻转动的角度都与跖屈和背屈的转动角度有关，而且这种关系容易被量化，因此选定跖屈和背屈的转动角度为自变量，设定跖屈时角度为负值，背屈时角度为正值。内旋、外旋、内翻和外翻的最大转动角度则被设定为因变量。通过角度测量设备可以获得每个跖屈或背屈角度所对应的内旋、外旋、内翻和外翻的最大转动角度。再利用数据拟合的方法得到它们的最大转动角度随跖屈和背屈转动角度变化的函数。
因为足踝运动空间的大小是通过结合各方向上运动的三维转动角度进行积分得到的，所以本文只需要考虑它在各个方向上的转动角度以及它们之间的相互约束关系，而不需要考虑它的形状。为了更加形象地把足踝运动空间表达和量化出来，本文把踝关节的转动部分(距骨滑车)看成一个半径为1的球体。如图1所示，当踝关节做跖屈和背屈运动时，该球绕X轴旋转；当踝关节做内旋和外旋运动时，该球绕Z轴旋转；当踝关节做内翻和外翻运动时，该球绕Y轴旋转。本研究在球的上半部分选择一个半圆作为基线(如图1中的黑色实线)，同
时，这条基线(半圆)在Y轴和Z轴构成的平面上。在每一个跖屈或背屈角度时，踝关节分别做最大内旋、外旋、内翻和外翻运动后基线所扫过的空间构成了一个不规则球面，而这个不规则球面所对应的那部分球体体积就可以代表踝关节在这个跖屈或背屈角度下，通过做四种运动所能达到的空间，称为子空间。由于本研究规定球的半径为1，所以得到的该子空间大小完全取决于某一个跖屈或背屈角度下内旋、外旋、内翻和外翻转动范围的大小以及它们之间的约束关系。再根据前面拟合的内旋、外旋、内翻和外翻的转动范围随跖屈和背屈的角度变化的方程，便可以得到子空间随跖屈和背屈角度变化的方程。最后对这个子空间进行关于跖屈和背屈的积分得到踝关节三维运动空间的大小。
其计算过程如下。
首先，假设内旋、外旋、内翻和外翻的最大转动角度随跖屈和背屈转动角度变化的函数分别为 f1(x)、f2(x)、 f3(x)、 f4(x)。
接下来计算每一个跖屈或背屈角度所对应的子空间的大小。只要自变量(跖屈和背屈的转动角度)确定，根据以上四个函数就可以确定相对应的内旋、外旋、内翻和外翻的最大转动角度，分别记为y1、y2、y3、y4。由于最大外翻与外旋、外翻与内旋、内翻和外旋、内翻与
内旋之间都分别存在重合部分，所以在计算某个跖屈或背屈角度所对应的子空间时不能直接使用公式(1)进行计算。应该分别求出外旋能达到而外翻不能达到的部分子空间v1；外旋能达到而内翻不能达到的部分子空间v2；内旋能达到而外翻不能达到的部分子空间v3；内旋能达到而内翻不能达到的部分子空间v4。然后通过公式(2)计算出通过内翻和外翻能达到的部分子空间v5。进而在某个跖屈或背屈角度下的运动子空间的大小可以通过公式(3)计算得到。
以计算外旋能达到而外翻不能达到的部分子空间为例。如图1所示，为了便于分析，假设球是相对静止而基线是转动的。当踝关节外翻时基线扫过球面旋转到图中的黑色虚线位置。当踝关节外旋时基线扫过球面旋转到图1中的白色实线位置。其中α作为一个辅助角，以便于几何分析，范围为[0,π/2]； β为踝关节外翻的最大转动角度，即为y4；θ为踝关节外旋的转动角度，其范围为[0,y2]。φ表示为球面上的任意一点与圆心构成的直线与Z轴的夹角。
基线的初始方程为：
踝关节最大外翻时，基线上点的坐标为：
结合公式(5)、(6)、(7)可得：
对公式(5)、(6)、(7)、(8)、(9)、(10)推导可得：
所以积分变量θ、φ、r的取值范围分别为：
故外旋能达到而外翻不能达到的部分子空间的大小v1为：
rumours
因此，v1为自变量(跖屈和背屈的转动角度)的变化函数。
同理，也可以得到外旋能达到而内翻不能达到的部分子空间的大小v2。公式(11)中的 β值是踝关节内翻的最大转动角度，即为y3。θ为踝关节外旋的转动角度，其范围为[0,y2]。
在计算内旋能达到而外翻不能达到的部分子空间的大小v3时，公式(11)中的β值是踝关节外翻的最大转动角度，即为y4。θ为踝关节内旋的转动角度，其范围为[0,y1]。
在计算内旋能达到而内翻不能达到的部分子空间的大小v4时，公式(11)中的β值是踝关节内翻的最大转动角度，即为y3。θ为踝关节内旋的转动角度，其范围为[0,y1]。
因此每个跖屈或背屈角度所对应的子空间大小为：
所以，v也是关于自变量(跖屈和背屈的转动角度)的变化函数。最后通过对v进行关于跖屈和背屈的积分就得到足踝三维运动空间的大小v3D：(xmin表示跖屈的最大角度，其值为负数；xmax表示背屈的最大角度，其值为正数)
月亮的诗3.1 实验器材
本研究采用光学式运动捕捉系统(NaturePoint公司)[19-20]。其基本原理是通过使用12个高分辨率红外摄像机跟踪表面涂有一层反光物质的小球，这些小球被称为标记点。在系统标定后，相机连续拍摄表演者的动作，并将图像序列保存下来，然后再进行分析和处理，识别其中的标记点，并计算其在每一瞬间的空间位置，进而得到其运动轨迹。这种光学式运动捕捉的优点是采样速率和精度较高；使用很方便；标记点的数量可根据实际情况任意添加；同时无电缆、机械装置的限制，表演者可以自由地表演。
3.2 试验方法
运动捕捉系统不能直接得到踝关节在各个方向上转动的角度，但能得到每个标记点在每一瞬间的空间位置。所以可以利用向量和三角函数将捕捉到的标记点的坐标值转化为踝关节的转动角度。此实验使用了三个标记点，分别记为mark1、mark2、mark3(图2)。
由于运动捕捉系统是通过标记点反光来得到它的三维坐标。在测量过程中，为了避免某个标记点可能被遮挡，以及脚的凹凸不平和表面皮肤的移动给测量结果带来的影响，本实验用绷带在人的脚底水平固定一个薄的长方形硬板(40×20 cm)，将三个标记点按照图2所描述的位置贴装在硬板上，贴装完成后mark1与mark2构成的直线与mark2与mark3构成的直线是垂直的。
为了避免测量过程中小腿随脚转动对测量数据带来影响，本实验选用了一个带有可调节铁架的医用床(图3a)，在测量过程中用绷带把小腿固定在铁架上，这样就能避免小腿随脚转动。而且可以根据被测者的身高调整铁架，避免铁架对数据测量带来影响。
为了测得足踝在各个方向上的旋转角度，需要有参考零点。如图3b所示，将人的小腿垂直于水平地面，脚底所在的平面和小腿轴线垂直，脚趾指向正前方的状态定义为初始状态，此时认为足踝各方向的旋转角度都为0°。由于小腿呈不规则圆锥状，所以当小腿被固定在铁架上的时候，小腿的轴线难以绝对垂直于地面。但在初始状态时要求脚底所在的平面垂直于小腿轴线，因此平行于脚底所在的平面固定的薄硬板与水平地面并非绝对平行，存在一定的初始误差。此误差大小与被测者的腿型和固定薄硬板时的精确度有关。经过统计分好像近义词>男女之间有纯洁的友谊吗
析，在初始状态下，不管相对于哪个转动方向，硬板与水平地面的夹角都小于4°。所以在后续的数据处理过程中，本研究将采集过程得到的角度减去在初始位置时的角度以消除初始误差带来的影响。
开始测量前，要求被测者坐在床上，将带有标记点的硬板固定在脚底，同时把小腿固定在铁架上，让受测者的脚保持在初始状态。开始测量后，要求被测者在改变跖屈或者背屈角度的时候分别做最大内旋、外旋、内翻和外翻运动，尽量得到每个跖屈或背屈角度所对应的最大内旋转动角度、最大外旋转动角度、最大内翻转动角度和最大外翻转动角度。如图3c所示，当测量右脚数据时候，人坐在铁架的右侧，这样当小腿束缚在铁架上时，踝关节更容易处于初始状态，使测量数据更加精确。同理，当测量左脚数据时让被测者坐在铁架的左侧，如图3d所示。每个人两只脚都要被测量。
