趙今朝,劉銘, 江秀全,史維友,婁宜泰,冷建材,徐舫舟△,馮超,楊清波,唐吉友,魯珊珊

(1．齊魯工業(yè)大學(xué)(山東省科學(xué)院) 光電工程國際化學(xué)院,濟南 250353;2．齊魯工業(yè)大學(xué)(山東省科學(xué)院)數(shù)學(xué)與統(tǒng)計學(xué)院,濟南 250353;3．山東第一醫(yī)科大學(xué)附屬第一醫(yī)院(山東省千佛山醫(yī)院) 神經(jīng)內(nèi)科,山東省免疫研究所,山東省風(fēng)濕免疫與轉(zhuǎn)化醫(yī)學(xué)重點實驗室,濟南 250014)

0 引言

睡眠呼吸暫停使測試者的睡眠結(jié)構(gòu)發(fā)生變化,從而影響睡眠分期的準確評估。為此,睡眠呼吸暫停和睡眠分期的準確檢測在睡眠質(zhì)量評估中至關(guān)重要[1]。腦電圖(electroencephalogram, EEG)、心電圖(electrocardiogram, ECG)、肌電圖(electromyography, EMG)和其他生理信號可用于睡眠階段劃分[2]。腦電信號包含更多的信息,高振幅的信號通常表示大腦處于更興奮的狀態(tài),通過觀察腦電信號的振幅變化,可了解大腦的興奮程度,而頻率變化則反映大腦神經(jīng)元之間的同步性。不同頻率范圍內(nèi)的腦電節(jié)律活動,如α、β、θ節(jié)律等,結(jié)合在一起形成了多種特定振蕩模式。這些節(jié)律活動對應(yīng)不同的大腦狀態(tài),如放松、專注、睡眠等。通過分析腦電信號的頻率和節(jié)律,可更好地了解大腦的整體功能和狀態(tài)。

許多學(xué)者利用腦電信號,運用機器學(xué)習(xí)和深度學(xué)習(xí)技術(shù)進行睡眠階段檢測和睡眠分期。Zhao等[3]將C3-A2和C4-A1的EEG信號分為五個子帶,提取了每個子帶的樣本熵和方差,采用鄰域組成分析(neighbourhood component analysis, NCA)方法進行特征選擇、支持向量機 (support vector machine,SVM)作為分類器,進行呼吸暫停事件的檢測。Diykh等[4]使用中值、峰度等特征進行睡眠分期。Li等[5]針對腦電信號的時域、頻域和非線性特征,采用隨機森林方法進行分類,并考慮了睡眠階段之間的轉(zhuǎn)移特征,選擇適當?shù)奶卣?將其輸入機器學(xué)習(xí)分類器,進行睡眠分期。隨著深度學(xué)習(xí)的發(fā)展,許多研究開始直接從生理信號中學(xué)習(xí)特征。例如,Zhang等[6]提出了一種高效的卷積神經(jīng)網(wǎng)絡(luò),通過在第一層的每個卷積核上附加可訓(xùn)練的系數(shù),以增強學(xué)習(xí)能力和靈活性。Zhuang等[7]則提出了一種基于多通道腦電信號的自動睡眠分期方法,結(jié)合長短期記憶(long short term memory, LSTM)作為分類網(wǎng)絡(luò)。

大多數(shù)研究并未關(guān)注不同通道之間的相互協(xié)同,難以從全局層面探索睡眠階段信息[8-9]。fMRI已經(jīng)證實,每個睡眠階段都與特定的功能連接模式相關(guān)聯(lián)[10-12]。An等[13]研究了不同年齡段的睡眠腦電圖的連接性,發(fā)現(xiàn)不同睡眠階段的連接性不同。相位耦合反映了不同腦區(qū)之間的相互作用和協(xié)調(diào)程度。因為鎖相值(phase locking value, PLV)對相位敏感,可提供腦電信號中不同腦區(qū)之間相位耦合的信息[14-16],所以被廣泛應(yīng)用于分析腦電圖的連接。di Biase等[14]通過構(gòu)建大腦連接,根據(jù)其連通性區(qū)分帕金森病患者和健康者。

