少通道下运动想象脑电特征的提取方法

摘要

本发明提供了一种少通道下运动想象脑电特征的提取方法，该方法利用每个通道时间序列的多个坐标延迟来同时优化空间滤波器与一个高阶多参数FIR滤波器，并用非参数化方法估计特征量与类别间信息熵的方法来选取最优个数的时间坐标延迟，基于多个时间坐标延迟的高阶FIR滤波器比一个参数的FIR滤波器有更好的幅频特性，该方法提供了在少通道情况下提取想象运动类别特征的有效途径。

权利要求

1.一种少通道下运动想象脑电特征的提取方法，其特征在于，包括步骤：步骤(1)脑电信号采集：根据不同运动想象实验范式，采集运动想象区的少通道脑电信号，所述少通道脑电信号为2～10路通道；

步骤(2)数据预处理：对采集到的脑电信号用滤波器做带通滤波，在运动想象实验范式中“提示”出现后的信号上做滑动窗记信号为矩阵Xt，时间间隔为一个采样点；

τ τ

步骤(3)信号延时：δ 表示信号的延时操作，即δ (Xt)＝Xt-τ；将采集到的每个通道的脑电信号做τ、2τ、…nτ的延时操作增加到采集到的信号中，将采集到的脑电信号扩展：i＝1,...,L，L为训练集中训练样本的总数目；

步骤(4)对扩展得到的高维信号进行CSP特征提取，分别由延时τ即n＝1、延时2τ即n＝2至延时Nτ构成的N个高维信号提取到N个特征集{f1(i)}i、{f2(i)}i…{fN(i)}

2*p×1

i，其中N是预定义的最大延时操作数，fn(i)∈R ,i＝1,...,L，L为训练集中训练样本的总数目，2*p为利用CSP算法降维后特征向量的维数；

步骤(5)使用非参数化方法估计特征集{fn(i)}i与类别间的互信息熵I(fn；C)＝H(C)-H(C|fn)，

其中I(fn；C)是特征集{fn(i)}i与类别间的互信息熵，H(C)是类别标签集合C的熵函数值，H(C|fn)是在给定特征集fn的情况下类别标签集合C的条件熵函数值；

步骤(6)选取 为最优特征集，其中 表示使互信息*

熵值取最大的特征集，n 为最优延时操作数；

*

步骤(7)对测试集脑电信号做τ、2τ…nτ的延时操作后，扩展成新的较高维数的脑电信号，用CSP算法提取特征；

步骤(8)输出分类结果。

2.如权利要求1所述的少通道下运动想象脑电特征的提取方法，其特征在于，步骤(1)中，不同运动想象实验范式为左右手运动想象实验范式。

3.如权利要求1所述的少通道下运动想象脑电特征的提取方法，其特征在于，步骤(1)中采集的采样频率为250Hz。

4.如权利要求1所述的少通道下运动想象脑电特征的提取方法，其特征在于，步骤(2)中在运动想象实验范式中提示出现后的信号上做2s的滑动窗记信号为矩阵Xt。

5.如权利要求1所述的少通道下运动想象脑电特征的提取方法，其特征在于，步骤(2)中滤波器为四阶巴特沃斯滤波器。

6.如权利要求1所述的少通道下运动想象脑电特征的提取方法，其特征在于，步骤(2)中带通滤波为通频带8～30Hz带通滤波。

7.如权利要求1所述的少通道下运动想象脑电特征的提取方法，其特征在于，步骤(2)中带通滤波为通频带7～32Hz带通滤波。

8.如权利要求1所述的少通道下运动想象脑电特征的提取方法，其特征在于，步骤(8)中，利用线性分类器LDA或是支持向量机做分类器输出分类结果。

说明书

少通道下运动想象脑电特征的提取方法

技术领域

[0001] 本发明涉及生物信息技术领域，尤其涉及一种少通道下运动想象脑电特征的提取方法。

背景技术

[0002] 目前，由于病变或意外造成的大脑与神经肌肉通路损伤的疾病在现有的医疗条件下还无法治愈。如在中国大约有六百万脑瘫患者，十二岁以下的儿童就占有一百八十万左右，这些患者失去肢体控制、自理等能力。而脑-机接口(Brain-Computer Interface，BCI)提供了制造治疗和改善这些患者疾病的医疗器械设备的技术基础，其主要机理是：训练患者进行相关的想象运动，同时采集其脑电信号，对信号进行相应的特征提取与识别，然后将识别结果作为控制量输入到相应的控制器以控制相应器械帮助患者实现日常的一些基本活动，从而提高他们的生活质量。

