EEGLAB使用流程
EEGLAB使用流程
MATLAB版本 2016a; EEGLAB版本 v2020.0
数据的准备
%matlab code clear all ;clc; filepath = strcat('200mM yuanshi\subject\subject3_',num2str(1),'.txt') sig=importdata(filepath); x=sig.data(2:1:end,2:2:end);%偶数项为脑电信号 event = zeros(1,length(x(:,1))); kaishi = 2.5;%这是应看到,5s一次刺激,后面改进使用算法寻找事件标记的开始 event(kaishi*256:kaishi*256 + 10) = 1; event((kaishi+5)*256:(kaishi+5)*256 + 10) = 1; event((kaishi+10)*256:(kaishi+10)*256 + 10) = 1; event((kaishi+15)*256:(kaishi+15)*256 + 10) = 1; event((kaishi+20)*256:(kaishi+20)*256 + 10) = 1; event((kaishi+25)*256:(kaishi+25)*256 + 10) = 1; x(:,end+1) = event;%事件标记通道 x = x';通道定位的文件
Standard-10-20-CCostmos.locs:
1 -18 0.51111 Fp1 2 18 0.51111 Fp2 3 -39 0.33333 F3 4 39 0.33333 F4 5 -90 0.25556 C3 6 90 0.25556 C4 7 -141 0.33333 P3 8 141 0.33333 P4 9 -162 0.51111 O1 10 162 0.51111 O2 11 -54 0.51111 F7 12 54 0.51111 F8 13 -90 0.51111 T3 14 90 0.51111 T4 15 -126 0.51111 T5 16 126 0.51111 T6 17 -43 0.5458 Sp1 18 43 0.5458 Sp2 19 0 0.25556 Fz 20 90 0 Cz 21 180 0.25556 Pz 22 180 0.51111 Oz 23 -90 0.63889 A1 24 90 0.63889 A2数据的导入(MATLAB array)
本数据采集的时候没有事件标记通道,是将采集的被试者的脑电信号后,通过算法计算出来的事件标记点,该数据一共有25个通道,24个通道是脑电通道,第25个通道是后期通过计算得出的事件标记点。
File->Import data->Using EEGLAB functions and plugins->From ASCII/float or Matlab array.MATLAB工作空间变量的名字是 x,采样率为256Hz(Data sampling rate),点击OK即可,也可以其他数据格式的导入。
输入数据的命名。点击OK
可以通过Plot->Channel data(scroll)把数据EEG波形图画出来观察一下。
现在点击REJECT就会把选中的数据丢失掉,就可以观察将不正常的数据剔除。只要是类似的界面,都可以实现数据的剔除
确定事件标记通道
File->Import event info->From data chanel
因为第25个通道是事件标记通道,所以输入25,也可以采取其他事件标记形式。
也可以在Edit -> Events values下面窗口进行事件标记的添加、命名和删除。
通道定位
Edit->Channels locations,下拉列表选择第二个参数,点击OK即可
弹出下面窗口
点击Read locations,选择通道定位文件,我的文件是Standard-10-20-CCostmos.locs(不同实验定位文件不同,需要根据实验修改文件),直接自动监测,点击OK。
数据的滤波
只说明最经典的FIR滤波器,Tools->Filter the data->Basic FIR filter(new,dault)
第一个参数是高通滤波,第二个是低通滤波,当两个参数都有输入的时候就是带通滤波,但是EEGLAB官网不推荐带通滤波,会出现错误,所以要高通低通滤波两次,这里参数选择0.5Hz-45Hz。将滤波的数据命名为subject_FIR。
运行ICA,去除伪迹
Tools->Decompose data by ICA,其他默认即可,点击OK,
ICA运行的时候是可以打断的,但是打断之后再运行是重新开始计算。
单纯查看ICA分量:Plot->Component maps等
查看并剔除ICA分量:Tools->Inspect/label component by maps
点击每个分量进行查看判断是否为伪迹
点击ACCEPT将其转化为REJECT就是将该分量舍弃。点击OK。
绘制ERP
Tools->Extract epochs
第一个参数(bengin):事件标记点
第二个参数(-0.5 1.5):数据选取事件标记点前面0.5秒和事件标记点1.5秒
第三个参数(subject_FIR_epochs):新数据组合的名字
点击OK,
数据的基准线是事件标记点前500ms。
之后就可以在Plot菜单栏中绘制你想要的ERP波形了
Plot->Sum/Compare comp.ERPs可以绘制多个数据的ERP的均值方差对比
第一个参数是选用哪个数据,直接填写数据的编号集
总结
以上是生活随笔为你收集整理的EEGLAB使用流程的全部内容,希望文章能够帮你解决所遇到的问题。
- 上一篇: openGL天空盒实现-立方体贴图
- 下一篇: eeglab导出图片