脑卒中运动功能评价模型建立方法、评价方法及系统与流程

1.本发明涉及康复工程、神经功能评估领域，本发明涉及脑卒中运动功能评价模型建立方法、评价方法及系统。


背景技术：

2.国际权威杂志《circulation》概述性报告指出，中国的脑卒中发病率呈逐年上升趋势，且势头迅猛。脑卒中仍为我国成年人致死和致残的首位原因，现有脑卒中患者1494万人，每年新发病例330万人。合理有效的康复评估策略对脑卒中患者功能恢复至关重要。研究表明，脑卒中致残的根本原因不是因为肢体本身受伤导致骨骼肌肉功能丧失，而是因为运动神经中枢系统受损伤没有得到有效治疗留下的肢体行为功能丧失，形成偏瘫症状。每年脑血管病死亡154万人，存活者中80％左右留有不同程度的功能障碍。
3.现有临床康复训练过程中缺乏对患者运动功能恢复及改善的个性化实时动态监测及状态评估，存在场地受限、形式单一、缺乏及时评价机制等问题，使得诸如虚拟康复、康复机器人等新型康复手段的应用普及受到严重制约和限制。
4.当前市场的评估系统主要从生理或运动信息的某一个层面进行运动功能评估，缺少从病理机制和多信息融合角度对脑卒中患者运动功能状态进行客观全面的数字评估，不能够满足日益增长的康复需求。


技术实现要素：

