重慶郵電大學(xué)光電工程學(xué)院 吳 優(yōu) 龐 宇 黃博強(qiáng) 彭良廣

穿戴式抗干擾血氧飽和度采集裝置設(shè)計

重慶郵電大學(xué)光電工程學(xué)院 吳 優(yōu) 龐 宇 黃博強(qiáng) 彭良廣

在消防營救和長時間有氧運(yùn)動中,由于缺乏對生理參數(shù)的及時獲取,時常造成嚴(yán)重的人員傷亡.為此,本文設(shè)計了一種以手套為載體的穿戴式血氧飽和度采集系統(tǒng),利用反射式探頭采集脈搏波信號,經(jīng)由STM32L151單片機(jī)進(jìn)行數(shù)據(jù)處理,并引入雙樹復(fù)小波變換和形態(tài)學(xué)濾波去除噪聲干擾和基線漂移,最后通過藍(lán)牙發(fā)送至手機(jī)端app獲取生理參數(shù).經(jīng)測試,裝置體積微小,佩戴方便,具備低功耗與抗運(yùn)動干擾能力,并且和標(biāo)準(zhǔn)監(jiān)護(hù)儀相比血氧飽和度平均誤差不超過2%.

穿戴式;血氧飽和度;雙樹復(fù)小波變換;形態(tài)學(xué)濾波;抗運(yùn)動干擾

0 引言

血氧飽和度作為一項(xiàng)重要的生理參數(shù),它間接反映出人體細(xì)胞的新陳代謝狀態(tài).

據(jù)統(tǒng)計,我國每年有近30名消防員犧牲,超過300人受傷致殘;在運(yùn)動領(lǐng)域,每年也有近150人發(fā)生心源性猝死.由于高危高壓的消防營救與過勞性的"有氧運(yùn)動"時常使身體處于缺氧狀態(tài),若長時間得不到恢復(fù)則會對生命造成嚴(yán)重危害.

出于對生命安全的考慮,消防和運(yùn)動領(lǐng)域?qū)τ谘躏柡投鹊膶?shí)時檢測尤為重要.目前,市面上存在的血氧飽和度檢測裝置主要運(yùn)用于醫(yī)院監(jiān)護(hù)和家庭監(jiān)護(hù),在靜態(tài)環(huán)境下具有較高的準(zhǔn)確性,但在消防施救和運(yùn)動過程中則存在很大局限.

考慮到裝置的便攜性、實(shí)時性、續(xù)航時間,以及抗運(yùn)動干擾能力,本文提出了一種新的解決方案:通過佩戴內(nèi)嵌微型裝置的手套即可完成血氧飽和度和脈率的采集,并利用藍(lán)牙將數(shù)據(jù)發(fā)送至手機(jī)端,再通過手機(jī)推送至后臺,從而達(dá)到檢測與反饋的目的.

1 總體系統(tǒng)設(shè)計

裝置采用透射式血氧探頭采集脈搏波信號,采集點(diǎn)為指尖部位.探頭內(nèi)置于一個能完全裹住手指第一關(guān)節(jié)的橡膠套中,由于橡膠材質(zhì)的伸縮性,因而可以使探頭很好的貼合于皮膚表面,而不會因?yàn)槭种傅幕顒又率固筋^移位.微型電路模塊內(nèi)藏于手套背部口袋,由于體積微小,佩戴手套后不會對消防和運(yùn)動人員造成任何不便和影響.系統(tǒng)總設(shè)計框圖如圖1所示.

圖1 系統(tǒng)總設(shè)計框圖

2 硬件設(shè)計

裝置硬件部分主要分為電源管理模塊、脈搏波信號采集模塊、單片機(jī)控制和處理模塊、藍(lán)牙模塊.硬件整體框圖如圖2所示.

2.1 單片機(jī)控制電路

本設(shè)計選擇低功耗微處理器STM32L151CBT6,其運(yùn)用了ST自主研發(fā)的130nm的超低漏電工藝技術(shù),具備先進(jìn)的超低功耗模式、優(yōu)化的動態(tài)運(yùn)行功耗以及特殊的安全性能,在高性能和超低功耗之間取得了很好的平衡.控制電路如圖3所示.

圖2 硬件整體框圖

圖3 單片機(jī)控制電路

