脑卒中运动功能评估模型构建方法和运动功能评估方法与流程

1.本发明涉及康复医学领域，特别涉及脑卒中患者的运动功能评估模型构建方法和运动功能评估方法。


背景技术：

2.脑卒中又称脑中风，是一种致残率非常高的疾病，超过半数的患者发病半年后仍有上肢功能障碍，80％的患者留有残疾。对脑卒中患者运动功能状态的准确评估，是制定有效脑卒中康复策略的重要依据。临床常用医学运动量表评估脑卒中患者的运动功能，基于医学运动量表的康复评估一般是通过康复医师的眼睛观察患者的直观动作，对患者的运动功能和康复效果进行评分，评分结果的客观性难以保障，评分结果的灵敏度和可靠性较低。除此之外，患者在康复训练过程中身体状态的微小变化难以被察觉，从而使得康复效果难以即时准确的评价。
3.为了提高评估的准确性，目前不少研究者以运动中的外周的生物力学、肌电等作为生理指标评估脑卒中患者的运动功能，这类型的评估方法要求患者完成指定的动作以获取相应数据，对于部分运动能力较弱的患者，某些必要动作无法完成，若采用代偿动作则存在无法真正评估患者恢复程度的问题。
4.可见，现有技术还有待改进和提高。


技术实现要素：

