重慶郵電大學(xué)光電工程學(xué)院 吳 優(yōu) 龐 宇 黃博強(qiáng) 彭良廣
穿戴式抗干擾血氧飽和度采集裝置設(shè)計(jì)
重慶郵電大學(xué)光電工程學(xué)院 吳 優(yōu) 龐 宇 黃博強(qiáng) 彭良廣
在消防營(yíng)救和長(zhǎng)時(shí)間有氧運(yùn)動(dòng)中,由于缺乏對(duì)生理參數(shù)的及時(shí)獲取,時(shí)常造成嚴(yán)重的人員傷亡.為此,本文設(shè)計(jì)了一種以手套為載體的穿戴式血氧飽和度采集系統(tǒng),利用反射式探頭采集脈搏波信號(hào),經(jīng)由STM32L151單片機(jī)進(jìn)行數(shù)據(jù)處理,并引入雙樹復(fù)小波變換和形態(tài)學(xué)濾波去除噪聲干擾和基線漂移,最后通過(guò)藍(lán)牙發(fā)送至手機(jī)端app獲取生理參數(shù).經(jīng)測(cè)試,裝置體積微小,佩戴方便,具備低功耗與抗運(yùn)動(dòng)干擾能力,并且和標(biāo)準(zhǔn)監(jiān)護(hù)儀相比血氧飽和度平均誤差不超過(guò)2%.
穿戴式;血氧飽和度;雙樹復(fù)小波變換;形態(tài)學(xué)濾波;抗運(yùn)動(dòng)干擾
血氧飽和度作為一項(xiàng)重要的生理參數(shù),它間接反映出人體細(xì)胞的新陳代謝狀態(tài).
據(jù)統(tǒng)計(jì),我國(guó)每年有近30名消防員犧牲,超過(guò)300人受傷致殘;在運(yùn)動(dòng)領(lǐng)域,每年也有近150人發(fā)生心源性猝死.由于高危高壓的消防營(yíng)救與過(guò)勞性的"有氧運(yùn)動(dòng)"時(shí)常使身體處于缺氧狀態(tài),若長(zhǎng)時(shí)間得不到恢復(fù)則會(huì)對(duì)生命造成嚴(yán)重危害.
出于對(duì)生命安全的考慮,消防和運(yùn)動(dòng)領(lǐng)域?qū)τ谘躏柡投鹊膶?shí)時(shí)檢測(cè)尤為重要.目前,市面上存在的血氧飽和度檢測(cè)裝置主要運(yùn)用于醫(yī)院監(jiān)護(hù)和家庭監(jiān)護(hù),在靜態(tài)環(huán)境下具有較高的準(zhǔn)確性,但在消防施救和運(yùn)動(dòng)過(guò)程中則存在很大局限.
考慮到裝置的便攜性、實(shí)時(shí)性、續(xù)航時(shí)間,以及抗運(yùn)動(dòng)干擾能力,本文提出了一種新的解決方案:通過(guò)佩戴內(nèi)嵌微型裝置的手套即可完成血氧飽和度和脈率的采集,并利用藍(lán)牙將數(shù)據(jù)發(fā)送至手機(jī)端,再通過(guò)手機(jī)推送至后臺(tái),從而達(dá)到檢測(cè)與反饋的目的.
裝置采用透射式血氧探頭采集脈搏波信號(hào),采集點(diǎn)為指尖部位.探頭內(nèi)置于一個(gè)能完全裹住手指第一關(guān)節(jié)的橡膠套中,由于橡膠材質(zhì)的伸縮性,因而可以使探頭很好的貼合于皮膚表面,而不會(huì)因?yàn)槭种傅幕顒?dòng)致使探頭移位.微型電路模塊內(nèi)藏于手套背部口袋,由于體積微小,佩戴手套后不會(huì)對(duì)消防和運(yùn)動(dòng)人員造成任何不便和影響.系統(tǒng)總設(shè)計(jì)框圖如圖1所示.
圖1 系統(tǒng)總設(shè)計(jì)框圖
裝置硬件部分主要分為電源管理模塊、脈搏波信號(hào)采集模塊、單片機(jī)控制和處理模塊、藍(lán)牙模塊.硬件整體框圖如圖2所示.
本設(shè)計(jì)選擇低功耗微處理器STM32L151CBT6,其運(yùn)用了ST自主研發(fā)的130nm的超低漏電工藝技術(shù),具備先進(jìn)的超低功耗模式、優(yōu)化的動(dòng)態(tài)運(yùn)行功耗以及特殊的安全性能,在高性能和超低功耗之間取得了很好的平衡.控制電路如圖3所示.
圖2 硬件整體框圖
圖3 單片機(jī)控制電路
單片機(jī)系統(tǒng)時(shí)鐘采用HSI,以16Mhz頻率運(yùn)行;程序下載方式為SWD;引腳red_on、ir_on、LEDir、LEDred與脈搏波信號(hào)采集電路相連,在通用定時(shí)器的時(shí)序下控制紅光和紅外光的交替產(chǎn)生;引腳OUT用于捕獲每一次光強(qiáng)探測(cè)器傳回的原始脈搏波信號(hào);RX與TX則為UART串口收發(fā)引腳,負(fù)責(zé)將處理完的數(shù)據(jù)傳輸給藍(lán)牙模組.
脈搏波采集電路由兩部分組成:一是由三極管和電阻構(gòu)成的驅(qū)動(dòng)電路,二是由二極管和光強(qiáng)探測(cè)器組成的信號(hào)采集電路.如圖4所示.
圖4 脈搏波信號(hào)采集電路
三極管的型號(hào)為S9012與S9013,起電流放大的作用,以驅(qū)動(dòng)二極管發(fā)光.
考慮到氧合血紅蛋白(oxyhemoglobin,HbO2)與脫氧血紅蛋白(deoxyhemoglobin,Hb)對(duì)不同波長(zhǎng)的光吸收度不同,而在波長(zhǎng)660nm與940nm處區(qū)別最大[1],因此選擇既能發(fā)出紅光(波長(zhǎng)660nm)又能發(fā)出紅外光(波長(zhǎng)940nm)的雙光源led,其由兩個(gè)極性相反的led組合封裝而成,具有兩個(gè)引腳,型號(hào)為EZDL-R660/IR905.
光強(qiáng)探測(cè)器為TSL235R型光頻轉(zhuǎn)換傳感器,其頻率響應(yīng)范圍為320nm至1050nm,可以將接收到的光轉(zhuǎn)換成一定頻率的數(shù)字信號(hào),由輸出管腳發(fā)送給單片機(jī)捕獲和處理.
軟件設(shè)計(jì)包括脈搏波信號(hào)的采集、傳輸、處理、計(jì)算和發(fā)送.其流程如圖5所示.
圖5 軟件流程圖
程序首先初始化各模塊,包括配置系統(tǒng)與各外設(shè)時(shí)鐘;配置I/O口的輸入、輸出、復(fù)用等模式;配置定時(shí)器及其中斷模式;配置ADC與DMA傳輸方式;初始化UART串口等.
led時(shí)序由定時(shí)器控制,以12ms為1個(gè)周期采集兩組脈搏波信號(hào),其中紅光與紅外光分別點(diǎn)亮2ms,其余時(shí)間處于熄滅狀態(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ù)代入公式計(jì)算出血氧飽和度和脈率.
UART與藍(lán)牙模組相連,當(dāng)單片機(jī)計(jì)算完成后即通過(guò)串口將數(shù)據(jù)傳輸給藍(lán)牙,藍(lán)牙再將數(shù)據(jù)發(fā)送給手機(jī)端顯示出血氧飽和度和脈率值,最后推送至后臺(tái).
電池電量檢測(cè)運(yùn)用到A/D轉(zhuǎn)換,當(dāng)檢測(cè)出電量低于預(yù)設(shè)值時(shí),利用串口向手機(jī)發(fā)送低電量提醒,一定時(shí)間后采取自動(dòng)關(guān)機(jī).
4.1.1 形態(tài)學(xué)濾波
正常情況下信號(hào)的中心線為直線,但實(shí)際波形里中心線往往隨時(shí)間上下擺動(dòng),稱其為基線漂移.基線漂移的實(shí)質(zhì)是信號(hào)直流分量不穩(wěn)定,主要由人體自身的呼吸運(yùn)動(dòng),以及探測(cè)器與皮膚之間接觸位移所致[2].采用形態(tài)學(xué)濾波可以消除基線漂移.
形態(tài)學(xué)濾波是一種非線性濾波技術(shù),算法核心為通過(guò)腐蝕、膨脹、形態(tài)開和形態(tài)閉等運(yùn)算[3],濾除脈搏波信號(hào)特征波形,再利用原始信號(hào)減去剩下的基線漂移信號(hào),從而得到直流分量穩(wěn)定的正常信號(hào).
腐蝕是一種使邊界向內(nèi)收縮的過(guò)程,而膨脹則使邊界向外擴(kuò)張.信號(hào)f關(guān)于結(jié)構(gòu)元素g的腐蝕、膨脹運(yùn)算分別如下:
若對(duì)信號(hào)先腐蝕,再膨脹,則為開運(yùn)算,公式如下:
若對(duì)信號(hào)先膨脹,再腐蝕,則為閉運(yùn)算,公式如下:
開運(yùn)算能夠去除毛刺與孤立的小點(diǎn),閉運(yùn)算則能填合孔洞與彌補(bǔ)縫隙.形態(tài)學(xué)濾波通過(guò)形態(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)信號(hào),左子樹L對(duì)應(yīng)實(shí)部,右子樹R對(duì)應(yīng)虛部,分別包含一組低通和高通濾波器,兩組濾波器可構(gòu)成Hilbert變換對(duì),使雙樹復(fù)小波變換近似解析[7].變換原理如圖7所示.
圖7 DTCWT原理圖
實(shí)樹與虛樹對(duì)應(yīng)濾波器之間的延遲相差一個(gè)采樣間隔,虛樹的采樣位置始終位于實(shí)樹的中間,以使雙樹抽取所得數(shù)據(jù)成互補(bǔ)關(guān)系,從而降低信息的丟失,實(shí)現(xiàn)近似平移不變性[8].
雙樹復(fù)小波變換通過(guò)分解脈搏波信號(hào)各層的低頻小波系數(shù),再利用閾值去除高頻小波系數(shù)中所含的噪聲成分,然后將分解的低頻小波系數(shù)和閾值處理后的高頻小波系數(shù)進(jìn)行信號(hào)重構(gòu),從而得到去除噪聲的光電容積脈搏波信號(hào).
血氧飽和度(SpO2),即血液中氧合血紅蛋白占全部可結(jié)合的血紅蛋白的百分比.其原始計(jì)算公式為:
利用朗伯比爾定律與光電轉(zhuǎn)換原理,可將血氧飽和度公式轉(zhuǎn)換為:
其中,A、B、C為二次曲線定標(biāo)擬合系數(shù),λ1和λ2分別為紅光和紅外光波長(zhǎng),IAC和IDC分別脈搏波信號(hào)的交流分量和直流分量.通過(guò)計(jì)算脈搏波信號(hào)上下包絡(luò)線的差值可得交流分量,計(jì)算脈搏波信號(hào)上下包絡(luò)線的均值可得直流分量.
脈率表示每分鐘脈搏的次數(shù).通過(guò)脈搏波信號(hào)找到兩個(gè)脈搏波的同一參考點(diǎn),由此計(jì)算出兩參考點(diǎn)間采樣點(diǎn)個(gè)數(shù)n和采樣率f,即可得一個(gè)完整脈搏波時(shí)間T:
則脈率R為:
采用光電轉(zhuǎn)換方式采集的脈搏波信號(hào)不可避免地會(huì)包含許多噪聲:包括人體自身呼吸運(yùn)動(dòng)和微動(dòng)脈搏動(dòng)引起的低頻干擾和高頻干擾[9],以及運(yùn)動(dòng)過(guò)程引入的頻率混疊干擾[10]等.
圖8為在運(yùn)動(dòng)狀態(tài)下采集的原始脈搏波信號(hào),可以看出存在較多毛刺,以及明顯的基線漂移,說(shuō)明信號(hào)中存在高頻干擾、直流干擾等噪聲,這會(huì)導(dǎo)致血氧飽和度的計(jì)算出現(xiàn)較大誤差.如圖9所示,當(dāng)采用形態(tài)學(xué)濾波和雙樹復(fù)小波變換對(duì)信號(hào)進(jìn)行處理后,波形變得較為平滑,并且基線漂移也得到抑除.
圖8 原始脈搏波信號(hào)
圖9 濾波后脈搏波信號(hào)
采集系統(tǒng)由手套、可嵌入手套的血氧采集裝置、手機(jī)端app組成.實(shí)物圖與app界面如圖10所示.
圖10 系統(tǒng)實(shí)物圖
同時(shí)用本裝置與標(biāo)準(zhǔn)監(jiān)護(hù)儀PC-304采集一名實(shí)驗(yàn)者的體征參數(shù).其數(shù)據(jù)對(duì)比如表1和表2所示.
表1 實(shí)驗(yàn)裝置與PC-304血氧飽和度對(duì)比
表2 實(shí)驗(yàn)裝置與PC-304脈率對(duì)比 (次/min)
通過(guò)上表可以看出,與監(jiān)護(hù)設(shè)備相比,實(shí)驗(yàn)裝置計(jì)算出的血氧飽和度誤差不超過(guò)3.5%,脈率差值不超過(guò)8次,因而可以滿足實(shí)際需求.
表3 實(shí)驗(yàn)裝置功耗測(cè)試
用3.7V電壓輸出的鋰電池給裝置供電,測(cè)試功耗如表3所示.
工作模式下,若用310mAh的鋰電池供電,裝置續(xù)航時(shí)間接近14小時(shí),能夠滿足運(yùn)動(dòng)或消防救險(xiǎn)的一次性供電.
目前市面上缺乏動(dòng)態(tài)環(huán)境下監(jiān)測(cè)血氧飽和度的完整系統(tǒng)方案,本文提供的設(shè)計(jì)具有較強(qiáng)的可移植性,能夠應(yīng)用于運(yùn)動(dòng)和消防領(lǐng)域,對(duì)于生命安全的維護(hù)具有重要意義.裝置在保證便攜性的同時(shí),能夠較為準(zhǔn)確地采集生理參數(shù),較低的功耗也滿足其長(zhǎng)時(shí)間地運(yùn)行.
[1]陳潔,蹇文淵,段俊國(guó).視網(wǎng)膜血氧飽和度測(cè)定技術(shù)研究及應(yīng)用進(jìn)展[J].西部醫(yī)學(xué),2013(11):1740-1741+1744.
[2]劉艷麗,趙為松,李海坤,汪方斌.基于形態(tài)濾波的脈搏波信號(hào)基線漂移消除方法研究[J].合肥工業(yè)大學(xué)學(xué)報(bào),2011(04):525-528.
[3]王際忠,王志武,顏國(guó)正,劉大生.血氧測(cè)量快速濾波與實(shí)時(shí)算法設(shè)計(jì)[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ù)小波變換的信號(hào)去噪算法[J].火力與指揮控制,2014(12):84-87.
[7]文苗,王成,白麗紅,等.基于DTCWT和cICA的光電容積脈搏波運(yùn)動(dòng)干擾消除算法[J].中國(guó)醫(yī)學(xué)物理雜志,2015,32(2):229-233.
[8]羅志增,周鎮(zhèn)定,周瑛,等.雙樹復(fù)小波特征在運(yùn)動(dòng)想象腦電識(shí)別中的應(yīng)用[J].傳感技術(shù)學(xué)報(bào),2014,27(5):575-580.
[9]韓慶陽(yáng),李丙玉,王曉東.一種同時(shí)消除脈搏波信號(hào)中呼吸基線漂移和高頻噪聲的方法[J].中國(guó)醫(yī)學(xué)物理學(xué)雜志,2014(02):4801-4805.
[10]龔渝順,吳寶明,高丹丹,蔣洪,閆慶廣.動(dòng)態(tài)環(huán)境血氧飽和度監(jiān)測(cè)的運(yùn)動(dòng)干擾分離自適應(yīng)對(duì)消方法[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
國(guó)家自然科學(xué)基金(61671091).本項(xiàng)目受重慶市高校創(chuàng)新團(tuán)隊(duì)(智慧醫(yī)療與系統(tǒng)核心技術(shù)),重慶市研究生科研創(chuàng)新項(xiàng)目(CYS17239)資助.
吳優(yōu)(1992-),男,碩士研究生,主要研究方向:體征信號(hào)處理與可穿戴設(shè)備研究.
黃博強(qiáng)(1992-),男,碩士研究生,主要研究方向:生命體征信號(hào)采集與處理.
彭良廣(1993-),男,碩士研究生,主要研究方向:生命體征信號(hào)采集與處理.