使用EEGLAB函数处理自己的脑电数据流程及代码

新建电极位置

新建电极位置文件

创建对应脑电数据结构

导入数据

创建事件

创建events事件

ica去除伪迹

ICA

分段

depoch=[-0.2,0.5];%事件-0.2s到0.5s的数据分段
outEEG = pop_epoch(EEG,{'erp'},depoch);%自定义的事件名称是erp，对这个事件分段

平均参考

%reref data
outEEG = pop_reref(EEG,[]);

去除基线

% remove baseline
dbase=[-200,0];%单位是ms
outEEG = pop_rmbase(EEG,dbase);

滤波

[outEEG, com, b] = pop_eegfiltnew(EEG, 'locutoff', 1,'hicutoff',10);%1-10Hz滤波

绘制所有通道的电压图

figure;
pop_plottopo(EEG,1:16,'erp',0);%这里的erp是给图取的名字

所有trial的erp图

figure
pop_erpimage(EEG,1);%在弹出的窗口选择通道

通道平均电压图和头皮图

dtimp=[-1000,8000];%绘制窗口
dtopo=[5000:500:7000];%绘制头皮图的时间点
figure
pop_timtopo(EEG,dtimp,dtopo,'erp','plotchans',[3 4]);%绘制底3 4两个通道的

时序头皮图

pop_topoplot(EEG,1,[-500:500:8000],'图',0,'plotchans',[1:1:5]);

ERSP

%spband为显示的频带范围
spband=[8,30];
%erspmax为色标范围
erspmax=[];%设置为空则为默认
pop_newtimef(EEG,1,3,[],0,'freqs',spband,'plotitc','off','erspmax',erspmax,'plottype','image','title','c3');%第三个参数为设定通道3，

PSD

figure;
ftopo=[10 12 20 25];%显示头皮图的频率值
fwin=[8,30];%显示的频率范围
pop_spectopo(EEG,1,[],'EEG','freq',ftopo,'plotchans',[1 2 3 4 5],'freqrange',[fwin(1) fwin(2)],'title','closePo');