針對多通道研究的不足,本研究通過探索呼吸暫停患者在不同睡眠階段下腦電信號中的連通性,嘗試以此為基礎(chǔ)闡述呼吸暫?；颊吆粑鼤和：退唠A段的變化規(guī)律,并將不同腦區(qū)之間的連通性作為特征,實現(xiàn)呼吸暫停檢測和睡眠分期。針對睡眠過程中的特殊波形會對腦電圖的連接產(chǎn)生影響[17],本研究設(shè)計了一個基于子時段呼吸暫停和睡眠階段的分類策略,主要包括:(1)使用鎖相值構(gòu)建不同睡眠階段的大腦網(wǎng)絡(luò),將鎖相值當作特征,尋找出最優(yōu)的頻段和通道數(shù)目。(2)利用提出的多子時段決策的方式,進行睡眠階段分類和睡眠呼吸暫停檢測。(3)根據(jù)構(gòu)建的大腦功能網(wǎng)絡(luò)分析不同睡眠階段及呼吸暫停健康者的大腦相互作用和連接機制。

1 材料和方法

1.1 材料

本研究使用循環(huán)交替模式(cyclic alternating pattern, CAP)睡眠數(shù)據(jù)庫[18-19]。CAP數(shù)據(jù)庫中有108個樣本,包括92名睡眠障礙患者和16名健康受試者,本研究選擇10名受試者用于分析。此外,使用ISRUC-SLEEP 數(shù)據(jù)庫,此數(shù)據(jù)庫的Data of Subgroup_1包含100名受試者的數(shù)據(jù),本研究去除其中患有其他疾病及信號質(zhì)量差的受試者,選擇其中27名受試者的數(shù)據(jù)。

1.2 方法

1.2.1數(shù)據(jù)預(yù)處理 首先,將腦電信號分為多個不同的頻段。由于睡眠頻率范圍一般集中在0～30 Hz,因此在濾波過程中,本研究將信號分別在δ波(0.5～4 Hz),θ波(4～8 Hz),α波(8～13 Hz),β波(13～30 Hz),γ(30～40 Hz)五個頻段進行濾波。將濾波后的數(shù)據(jù)分別劃分為5段(6 s)、10段(3 s)、15段(2 s)、20段(1.5 s)和30段(1 s)方法。設(shè)L為30 s睡眠腦電數(shù)據(jù)樣本。當子時段的數(shù)量被設(shè)置為長Ls時,被劃分的子時段的長度滿足以下約束:

L=Ls×Ns

(1)

Ns是子時段數(shù)目。CAP數(shù)據(jù)庫中的采樣頻率為512 Hz,數(shù)據(jù)樣本持續(xù)時間為30 s,則子時段長度的計算如下:

Ls=L/Ns=(30×512)/Ns=15 360/Ns

(2)

由式(2)可知,子時段數(shù)Ns在實驗中起著重要作用。通過設(shè)置適當數(shù)量的子時段,可從中獲得有效的分類特征。然后,將數(shù)據(jù)在δ(0.5～4 Hz)、θ(4～8 Hz)、α(8～13 Hz)、β(13～30 Hz)和γ(30～40 Hz)五個頻帶中濾波。

1.2.2構(gòu)建功能網(wǎng)絡(luò) 本研究采用PLV評估大腦的功能連接[20],計算如下:

(3)

其中,N代表睡眠信號的樣本數(shù),t代表采樣周期,兩個時間序列分別為x(t)和y(t),兩個瞬時相位是φx(t)和φy(t)。實驗得到PLV對稱矩陣,矩陣中的每個值代表一對通道的耦合關(guān)系。

