劉學(xué)森，張俊仕

(1.河南省駐馬店市中心醫(yī)院心內(nèi)科，河南駐馬店 463000；2.新疆醫(yī)科大學(xué)，新疆烏魯木齊 830011)

1 引 言

心房顫動(dòng)簡(jiǎn)稱為房顫[1],是一種常見的心血管疾病，主要特征為快速性心律失常。房顫容易誘發(fā)腦中風(fēng)、血栓、心力衰竭等疾病，情況嚴(yán)重者會(huì)危及生命?，F(xiàn)今，我國(guó)房顫患者高達(dá)1 000多萬(wàn)，隨著居民生活質(zhì)量的不斷提升，房顫患者出現(xiàn)了明顯的年輕化趨勢(shì)，引起醫(yī)學(xué)領(lǐng)域的極大關(guān)注[2]。

依據(jù)房顫發(fā)作持續(xù)時(shí)間，可以將其劃分為三個(gè)階段，分別為陣發(fā)性房顫、持續(xù)性房顫與永久性房顫。由于陣發(fā)性房顫具有發(fā)作持續(xù)時(shí)間短、復(fù)發(fā)率高、心電圖難捕捉等特點(diǎn)，隨著患者病程的延長(zhǎng)，會(huì)逐步演化成持續(xù)性房顫或永久性房顫，威脅患者的生命[3]。故如何檢測(cè)陣發(fā)性房顫被稱為醫(yī)學(xué)領(lǐng)域重點(diǎn)研究課題之一。心電圖是常用的心血管疾病檢測(cè)手段，主要是記錄患者心臟活動(dòng)發(fā)生電位變化曲線，以此來(lái)反映患者心臟情況。由于現(xiàn)有陣發(fā)性房顫自動(dòng)檢測(cè)算法提取心電特征過(guò)于單一，無(wú)法滿足患者的醫(yī)學(xué)檢測(cè)需求，故本研究提出基于多特征融合的陣發(fā)性房顫自動(dòng)檢測(cè)算法研究。

2 陣發(fā)性房顫自動(dòng)檢測(cè)算法設(shè)計(jì)

2.1 心電信號(hào)降噪方法設(shè)計(jì)

心電信號(hào)采集過(guò)程中，由于受多種背景噪聲的影響難以檢測(cè)較微弱的心電信號(hào)，且嚴(yán)重的噪聲環(huán)境會(huì)掩蓋大量有效的疾病信息。心電信號(hào)中噪聲通常包含工頻干擾、肌電尾跡、電級(jí)接觸噪聲、內(nèi)部噪聲與隨機(jī)噪聲等[4]。為了更加精確地檢測(cè)陣發(fā)性房顫，基于小波包分解算法對(duì)心電信號(hào)進(jìn)行降噪處理[5]，步驟如下：

(1)

ρi表示的是能量比率，由子帶信號(hào)經(jīng)過(guò)計(jì)算而得。

步驟三：更新全部子帶信號(hào)相應(yīng)小波包系數(shù)，公式為：

(2)

(3)

α表示的是調(diào)節(jié)參數(shù)；λ表示的是固定閾值。

步驟四：以步驟三更新小波包系數(shù)為基礎(chǔ)，重構(gòu)心電信號(hào)，獲得降噪后心電信號(hào)S。

2.2 R峰自動(dòng)檢測(cè)

以降噪后心電信號(hào)為基礎(chǔ)，基于小波變換技術(shù)自動(dòng)檢測(cè)R波波峰，為下述心電信號(hào)特征提取做準(zhǔn)備[6]。

利用寬度為T秒的滑動(dòng)窗對(duì)降噪后心電信號(hào)S進(jìn)行分段，心電片段記為S=(S1,S2,…,SN)，N表示的是片段總數(shù)量?；谛〔ㄗ儞Q對(duì)心電片段Si進(jìn)行離散化處理，將其分解為K個(gè)子帶信號(hào)，記為Si(CK)。

依據(jù)心電信號(hào)QRS波群的頻率分布范圍，選取適當(dāng)?shù)淖訋盘?hào)作為QRS占優(yōu)子帶[7]，記為Si(CK*)，以此為基礎(chǔ)，檢測(cè)并確定候選R峰位置[8]。首先，定義Si(CK*)上轉(zhuǎn)折點(diǎn)集合為TURi，動(dòng)態(tài)閾值為THRi，在集合TURi中篩選出全部大于動(dòng)態(tài)閾值的點(diǎn)，表達(dá)式為：

(4)

將式(4)結(jié)果在QRS占優(yōu)子帶Si(CK*)上相應(yīng)的位置序列記為ri，即為候選R峰位置序列。

