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导出图片