1.2.3選擇最優(yōu)頻段和通道數(shù)目 為降低算法復(fù)雜度,提高運行效率,減少計算負擔,本研究將CAP數(shù)據(jù)庫中12個通道的數(shù)據(jù)逐一連接,在矩陣中設(shè)置不同的閾值,以尋找連通性佳的通道。為進行對比,本研究選擇五個不同數(shù)量的EEG通道,包括12、10、8、6和4通道,在2.1章劃分的5個頻段中使用,選出最優(yōu)頻段,并采用10折交叉驗證來驗證所提算法的有效性。

1.2.4分類器及分類策略 LibSVM使用一對一法(one-versus-one,OVO)進行多分類。該方法通過設(shè)計一個SVM來處理任意兩類樣本之間的分類問題,因此,對于k個類別的樣本,需要設(shè)計k(k-1)/2個SVM。當對未知樣本進行分類時,根據(jù)得票最多的類別確定該未知樣本的類別。

為評估不同頻帶和不同數(shù)量的子時段中的腦功能連接,本研究使用三種特征處理策略對睡眠階段進行分類。(1)從30 s腦電數(shù)據(jù)中直接提取功能連接特征;(2)融合多個子時段提取的功能連接性特征,然后將特征作為一個整體輸入到分類器;(3)直接分析子時段的功能連通性特征并進行分類,然后將分類結(jié)果作為整個30 s片段的結(jié)果。具體描述見圖1—圖3。

圖1 用子時段數(shù)據(jù)構(gòu)建PLV矩陣

圖2 構(gòu)建PLV矩陣并特征融合

圖3 構(gòu)建PLV矩陣并經(jīng)過分類決策

2 結(jié)果

2.1 頻段和通道選擇結(jié)果

由表1可知,選擇六個通道時可獲得最佳性能。本研究分別測試了δ、θ、α、β和γ五個頻率范圍的5類睡眠階段任務(wù),具體分類結(jié)果見表2。其中,六通道獲得了最好的睡眠分期效果,且α頻段的指標優(yōu)于其他頻段,分類精度為69.4%。

表1 選擇不同通道數(shù)目及每百次事件運行時間

表2 從5種頻段提取特征對5種睡眠階段分類結(jié)果

2.2 分類策略對睡眠階段分類的結(jié)果

由表3可知,多段特征融合方法可獲得更高的分類性能。其中,30段法準確率最高,α頻帶睡眠分期準確率達到76.89%。多子段策略是將子時段提取的特征直接分類,分類結(jié)果是子時段的結(jié)果,不能代表原始腦電的類別,因此將分割的數(shù)據(jù)按時間順序重新組合。在組合的分類結(jié)果中,找到類別最多的分類結(jié)果。由表4可知,30個子時段在α頻段仍具有最佳分類性能,準確率達到89.73%。

表3 子時段特征融合方法下,5頻段5子時段數(shù)目對睡眠階段分類的準確率

表4 子時段分類決策方法下,5個頻段和5個子時段數(shù)目對睡眠階段分類準確率

2.3 分類策略對呼吸睡眠暫停分類的結(jié)果

在ISRUC-SLEEP數(shù)據(jù)庫Data of Subgroup_1的27個受試者數(shù)據(jù)中,找出相同數(shù)量的睡眠呼吸暫停和正常呼吸的數(shù)據(jù)。使用2.2所提方法進行呼吸睡眠暫停信號和正常呼吸信號分類,準確率達到了94.2%。

2.4 PLV的大腦網(wǎng)絡(luò)分析