[0003] 由于大脑皮层对躯体运动具有交叉支配的特点，即一侧的皮质运动区支配对侧的躯体骨骼肌运动。人脑在有运动的意愿后，无论是实际输出运动还是仅仅想象运动(例如左手做抓、握运动或是右手做抓、握运动)，都会导致对侧大脑运动感觉区mu和beta节律的脑电信号失去同步的节律性活动而导致能量减弱，而其同侧运动感觉区的mu和beta节律的脑电信号节律性活动加强。这种现象称为事件相关去同步(Event-Related Desynchronization,ERD)和事件相关同步(Event-Related Synchronization,ERS)的生理现象。这种现象构成区分运动想象脑电信号的生理基础。

[0004] 在研究大脑生理现象时，现有的仪器如Neuroscan等商用脑电采集仪具有32导、64导和128导等多通道采集信号的能力。较多的通道数在一定程度上可以提高脑电信号的空间分辨率，但存在成本高、不易便携等不足。而BCI在实际应用中，易用性、便携性、成本低、抗干扰能力强等是实用化主要关注的问题。本领域的技术人员考虑到如果仅采集几个通道，甚至是2-3个通道(除参考极、参考地外)的脑电信号，采集设备的通道数就大大降低，从而使得设备简单、经济，同时操作方便、鲁棒性也增强。

[0005] 共空域模式(Common Spatial Pattern，CSP)方法在通道数比较多的情况下是基于运动想象脑-机接口中非常有效的一种空间滤波方法，而在只有几个通道EEG(脑电图描记器)信号的情况下，空间域信息量本身有限，直接应用CSP的效果甚至不如直接用带通能量法(Band Power，BP)。现有技术中存在一种想象单侧肢体运动的脑电特征的提取方法，该方法基于共空域模式CSP方法的脑电特征提取方法。这种方法在少通道脑电信号的情况下就存在不足，无法提供较高的分类正确率。

[0006] 最近的一些以减少脑电通道数的研究，如Lal、Arvaneh等学者通过研究通道优化选择问题，优化结果表明对于运动想象的脑-机接口，并非通道数目越多越好，在他们各自提出的优化目标下不同的受试者有不同数目和位置的最优通道布局。利用最优通道布局可以在用较少通道的脑电信号下有较高的分类正确率。然而，针对少通道数目脑电信号情况下，如何在保持分类正确率最高的情况下提取运动想象相关的特征这一关键问题并没有得到关注和解决。

[0007] 因此，本领域的技术人员致力于开发一种脑-机接口(Brain-Computer Interface，BCI)系统采集少量通道数目脑电信号(Electroencephalogram，EEG)情况下，运动想象脑电特征的提取方法。

发明内容

[0008] 为了解决上述技术问题，本发明提供了一种少通道下运动想象脑电特征的提取方法，该方法针对采集少通道脑电信号的情况，提出一种基于信息熵优化选取多个坐标延迟的共空域频谱模式的特征提取方法，以提高基于运动想象的脑电特征分类正确率。

[0009] 本发明的基本思想是：由于大脑是公认的复杂的非线性动力学系统，将时间坐标延迟利用低维时间序列重构多维相空间是重构复杂动力学系统的等价状态空间的常用方法之一。最初，共空域频谱模型的提出，是基于每个通道时间序列的一个坐标延迟提出的。该方法等价于同时优化空间滤波器与一个只有单个可变参数的有限冲击响应(Finite Impulse Response，FIR)滤波器的原理。由于一个参数的FIR滤波器的性能有限，本发明利用每个通道时间序列的多个坐标延迟来同时优化空间滤波器与一个高阶多参数FIR滤波器，并用非参数化方法估计特征量与类别间信息熵的方法来选取最优个数的时间坐标延迟。基于多个时间坐标延迟的高阶FIR滤波器比一个参数的FIR滤波器有更好的幅频特性，该方法提供了在少通道情况下提取想象运动类别特征的有效途径。

[0010] 为实现上述目的，本发明提供的少通道下运动想象脑电特征的提取方法包括以下步骤：

[0011] 步骤(1)脑电信号采集：根据不同运动想象实验范式，采集运动想象区的少通道脑电信号，少通道脑电信号为2～10路通道。

[0012] 优选地，步骤(1)中的所述不同运动想象实验范式为左右手运动想象实验范式。

[0013] 优选地，步骤(1)中所述采集的采样频率为250Hz。

[0014] 步骤(2)数据预处理：对采集到的脑电信号用滤波器做带通滤波，在运动想象实验范式中“提示”出现后的信号上做滑动窗记信号为矩阵Xt，时间间隔为一个采样点。

[0015] 优选地，步骤(2)中在运动想象实验范式中提示出现后的信号上做2秒钟的滑动窗记信号为矩阵Xt。

