• 沒有找到結果。

數學認知研究腦電波訊號前處理

第四章 實驗方法

4.1 數學認知研究腦電波訊號前處理

我們在本研究中,腦電波擷取所設定的取樣頻率是1000Hz,也就是每 秒擷取1000個腦電波資料點。經過實驗得到的腦電波訊號資料中含有許多 的雜訊,其中以眼動訊號為影響最為顯著的雜訊,所以首先我們必需先將 眼動雜訊去除才可以進行往後的分波析作業。以往眼動訊號只採用一個 VEOU電極進行量測,之後再將原始腦電波被眼動雜訊干擾的部分濾除,

此種採用單電極量測眼動訊號的方法雖然簡便,但因為眼動電極所量測到 的眼動訊號中,實際上也含有腦電波訊號,以致將原始腦電波濾除眼動雜 訊時,也將一部分不屬於眼動雜訊的腦電波訊號濾除掉了,可能造成腦電 波訊號的不完整性。

本研究為了避免腦電波訊號的不完整性,將傳統量測眼動訊號的電極 增加為四個電極VEOU、VEOL、HEOR、HROL,VEOU電極置於左眼眉 毛上方,VEOU電極的位置垂直下方經過眼球的左眼下眼瞼位置,則是 VEOL電極的擺置位置,而HEOL電極則擺置於左眼與太陽穴中間,HEOR 電極則是擺放在右眼和右邊太陽穴中間,四個電極擺置位置如圖4.1所示。

圖4.1 眼動訊號電極擺置示意圖

眼動訊號的產生是由於眨眼與眼球移動所造成的,眨眼只有垂直方向 的運動,而眼球除了垂直方向外還會有水平方向的運動,所以我們另外增 加一組水平方向的電極HEOR、HROL,另外垂直方向也增加了VEOL電 極。因為每一個眼動電極所量測到的眼動訊號都包含了腦電波,首先我們 先將垂直方向的眼動電極組VEOU、VEOL的眼動訊號相減,眼動訊號相減 後已將腦電波訊號扣除,剩餘的眼動訊後我們視為真正垂直方向的眼動雜 訊,定義為VEOG眼動雜訊。水平方向的電極組HEOR、HROL也使用相減 的方法,得到HEOG眼動雜訊,我們將之視為真正的水平眼動雜訊。

我們量測處理眼動雜訊的方法稱為線性誘導(Linear derivation),此方法 就是將水平和垂直的眼動訊號進行兩兩相減,去除所含的腦電波訊號後只 剩下VEOG和HEOG訊號,分別代表垂直和水平的眼動雜訊,如圖4.2。

如前段所提到的,腦電波是相當微弱的訊號,所以如果我們沒將眼動訊號 濾除乾淨,這樣我們之後所分析的腦電波資料將包含極高成分的眼動雜 訊,有鑑於此,本研究將以往研究用來量測眼動雜訊的電極,由一個增加 到四個電極,並以線性誘導處理,以利我們之後的腦電波分析的準確性。

經由線性誘導(Linear derivation)後我們得到垂直與水平的眼動雜訊,腦 電波訊號是很微小的,所以眼動對腦電波所產生的雜訊相對是很大的,因 此在進行所有腦電波訊號分析,一定要先做去眼動雜訊的工作。我們使用 NeuroScan 4.3中Transforms中的Ocular artifact reduction指令,接著在出現如 圖3.10的指令框中的Trigger選項勾選Negative,此指令是依照眼動訊號發生 時的峰值為正或為負來決定,而在此指令框的右邊有一名為Blink Value的 欄位,需選取channel中的VEOG的電極選項,因為VEOG是實驗中經過線 性誘導後的眼動雜訊,選取後才可以讓軟體依此點的訊號來做眼動訊號的 消除,如圖4.3所示。在圖4.3指令框中的Output選項中,選取LDR+CNT選 項,在Review欄位中則都不勾選,如此即設定好去眼動雜訊的流程。

圖 4.3 去眼動指令框

我們以其中一位受測者的腦電波訊號圖來看,上圖為尚未去除眼動雜 訊的腦電波圖,在圖上最上方為眼動訊號,可以發現每當產生眼動訊號

時,會影響到其他正常腦電波的量測,但在經過去眼動訊號處理後,可以 發現本來受干擾的腦電波訊後都變得更加平順清楚,如圖4.4 所示。

圖 4.4 去眼動前、後腦電波圖

色的背景,代表已成功濾除。

經過以上步驟後,我們進行腦電波的分段工作,使用的則是NeuroScan 4.3 軟體中 Epoch file 的功能,在 Epoching properties 指令框中,我們在 Mode 選項中選取no trigger,並且在腦電波時間段落中設定為 1024 個資料點,

因為當初本研究設定的腦電波取樣頻率為1000Hz,所以分段後的每段腦電 波時間段落約為1 秒鐘,而經分段後的 1 秒鐘腦電波才是適合我們分析的 腦電波長度且又同時包含足夠的資訊,如圖4.5 所示。

圖 4.5 各電極分段後腦電波圖

圖 4.6 單一電極分段後腦電波圖

我們經由量測所得到的腦電波資料中,除了眼動雜訊外,還包含著許 多我們不需要的生理雜訊和室電訊號,這些生理雜訊與室電訊號我們使用

Artifact reduction 的功能來濾除,在 Operation 選項中選取 Reject on

criteria,而在 Criteria 選項中選取 Time domain,並將區將設定成符合腦電 波範圍的-85 至 85 之間,如圖 4.7 所示,最後選取需要濾除生理訊號的電 極,在這裡我們選取所有的電極,除了量測眼動訊號的電極VEOG 和 HEOG 除外。

圖 4.7 Artifact reduction 指令框

將分段後的腦電波濾除生理雜訊後,我們必須將每個電極的每段腦電 波進行平均的工作,因為我們分析著重於腦電波頻譜上的特徵,在

Averaging 指示框中我們就選取 Frequency domain,並且選擇我們所要平均 的Scaling 為 Amplitude,因為腦電波資料的振幅代表的電壓的大小值才是 我們分析的重點,最後在Window 選項中必須選取 Hanning 視窗來進行平 均腦電波的動作,如圖4.8 所示。經過平均後的每個電極的腦電波圖如圖

圖 4.8 Averaging 指示框

圖 4.9 各電極平均後腦電波圖

圖 4.10 單一電極平均後腦電波圖

經過以上腦電波的資料前處理工作後,我們大致上可以開始進行腦電 波的分析作業,本研究部分的目的為探討進行數學四則運算時,大腦的活 動特徵,為了明顯且有效率地達成此目的,我們必須先加以分析腦電波每 個主要頻帶,我們將在大腦認知活動中經常被各研究團隊所探討的Delta、

Theta、Alpha、Beta 頻段,分別作平均處理得到每個電極的腦電波位於各 個頻段的能量大小,並輸出成文件檔以利於我們進行之後的腦電波資料分 析,如圖4.11 所示。

圖 4.11 輸出各頻帶腦電波指令框