單片機(jī)系統(tǒng)時鐘采用HSI,以16Mhz頻率運(yùn)行;程序下載方式為SWD;引腳red_on、ir_on、LEDir、LEDred與脈搏波信號采集電路相連,在通用定時器的時序下控制紅光和紅外光的交替產(chǎn)生;引腳OUT用于捕獲每一次光強(qiáng)探測器傳回的原始脈搏波信號;RX與TX則為UART串口收發(fā)引腳,負(fù)責(zé)將處理完的數(shù)據(jù)傳輸給藍(lán)牙模組.

2.2 脈搏波信號采集電路

脈搏波采集電路由兩部分組成:一是由三極管和電阻構(gòu)成的驅(qū)動電路,二是由二極管和光強(qiáng)探測器組成的信號采集電路.如圖4所示.

圖4 脈搏波信號采集電路

三極管的型號為S9012與S9013,起電流放大的作用,以驅(qū)動二極管發(fā)光.

考慮到氧合血紅蛋白(oxyhemoglobin,HbO2)與脫氧血紅蛋白(deoxyhemoglobin,Hb)對不同波長的光吸收度不同,而在波長660nm與940nm處區(qū)別最大[1],因此選擇既能發(fā)出紅光(波長660nm)又能發(fā)出紅外光(波長940nm)的雙光源led,其由兩個極性相反的led組合封裝而成,具有兩個引腳,型號為EZDL-R660/IR905.

光強(qiáng)探測器為TSL235R型光頻轉(zhuǎn)換傳感器,其頻率響應(yīng)范圍為320nm至1050nm,可以將接收到的光轉(zhuǎn)換成一定頻率的數(shù)字信號,由輸出管腳發(fā)送給單片機(jī)捕獲和處理.

3 軟件設(shè)計

軟件設(shè)計包括脈搏波信號的采集、傳輸、處理、計算和發(fā)送.其流程如圖5所示.

圖5 軟件流程圖

程序首先初始化各模塊,包括配置系統(tǒng)與各外設(shè)時鐘;配置I/O口的輸入、輸出、復(fù)用等模式;配置定時器及其中斷模式;配置ADC與DMA傳輸方式;初始化UART串口等.

led時序由定時器控制,以12ms為1個周期采集兩組脈搏波信號,其中紅光與紅外光分別點(diǎn)亮2ms,其余時間處于熄滅狀態(tài),以降低功耗;當(dāng)led點(diǎn)亮?xí)r,單片機(jī)的TIM4_CH2通道也開始捕獲數(shù)據(jù),如此循環(huán)采集2000組數(shù)據(jù)后統(tǒng)一進(jìn)行處理,最后將處理完的數(shù)據(jù)代入公式計算出血氧飽和度和脈率.

UART與藍(lán)牙模組相連,當(dāng)單片機(jī)計算完成后即通過串口將數(shù)據(jù)傳輸給藍(lán)牙,藍(lán)牙再將數(shù)據(jù)發(fā)送給手機(jī)端顯示出血氧飽和度和脈率值,最后推送至后臺.

電池電量檢測運(yùn)用到A/D轉(zhuǎn)換,當(dāng)檢測出電量低于預(yù)設(shè)值時,利用串口向手機(jī)發(fā)送低電量提醒,一定時間后采取自動關(guān)機(jī).

4 數(shù)據(jù)處理與計算

4.1 信號處理

4.1.1 形態(tài)學(xué)濾波

正常情況下信號的中心線為直線,但實(shí)際波形里中心線往往隨時間上下擺動,稱其為基線漂移.基線漂移的實(shí)質(zhì)是信號直流分量不穩(wěn)定,主要由人體自身的呼吸運(yùn)動,以及探測器與皮膚之間接觸位移所致[2].采用形態(tài)學(xué)濾波可以消除基線漂移.

形態(tài)學(xué)濾波是一種非線性濾波技術(shù),算法核心為通過腐蝕、膨脹、形態(tài)開和形態(tài)閉等運(yùn)算[3],濾除脈搏波信號特征波形,再利用原始信號減去剩下的基線漂移信號,從而得到直流分量穩(wěn)定的正常信號.

腐蝕是一種使邊界向內(nèi)收縮的過程,而膨脹則使邊界向外擴(kuò)張.信號f關(guān)于結(jié)構(gòu)元素g的腐蝕、膨脹運(yùn)算分別如下:

若對信號先腐蝕,再膨脹,則為開運(yùn)算,公式如下:

若對信號先膨脹,再腐蝕,則為閉運(yùn)算,公式如下:

開運(yùn)算能夠去除毛刺與孤立的小點(diǎn),閉運(yùn)算則能填合孔洞與彌補(bǔ)縫隙.形態(tài)學(xué)濾波通過形態(tài)開閉和形態(tài)閉開運(yùn)算的組合達(dá)到濾波的效果,其流程如圖6所示.

