陳偉偉 大連理工大學(xué)信息與通信工程學(xué)院 116024

基于K均值的阻塞性睡眠呼吸暫停低通氣綜合癥診斷

陳偉偉 大連理工大學(xué)信息與通信工程學(xué)院 116024

引言

阻塞性睡眠呼吸暫停低通氣綜合癥[1]（Obstructive Sleep Apnea Hypopnea Syndrome, OSAHS）是一種具有嚴(yán)重危害及潛在危險的高發(fā)性疾病，由于睡眠時發(fā)生低氧及高碳酸血癥，最終引起全身多系統(tǒng)、多器官的漸進(jìn)性危害。OSAHS最顯著地臨床表現(xiàn)就是睡眠打鼾，表現(xiàn)為鼾聲時高時低，并可以完全中斷，嚴(yán)重患者可以憋醒。多導(dǎo)睡眠圖（Polysomnography, PSG）是公認(rèn)的監(jiān)測OSAHS的金標(biāo)準(zhǔn)。通過監(jiān)測患者在安靜入睡的狀態(tài)下的連續(xù)6～8小時的腦電圖、心電圖、肌電圖、血氧飽和度、鼾聲、呼吸動度、腿動等十多項指標(biāo)的改變診斷OSAHS， PSG[2]以AHI為標(biāo)準(zhǔn)對OSAHS病情程度評判。呼吸暫停指數(shù)（Apnea Hypopnea Index ，AHI）指數(shù),也稱為呼吸紊亂指數(shù)，是指每小時呼吸暫停加上低通氣的次數(shù)。呼吸暫停是指睡眠過程中口鼻呼吸氣流完全停止10秒以上；低通氣是指睡眠過程中呼吸氣流強(qiáng)度（幅度）較基礎(chǔ)水平降低50%以上，并伴有血氧飽和度較基礎(chǔ)水平下降大于等于4%。因為PSG監(jiān)測價錢偏高而且需要檢測者在醫(yī)院的檢查室中進(jìn)行整晚大約8個小時的監(jiān)測。調(diào)查顯示，大約93%的女性和82%的男性患者沒有得到診斷，原因就在于目前尚未找到一種簡捷，方便，低費用的檢測方法來實現(xiàn)OSAHS的篩查。

鼾聲檢測主要針對患者的鼾聲特性進(jìn)行研究與判斷，鼾聲的獲取只需一組無接觸式的麥克風(fēng)和數(shù)據(jù)記錄裝置。目前多數(shù)研究者主要從鼾聲的時域和頻域參數(shù)著手。時域參數(shù)多數(shù)以單位時間內(nèi)的鼾聲數(shù)，鼾聲持續(xù)時間和間隔，間斷鼾聲數(shù)，最后以鼾聲數(shù)和間斷鼾聲數(shù)為指標(biāo)進(jìn)行判斷[3]?；邝暀z測的OSAHS診斷研究在國內(nèi)外均尚屬起步階段。而且研究涉及廣大人民的身體健康問題，研究意義重大。本文主要應(yīng)用端點檢測的方法模擬呼吸暫停低通氣指數(shù)，并基于K均值分類的方法實現(xiàn)OSAHS的診斷。

1 鼾聲端點檢測

鼾聲實際上也是一種聲音。通過鼾聲檢測患者有病無病，可以廣義的看作為語音檢測、識別問題。不過，相對于常見的語音識別問題，鼾聲檢測從某種意義來說更簡單。因為它只需要檢測兩個狀態(tài)，即患者“有”或者“沒有”患病。因此語音信號處理中的理論、方法是研究該問題的主要研究方法。

現(xiàn)階段實驗用的鼾聲數(shù)據(jù)多數(shù)是與PSG檢測同時錄制的。由麥克風(fēng)采集到的鼾聲經(jīng)低通濾波和A/D轉(zhuǎn)換后得到鼾聲數(shù)字信號。鼾聲數(shù)據(jù)一般較長，基于音頻信號的短時性分析技術(shù)首先將鼾聲信號分幀和預(yù)加重。端點檢測[4]的目的就是區(qū)分鼾聲段和靜音段，有效地計算鼾聲片段的時間間隔，并可以將真正的鼾聲數(shù)據(jù)分離出進(jìn)行后續(xù)的處理。因為背景噪聲的不穩(wěn)定性依據(jù)固定閾值進(jìn)行端點檢測的正確率往往受限制，本文采用基于短時能量的自適應(yīng)端點檢測，并將鼾聲最短持續(xù)時間作為判斷條件加入來提高正確率。

1.1 端點檢測

將鼾聲信號分幀處理，因為鼾聲的時變性相對于語音信號的時變性較弱，本文取分幀幀長為128ms,幀移取64ms。應(yīng)用COOLEDIT軟件觀察鼾聲數(shù)據(jù)發(fā)現(xiàn)鼾聲的持續(xù)時間一般大于0.6s，故在端點檢測時將鼾聲的最短持續(xù)時間大于0.6s作為區(qū)分鼾聲的初始條件。這樣可以把持續(xù)時間稍長的噪音片段隔離掉，保證有效地提取真正的鼾聲數(shù)據(jù)。