由2.2和2.3小節(jié)得到的結(jié)果可見,α頻段具有分類的效果。圖4為在α頻段不同睡眠階段的連通性,紅線代表連通性高,藍線則相反。可見,N1階段的整體連接性大于快速眼動階段;N2階段的整體連接性大于N3階段;N3和W階段的連接性一般小于快速眼動階段。N1階段左右腦之間的連通性比N2階段強,但枕葉、頂葉和額葉之間的連通性比N2階段弱。與N3階段相比,N1階段的左右腦較弱,但連通性比N3階段強,且不同腦區(qū)的額葉和頂葉之間的連通性也比N3階段強。比較N2、W和REM階段,N2階段左右腦之間的連通性比REM階段弱,而在W階段,不同腦區(qū)之間的連通性更強。

圖4 α頻段不同睡眠階段大腦網(wǎng)絡(luò)連接性的比較

研究發(fā)現(xiàn),睡眠呼吸暫停測試者的各睡眠階段的數(shù)量會發(fā)生變化,其中深度睡眠階段N3的數(shù)量減少,過渡階段N1的數(shù)量增加。由圖5可知,睡眠呼吸暫停時的連通性明顯大于正常呼吸時,與N1階段和REM階段的連通性非常相似。

圖5 睡眠呼吸暫停和正常呼吸大腦網(wǎng)絡(luò)連接性的比較

2.5 交叉驗證

根據(jù)以上實驗,本研究通過PLV特征提取、多子時段特征融合和多子時段決策方法,最終發(fā)現(xiàn)30個子時段特征對決策具有最佳分類效果。為驗證該方法的有效性,在被試之間進行了10次交叉驗證實驗。睡眠階段分類的最終交叉驗證結(jié)果見圖6,結(jié)果是88.87%。正常呼吸和呼吸暫停分類的最終交叉驗證結(jié)果見圖7,結(jié)果是93.64%。

圖6 睡眠階段分類的結(jié)果

圖7 正常呼吸和呼吸暫停分類的結(jié)果

3 討論

本研究算法在N1階段的左右腦區(qū)的協(xié)同性大于REM階段,呈現(xiàn)出明顯的左右腦區(qū)對稱,N2階段全腦的活躍度明顯大于N3階段,因此,N2階段相對N3階段全腦連接更加活躍,N3階段和W階段的連通性小于REM階段,原因在于N3處于深度睡眠階段,腦區(qū)之間的活動減少,而W階段呈現(xiàn)出不規(guī)則的變化,相對REM階段的連接較弱。對于N2和REM階段,N2階段左右腦區(qū)的協(xié)同性比REM階段小,但同側(cè)活躍性大于REM階段。呼吸暫停期間的腦區(qū)更加活躍,明顯大于正常呼吸時,原因在于從深度睡眠逐步到REM階段和W期間過渡,因此腦區(qū)活動增加呈現(xiàn)出更加活躍的趨勢。不同數(shù)目的子時段對區(qū)分睡眠階段有不同的表現(xiàn)。30-Sub的情況顯示了良好的性能,但在α頻段使用30-Sub在N1和W級之間具有更高的錯誤率。本研究后續(xù)計劃將CAP數(shù)據(jù)庫中的多通道EEG信號用于不同的睡眠障礙的分類,如失眠和REM運動障礙。

4 結(jié)論

本研究通過探討多通道腦電信號在睡眠分期上的表現(xiàn),利用多個子時段獲得有效的融合特征,提出分類優(yōu)化策略,利用腦功能網(wǎng)絡(luò)分析睡眠分期的生理現(xiàn)象。不同的子時段對特征學(xué)習(xí)有不同的分類效果,在20-Sub的情況下,REM階段的分類比其他子時段的劃分更好,表明在20-Sub中提取的特征更適合REM階段的分類。最優(yōu)時段劃分數(shù)目與特征之間的關(guān)系將在后續(xù)工作中開展進一步研究。本研究探索了子時段腦電特征學(xué)習(xí)在多通道腦電睡眠分期中的應(yīng)用,實驗結(jié)果證明,本研究的多通道EEG睡眠階段算法的有效性和魯棒性較強,對提高睡眠分期的應(yīng)用具有重要潛力。