5.本发明提供了脑卒中运动功能评价模型建立方法、评价方法及系统，目的在于实现脑卒中后运动功能康复的社区化、普适化、个性化，实现对脑卒中后患者运动功能状态的有效评估。
6.为实现上述目的，本发明提供如下技术方案：
7.一方面，本发明提供了一种脑卒中运动功能评价模型建立方法，该方法包括特征指标筛选阶段和可视化模型构建阶段；
8.其中，特征指标筛选阶段包括：通过皮尔逊相关性分析进行上田敏临床定性特征指标与定量特征指标相关性分析确定其一致性；对不同的特性采用不同的方法进行筛选，得到最显著性最强的频段来进行评价指标计算；所述定量特征指标为体现神经肌肉功能耦合状态的指标，包括：局部频带特性指标、异频耦合特性指标、时延信息传递指数指标、复杂多尺度一致性指标；所述异频耦合特性是指从复杂性的角度刻画大脑与肌肉间潜在的非线性耦合机制，所述局部频带特性是指在特定频段内脑电信号和肌电信号间不同信息传递方向的非线性复杂性度量，所述时延信息传递指数能够表现信息交互的方向性和两信号间的耦合强度，所述复杂多尺度一致性能够描述两耦合信号之间的多尺度同步关系；
9.可视化模型构建阶段包括：基于所述特征指标筛选阶段筛选出的定量特征指标构建可视化运动功能评价模型，所述可视化运动功能评价模型通过ahp
‑
fce算法进行定性
‑
定量评价。
10.进一步地，对不同的特性采用不同的方法进行筛选，包括：
11.局部频带特性是指通过二维傅里叶变换计算多维拓展空间向量的频域表示，在此基础上对各频域点的维数向量进行传递熵计算，以此实现脑电信号、肌电信号不同频率值上的信息交互定量描述；局部频带特性指标采用的筛选方法包括：
12.分别计算出脑卒中患者和健康对照组的传递谱熵值，计算出emg到eeg和eeg到emg方向alpha、beta1、beta2、gamma1、gamma2及piper频段的tse显著面积指标，并对上述指标进行统计性分析，选出具有显著性差异的频段；通过皮尔逊相关性分析上述频段显著面积与上田敏临床评定量表等级的相关性，选出具有显著相关性关系的频段；
13.异频耦合特性是指通过小波包分解后获取多组不同时域尺度的子带信号，对各频率交叉的子带信号进行传递熵计算，以此实现不同时频尺度上两信号间信息交互的计算，通过基于异频耦合模型产生的异频耦合关系的仿真数据的分析，从仿真层面为神经肌肉功能耦合局部频带特性提供理论依据；异频耦合特性指标采用的筛选方法包括：
14.分别求取不同耦合比下基于小波包传递熵的delta、theta、alpha、beta及gamma频段显著面积，对每个耦合比下的显著面积进行统计性分析，选出具有显著性差异的频段；通过皮尔逊相关性分析不同耦合比下该频段的显著面积与上田敏临床评定量表等级的相关性，选出具有显著相关性关系的耦合比；
15.时间延迟是指神经肌肉系统包含的若干个子系统的信息交互和控制反馈之间存在的时间延迟；时延信息传递指数是指为了能够更加准确的量化描述两序列间的延迟时间而提出的一种通过对耦合信号进行排序重构以获取其排序模式，对排序模型的信号传递指数分析，从而实现对延迟时间的有效估计；
16.时延特性指标筛选方法包括：
17.经过对脑卒中患者和健康对照组emg到eeg和eeg到emg方向的延迟时间进行统计性分析，筛选出具有显著性差异的方向；通过皮尔逊相关性分析该方向的延迟时间与上田敏临床评定量表等级的相关性，出具有显著相关性关系的方向；
18.复杂多尺度一致性是指为了描述两耦合信号间的多尺度同步关系而提出的一种基于粗粒化过程进行的多尺度处理，通过多尺度一致性模型的显著面积指标，用以量化两耦合信号间不同尺度上的同步耦合关系；复杂多尺度一致性特性指标筛选方法包括：
19.分别计算出脑卒中患者和健康对照组1
‑
20尺度上的显著0
‑
60hz频段范围内的显著面积指标，并通过其变化趋势图分析，筛选出显著面积指标数值显著低于健康对照组的尺度；通过该尺度下显著面积指标与上田敏临床评定量表等级进行皮尔逊相关性分析，筛选处存在显著的正相关关系的尺度。
20.进一步地，选用eeg到emg方向beta2和gamma2频段的tse显著面积指标来体现局部频带特性；
21.选用耦合比为3：2及2：1条件下eeg到emg方向gamma频段的wpte显著面积指标来体现异频耦合特性；
22.选用eeg到emg方向的延迟时间来体现时间延迟特性；
23.选用尺度为9、10、12下的显著面积指标来体现复杂多尺度一致性特性。
24.进一步地，ahp
‑
fce算法中，确定因素集u为定量评估指标，即脑电特征值和肌电特征值，满足u＝{u1,u2}；
25.确定权重集w为定性评估指标，即根据上田敏量表得到的阶段值，满足w＝{w1,w2,w3,w4,w5}；其中，u1表示脑电特征值矩阵，u2表示肌电特征值矩阵；
26.综合评价结果表示为：
[0027][0028]
依据肌电评价指标u1表示为：
[0029][0030]
u2＝r6。
[0031]
又一方面，本发明还提供了一种基于上述脑卒中运动功能评价模型建立方法建立的脑卒中运动功能评价模型进行脑卒中运动功能评价的方法，该方法包括：
[0032]
对用户的脑电信号和肌电信号进行数据采集；
[0033]
对采集到的脑电信号和肌电信号进行信号预处理；
[0034]
对经过预处理的信号进行特征提取；特征提取为神经肌肉耦合指标提取，包括：局部频带特性、异频耦合特性、时延信息传递指数、复杂多尺度一致性四项特征值；所述异频耦合特性是指从复杂性的角度刻画大脑与肌肉间潜在的非线性耦合机制，所述局部频带特性是指在特定频段内脑电信号和肌电信号间不同信息传递方向的非线性复杂性度量，所述时延信息传递指数能够表现信息交互的方向性和两信号间的耦合强度，所述复杂多尺度一致性能够描述两耦合信号之间的多尺度同步关系；
[0035]
对提取出的特征进行特征层融合；
[0036]
将得到的融合特征输入到建立的可视化运动功能评价模型，得到用户的可视化康复评价结果。
[0037]
进一步地，采用64通道neuracle采集系统、trigno无线肌电系统和无线同步脉冲触发器进行数据采集。
[0038]
进一步地，对采集到的脑电信号进行预处理，包括：基于高通滤波器去除基线漂移，通过自适应滤波器去除50hz工频干扰；基于独立分量分析去除心电和眼电信号；基于典型相关性分析去除脑电信号中的肌电信号，筛选出cp3、c3、fc3三个电极的信号；
[0039]
对采集到的肌电信号进行预处理，包括：通过自适应陷波器去除50hz工频信号，基于0.5
‑
150hz的带通滤波器去除直流干扰，最后对肌电信号进行整流。
[0040]
进一步地，提取异频耦合特性包括：使用传递谱熵的分析方法，基于相空间重构方法将一维脑电信号构造成多维矩阵，对其进行二维傅里叶变化，在其频域上进行传递熵计算；
[0041]
提取局部频带特性包括：
[0042]
使用小波包传递熵的方法，基于代理数据对其进行统计分析，给出小波包传递熵方法的显著性指标用以量化两耦合信号间的异频耦合特性；
[0043]
提取时延信息传递指数，包括：
[0044]
序列x到序列y的时延信息传递指数为：
[0045][0046]
式中，x
′
和y
′
分别为x和y序列在t
‑
u时刻的嵌入维数为m
x
′
和m
y
′
的空间向量，空间向量x
′
和y
′
的第ψ
x
′
，1≤ψ
x
′
≤m
x
′
！和ψ
y
′
，1≤ψ
y
′
≤m
y
′
！种模式出现的次数分别为c
ψx
′
，1≤c
ψx
′
≤m
x
′
！和c
ψy
′
，1≤c
ψy
′
≤m
y
′
！，和分别为空间向量x
′
和y
′
的加权概率，为y在y
′
条件上的加权排序条件概率，为y在x
′
及y
′
共同条件上的加权排序条件概率；
[0047]
提取复杂多尺度一致性，包括：
[0048]
频率点f处的多尺度一致性和复杂多尺度一致性可以表示如下：
[0049][0050][0051]
式中，和分别为尺度s的时间序列和在频率点f处的自谱；和则分别
为时间序列和频率点f处的自谱；是时间序列和的频率点f处的互谱；是时间序列和频率点f处的互谱。
[0052]
进一步地，还包括：绘制不同康复阶段用户各项指标细化图。
[0053]
又一方面，本发明还提供了一种基于上述脑卒中运动功能评价模型建立方法建立的脑卒中运动功能评价模型进行脑卒中运动功能评价的系统，所述系统包括：
[0054]
数据采集模块，用于对用户的脑电信号和肌电信号进行数据采集；
[0055]
数据预处理模块，用于对所述数据采集模块采集到的脑电信号和肌电信号进行信号预处理；
[0056]
特征提取模块，用于对经过所述数据预处理模块预处理的信号进行特征提取；特征提取为神经肌肉耦合指标提取，包括：局部频带特性、异频耦合特性、时延信息传递指数、复杂多尺度一致性四项特征值；所述异频耦合特性是指从复杂性的角度刻画大脑与肌肉间潜在的非线性耦合机制，所述局部频带特性是指在特定频段内脑电信号和肌电信号间不同信息传递方向的非线性复杂性度量，所述时延信息传递指数能够表现信息交互的方向性和两信号间的耦合强度，所述复杂多尺度一致性能够描述两耦合信号之间的多尺度同步关系；
[0057]
特征融合模块，用于对所述特征提取模块提取出的特征进行特征层融合；
[0058]
评价模块，用于将所述特征融合模块得到的融合特征输入到建立的可视化运动功能评价模型，得到用户的可视化康复评价报告。
[0059]
本发明技术方案产生的积极效果如下：
[0060]
1、本发明从脑电、肌电、脑肌间三方面采集到的信息进行处理，从定性和定量综合得到患者的神经功能恢复程度，通过不同的特征指标分析处理，得到评估结果。从病理机制和多信息融合角度对脑卒中患者运动功能状态进行客观全面的数字评估，使得到的数据更科学准确。
[0061]
2、本发明能够实现康复评估的家庭化、社区化，通过将医生指导数字化，实现患者自主评估，降低了医生主观判断失误带来的评估结果差异。
附图说明
[0062]
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图做以简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
[0063]
图1为本发明实施例中一种脑卒中运动功能评价模型建立方法的流程图；
[0064]
图2为本发明实施例中一种脑卒中运动功能评价方法的流程图；
[0065]
图3为本发明实施例中不同康复阶段患者各项指标细化图(左)、患者康复训练效果评估图(右)；
[0066]
图4为本发明实施例中康复周期5和康复周期6的数据分析图；
[0067]
图5为本发明实施例中一种脑卒中运动功能评价系统的结构框图。
具体实施方式
[0068]
为了使本技术领域的人员更好地理解本发明方案，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分的实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本发明保护的范围。
[0069]
需要说明的是，本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换，以便这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外，术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
[0070]
参见图1，其示出了本发明实施例中一种脑卒中运动功能评价模型建立方法的流程图，该方法包括特征指标筛选阶段和可视化模型构建阶段；
[0071]
其中，特征指标筛选阶段包括：多指标融合筛选，得到最显著性最强的频段来进行评价指标计算。
[0072]
多指标融合筛选是指首先通过皮尔逊相关性分析进行上田敏临床定性特征指标与定量特征指标相关性分析确定其一致性，其次对不同的特性采用了不同的方法进行筛选，得到最显著性最强的频段来进行评价指标计算
[0073]
本发明实施例中选取分别体现时延特性、多尺度耦合特性、异频耦合特性、局部频带特性的时延信息传递指数、复杂多尺度一致性、小波包传递熵、传递谱熵的多维特征融合矩阵作为定量评价指标。
[0074]
其中，局部频带特性是指通过二维傅里叶变换计算多维拓展空间向量的频域表示，在此基础上对各频域点的维数向量进行传递熵计算，以此实现两信号即脑电信号、肌电信号不同频率值上的信息交互定量描述。局部频带特性指标采用的筛选方法包括：首先分别计算出脑卒中患者和健康对照组的传递谱熵值(tse)，计算出emg到eeg和eeg到emg方向alpha(8
‑
15hz)、beta1(15
‑
25hz)、beta2(25
‑
35hz)、gamma1(35
‑
45hz)、gamma2(45
‑
60hz)及piper(>60hz)频段的tse显著面积指标，并对上述指标进行统计性分析，选出具有显著性差异的频段；其次通过皮尔逊相关性分析上述频段显著面积与上田敏临床评定量表等级的相关性，选出具有显著相关性关系的频段。筛选结果为：选用eeg到emg方向beta2和gamma2频段的tse显著面积指标来体现局部频带特性。
[0075]
异频耦合特性是指通过小波包分解后获取多组不同时域尺度的子带信号，对各频率交叉的子带信号进行传递熵计算，以此实现不同时频尺度上两信号间信息交互的计算，通过基于异频耦合模型产生的异频耦合关系的仿真数据的分析，从仿真层面为神经肌肉功能耦合局部频带特性提供理论依据。异频耦合特性指标采用的筛选方法包括：首先分别求
取不同耦合比下基于小波包传递熵(wpte)的delta、theta、alpha、beta及gamma频段显著面积，对每个耦合比下的显著面积进行统计性分析，选出具有显著性差异的频段；其次通过皮尔逊相关性分析不同耦合比下该频段的显著面积与上田敏临床评定量表等级的相关性，选出具有显著相关性关系的耦合比。筛选结果为：选用耦合比为3：2及2：1条件下eeg到emg方向gamma频段的wpte显著面积指标来体现异频耦合特性。
[0076]
时间延迟是指神经肌肉系统包含的若干个子系统的信息交互和控制反馈之间存在的时间延迟。时延信息传递指数是指为了能够更加准确的量化描述两序列间的延迟时间而提出的一种通过对耦合信号进行排序重构以获取其排序模式，对排序模型的信号传递指数分析，从而实现对延迟时间的有效估计。时延特性指标筛选方法包括：经过对脑卒中患者和健康对照组emg到eeg和eeg到emg方向的延迟时间进行统计性分析，筛选出具有显著性差异的方向；其次通过皮尔逊相关性分析该方向的延迟时间与上田敏临床评定量表等级的相关性，出具有显著相关性关系的方向。筛选结果为：选用eeg到emg方向的延迟时间来体现时间延迟特性。
[0077]
复杂多尺度一致性是指为了描述两耦合信号间的多尺度同步关系而提出的一种基于粗粒化过程进行的多尺度处理，通过多尺度一致性模型的显著面积指标，用以量化两耦合信号间不同尺度上的同步耦合关系。复杂多尺度一致性特性指标筛选方法包括：分别计算出脑卒中患者和健康对照组1
‑
20尺度上的显著0
‑
60hz频段范围内的显著面积指标，并通过其变化趋势图分析，筛选出显著面积指标数值显著低于健康对照组的尺度；其次通过该尺度下显著面积指标与上田敏临床评定量表等级进行皮尔逊相关性分析，筛选处存在显著的正相关关系的尺度。筛选结果为：选用尺度为9、10、12下的显著面积指标来体现复杂多尺度一致性特性。
[0078]
可视化模型构建阶段是指基于特征指标筛选阶段筛选出的定量特征指标构建可视化运动功能评价模型。
[0079]
可视化运动功能评价模型中，通过ahp
‑
fce法进行定性
‑
定量评价，ahp
‑
fce法是一种定性和定量结合的评价方法，ahp法确定指标权重，fce法结合ahp法的多层次指标权重来决定评价结果，能全面准确地评价目标问题实际状况，能够解决传统临床评价方法受医师主观因素影响、评价精准度不高等缺点。
[0080]
ahp
‑
fce算法中，确定因素集u为定量评估指标，即脑电特征值和肌电特征值,满足u＝{u1,u2}；确定权重集w为定性评估指标，即根据上田敏量表得到的阶段值，满足w＝{w1,w2,w3,w4,w5}。其中，u1表示脑电特征值矩阵，u2表示肌电特征值矩阵；
[0081]
综合评价结果表示为：
[0082][0083]
依据肌电评价指标u1可表示为：
[0084][0085]
u2＝r6。
[0086]
上述实施例中，基于上田敏临床定性特征指标与神经肌肉耦合状态定量特征指标建立的可视化运动功能评价模型，综合描述患者在不同临床康复阶段下的神经肌肉功能耦合特征指标及其变化规律，实现了对患者运动功能的多层次定性定量评价的可视化综合体现。
[0087]
如图2所示，其示出了本实施例中利用上述实施例中建立的可视化运动功能评价模型进行脑卒中运动功能评价方法的流程图，该方法包括以下三个部分：
[0088]
s101、数据采集；
[0089]
数据采集包括采用64通道neuracle采集系统、trigno无线肌电系统和无线同步脉冲触发器进行脑电信号和肌电信号采集。
[0090]
s102、数据处理；
[0091]
数据处理包括三个步骤：
[0092]
s201、对采集到的脑电信号和肌电信号进行信号预处理；
[0093]
其中，信号包括脑电信号和肌电信号。
[0094]
脑电信号预处理包括：首先基于高通滤波器去除基线漂移，通过自适应滤波器去除50hz工频干扰，其次基于独立分量分析去除心电和眼电信号，最后基于典型相关性分析去除脑电信号中的肌电信号，筛选出cp3、c3、fc3三个电极的信号，因为该通道处的神经肌肉功能耦合最显著。
[0095]
肌电信号预处理包括：通过自适应陷波器去除50hz工频信号，基于0.5
‑
150hz的带通滤波器去除直流干扰，最后对肌电信号进行整流。
[0096]
s202、对经过信号预处理的信号进行特征提取；
[0097]
特征提取为神经肌肉耦合指标提取，包括：局部频带特性、异频耦合特性、时延信息传递指数、复杂多尺度一致性四项特征值。
[0098]
其中，异频耦合特性是指从复杂性的角度刻画大脑与肌肉间潜在的非线性耦合机制。本发明实施例中使用传递谱熵的分析方法，基于相空间重构方法将一维脑电信号构造成多维矩阵，然后对其进行二维傅里叶变化，在其频域上进行传递熵计算。将原始序列x和t拓展到m维的空间，得到维数为m
×
n的矩阵s和v；
[0099]
[0100][0101]
式中，m及n均为整数且1≤m≤m,1≤n≤n；
[0102]
对空间向量s和v进行二维傅里叶变换，设定其采样频率为f
s
，则其变换为：
[0103][0104][0105]
式中，参数k为整数且k＝1,2,
…
,m，表示向量空间s和v中的不同方向点，即维数点；
[0106]
二维向量矩阵h(k,f)和g(k,f)表示为：
[0107]
w(f)＝{h(1,f),h(2,f),
…
,h(k,f),
…
,h(m,f)}；
[0108]
v(f)＝{g(1,f),g(2,f),
…
,g(k,f),
…
,g(m,f)}；
[0109]
则传递熵公式表示为：
[0110][0111]
式中，f是离散频率值，τ是频率预测值，和分别是向量w(f)和v(f)在频率点f处的φ和维的延时向量，p表示概率计算。
[0112]
局部频带特性是指通过对耦合信号进行小波包分解获取多个时频局部频带信号，并对每个局部频带信号进行传递熵的计算，以此探究在特定频段内脑电信号和肌电信号间不同信息传递方向的非线性复杂性度量。本发明实施例中使用小波包传递熵的方法，基于代理数据对其进行统计分析，给出小波包传递熵方法的显著性指标用以量化两耦合信号间的异频耦合特性。小波包传递熵用下式表示：
[0113][0114]
式中t是离散时间变量，u是离散的预测时间变量；和分别是序列和的m
‑
和n
‑
维的时间延迟空间向量，和分别为信号x在第i层的第n个子空间的子频带信号和y信号在第j层的第m个子空间的子频带信号，参数i,j∈z；
[0115]
特定频带内的显著面积指标为：
[0116][0117]
式中，δf1和δf2为分辨率，分辨率越小则表明wpte所描述的fcmc特性的细节信息越多，和分别是x信号频率范围的最大值和最小值，且参数r∈z，且
[0118]
时延信息传递指数不仅能够用于两信号间信息交互的方向性，更能描述两信号间的耦合强化。序列x到序列y的时延信息传递指数为：
[0119][0120]
式中，x
′
和y
′
分别为x和y序列在t
‑
u时刻的嵌入维数为m
x
′
和m
y
′
的空间向量，空间向量x
′
和y
′
的第ψ
x
′
，1≤ψ
x
′
≤m
x
′
！和ψ
y
′
，1≤ψ
y
′
≤m
y
′
！种模式出现的次数分别为c
ψx
′
，1≤c
ψx
′
≤m
x
′
！和c
ψy
′
，1≤c
ψy
′
≤m
y
′
！，和分别为空间向量x
′
和y
′
的加权概率，为y在y
′
条件上的加权排序条件概率，为y在x
′
及y
′
共同条件上的加权排序条件概率。
[0121]
复杂多尺度一致性是指通过引入k阶粗粒化处理，用于一致性方法分析中，来描述两耦合信号之间的多尺度同步关系。频率点f处的多尺度一致性和复杂多尺度一致性可以表示如下：
[0122][0123][0124]
式中，和分别为尺度s的时间序列和在频率点f处的自谱；和则分别为时间序列和频率点f处的自谱；是时间序列
和的频率点f处的互谱；是时间序列和频率点f处的互谱。
[0125]
s203、对提取出的特征进行特征层融合；
[0126]
特征层融合是指一种在信号特征层面进行融合的方法，先对各类信号源分别完成特征提取工作，然后把各类特征进行融合，对融合后的数据进行模式识别。特征融合是一种次级融合方式，相对于数据层融合方法对原始数据进行了特征提取，使得数据特征更加明显且降低了数据量提高处理速度。在进行特征融合时对于原始信号而言保留了信号的关键信息同时去掉了冗余信息，减少了信号处理的复杂度又最大程度保证了信号的真实度。
[0127]
s103、可视化显示；
[0128]
可视化显示是指将得到的患者融合特征输入到上述实施例中建立的可视化运动功能评价模型，得到患者的可视化康复评价结果。
[0129]
绘制不同康复阶段患者各项指标细化图能够直接准确的展示了各阶段患者的详细指标数值，包含筛选出的tse、wpte、tdit、cmsc的特定频段指标值。图3(左)能够直观看出，随着上田敏评定等级增加，各项指标呈不同程度的线性关系。
[0130]
为直观显示评价结果和康复训练的效果，建立基于康复阶段、康复周期、评价指标的患者康复训练效果评价图，构建多元可视化定性定量融合评估模型，综合描述患者在一段康复训练阶段内的康复效果。康复阶段是指上田敏临床评定量表等级，康复周期是指患者经历的康复训练时间，本专利中一周期指代十天。图3(右)中，根据康复阶段图例显示不同的康复阶段对应不同深度的颜色，即颜色越浅康复阶段值越大，对应的上田敏临床评定量表的等级越高，患者的康复效果就越好。该图能够直观的显示患者康复训练的效果对比。
[0131]
如图4所示，根据患者在康复周期5和康复周期6的数据进行分析，数据的间隔为一个周期(十天)，绘制条形图进行比对可以看出各项指标均有明显的差异，随患者康复训练过程的进行，gamma2频段的tse值和时延指标数值减小，其他指标均增加，与不同康复阶段患者各项指标细化图所表现的整体趋势保持一致。
[0132]
上述实施例中的脑卒中运动功能评价方法，从脑电、肌电、脑肌间三方面采集到的信息进行处理，从定性和定量综合得到患者的神经功能恢复程度，通过不同的特征指标分析处理，得到评估结果。从病理机制和多信息融合角度对脑卒中患者运动功能状态进行客观全面的数字评估，使得到的数据更科学准确。同时，本发明实施例中的脑卒中运动功能评价方法能够实现康复评估的家庭化、社区化，通过将医生指导数字化，实现患者自主评估，降低了医生主观判断失误带来的评估结果差异。
[0133]
对应于上述脑卒中运动功能评价方法，本发明实施例中还提供了一种脑卒中运动功能评价系统，如图5所示，其示出了本实施例中利用上述实施例中建立的可视化运动功能评价模型进行脑卒中运动功能评价的系统框图，该系统包括：
[0134]
数据采集模块301，用于对用户的脑电信号和肌电信号进行数据采集；
[0135]
数据预处理模块302，用于对所述数据采集模块采集到的脑电信号和肌电信号进行信号预处理；
[0136]
特征提取模块303，用于对经过所述数据预处理模块预处理的信号进行特征提取；脑电信号特征提取为神经肌肉耦合指标提取，包括：局部频带特性、异频耦合特性、时延信
息传递指数、复杂多尺度一致性四项特征值；所述异频耦合特性是指从复杂性的角度刻画大脑与肌肉间潜在的非线性耦合机制，所述局部频带特性是指在特定频段内脑电信号和肌电信号间不同信息传递方向的非线性复杂性度量，所述时延信息传递指数能够表现信息交互的方向性和两信号间的耦合强度，所述复杂多尺度一致性能够描述两耦合信号之间的多尺度同步关系；
[0137]
特征融合模块304，用于对所述特征提取模块提取出的特征进行特征层融合；
[0138]
评价模块305，用于将所述特征融合模块得到的融合特征输入到建立的可视化运动功能评价模型，得到用户的可视化康复评价报告。
[0139]
在所有评估结束后，系统将自动生成康复评估报报告，自行选择下载电子版或打印，所有流程结束，关闭系统。
[0140]
对于本发明实施例的脑卒中运动功能评价系统而言，由于其与上面实施例中的脑卒中运动功能评价方法相对应，所以描述的比较简单，相关相似之处请参见上面实施例中部分的说明即可，此处不再详述。
[0141]
在本技术所提供的几个实施例中，应该理解到，所揭露的技术内容，可通过其它的方式实现。其中，以上所描述的装置实施例仅仅是示意性的，例如所述单元的划分，可以为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，单元或模块的间接耦合或通信连接，可以是电性或其它的形式。
[0142]
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
[0143]
另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。
[0144]
所述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可为个人计算机、服务器或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：u盘、只读存储器(rom，read
‑
only memory)、随机存取存储器(ram，random access memory)、移动硬盘、磁碟或者光盘等各种可以存储程序代码的介质。
[0145]
最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。