第i幀的短時能量表示為：

其中l(wèi)en表示幀長。因為鼾聲數(shù)據(jù)的開端往往是以靜音開始的，本文將預(yù)讀入的5幀鼾聲數(shù)據(jù)的平均能量值為參考設(shè)置初始的閾值，自適應(yīng)的過程表示為:

其中LT和HT為經(jīng)驗數(shù)據(jù)，EL為低閾值，EH為高閾值，E為每幀能量值。自適應(yīng)的端點檢測算法能夠根據(jù)實時的噪聲變化調(diào)整門限。

1.2 鼾聲事件的處理

經(jīng)過端點檢測后得到每個鼾聲的開始時間和結(jié)束時間。鼾聲的持續(xù)時間定義為一次鼾聲的開始時間到結(jié)束時間，間隔時間定義為本次鼾聲的結(jié)束時間到下次鼾聲的開始時間。

觀察所有鼾聲的得出如下結(jié)論：98. 5%正常鼾聲的間隔時間的分布在1.4s到4.0s。而有的鼾聲因為個人呼吸頻率的不同，往往吸氣和呼氣的間隔時間相對他人的時間較長，一個完整鼾聲的中間就會出現(xiàn)小間隔。而這個間隔時間往往較小，所以我們對端點檢測后的數(shù)據(jù)做以下處理：將間隔時間小于0.1s的兩個鼾聲歸并成一個鼾聲，完成鼾聲事件的整合。臨床上的呼吸暫停是指睡眠過程中口鼻呼吸氣流完全停止10s以上，又因為打鼾者打鼾有時并不是連續(xù)的而是存在間斷時間，所以我們將間隔時間大于10s小于60s的鼾聲事件定義為呼吸暫停。

2 基于K均值的OSAHS診斷

臨床上規(guī)定每晚7小時睡眠中，呼吸暫停反復(fù)發(fā)作30次以上或AHI大于等于5次/小時以上為患有OSAHS。根據(jù)AHI指數(shù)的定義將整晚發(fā)生呼吸暫停的次數(shù)取平均數(shù)并以一定的閾值做判斷的方法沒有考慮到個體差異帶來的影響。因為人與人的呼吸道是不同的，打鼾的頻率也往往會有差異，有的人本身打鼾的頻率就比較低，間隔時間也就比較長，對診斷的結(jié)果就會產(chǎn)生影響。為了減小因個體差異帶來的影響首先應(yīng)用K均值分類的方法對鼾聲間隔時間分為兩類，并以兩類中心點的差值和發(fā)生暫停的次數(shù)為診斷依據(jù)。

2.1 K均值分類原理

K均值聚類算法[5]是將給定的數(shù)據(jù)集合分成確定的若干組。定義k個中心點,每組一個,由于不同的初始中心位置產(chǎn)生不同的聚類結(jié)果,所以選取適當(dāng)中心點是聚類的關(guān)鍵。 通過使它們之間的距離盡可能大,使給定的所有數(shù)據(jù)點結(jié)合到離它最近的中心點的聚類中。 當(dāng)所有的數(shù)據(jù)點都分配到中心點的范圍內(nèi)后,初始的聚類形成。通過先前階段的聚類中心重新計算k個新中心點,再將給定的數(shù)據(jù)重新分配到離它最近的新中心點。不斷進(jìn)行循環(huán),由循環(huán)的結(jié)果得知, k個中心點逐步地改變直到它們的位置不再變化為止,即聚類中心不再移動.則標(biāo)準(zhǔn)函數(shù)的最小值函數(shù)的定義為:

式中：為聚類組的數(shù)量，為第i個聚類的均值向量；z為得分向量；為聚類中的第i個聚類。K均值算法的具體步驟是：

1）為每個聚類確定一個初始的聚類中心，這樣k個聚類存在k個聚類中心；

2）將樣本集中的每一個樣本按照最小距離原則分配到k個聚類中某一個；

3）使用每個聚類中的所有樣本的均值作為新的聚類中心；

4）如果聚類中心有變化則重復(fù)2），3）步直到聚類中心不再變化為止；

5）最后得到的k個聚類中心就是聚類的結(jié)果。

2.2 OSAHS診斷

本文對鼾聲中存在的兩類間隔時間的定義如下：第一類，若間隔時間小于10s定義為正常的鼾聲間隔時間；第二類，若間隔時間大于10s小于60s定義為發(fā)生呼吸暫停的鼾聲間隔時間。首先確定K均值分類方法的初始中心點的取值，分別統(tǒng)計第一類間隔時間和第二類間隔時間的平均值作為初始中心點。

OSAHS診斷的依據(jù)：

1）當(dāng)只存在第一類的鼾聲間隔時間：即沒有發(fā)生呼吸暫停的，診斷未患OSAHS。