圖6 形態(tài)學(xué)濾波流程圖

4.1.2 雙樹復(fù)小波變換

雙樹復(fù)小波變換(dual-tree complex wavelet transform,DTCWT)最初由Kingsbury提出[4],后又經(jīng)Selesnick等人不斷完善[5].

DTCWT基于離散小波變換,但克服了其平移敏感性與易造成頻率混疊的缺陷[6].具體形式為采用二叉樹結(jié)構(gòu)的雙路離散小波變換并行分解和重構(gòu)信號,左子樹L對應(yīng)實(shí)部,右子樹R對應(yīng)虛部,分別包含一組低通和高通濾波器,兩組濾波器可構(gòu)成Hilbert變換對,使雙樹復(fù)小波變換近似解析[7].變換原理如圖7所示.

圖7 DTCWT原理圖

實(shí)樹與虛樹對應(yīng)濾波器之間的延遲相差一個采樣間隔,虛樹的采樣位置始終位于實(shí)樹的中間,以使雙樹抽取所得數(shù)據(jù)成互補(bǔ)關(guān)系,從而降低信息的丟失,實(shí)現(xiàn)近似平移不變性[8].

雙樹復(fù)小波變換通過分解脈搏波信號各層的低頻小波系數(shù),再利用閾值去除高頻小波系數(shù)中所含的噪聲成分,然后將分解的低頻小波系數(shù)和閾值處理后的高頻小波系數(shù)進(jìn)行信號重構(gòu),從而得到去除噪聲的光電容積脈搏波信號.

4.2 數(shù)據(jù)計算

血氧飽和度(SpO2),即血液中氧合血紅蛋白占全部可結(jié)合的血紅蛋白的百分比.其原始計算公式為:

利用朗伯比爾定律與光電轉(zhuǎn)換原理,可將血氧飽和度公式轉(zhuǎn)換為:

其中,A、B、C為二次曲線定標(biāo)擬合系數(shù),λ1和λ2分別為紅光和紅外光波長,IAC和IDC分別脈搏波信號的交流分量和直流分量.通過計算脈搏波信號上下包絡(luò)線的差值可得交流分量,計算脈搏波信號上下包絡(luò)線的均值可得直流分量.

脈率表示每分鐘脈搏的次數(shù).通過脈搏波信號找到兩個脈搏波的同一參考點(diǎn),由此計算出兩參考點(diǎn)間采樣點(diǎn)個數(shù)n和采樣率f,即可得一個完整脈搏波時間T:

則脈率R為:

5 實(shí)驗(yàn)結(jié)果

采用光電轉(zhuǎn)換方式采集的脈搏波信號不可避免地會包含許多噪聲:包括人體自身呼吸運(yùn)動和微動脈搏動引起的低頻干擾和高頻干擾[9],以及運(yùn)動過程引入的頻率混疊干擾[10]等.

圖8為在運(yùn)動狀態(tài)下采集的原始脈搏波信號,可以看出存在較多毛刺,以及明顯的基線漂移,說明信號中存在高頻干擾、直流干擾等噪聲,這會導(dǎo)致血氧飽和度的計算出現(xiàn)較大誤差.如圖9所示,當(dāng)采用形態(tài)學(xué)濾波和雙樹復(fù)小波變換對信號進(jìn)行處理后,波形變得較為平滑,并且基線漂移也得到抑除.

圖8 原始脈搏波信號

圖9 濾波后脈搏波信號

6 測試分析

6.1 系統(tǒng)測試

采集系統(tǒng)由手套、可嵌入手套的血氧采集裝置、手機(jī)端app組成.實(shí)物圖與app界面如圖10所示.

圖10 系統(tǒng)實(shí)物圖

同時用本裝置與標(biāo)準(zhǔn)監(jiān)護(hù)儀PC-304采集一名實(shí)驗(yàn)者的體征參數(shù).其數(shù)據(jù)對比如表1和表2所示.

表1 實(shí)驗(yàn)裝置與PC-304血氧飽和度對比

表2 實(shí)驗(yàn)裝置與PC-304脈率對比 (次/min)

通過上表可以看出,與監(jiān)護(hù)設(shè)備相比,實(shí)驗(yàn)裝置計算出的血氧飽和度誤差不超過3.5%,脈率差值不超過8次,因而可以滿足實(shí)際需求.

表3 實(shí)驗(yàn)裝置功耗測試

6.2 功耗測試

用3.7V電壓輸出的鋰電池給裝置供電,測試功耗如表3所示.

