本发明专利技术提出了一种基于改进的微状态分析的运动想象脑电信号分类方法，考虑到运动想象这一特定任务状态下微状态分析步骤是否可以优化的问题。首先计算预处理后的多通道脑电信号全局场功率，去除数值较高和较低的点，将剩余局部峰值点时刻对应的电势地形图进行聚类分析；将聚类得到的四种微状态分别与原始脑电信号的电势地形图进行空间相关性计算，构建出微状态时间序列，并计算出每个微状态对应的持续时间、出现频率等参数；将训练集数据的微状态时间序列参数作为特征输入到SVM分类器中训练出分类模型，对测试集特征进行分类得到测试结果。本发明专利技术创新性地将微状态分析法应用在运动想象脑电信号上以达到分类的目的，并取得了较好的分类结果。较好的分类结果。
【技术实现步骤摘要】
基于改进的微状态分析的运动想象脑电信号分类方法
[0001]本专利技术属于脑科学、机器学习交叉研究领域。
技术介绍
[0002]脑
‑
机接口(Brain Computer Interface，BCI)技术在人脑思维与外部设备之间构建了直接的通路，使得大脑可以在无需神经、肌肉等的参与下便可直接与外界进行通讯，对残疾人、老年人的生活质量的提高有重大意义。脑电图(Electroencephalogram，EEG)是目前检测大脑功能活动的常用手段，具有操作简便、时间分辨率高、记录无损等优点。由运动想象任务产生的脑电信号无需外界刺激，不易引起被试疲劳，容易被残障人群所接受，因此如何识别出运动想象脑电信号的特征对BCI的实现尤为重要。
[0003]1987年，Lehmann等人证明了多通道静息状态脑电信号的α频带可以解析为有限数量的不同的准稳态，这种准稳态就是微状态。微状态由记录在头皮上的多通道阵列中的电势地形拓扑图定义，具有以下两个特点。第一，一个微观状态在快速转变为另一个不同微观状态之前大约保持80
‑
120ms的稳定状态。第二，尽管在多通道脑电信号中存在许多类型的电势地形图拓扑图，但是大部分信号可以由少量电势地形图来表示。值得注意的是，通过交叉验证准则，得出微状态的最佳数目为四，且四种微状态可以代表所有数据70％～80％的特征。
[0004]提取脑电信号在不同运动想象任务中存在的差异，并将这些差异作为脑电分类的依据便是脑功能特征提取中关键的一步。传统的脑电特征提取方法往往是基于少数电极的脑电数据进行的，存在忽略时空信息、分类准确率低等问题，而脑电信号微状态分析法是通过考虑来自所有电极的信号来建立功能状态的全局表示，解决了传统方法很难处理多通道脑电信号的问题，并且在保留时空信息的基础上对脑电信号进行分析处理。目前已经存在将脑电信号微状态分析法应用到运动想象任务中的例子，但只是按照经典的微状态分析算法步骤进行处理，并没有考虑到运动想象这一特定任务状态下EEG微状态分析步骤是否可以优化的问题。
技术实现思路
[0005]本专利技术提供了一种基于改进的微状态分析的运动想象脑电信号分类方法，对于微状态分析步骤的改进主要体现在需要进行聚类的电势地形图的选取方式上。将数值较高和较低的全局场功率(Global field power,GFP)对应的采样点分别称为“高GFP”和“低GFP”，而分析表明，“高GFP”和“低GFP”与四类微状态的空间相关度是相似的，也就是说“高GFP”和“低GFP”的空间相关度的可分性较差，这些采样点对于后续步骤中提取特征的贡献不大。因此，本专利技术在电势地形图的选取方式上进行改进，先将“高GFP”和“低GFP”进行去除再将剩余GFP局部峰值点时刻所对应的电势地形图进行聚类分析。
技术实现思路
整体框图如附图1所示。
[0006]本专利技术包括以下步骤：
[0007](1)计算脑电信号全局场功率(Global field power,GFP)，去除数值较高和较低的点，将剩余GFP局部峰值点时刻对应的电势地形图进行聚类分析，获得四种微状态；
[0008](2)将聚类得到的四种微状态分别与原始脑电信号的电势地形图进行空间相关度计算，构建出微状态时间序列，并计算出每个微状态对应的持续时间、出现频率、所占时间比率、转换概率等参数；
[0009](3)将训练集数据的微状态时间序列参数作为特征输入到支持向量机(Support Vector Machine，SVM)分类器中训练出分类模型，对测试集特征进行分类得到测试结果。
[0010]本专利技术与现有的运动想象脑电信号分类方法相比，具有以下创新点：
[0011](1)将脑电信号微状态分析法运用到运动想象脑电信号分类中，并对运动想象这一特定任务状态下EEG微状态分析步骤进行了改进。具体改进步骤为去除掉在空间相关度上的可分性较差“高GFP”和“低GFP”。
[0012](2)将“高GFP”和“低GFP”对应的电势地形图与四种微状态电势地形图的空间相关度进行统计分析，通过经验统计和Wilcoxon秩和检验分析出可分性，为改进步骤提供依据。
附图说明
[0013]附图1为
技术实现思路
整体框图
[0014]附图2为一个被试四种微状态电势地形图
[0015]具体实现方式
[0016]下面结合附图对本专利技术做详细叙述。以下实施例用于说明本专利技术，但不用来限制本专利技术的范围。
[0017]本专利技术是一种基于微状态分析的运动想象脑电信号分类方法，该方法包括以下步骤：
[0018](1)数据包含了九个被试四种运动想象思维任务，即左手想象运动、右手想象运动、舌头想象运动和双脚想象运动。每个被试在两个不同日期进行了两次运动想象脑电数据采集，每次采集可以分成若干个时间段，每个时间段进行4N次单次的运动想象实验任务，其中每个类别各N次。所有实验均使用22个电极记录，电极位置采用国际10/20系统。本实验选用其中的左手和右手的EEG信号进行二分类任务，截取了每位被试在视觉提示后的2s时间窗内运动想象脑电信号数据。
[0019](2)计算脑电信号全局场功率(Global field power,GFP)局部峰值点，其计算公式如下：
[0020][0021]其中v
i
(t)是在时间t的电极i处的电压，是在时间t的所以电极上的平均电压，N表示电极的数目。
[0022]对GFP进行统计分析，按数值大小排列，去除前15％和后15％的点，将剩余GFP局部峰值点时刻对应的电势地形图通过Modified k
‑
means算法进行聚类分析，获得四种微状态电势地形图。
[0023](3)全局电势图差异(Global Map Dissimilarity，GMD)是一种空间距离度量值，它不受信号强度的影响，而只考虑电势地形图的相似程度，其计算公式如下：
[0024][0025]其中，x
n
和x
n'
分别代表两个时刻的脑电数据，GFP
n
和GFP
n'
是对应时刻的全局场功率，N表示电极的数目。因此用GMD来表示四种微状态分别与原始脑电信号的电势地形图空间相关度，取每个时刻最高空间相关度的微状态电势地形图作为这个时刻的微状态，即标记此时刻的脑电为该微状态的序号(即1、2、3和4)，那么脑电信号就被表示成四种微状态交替出现的时间序列，也就是微状态时间序列。
[0026]计算出每个微状态对应微状态时间参数，包括：
[0027]1)平均持续时间：每种微状态稳定时，所占时间长度的均值。
[0028]2)每秒出现的频率：每种微状态每秒出现的次数。
[0029]3)所占时间比：每种微状态持续时间总和与脑电信号持续时间的比值。
[0030]4)不同微状态之间的转换概率：从当前微状态转换成另一种微状态的概率。
[0031](4)采用十折交叉验证的方本文档来自技高网...
【技术保护点】
【技术特征摘要】
1.基于改进的微状态分析的运动想象脑电信号分类方法，其特征在于计算脑电信号全局场功率(Global field power,GFP)，去除数值较高和较低的点，将剩余GFP局部峰值点时刻对应的电势地形图进行聚类分析获取微状态，将微状态分别与原始脑电信号的电势地形图进行空间相关性计算，构建出微状态时间序列并计算出每个微状态对应的时间参数，将训练集数据的微状态时间序列参数作为特征输入到SVM分类器中训练出分类模型，对测试集特征进行分类得到测试结果。2.根据权利要求1所述的基于改进的微状态分析的运动想象脑电信号分类方法，其特征在于改进的GFP的选取方式，具体包括：计算脑电信号全局场功率局部峰值点，对其进行统计分析，按数值大小排列，去除前15％和后15％的点，将剩余GFP局部峰值点时刻对应的电势地形图通过Modified k
‑
mean...
【专利技术属性】
技术研发人员：谢松云，付荧鑫，谢辛舟，王珍珍，段绪，
申请(专利权)人：西北工业大学，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。