预处理脑电数据一般包括以下步骤,本文不涉及未添加链接的步骤:
顺序 | 步骤 | 附注 |
---|---|---|
1 | 导入数据 | ·电极信息和事件手动导入电极信息和事件 |
2 | 定位电极 | ·导入电极位置相关文件 |
3 | 删除无用电极 | ·如剔除HEOG、VEOG、M1、M2 等电极点 |
4 | 重参考 | ·计算总平均参考或选择特定电极作为参考 |
5 | 滤波 | ·低通滤波 (设置范围为 30-100Hz) ·高通滤波(设定范围为0.1-1Hz) |
6-1 | 截取正确的反应期 | ·提取正确反应期[-3-3.5s,0s] |
6-2 | 截取刺激相关时段 | ·事件相关时段提取[-0.5s,-0.1s] |
7 | 基线矫正 | ·如事件前[-0.5s,-0.1s] |
8 | 剔除坏段 | ·消除极值(信号超过)±100uV) ·消除不可信的试验(如6个单电极标准或5个所有电极标准)·除去异常分布时间(如平均维度的5个标准差) |
9 | 去除伪迹 | ·独立成分分析 ·手动去除成分(如与眼电伪影相关的典型成分) ·基于ICA的分解 |
基于数据预处理的使用脑电数据采集,它是一种稳态视觉诱发电位,属于一种稳态诱发电位。稳态视觉诱发电位的响应通常在视觉皮层中产生,可用于视觉注意、工作记忆等认知评价。本文并不全部涉及数据预处理。
首先将Curry导入采集的脑电信号eeglab,如图所示:数据: 我的eeglab不能直接导入CDT这里我们打开类型文件eeglab导入curry插件,选择File->Manage EEG extensions ,在搜索框中搜索loadcurry并安装 选择安装成功后File->Import data->Using EEG functions and plugins->From Neuroscan Curry files导入CDT文件结果如下图所示:
需要知道的是,在用Curry收集脑电数据后,保存了相关电极位置信息.dpa在文件中,如果此时直接使用,EEGLAB读取数据文件.cdt和配置文件.dpa,电极位置等信息可自动读入eeglab,我们不再需要手动添加电极信息。选择Plot->Channel data 显示未处理的数据和选择Plot->Channel locations显示电极位置:
在这里,我们选择保留通道C3、Cz、C4、P7、P3、Pz、P4、P8、O1、Oz、O2,应该删除HEOG、FP1、FP2、VEOG、F7、F3、Fz、F4、F8、T7、T8、M1、M2、CP3、CPz、CP4这些电极。我在这里没有保留T7、T8.这两个电极可以作为平均参考,我在这里使用Cz详细操作见下一步。选择Edit->Selet data获得以下界面: 在Channel range处勾选 on ->remove these表示删除(不检查表示保留),选择要删除的电极,点击OK: 在保留文件界面,勾选Overwrie点击覆盖原始文件OK,发现文件大小变小,删除成功,如下图所示:
勾选Re-reference data to channels 选择我们的参考电极,如果选择两个平均电位值作为参考电极(个人理解)。选择后点击OK,再次Overwrite先保存更改后的数据,获取以下界面Reference处出现Cz
两个一起填代表进行带通滤波,将处理好的结果Overwrite覆盖原始数据。
点击OK然后继续覆盖保存数据,然后出现界面: 默认点击OK覆盖保存。在这里,我们可以观察脑电图1 2 3 4 5 …基线已包括:
采用去除伪迹ICA分离伪迹将分离的独立成分与神经活动相关,然后去除标记为伪迹的成分,重组获得真实相关的脑电图数据。常见伪迹异常成分的地形图如下: 1.地形图上的能量只集中在额叶上(眼电伪) 2.地形图不连续(噪声伪)) 3.地形图受限于单个电极(电极伪迹) 选择Tools->Decompose data by ICA得到如下界面,默认设置点击OK
进行ICA分解,等待一段时间。 分析完成后可以看到 ICA weigths完成,如下图: 展示一下除去伪迹后的脑电光谱图,如下图:
最后将我们简单处理好的脑电数据导出保存,选择File->save current dataset as保存为本地.set文件。