利用小波變換[9]將候選R峰位置序列ri映射到心電片段上，獲取心電片段上的候選R峰位置序列，記為Ri。由于Ri存在著一定的誤差，利用半徑為v的加窗微調(diào)心電片段候選R峰位置序列，并仍將其記為Ri。

觀察心電信號(hào)RR間期局部變化趨勢(shì)，依據(jù)RR間期序列的左、右鄰平均值篩查R峰，降低個(gè)體差異，避免誤檢情況發(fā)生，并能夠避免房顫心電片段中的R峰被去除，提升陣發(fā)性房顫自動(dòng)檢測(cè)的精度[10]。其中，RR間期序列的左、右鄰平均值計(jì)算公式為：

(5)

Mleft(p)與Mright(p)分別表示的是RR間期序列的左、右鄰平均值；RRi(j)表示的是RR間期序列。

2.3 心電信號(hào)特征提取

基于上述R峰自動(dòng)檢測(cè)結(jié)果，依據(jù)熵理念，分別在P波缺失與RR間期不規(guī)則兩個(gè)角度提取心電信號(hào)特征，為下述陣發(fā)性房顫自動(dòng)檢測(cè)提供可靠的、多樣的數(shù)據(jù)支撐。

2.3.1基于P波的特征提取 正常心電信號(hào)中，每個(gè)心動(dòng)周期內(nèi)均包含一個(gè)圓鈍并平滑的P波[11]。但若是人體出現(xiàn)房顫，正常竇性P波缺失，取而代之的為f波，其具備大小不一、間隔不均、快速振蕩等特性。故本研究通過(guò)能量變化指數(shù)來(lái)描述P波與f波之間的不同。

根據(jù)已有研究可知，P波頻率分布范圍約為0～20 Hz，故該頻率范圍內(nèi)的子帶信號(hào)對(duì)陣發(fā)性房顫?rùn)z測(cè)更加有效。將0～20 Hz頻率范圍內(nèi)的子帶信號(hào)稱為P波占優(yōu)子帶，記為Si(Cj)。以此為基礎(chǔ)，通過(guò)短時(shí)傅里葉變換獲取心電片段在子帶信號(hào)Si(Cj)的頻譜圖，見圖1。

由圖1可知，陣發(fā)性房顫心電功率譜在大部分頻率范圍內(nèi)高于正常心電功率譜。當(dāng)時(shí)間t固定時(shí)，陣發(fā)性房顫心電功率譜值隨著頻率變化產(chǎn)生明顯改變，而正常心電功率譜值無(wú)明顯變化，并保持基本穩(wěn)定[12]。

圖1 心電片段頻譜圖

將頻譜圖按照時(shí)間軸方向切割為多個(gè)薄片，記為“SP切片”。則Si(Cj)能量變化指數(shù)計(jì)算公式為：

(6)

則對(duì)于K*個(gè)可用子帶信號(hào)來(lái)說(shuō)，心電信號(hào)相應(yīng)P波特征為Fp={EVINK*}N×K*。

2.3.2基于RR間期的特征提取 通常心房率與心室率保持一致，約為60～100次/min，RR間期規(guī)則劃一；當(dāng)人體發(fā)生陣發(fā)性房顫時(shí)，心房率與心室率發(fā)生變化，心房率可以達(dá)到心室率的2～3倍，約為350～600次/min，RR間期呈現(xiàn)出絕對(duì)不規(guī)則的特征[13]。故本研究通過(guò)RR間期振蕩指數(shù)來(lái)描述正常心電與陣發(fā)性房顫心電之間的不同。

圖2 心電片段RR間期分布圖

正常心電信號(hào)中RR間期長(zhǎng)度幾乎保持一致，而陣發(fā)性房顫心電信號(hào)中RR間期分布較為散亂，長(zhǎng)度變化區(qū)間較大[15]。

心電信號(hào)RR間期振蕩指數(shù)計(jì)算公式為：

(7)

2.4 陣發(fā)性房顫自動(dòng)檢測(cè)

為了提升陣發(fā)性房顫自動(dòng)檢測(cè)精度，利用決策級(jí)融合算法將上節(jié)提取的心電信號(hào)特征進(jìn)行融合，通過(guò)支撐向量機(jī)自動(dòng)檢測(cè)陣發(fā)性房顫[16]。

利用決策級(jí)融合算法獲得陣發(fā)性房顫心電融合特征：

F={Fp,FRR}N×(K*+4)*

(8)

支撐向量機(jī)實(shí)質(zhì)上是一個(gè)二分類模型，可以總結(jié)為特征空間上的最大間隔超平面，利用數(shù)學(xué)方法對(duì)特定問(wèn)題進(jìn)行求解。支撐向量機(jī)最關(guān)鍵部分為分類決策函數(shù)，表達(dá)式為：

