脑电波(Electroencephalogram,EEG)是一种使用电生理指标记录大脑活动的方法,大脑在活动时,大量神经元同步发生的突触后电位经总和后形成的。它记录大脑活动时的电波变化,是脑神经细胞的电生理活动在大脑皮层或头皮表面的总体反映。脑机接 口(BCI),是一种通过检测中枢神经系统活动并将其转化为人工输出,来替代、修复、增强、补充或者改善中枢神经系统的正常输出,从而改变中枢神经系统与内外环境之间的交互作用的技术。
将EEG数据加载进eeglab之后,
eeglab
能够读取到的仅仅是每一个通道的名称和数值,但是它并不知道这个通道是位于头皮上的什么地方,因此我们需要加载进一个跟记录数据时相匹配的通道位置信息eeglab中默认的文件是 standard-10-5-cap385(国际10-5系统)
在EEGLAB中预处理数据的时候,使用
ICA
方式去除眼电成分,这种方式可以不用眼电通道数据的参与就可以进行。
滤波分为
高通滤波
,低通滤波
,带通滤波
和凹陷滤波
。
high-pass
,低通low-pass
,英文看起来就更为直观了。数据打上mark,将数据切分为一段一段,在这些范围以外的的数据,直接去除。 一般的ERP分段在事件前200毫秒到事件后1000毫秒这个时间段就足够了,因为几乎所有的ERP成分,都是在1秒以内产生的。但是如果后期还要做时频分析的话,则需要将分段时间拉长,扩大到事件前1秒到事件后2秒这个时间段内,因为时频分析的算法,要求数据要有一定的长度。
被试被认为是处于一个相对平静的状态,此时的脑电活动代表了一个平静状态下的脑电活动,将这段时间内的脑电活动当做成一个基线。 一般来说,都是以分段中,事件之前这段时间来作为基线的,比如分段为-200毫秒到1000毫秒,那就将所有的数据,减去前200毫米以内数据的均值。 基线校正的第二个作用是防止数据漂移带来的影响。有时候因为某些原因,我们采集到的数据,会呈现一个逐渐缓慢向上漂移的趋势,而在分完段之后,由于向上飘的影响,每一段数据的起点将不在同一个地方,也会使得这段数据的绝对波幅变高。而基线校正就可以校正这种漂移带来的影响,让每一段数据都拥有一个差不多的起点。
横向操作
和纵向操作
,横向即对通道进行校正操作,对数据不好的导联进行插值处理,纵向是指挑出数据不好的trials,删除掉。建议是先横向处理后再纵向剔除,尽量保留下更多的trials数。这是一个很有用但也很难的方法,我并不是十分清楚它的原理,只能大概理解,也只能大概描述一下。前面说过,我们记录到的是头皮脑电,脑电帽上划出了一个个的点来表示位置,但是我们在
FPz
点记录到的数据,就真的是FPz
点这个位置头皮下方的区域所产生的电活动?不一定。做过实验的都知道,眨眼会产生影响,左右看会产生影响,帽子戴得太紧导致肌肉紧张也会产生影响,还有左右的FP1,FP2
处产生的电活动也会影响到FPz
,甚至在离它最远的Oz点下方产生的电活动,也有可能对前方FPz
点产生微弱的影响。因为真正的电活动产生于头皮下方颅骨内部,它经过了这么一层又一层的传播之后,不同源的电活动肯定会相互影响,从而导致某个记录点记录到的数据,混合了很多不同成分的电活动。而独立主成分分析,就是要对记录点的数据做一个逆运算,把每个记录点的数据, 分解成一个又一个的成分组成。然后我们再从中剔除掉伪迹成分,比如眨眼的成分,肌肉紧张的成分等,从而得到一个相对干净的数据。
ICA是一种分离线性混合源的技术。例如,我们可以尝试混合然后分开两个源。
基于ICA的伪影校正可以通过线性分解从EEG数据中分离并去除多种伪影。ICA方法基于以下假设:
共空间模式(Common Spatial Pattern, CSP)
是一种对两分类任务下的空域滤波特征提取算法,能够从多通道的脑机接口数据里面提取出每一类的空间分布成分。公共空间模式算法的基本原理是利用矩阵的对角化,找到一组最优空间滤波器进行投影,使得两类信号的方差值差异最大化,从而得到具有较高区分度的特征向量。
由于不同特征的信息包含在信号的不同频段中,因此考虑EEG信号,可以将其分离到不同的特定频段进行分析。在当前研究的特征提取和特征分类算法中,
CSP
是性能最好和应用最广泛的特征提取方法。CSP
将当前数据转换到另一个新的空间,使得类间数据的方差最大化,类内的方差最小化。通常,在运动想象EEG信号分类中,
CSP
算法往往选择宽频段(4-40HZ)来进行信号处理。研究表明,这个宽频段包含了大多数的运动想象相关的信息,同时,也包含着其他冗余信息,所以不少研究工作重点考虑这个宽频段的选择最佳问题,以此来改善CSP的分类性能。在提出的共空间频谱模式(CSSP)
中,有限脉冲响应(FIR)滤波器
在CSP滤波器
的基础上进行了优化,通过插入时间延迟τ来实现,该时间延迟τ允许单独调谐频率滤波器,评估若干样本点的时间延迟τ值,并且使用10倍交叉验证来选择最佳值。使用带通滤波器对信号进行滤波,然后进行空间滤波和特征提取。实验结果显示CSSP改进了性能。提出的共稀疏频谱空间模式(CSSSP)
进一步改进了CSSP方法,其找到所有信道共用的频谱模式,而非CSSP为每个信道找到不同的频谱模式。已经提出的
子带公空间模式(SBCSP)
,其中在多个子带处对运动图像EEG信号进行滤波,并且从每个子带提取CSP特征。为了减小子带的维数,线性判别分析(LDA)
分别应用于每个子带的特征,并融合在一起以进行分类。 SBCSP实现了比CSP,CSSP和CSSSP更高的分类精度。然而,SBCSP忽略了从不同子带获得的CSP特征的可能关联,因此提出了滤波器组CSP(FBCSP)
来解决该问题。 FBCSP计算来自多个子带的CSP特征的互信息,以便选择最具辨识性的特征,所选特征使用支持向量机(SVM)分类器进行分类。 FBCSP的表现优于SBCSP,然而,它仍然使用了若干子带,这些子带依然会带来不少的计算成本的增加。所以提出DFBCSP
来解决这个问题,DFBCSP利用单信道(C3,C4或Cz)频带功率的Fisher比(FR)来从多个重叠子频带中选择最具判别性的子频带。然后为每个子带提取CSP特征,并使用SVM分类器进行分类。与SBCSP和FBCSP相比,DFBCSP实现了更高的分类准确度和更低的计算成本。
考虑
DFBCSP
的优劣,因此提出了一种改进的DFBCSP
方法来解决这些问题。首先,使用从使用所有通道数据生成的特征计算的互信息,而不是像在DFBCSP-FR
中那样使用单通道频带功率的FR,来选择能够产生最佳结果的频带。而不使用仅使用单个信道频带功率和FR作为选择子频带(DFBCSP-FR
)的方法。这是因为EEG信号主要受噪声污染。因此,如果用于计算FR的单个信道被噪声破坏,则选择出来的频段将会影响分类器的分类精度。这导致选择的子带往往不是最有效的子带,因为可能选择了具有冗余信息的子带。因此,文章建议利用所有可用信道数据通过利用互信息来选择最具判别性的子带,以便获得最佳结果。与使用用于频带选择的单信道信息相比,使用所有信道数据进行频带选择减少了选择具有冗余信息的子频带的机会。 其次,文章不再像在DFBCSP-FR中那样仅使用来自重叠子带的CSP特征,而是结合利用宽带CSP和CSSP的特征信息,文章引入了7-30Hz的单个宽带以及在4-30Hz范围内的12个重叠子带,带宽为4Hz并且重叠2Hz。
改进DFBCSP算法的详细流程框图
EEG信号特征提取就是以脑电信号作为源信号,确定各种参数并以此为向量组成表征信号特征的特征向量。
特征參数主要包括
时域信号(如幅值)
和频域信号(如频率)
两大类,相应的特征提取方法也分为时域法
、频域法
和时-频域
方法。
频域分析方法主要是基于EEG信号各频段功率、相干等。
这类方法是建立在假设EEG信号具有平稳特性的基础上,同时只是考虑信号的频域信息,忽略信号在时间上的分辨率。在EEG信号研究中,常用的频域分析方法包括
功率谱估计
(直接法和间接法
)。 运动想象EEG信号的ERD
和ERS
现象只出现在特定频率范围,比如8-12Hz 的Mu波
和18-26Hz 的Beta 波
。因此,通过功率谱等谱分析方法,也可以有效地从EEG提取中特征。其中比较出名的方法有自回归功率谱分析、双谱分析法等。
时域分析方法则主要分析EEG,如
幅度
、均值
、方差
、偏歪度
、峭度
等。这类方法主要是利用EEG波形的性质,如波幅、均值、方差、偏歪度和峭度等对临床EEG记录进行观察分析。在EEG信号研究中,常用的时域分析方法有:
过零点分析
、直方图分析
、方差分析
、相关分析
、峰值检测
、波形参数分析
和波形识别
等。
这是近年运动想象领域比较通用的特征提取方法,主要通过设计空域滤波器对EEG 的多通道空间分布进行处理,提取可分的特征。包括:
共空域模式法(Common spatialpattern,CSP)
,以及基于CSP改进的方法。
基于CSP改进的方法包括:
鉴于脑电信号的非线性特性和运动想象时的节律特性,提出了小波模糊熵的特征提取方法,利用小波变换将EEG信号进行小波分解,得到对应运动想象EEG信号的alpha和beta节律,然后采用模糊熵方法提取特征。
快速傅里叶变换(Fast Fourier transform, FFT)
经典的FFT在分析确定信号和平稳信号时很有效,但在分析突变信号的频谱时具有一定的局限性。对于alpha波、beta波等脑电基本节律,功率谱分析或FFT等方法可以从频域中直接提取出特征。
稀疏傅里叶变换(Sparse Fourier transform, SFT)
SFT是时间窗内信号特征的平均,时间窗里信号越短,获得的时间分辨率就越高。该方法在一定程度上克服了FFT不具有的局部分析能力的缺陷,在某些信号处理中有一定的作用。
自适应自回归(Adaptive Auto Regressive ,AAR)
AAR模型参数方法随每一样本点的输入而改变,能够很好地反映大脑的状态。不过,该方法更适合分析平稳的脑电信号,对包含高度非平稳信号的运动想象EEG,该模型是有局限性的。
小波变换(Wavelettransform, WT)
小波变换在低频时的时间分辨率较低,而频率分辨率较高;在高频时的时间分辨率较高,而频率分辨率较低,正符合低频信号变化緩慢而高频信号变化较快的特点。因此小波分析比傅里叶变换和短时傅里叶变换更适合分析非平稳的脑电信号。