[0016] 优选地，步骤(2)中所述滤波器为四阶巴特沃斯滤波器。

[0017] 优选地，步骤(2)中所述带通滤波为通频带8～30Hz带通滤波。

[0018] 优选地，步骤(2)中所述带通滤波为通频带7～32Hz带通滤波

[0019] 步骤(3)信号延时：δτ表示信号的延时操作，即δτ(Xt)＝Xt-τ；将采集到的每个通道的脑电信号做τ、2τ、…nτ的延时操作增加到原始输入信号中，这样少通道的原始脑电信号扩展为多通道的较高维数的脑电信号。也即原信号延时nτ后并入原信号成为一个新的高维信号：

[0020]

[0021] 步骤(4)对新的高维信号进行CSP特征提取，分别由延时τ(也即n＝1)、延时2τ(也即n＝2)至延时Nτ构成的N个高维信号提取到N个特征集{f1(i)}i、{f2(i)}i…
2*p×1
{fN(i)}i，其中fn(i)∈R ,i＝1,...,L，L为训练集中训练样本的总数目，2*p为利用CSP算法降维后特征向量的维数。

[0022] 步骤(5)利用非参数化方法估计特征集{fn(i)}i与类别间的互信息熵[0023] I(fn；C)＝H(C)-H(C|fn)

[0024] 其中I(fn；C)是特征集{fn(i)}i与类别间的互信息熵，H(C)是类别标签集合C的熵函数值，H(C|fn)是在给定特征集fn的情况下类别标签集合C的条件熵函数值。

[0025] 步骤(6)选取 为最优特征集，n*为最优延时操作数。

[0026] 步骤(7)对测试集脑电信号做τ、2τ…n*τ的延时操作后，扩展成新的较高维数的脑电信号，用CSP算法提取特征。

[0027] 步骤(8)输出分类结果。

[0028] 优选地，利用线性分类器LDA或是支持向量机做分类器输出分类结果。

[0029] 在上述描述中，训练集是指运动想象实验中按照提示的类别(如左、右手)进行运动想象过程中采集到的脑电信号集合。测试集是指运动想象实验中想象类别未知的过程中采集到的脑电信号集合。

[0030] 大脑是公认的复杂的非线性动力学系统，将时间坐标延迟利用低维时间序列重构多维相空间是重构复杂动力学系统的等价状态空间的经典方法。最初的共空域频谱模型是基于每个通道时间序列的一个坐标延迟提出的，其等价于同时优化空间滤波器与一个只有单个可变参数的有限冲击响应(Finite Impulse Response，FIR)滤波器。由于一个参数的FIR滤波器的性能有限，本发明利用每个通道时间序列的多个坐标延迟来同时优化空间滤波器与一个高阶多参数FIR滤波器，并用非参数化方法估计特征量与类别间信息熵的方法来选取最优个数的时间坐标延迟。基于多个时间坐标延迟的高阶FIR滤波器比一个参数的FIR滤波器有更好的幅频特性，该方法提供了在少通道情况下提取想象运动类别特征的有效途径。

[0031] 本发明提供的少通道下运动想象脑电特征的提取方法具有如下优点：本方法利用原始信号的多个延时组成新的高维信号以克服原始信号空间域信息量有限的不足，同时所求得的空频滤波器等价于空域滤波器中嵌入一个高阶FIR滤波器，在优化空间域信息的同时优化频域信息；其次，最优延时因子个数由特征集与类别间的互信息熵来优化，在最优延时因子确定后，特征提取只需要将原始信号作相应个数延时，然后作线性降维求特征即可，除初始需要滤波，如8-30Hz带通滤波外，不再需要其他滤波，计算量很小，在线实现方便。

附图说明

[0032] 图1为本发明的较佳实施例中的少通道运动想象脑电特征提取方法的原理示意图。

具体实施方式

[0033] 以下结合附图对本发明较佳实施例的技术方案作进一步描述。

[0034] 图1为本发明的较佳实施例中的少通道运动想象脑电特征提取方法的原理示意图。在本实施例中，利用覆盖运动感觉区的少通道电极帽(如覆盖C3、Cz、C4)采集脑电信号，根据左右手运动想象实验范式进行运动想象输出。

[0035] 受试者按要求佩戴好脑电帽后坐在椅子上，显示器距受试者双眼约1米左右。脑-机接口系统中脑电信号的采样频率设为250Hz，设置50Hz陷波器去除工频干扰。