(9)

sign(·)表示的是分類決策函數(shù)；b*表示的是分類決策函數(shù)參量。

通過(guò)上述過(guò)程實(shí)現(xiàn)的陣發(fā)性房顫自動(dòng)檢測(cè)，不僅為陣發(fā)性房顫患者提供更加有效的診斷手段，也為患者爭(zhēng)取到較多的治療時(shí)間，防止病情惡化為持續(xù)性房顫或永久性房顫。

3 數(shù)值實(shí)驗(yàn)與結(jié)果分析

3.1 實(shí)驗(yàn)數(shù)據(jù)庫(kù)選取

為驗(yàn)證本研究算法性能，選取MIT-BIH數(shù)據(jù)庫(kù)作為實(shí)驗(yàn)數(shù)據(jù)庫(kù)。本研究采用MIT-BIH數(shù)據(jù)庫(kù)中房顫數(shù)據(jù)集合，該集合中包含25條心電片段，需要注意的是，其中23條為陣發(fā)性房顫心電片段，另外兩條心電片段標(biāo)記未完全。陣發(fā)性房顫數(shù)據(jù)參數(shù)見表1。

表1 陣發(fā)性房顫數(shù)據(jù)參數(shù)表

在房顫數(shù)據(jù)集合中，選取04015、04043、04126、04908、04936、07910、08378等記錄中的30 min心電片段作為實(shí)驗(yàn)數(shù)據(jù)集。

由于篇幅的限制，本研究只展示記錄04043的心電片段，見圖3。

3.2 實(shí)驗(yàn)評(píng)價(jià)指標(biāo)選取

為了直觀顯示本研究算法的檢測(cè)性能，選取敏感度、特異度與準(zhǔn)確率作為算法檢測(cè)性能實(shí)驗(yàn)評(píng)價(jià)指標(biāo)，計(jì)算公式為：

(10)

圖3 記錄04043心電片段圖

sensitivity、specificity與accuracy分別表示的是敏感度、特異度與準(zhǔn)確率；TP表示的是真陽(yáng)性，指的是陣發(fā)性房顫心電被正確檢測(cè)；FN表示的是假陰性，指的是陣發(fā)性房顫心電被誤檢；TN表示的是真陰性，指的是正常心電被正確檢測(cè)；FP表示的是假陽(yáng)性，指的是正常心電被誤檢。

3.3 實(shí)驗(yàn)結(jié)果分析

依據(jù)選取的實(shí)驗(yàn)數(shù)據(jù)進(jìn)行數(shù)值實(shí)驗(yàn)，通過(guò)六次不重復(fù)實(shí)驗(yàn)并計(jì)算平均數(shù)值作為實(shí)驗(yàn)結(jié)果。通過(guò)個(gè)人獨(dú)立檢驗(yàn)與交叉檢驗(yàn)兩種方式驗(yàn)證該算法的檢測(cè)性能。通過(guò)個(gè)人獨(dú)立檢驗(yàn)得到實(shí)驗(yàn)評(píng)價(jià)指標(biāo)數(shù)據(jù)見表2。

表2 個(gè)人獨(dú)立檢驗(yàn)實(shí)驗(yàn)評(píng)價(jià)指標(biāo)數(shù)據(jù)表

依據(jù)表2數(shù)據(jù)進(jìn)行交叉檢驗(yàn)，則得到實(shí)驗(yàn)評(píng)價(jià)指標(biāo)數(shù)據(jù)見表3。

由表2可知，個(gè)人獨(dú)立檢驗(yàn)實(shí)驗(yàn)評(píng)價(jià)指標(biāo)數(shù)據(jù)較好，表明本研究算法具備有效性。由表3可知交叉檢驗(yàn)中，檢測(cè)效果也優(yōu)良，并具備普適性，說(shuō)明本研究算法在臨床應(yīng)用中具備較高的輔助作用。兩種檢驗(yàn)方式結(jié)果存在一定的差距，主要由患者間個(gè)體差異性造成，并且該影響因素?zé)o法避免。

上述數(shù)值實(shí)驗(yàn)結(jié)果顯示，本研究算法敏感度為94.6%，特異度為93.7%，準(zhǔn)確率為94.0%，滿足現(xiàn)今陣發(fā)性房顫的檢測(cè)需求，具有可行性。

4 結(jié)束語(yǔ)

本研究基于多特征融合技術(shù)設(shè)計(jì)了陣發(fā)性房顫自動(dòng)檢測(cè)算法，極大地提升了算法的敏感度、特異度與準(zhǔn)確率，為陣發(fā)性房顫臨床檢測(cè)提供了更加有效的手段，為陣發(fā)性房顫醫(yī)學(xué)研究提供一定的參考。