5.鉴于上述现有技术的不足之处，本发明的目的在于提供一种脑卒中患者的运动功能评估模型构建方法和运动功能评估方法，旨在解决现有技术中的评估方法实用性有限的问题。
6.为了达到上述目的，本发明采取了以下技术方案：
7.一种脑卒中患者的运动功能评估模型的构建方法，其中，包括以下步骤：
8.步骤s101.构建数据库：通过多通道检测设备采集多名脑卒中患者额叶或运动区的静息态脑血氧信号数据，组成基础数据库；
9.步骤s102.数据预处理：对基础数据库中采集的静息态脑血氧信号数据进行带通滤波去除漂移噪声和生理干扰；再采用样条插值的方法对伪迹进行矫正；
10.步骤s103.特征提取：采用pearson相关性计算公式计算测量区域内两两通道之间的脑功能连接，然后对得到的相关性矩阵进行fisher’s z变换；采用lasso回归的方式在fisher’s z分数矩阵中选择最优的分类特征子集；
11.步骤s104.构建支持向量机模型：根据步骤s103中获得的最优的分类特征子集，构建基于支持向量机的脑卒中患者的运动功能评估模型。
12.所述的脑卒中患者的运动功能评估模型的构建方法，其中，所述步骤s101中通过近红外脑功能检测设备进行多通道采集脑卒中患者额叶或运动区的静息态脑血氧信号数据。
13.所述的脑卒中患者的运动功能评估模型的构建方法，其中，所述近红外脑功能检测设备包括多个光源和多个探测器；多个所述光源和探测器交错设置，每个所述光源和探测器之间形成一个测量通道。
14.所述的脑卒中患者的运动功能评估模型的构建方法，其中，所述静息态脑血氧信号数据为脑卒中患者保持休息状态，闭上眼睛并保持清醒时采集得到的脑血氧信号数据。
15.所述的脑卒中患者的运动功能评估模型的构建方法，其中，所述步骤s102还包括空间滤波，所述步骤s102中基础数据库中采集的静息态脑血氧信号数据先通过空间滤波以去除表层干扰，然后通过带通滤波去除漂移噪声和生理干扰；再采用样条插值的方法对伪迹进行矫正。
16.所述的脑卒中患者的运动功能评估模型的构建方法，其中，所述空间滤波的操作方法为：将采集的静息态脑血氧信号数据减去所有测量通道的平均值，得到去除表层干扰后的静息态脑血氧信号。
17.所述的脑卒中患者的运动功能评估模型的构建方法，其中，所述步骤s102中所述带通滤波的频率范围为0.01
‑
0.1hz。
18.所述的脑卒中患者的运动功能评估模型的构建方法，其中，所述步骤s102中采用三次样条插值的方法对伪迹进行矫正，包括以下过程：
19.以基线作为判断标准，以基线的三倍标准差作为阀值；
20.检测得到每个伪迹段的起点和终点，得出伪迹段的数量，遍历所有伪迹段并对每个伪迹段进行三次平滑样条插值，然后用每个伪迹段减去每个伪迹段插值后的信号；
21.对处理后的信号进行整体重构，获得重构的信号序列。
22.所述的脑卒中患者的运动功能评估模型的构建方法，其中，所述步骤s104中，所述支持向量模型的惩罚系数和核函数采用网格搜索分方式进行优化。
23.一种脑卒中患者的运动功能评估方法，包括以下步骤：
24.步骤s201.数据采集：使用近红外脑功能检测设备进行多通道采集脑卒中患者的静息态脑血氧信号数据；
25.步骤s202.数据预处理：对采集的静息态脑血氧信号数据通过空间滤波去除表层干扰，然后通过带通滤波去除漂移噪声和生理干扰；最后采用样条插值的方法对伪迹进行矫正；
26.步骤s203.特征提取：采用pearson相关性计算公式计算测量区域内两两通道之间的脑功能连接，然后对得到的相关性矩阵进行fisher’s z变换；采用lasso回归的方式在fisher’s z分数矩阵中选择最优的分类特征子集；
27.步骤s204.运动功能等级评估：采用上述方法构建的脑卒中患者的运动功能评估模型对脑卒中患者进行运动功能等级评估。
28.有益效果：
29.本发明提供了一种脑卒中患者的运动功能评估模型的构建方法。模型的构建先采集多名脑卒中患者额叶或运动区的静息态脑血氧信号数据构建数据库，然后将数据经过预处理后，通过回归的方法选取脑卒中患者静息状态下脑连接信息作为特征，构建机器学习的模型，可以对脑卒中患者不同运动功能状态进行分类预测。另外，本发明还提供基于上述模型的运动功能评估方法，对于运动能力较差的患者同样可以适用，也可以用于对患者康
复训练后即时效果的评估。本发明评估方法的检测操作简单，可以排除主观因素影响，且干扰因素少，准确性高，可反复多次测量。
附图说明
30.图1为构建脑卒中患者的运动功能评估模型的流程图。
31.图2为步骤s103特征提取的流程图。
32.图3为一种改进方式中步骤s102数据预处理的流程图。
33.图4为脑卒中患者的运动功能评估的流程图。
34.图5为近红外脑功能检测设备的光源和探测器的分布示意图。
具体实施方式
35.本发明提供一种运动功能评估模型构建方法和运动功能评估方法，为使本发明的目的、技术方案及效果更加清楚、明确，以下参照附图并举实施例对本发明进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
36.本发明的描述中，除非另有说明，“多个”的含义是两个或两个以上。
37.请参阅图1、图2，脑卒中患者的运动功能评估模型的构建方法，包括以下步骤：
38.步骤s101.构建数据库：通过多通道检测设备采集多个脑卒中患者额叶或运动区的静息态脑血氧信号数据，组成基础数据库；
39.步骤s102.数据预处理：对基础数据库中采集的静息态脑血氧信号数据进行带通滤波去除漂移噪声和生理干扰；再采用样条插值的方法对伪迹进行矫正，减少伪迹对计算静态功能连接的影响；
40.步骤s103.特征提取：s103a.采用pearson相关性计算公式计算测量区域内两两通道之间的脑功能连接，然后对得到的相关性矩阵进行fisher’sz变换，以改善数据的正态性，变换公式为z＝0.5
×
(ln(r 1)
‑
ln(r
‑
1))；s103b.采用lasso回归的方式在fisher’s z分数矩阵中选择最优的分类特征子集；
41.具体的，由于fisher’s z分数矩阵呈对称，因此，可以将fisher’sz分数矩阵的下三角或者上三角部分作为分类特征，再通过lasso回归的方式选择最优的分类特征子集；对于采用lasso进行特征的选择，可以根据最终lasso回归的系数的绝对值的大小作为衡量特征的贡献，剔除系数为0的特征和贡献度较小的特征；
42.步骤s104.构建支持向量机模型(svm)：根据步骤s103中获得的最优的分类特征子集，构建基于支持向量机的脑卒中患者的运动功能评估模型。
43.本实施例中，通过近红外脑功能检测设备进行多通道采集脑卒中患者额叶区域的静息态脑血氧信号数据。近红外脑功能检测设备可通过氧合血红蛋白和脱氧血红蛋白对不同波长的近红外光吸收系数的变化来测量两种的浓度变化。而本发明利用脑卒中患者的氧合血红蛋白浓度变化来构建静态功能连接。
44.本发明选择额叶或运动区作为检测区域。其中大脑的额叶主要与人体的记忆、思维、语言和运动等功能相关，额叶的损伤会导致无法集中精力观察和思考，更无法进行仔细的逻辑推理，对突发事件无能为力，健忘，行为反应迟钝，出现情绪波动，因此，以额叶作为检测区得到的数据可以反映患者的反应能力和运动能力。运动区为大脑运动中枢，与躯体
各部位的肌肉调节支配相关，检测得到的数据与运动功能的相关性更强。
45.所述近红外脑功能检测设备包括多个光源和多个探测器；多个所述光源和探测器交错设置，每个所述光源和探测器之间形成一个测量通道。请参阅图5，本实施例中，所述近红外脑功能检测设备包括8个光源和7个探测器，光源和探测器放置在患者的前额叶区，并按图5进行分布，图中圆形表示光源，矩形表示探测器，光源和检测器之间共形成22个测量通道。具体的每个光源和探测器之间的距离为3cm。近红外脑功能检测设备的采样率为10hz，采样时间为6分钟。在静息态脑血氧信号数据采集期间，脑卒中患者保持休息状态坐在椅子上，闭上眼睛并保持清醒，在检测期间尽量减少身体运动，以保证检测数据的稳定性和可靠性。
46.请参阅图3，作为一种改进，所述步骤s102还包括空间滤波，所述步骤s102中基础数据库中采集的静息态脑血氧信号数据预处理包括：s102a.先通过空间滤波以去除表层干扰；s102b.然后通过带通滤波去除漂移噪声和生理干扰；s102c.再采用样条插值的方法对伪迹进行矫正。
47.具体的，所述空间滤波的操作方法为：将采集的静息态脑血氧信号数据减去所有测量通道的平均值，得到去除表层干扰后的静息态脑血氧信号。
48.作为一种改进，所述步骤s102中所述带通滤波的频率范围为0.01
‑
0.1hz。静息状态是大脑的一种有组织的基线默认模式，在静息态中存在自发的特定低频震荡，这些自发的低频震荡的频率范围在0.01
‑
0.1hz之间。
49.作为一种改进，所述步骤s102中所述样条插值为三次样条插值。具体的，离散时间序列x(t)＝{x(t
i
)}，t
i
＝iδt，i＝1,2，...，n，其中δt为采样周期，n为采样点数。
50.三次样条插值的操作流程如下：
51.1.采用移动标准差的方式检测运动伪迹
52.1.1以基线作为判断标准，以基线的三倍标准差作为阀值t；
53.1.2滑动窗口的大小为w，w的取值为201(其中w可以根据实际的需求进行调整)，以1个单位的点数作为滑动步长，计算每个时间窗的标准差；把存储每个时间窗标准差的一维矩阵d的前k＝(w
‑
1)/2个元素置0，使每个时间窗内的标准差放在时间窗的中点的位置；
54.1.3创建一个长度与一维矩阵d一致的一维矩阵r，将一维矩阵d中大于阀值t的标准差的位置对应在一维矩阵r中置1，否则置0；
55.2.检测每个伪迹段的起点和终点
56.2.1遍历矩阵r并将其前后相邻两元素相减得出差值(后一个元素减前一个元素)，如果前后相邻两元素的差值为1则为起点，如果差值为
‑
1则为终点。然后分别用一个一维矩阵lms存储起点的位置，用一个一维矩阵lmf存储终点的位置。如果终点为空则将数据的最后一个点的位置作为终点，如果起点为空，则数据的第一个点作为起点。如果出现第一个起点位置大于第一个终点的位置则将1插入lms的第一个点的前一个位置即[1,lms]。如果最后一个起点大于最后一个终点则将数据的最后一个点的后一个位置插入1，即[lmf,1]；
[0057]
2.2根据lms和lmf可以确定信号存在的伪迹段的数量，遍历所有伪迹段并对每个伪迹段进行三次平滑样条插值，然后用每个伪迹段减去每个伪迹段插值后的信号；
[0058]
3.对处理后的信号进行整体重构，获得重构的信号序列。
[0059]
一种实施方式中，所述步骤s104中，所述支持向量模型的惩罚系数和核函数采用
网格搜索分方式进行优化。
[0060]
实施例
[0061]
随机选择21名脑卒中患者采用本发明构建的模型进行运动功能评估。采用本发明的评估方法前，这21名患者根据现有的运动评估量表评估后分为三组，分别为：运动能力较差组、运动能力中等组和运动能力较好组，该评估结果用于检验本发明的脑卒中患者的运动功能评估方法的正确率。
[0062]
请参阅图4，本实施例的脑卒中患者的运动功能评估方法，包括以下步骤：
[0063]
步骤s201.数据采集：上述21名脑卒中患者分别采用的近红外脑功能检测设备进行多通道采集；采集的区域为患者的前额叶区，近红外脑功能检测设备的光源和探测器按照图5所示进行分布；每个光源和探测器之间的距离为3cm。近红外脑功能检测设备的采样率为10hz，采样时间为6分钟。静息态脑血氧信号数据采集期间，脑卒中患者保持休息状态坐在椅子上，闭上眼睛并保持清醒，在检测期间尽量减少身体运动。
[0064]
步骤s202.数据预处理：对采集的静息态脑血氧信号数据通过空间滤波去除表层干扰；然后通过带通滤波去除漂移噪声和生理干扰，带通滤波的频率范围为0.01
‑
0.1hz；最后采用三次样条插值的方法对伪迹进行矫正；
[0065]
步骤s203.特征提取：采用pearson相关性计算公式计算测量区域内两两通道之间的脑功能连接，然后对得到的相关性矩阵进行fisher’s z变换；采用lasso回归的方式在fisher’s z分数矩阵中选择最优的分类特征子集；
[0066]
步骤s204.运动功能等级评估：采用本发明构建的脑卒中患者的运动功能评估模型对脑卒中患者进行运动功能等级评估。
[0067]
下表是本发明的运动功能评估方法对脑卒中患者的运动功能等级评估分类的正确率结果：
[0068]
[0069][0070]
从上述结果可以看出，本发明构建的脑卒中患者的运动功能评估方法正确率高达95.60％，可用于脑卒中患者的运动功能评估。相比现有的评估方法，本发明的运动功能评估方法不受主观因素影响，适用范围广，无需患者完成特定动作，适用于行动不便或因其他原因导致有行动障碍的患者。所采集的数据统一性和可比较性强，更有利于准确客观地评估脑卒中患者的运动功能，评估方法中所选取的特征也可以用于评估患者康复训练中或者训练后即时效果。此外，随着数据库的不断扩充，本发明所构建的模型准确性会进一步地得到提升。
[0071]
本发明中，对于评估脑卒中患者的运动功能等级的分类标准还可以根据实际需要作进一步替换或改变，如对脑卒中患者的运动功能等级可以根据fugl
‑
meyer评定量表的分值进行进一步细化划分或者是根据评定量表brunnstrom的功能分期来划分脑卒中患者的运动功能水平。
[0072]
可以理解的是，对本领域普通技术人员来说，可以根据本发明的技术方案及其发明构思加以等同替换或改变，而所有这些改变或替换都应属于本发明所附的权利要求的保护范围。