[0036] 在实验流程中，每个训练样本开始时屏幕全黑，同时出现白色十字，2s钟后‘哔’一声提示受试者集中精神准备开始。在第3s钟开始时有黑色条块叠加在白色十字上提示运动想象方式。当黑色条块在左侧时代表想象左手运动(如抓或握)，当黑色条块在右侧时代表想象右手运动，提示持续1.25s。要求受试者根据提示做相应的运动想象内容，持续想象3s以上，系统在第7s给出判断结果。在8.5s钟时十字消失，一个训练样本的实验结束，受试者可以休息，休息时间为1.5～2.5s，随机产生。左、右手运动想象各150次，顺序随机出现。取左、右运动想象各75次做训练集，其他训练样本做测试集。

[0037] 对采集到的脑电信号用四阶巴特沃斯滤波器做8-30Hz带通滤波，在运动想象实验范式中提示出现后的信号上做2s钟的滑动窗，时间间隔为一个采样点，为特征提取做准备。

[0038] 将每个通道原始信号做τ、2τ、…nτ的延时操作增加到原始输入信号中，将少通道的原始信号扩展为多通道的高维信号：

[0039] i＝1,...,L，L为训练集中训练样本的总数目。

[0040] 对新的高维信号进行CSP特征提取，具体按下列步骤完成：

[0041] 第一步：将截取后的滑动窗内由原始信号扩展后的高维数字信号组成的数字矩阵按提示信号的左XLi、右XRi类别进行分类，对每个动作类别的信号矩阵求协方差矩阵，并对所有动作求平均。类别平均协方差矩阵由下式给出：

[0042] CL＝Exp(XLi*XLiT)

[0043] CR＝Exp(XRi*XRiT)

[0044] 其中，Exp()为期望，下标L和R分别代表左、右两个类别，CL和CR为左、右手运动想象的协方差矩阵。XLi表示第i个训练样本是属于左手运动想象类别的脑电信号，T是转置操作符。

[0045] 第二步：将求得的两类个别平均协方差矩阵同时对角化，该目标由下述广义特征值分解实现：

[0046] CL·w＝λCR·w

[0047] 其中w是所求的空域滤波器，λ为对应的广义特征值。

[0048] 第三步：将广义特征值分解中对应最大的3个特征值以及最小的3个特征值的特征向量构成CSP降维矩阵(也即空域滤波矩阵)

[0049] 第四步：利用所求得的降维矩阵 构造特征向量

[0050]

[0051] 式中w1、w2…为对应最大、第二大…等广义特征值的特征向量。

[0052] 通过上述四个步骤完成CSP降维矩阵求解，以及求得每个训练样本的CSP特征向量。这样的特征向量放在一起形成特征向量集{fn(i)}i。

[0053] 改变延时操作nτ的数值n＝1,2,...N，求得N组CSP特征向量集{f1(i)}i、{f2(i)}i、…{fN(i)}i，i＝1,...,L，其中L为训练集中训练样本的总数目。

[0054] 利用非参数化方法估计特征量集{fn(i)}i与类别间的互信息熵，具体由下式给出：

[0055] I(fn；C)＝H(C)-H(C|fn)

[0056] 式中 为类别的先验概率；给定特征变量fn下的条件熵为：

[0057]

[0058] 其中fn(i)为延时因子个数为n的第i个训练样本的特征向量。P(ω|fn(i))为给定特征向量fn(i)时类别ω的条件概率，根据贝叶斯定理可以求得：

[0059] 特征fn(i)的条件概率密度函数p(fn(i)|ω)用基于Parzen窗的非参数化方法来估计，具体由下式给出。

[0060]

[0061] 其中 y(j)∈Iω表示第j个训练样本的特征向量的类别标签属于ω的指标集。在高维数据中，Ψ常取对角阵形式，其对角元素为：

[0062] Ψkk＝hσ2，k＝1,...,2*p

[0063] σ为对应特征分量的方差，h为平滑参数取为 用来控制窗函数的宽度，2*p为特征向量的维数。于是，由上述式可以求得不同延时因子个数下的I(fn；C)。

[0064] 接下来，选取 为最优特征集，n*为最优延时操作数。对测试集中脑电信号做τ、2τ…n*τ的延时操作后，得到测试集脑电信号的高维扩展矩阵。然后再用训练集得到的CSP降维矩阵对测试集高维扩展矩阵提取特征。

[0065] 最后，利用线性分类器LDA或是支持向量机做分类器判别测试集中特征向量所属类别，输出分类结果。

[0066] 以上详细描述了本发明的较佳具体实施例。应当理解，本领域的普通技术无需创造性劳动就可以根据本发明的构思做出诸多修改和变化。因此，凡本技术领域中技术人员依本发明的构思在现有技术的基础上通过逻辑分析、推理或者有限的实验可以得到的技术方案，皆应在由权利要求书所确定的保护范围内。