脑电数据预处理_脑电波处理源代码

记录一个数据洁癖在脑电数据预处理阶段的经验和踩过的坑，若有不对的地方请批评指正。

我的机器是BV64导，实验是类似于静息态的任务态。

预处理顺序：

1、定位电极

2、删除无用电极，我这里是IO活动电极

3、高低通、NOTCH滤波

到这里是可以批处理的，代码如下：（路径和文件名需要换成各位自己的）

clc;clear all;close all
for i = 1:22
EEG = pop_loadbv('D:\00data\eegdata\', strcat('subject', num2str(i), '.vhdr'),[], []);
EEG=pop_chanedit(EEG, 'lookup','D:\\MATLAB_2013b_install\\toolbox\\eeglab\\plugins\\dipfit2.2\\standard_BESA\\standard-10-5-cap385.elp');
EEG = pop_select( EEG,'nochannel',{'IO'});
EEG = pop_eegfiltnew(EEG, [], 0.5, 6600, true, [], 1);
EEG = pop_eegfiltnew(EEG, [], 45, 294, 0, [], 1);
EEG = pop_eegfiltnew(EEG, 48, 52, 1650, 1, [], 1);
EEG = pop_eegfiltnew(EEG, 98, 102, 1650, 1, [], 1);
EEG = pop_saveset( EEG, 'filename',strcat(num2str(i), '.set'),'filepath','D:\\00data\\1filter\\');
end

4、select data using event把不同任务数据成段提取出来

clc;clear all;close all
for i = 1:22
EEG = pop_loadset('filename',strcat(num2str(i),'.set'),'filepath','D:\00data\1filter\\');
EEG = eeg_checkset( EEG );
EEG = pop_rmdat( EEG, {'S 10'},[0 180] ,0);
EEG = eeg_checkset( EEG );
EEG = pop_saveset( EEG, 'filename',strcat(num2str(i),'.set'),'filepath','D:\\00data\\2intention\\rest\\');
EEG = eeg_checkset( EEG );
end

5、分段，用代码分段。

clc;clear all;close all
for i = 1:22
EEG = pop_loadset('filename',strcat(num2str(i),'.set'),'filepath','D:\\00data\\2intention\\rest\\');
EEG = eeg_regepochs(EEG, 'recurrence', 1, 'limits',[0 1], 'rmbase',NaN);
EEG = pop_saveset( EEG, 'filename',strcat(num2str(i), '.set'),'filepath','D:\\00data\\3epoch\\rest\\');
end

6、去坏段，别去太狠了，否则数据损失太多，而且ICA也是可以修正一些的，我一开始去的狠，导致都没有眼漂成分。

7、接下来重参考。64导不算密集，能不能用平均参考？目前我读的领域内文献还是基本双侧乳突，因此，我这里是用TP9和10作参考电极。（100以上通道再考虑平均参考吧，完全静息态也可以平均参考。如果是平均参考，可以在预处理最后一步做，在保证数据干净后做）。

8、ICA

先通过map删成分。经验是，前两个一般是眨眼和眼漂，必选。

后面如果有全头绿，只有耳朵附近的一点高能量，则是坏成分，删。但！是如果不是全头绿，而是有激发大范围的高能，同时耳朵处一点红或蓝，说明这个成分里有有用的成分，就不能当做坏成分删掉。

肌电在15hz以上，如果数据分析后期要分析高频，那就得删删肌电了。一般肌电在滤波时删掉一些了，但难免也会在ica之后也出来一些，也是要注意是完全坏成分再删掉，表现为只有一两个trial是能量高的。还有常见的成分是均匀集中在后半部分的心电。以及坏导成分。

坏成分识别的方法：

（1）成分的折线图10hz突然高了，正常应该是逐渐下降或持平。

（2）全头绿，只有一点高能

不确定的就不要删了，免得删掉有用信息。

这个ICA还是要看个人感觉，比较主观，千万别纠结。

9、去除极端值

100 -100

10、手动检查

删除明显的伪迹。

这步也比较纠结，有的人说个别波形凸起对后面分析不影响的，但我记得有人说过别让这种明显的伪迹混入后续分析。我个人还是习惯把一些非常明显眼电等删掉，删的幅度还是比较大的。