工作模式下,若用310mAh的鋰電池供電,裝置續(xù)航時間接近14小時,能夠滿足運(yùn)動或消防救險的一次性供電.

7 結(jié)論

目前市面上缺乏動態(tài)環(huán)境下監(jiān)測血氧飽和度的完整系統(tǒng)方案,本文提供的設(shè)計具有較強(qiáng)的可移植性,能夠應(yīng)用于運(yùn)動和消防領(lǐng)域,對于生命安全的維護(hù)具有重要意義.裝置在保證便攜性的同時,能夠較為準(zhǔn)確地采集生理參數(shù),較低的功耗也滿足其長時間地運(yùn)行.

[1]陳潔,蹇文淵,段俊國.視網(wǎng)膜血氧飽和度測定技術(shù)研究及應(yīng)用進(jìn)展[J].西部醫(yī)學(xué),2013(11):1740-1741+1744.

[2]劉艷麗,趙為松,李海坤,汪方斌.基于形態(tài)濾波的脈搏波信號基線漂移消除方法研究[J].合肥工業(yè)大學(xué)學(xué)報,2011(04):525-528.

[3]王際忠,王志武,顏國正,劉大生.血氧測量快速濾波與實(shí)時算法設(shè)計[J].北京生物醫(yī)學(xué)工程,2015(06):595-599.

[4]KINGSBURY N G.The dual-tree complex wavelet transform:A New technique for shift invariance and directional filter[J].IEEE Digital Signal Processing Workshop,1998,98(1):2-5.

[5]Selesnick I W,Baraniuk R G,Kingsbury N C.The dual-tree complex wavelet transform[J].IEEE signal processing magazine,2005,22(6):123-151.

[6]劉文濤,陳紅,蔡曉霞,劉俊彤.基于雙樹復(fù)小波變換的信號去噪算法[J].火力與指揮控制,2014(12):84-87.

[7]文苗,王成,白麗紅,等.基于DTCWT和cICA的光電容積脈搏波運(yùn)動干擾消除算法[J].中國醫(yī)學(xué)物理雜志,2015,32(2):229-233.

[8]羅志增,周鎮(zhèn)定,周瑛,等.雙樹復(fù)小波特征在運(yùn)動想象腦電識別中的應(yīng)用[J].傳感技術(shù)學(xué)報,2014,27(5):575-580.

[9]韓慶陽,李丙玉,王曉東.一種同時消除脈搏波信號中呼吸基線漂移和高頻噪聲的方法[J].中國醫(yī)學(xué)物理學(xué)雜志,2014(02):4801-4805.

[10]龔渝順,吳寶明,高丹丹,蔣洪,閆慶廣.動態(tài)環(huán)境血氧飽和度監(jiān)測的運(yùn)動干擾分離自適應(yīng)對消方法[J].航天醫(yī)學(xué)與醫(yī)學(xué)工程,2012(04):266-270.

Design of Wearable Anti-interference Oxygen Saturation Acquisition Device

Wu You,Pang Yu,Huang Boqiang,Peng Liangguang
(Chongqing University of Posts and Telecommunications Academy of Opto Engineering,Chongqing 400065)

Due to the lack of timely acquisition of physiological parameters in the fire rescue operation and long-time aerobic exercise,many people lost their life.For the reason,a wearable blood oxygen saturation acquisition system based on glove was designed.The device used reflection probe to acquire photoplethysmography signal and STM32L151 microprocessor to process data,and used the dual-tree complex wavelet transform(DTCWT)and morphological filtering algorithm to remove the noise and baseline drift,and finally sent the data to mobile phone app through Bluetooth.By the test,the device is small and easy to wear,and has low power consumption and anti-motion-interference ability,and compared with standard monitor,the error of blood oxygen saturation is within 3.5%.

wearable;blood oxygen saturation;DTCWT;morphological filtering;anti-motion-interference

國家自然科學(xué)基金(61671091).本項(xiàng)目受重慶市高校創(chuàng)新團(tuán)隊(duì)(智慧醫(yī)療與系統(tǒng)核心技術(shù)),重慶市研究生科研創(chuàng)新項(xiàng)目(CYS17239)資助.

吳優(yōu)(1992-),男,碩士研究生,主要研究方向:體征信號處理與可穿戴設(shè)備研究.

黃博強(qiáng)(1992-),男,碩士研究生,主要研究方向:生命體征信號采集與處理.

彭良廣(1993-),男,碩士研究生,主要研究方向:生命體征信號采集與處理.