2）當(dāng)間隔時間中既包含第一類鼾聲間隔時間又包含第二類鼾聲間隔時間：若兩中心點的距離小于10s，則診斷為未患有OSAHS；若中心點的距離大于30s，則診斷為患有OSAHS，而且可以確定為比較嚴(yán)重的阻塞暫停；若中心點的距離大于10s小于30s則計算第二類間隔時間每小時的平均發(fā)生次數(shù)。

3 實驗結(jié)果與分析

本文實驗用的鼾聲數(shù)據(jù)是由大連醫(yī)科大學(xué)提供。表1是對端點檢測方法的評估。截取不同人的大約12～28分鐘的數(shù)據(jù)段統(tǒng)計鼾聲總數(shù)，并與兩個不同人耳聽到的鼾聲總數(shù)做對比，正確率在86%以上。不同的人耳聽到的鼾聲個數(shù)也會存在一定的誤差。

表2 K均值分類方法的評估Tab2. Evaluation of the K means

表3 預(yù)測AHI指數(shù)Tab3. Predicted AHI

表1 端點檢測方法的評估Tab1. Evaluation of the endpoint detection

表2中編號為5和7的鼾聲分別取于同一OSAHS患者的患病段和不患病段，對比發(fā)現(xiàn)患者患病時和不患病時的間隔時間存在很大的差異。2號數(shù)據(jù)是明顯的打鼾頻率較慢，若只采用暫停次數(shù)做判斷就會誤判為OSAHS，基于新算法正確判斷為正常的鼾聲。

表3中數(shù)據(jù)不同于表1和表2的數(shù)據(jù)，為整晚約7個小時的鼾聲。AHI指數(shù)范圍從2.6～72不等，包含了未患病的，中輕度OSAHS患者和重度OSAHS患者，預(yù)測AHI和實際AHI存在一定的差值。差值的存在可能是因為背景噪聲的影響，但不影響OSAHS的篩查。

4 結(jié) 論

通過仿真實驗得到以下結(jié)論：分析鼾聲的間隔時間來診斷OSAHS確實是行之有效的方法，而基于K均值的鼾聲間隔時間分類減小了因不同人的打鼾的頻率造成的分析誤差。預(yù)測AHI指數(shù)基本能夠體現(xiàn)患病的程度，實現(xiàn)OSAHS的篩查。

[1]殷善開, 易紅良, 曹振宇. 阻塞性睡眠呼吸暫停低通氣綜合癥[M]. 北京: 科學(xué)技術(shù)文獻(xiàn)出版社.2006.

[2]齊志勇, 張治平, 呼和牧仁等.多導(dǎo)睡眠儀(PSG)在阻塞性睡眠呼吸暫停低通氣綜合征中的臨床應(yīng)用價值[J]. 當(dāng)代醫(yī)學(xué). 2007（118）.

[3]Yeh Liang Hsu, Ming Chou Chen, Chih Ming Cheng, et al. Development of a portable device for home monitoring of snoring[C].2005 IEEE International Conference on Systems, Man and Cybernetics, Hilton Waikoloa Village, Hawaii, USA , October, 2005.

[4]鮑長春. 數(shù)字語音編碼原理[M]. 西安:西安電子科技大學(xué)出版社. 2007.

[5]邊肇慶, 張學(xué)工. 模式識別[M]. 北京:清華大學(xué)出版社. 2000.

Snore Detection for Obstructive Sleep Apnea Hypopnea Syndrome Diagnosis

Chen Weiwei School of Information and Communication Engineering, Dalian University of Technology, Dalian 116024,China

阻塞性睡眠呼吸暫停低通氣綜合癥最常見的癥狀是睡眠打鼾，本文依據(jù)語音信號處理的方法對鼾聲信號分析，實現(xiàn)對阻塞性睡眠低通氣綜合癥的篩查，以減小多導(dǎo)睡眠監(jiān)測的檢查負(fù)荷?；诙虝r能量的自適應(yīng)端點檢測記錄鼾聲的間隔時間，并以K均值的方法對鼾聲間隔時間分兩類，以兩類中心點的距離和發(fā)生暫停的次數(shù)為依據(jù)實現(xiàn)OSAHS的診斷。本文算法旨在減小個體差異對判斷結(jié)果的影響。

呼吸暫停；低通氣；多導(dǎo)睡眠圖；鼾聲；K均值

Snore is the typical symptom of the obstructive sleep apnea hypopnea syndrome. Based on the method of speech signal process realized the cheap and portable method of screening OSAHS, where reduce the load of PSG. Based on the adaptive short-time energy endpoint detection record the time of interval of snores, and the method of K mean classified the time of internal to two types, the distance between the center and the occurrence of sleep apnea to diagnosis OSAHS. The new method is designed to reduce the individual differences in judging the results.

OSAHS; AHI; PSG; snore; K means

TN912 文獻(xiàn)標(biāo)識碼：A

10.3969/j.issn.1001-8972.2010.22.092

陳偉偉（1986-）女，山東泰安人，大連理工大學(xué)信息與通信工程院碩士研究生，研究方向為語言信號處理?！?/p>