脑电数据预处理,eeglab预处理采集的SSVEP数据
腦電數(shù)據(jù)預處理一般包括以下步驟,其中未添加鏈接的步驟本文未涉及:
| 1 | 導入數(shù)據(jù) | ·導入腦電數(shù)據(jù)以及手動導入電極信息和事件 |
| 2 | 定位電極 | ·導入電極位置相關文件 |
| 3 | 刪除無用電極 | ·如剔除HEOG、VEOG、M1、M2 等電極點 |
| 4 | 重參考 | ·計算總平均參考或者選擇特定電極作為參考 |
| 5 | 濾波 | ·低通濾波 (設置范圍為 30-100Hz) ·高通濾波(設置范圍為0.1-1Hz) |
| 6-1 | 截取正確反應時段 | ·正確反應時段提取[-3.5s,0s] |
| 6-2 | 截取刺激相關時段 | ·事件相關時段提取[-0.5s,-0.1s] |
| 7 | 基線矯正 | ·減去基線均值,如事件前[-0.5s,-0.1s] |
| 8 | 剔除壞段 | ·剔除極值時段(信號超過±100uV) ·剔除不可信試次(如單電極6個標準差外或所有電極5個標準差外) ·剔除異常分布時段(如維度在平均維度5個標準差外) |
| 9 | 去除偽跡 | ·獨立成分分析 ·手動剔除成分(如眼電偽跡相關的典型成分) ·基于ICA的分解 |
本次數(shù)據(jù)預處理使用的為基于SSVEP采集的腦電數(shù)據(jù),SSVEP是穩(wěn)態(tài)視覺誘發(fā)電位,屬于穩(wěn)態(tài)誘發(fā)電位的一種。穩(wěn)態(tài)視覺誘發(fā)電位的響應通常在視覺皮層的V1中產(chǎn)生,SSVEP可用于視覺注意和工作記憶等認知評估。在數(shù)據(jù)預處理方面,以上操作本文未全部涉及。
1.導入數(shù)據(jù)首先將Curry采集的腦電信號導入eeglab,數(shù)據(jù)如圖所示:
我的eeglab不能直接導入CDT類型的文件,這里我們打開eeglab導入curry插件,選擇File->Manage EEG extensions ,在搜索框內(nèi)搜索loadcurry并安裝
安裝成功后選擇File->Import data->Using EEG functions and plugins->From Neuroscan Curry files導入CDT文件。結(jié)果如下圖:
需要知道的是,在用Curry8采集腦電數(shù)據(jù)后,相關電極位置信息已經(jīng)保存在.dpa文件內(nèi)了,此時如果直接采用EEGLAB讀入該數(shù)據(jù)文件.cdt和配置文件.dpa,則有關電極位置等信息就可自動讀入eeglab,我們也就不再需要手動為其添加電極信息了。選擇Plot->Channel data 展示未處理之前的數(shù)據(jù)以及選擇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,發(fā)現(xiàn)文件大小已經(jīng)變小,刪除成功,如下圖:
勾選Re-reference data to channels 選擇我們的參考電極,如果選兩個即計算二者的平均電位值作為參考電極(個人理解)。選擇好之后點擊OK,再次Overwrite保存更改后的數(shù)據(jù),得到如下界面,發(fā)先在Reference處出現(xiàn)Cz
兩個一起填代表進行帶通濾波,將處理好的結(jié)果Overwrite覆蓋原有數(shù)據(jù)。
7.基線矯正根據(jù)實驗事件中的事件編碼,將腦電數(shù)據(jù)劃分為多個數(shù)據(jù)段,如刺激前100ms到刺激后600ms,在本腦電數(shù)據(jù)中包含11個類別,這里我都選上選擇Tools->Extract Epoch 如下圖:
點擊OK后繼續(xù)覆蓋保存數(shù)據(jù),之后出現(xiàn)界面:
默認點擊OK覆蓋保存即可。這里我們觀察下腦電圖,可以看到1 2 3 4 5 …基線已經(jīng)包含:
去除偽跡,采用ICA分離偽跡將分離的獨立成分偽跡相關成分和神經(jīng)活動相關成分,再將標記為偽跡的成分剔除,重組得到真實相關的腦電數(shù)據(jù)。常見的偽跡異常成分地形圖表現(xiàn)為:
1.在地形圖中能量只集中在額葉(眼電偽跡)
2.地形圖不連續(xù)(噪聲偽跡)
3.地形圖受限于單個電極(電極偽跡)
選擇Tools->Decompose data by ICA得到如下界面,默認設置點擊OK
進行ICA分解,等待一段時間。
分析完成后可以看到 ICA weigths完成,如下圖:
展示一下除去偽跡后的腦電光譜圖,如下圖:
最后將我們簡單處理好的腦電數(shù)據(jù)導出保存,選擇File->save current dataset as保存為本地.set文件。
總結(jié)
以上是生活随笔為你收集整理的脑电数据预处理,eeglab预处理采集的SSVEP数据的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: php用pkcs1还是pkcs8,PKC
- 下一篇: vdbench数据校验翻译