包 穎, 單新治, 王冠學(xué), 洪瑞金, 張大偉, 高秀敏

(上海理工大學(xué) 光電信息與計算機(jī)工程學(xué)院, 上海 200093)

一種多信息融合無創(chuàng)血糖檢測數(shù)據(jù)處理方法研究

包 穎, 單新治, 王冠學(xué), 洪瑞金, 張大偉, 高秀敏

(上海理工大學(xué) 光電信息與計算機(jī)工程學(xué)院, 上海 200093)

以生物特征檢測技術(shù)為基礎(chǔ),基于已經(jīng)構(gòu)建的無創(chuàng)血糖檢測系統(tǒng)進(jìn)行血糖檢測數(shù)據(jù)分析。以光電檢測技術(shù)為主,結(jié)合其他人體特征信息檢測,對人體進(jìn)行無創(chuàng)血糖測試得到大量數(shù)據(jù),建立數(shù)學(xué)模型并利用算法來計算血糖濃度。與一種常用有創(chuàng)血糖儀檢測數(shù)據(jù)進(jìn)行對比,有創(chuàng)與無創(chuàng)檢測結(jié)果的相關(guān)系數(shù)達(dá)到了85.4%,為多信息融合無創(chuàng)血糖檢測系統(tǒng)提供了一種有效的數(shù)據(jù)處理方法。

多信息融合; 無創(chuàng); 血糖檢測

引 言

從1980年至今,我國糖尿病患者呈現(xiàn)逐年遞增趨勢,到2010年全國共有糖尿病患者約1.13億人,預(yù)計到2020年將超過2億人,目前我國糖尿病患者診斷率不足40%。由這些數(shù)據(jù)可知,血糖值的檢測對糖尿病患者和處于亞健康的群體有著重要的意義。目前,血糖檢測的主要方法還是傳統(tǒng)的針扎式取血測量,但傳統(tǒng)扎針檢測法具有較多缺點:讓受測者有疼痛感,體驗較差;每次都需要一條試紙,儀器更換電池頻繁,測量成本較高;需要有扎針、采血、測量等繁瑣步驟,操作不便;試紙保存需要特定條件,要與測試機(jī)器匹配。

目前自我血糖濃度監(jiān)測一般是采用快速血糖儀進(jìn)行測量,由于該方法需要采集血樣,容易引起患者的生理疼痛并伴有感染的危險,因此這種有創(chuàng)的測量方法在很大程度上限制了血糖測量的頻率。

國內(nèi)外無創(chuàng)血糖測量方法主要有以下幾類:測量皮下滲出組織液中的血糖濃度、微波無創(chuàng)血糖檢測、皮下植入傳感器、人體的射頻阻抗無創(chuàng)檢測血糖值、利用唾液進(jìn)行無創(chuàng)血糖檢測、超聲波血糖檢測儀、聚光斷層攝影、旋光法、激光拉曼光譜法、紅外光譜法[1]。其中光學(xué)和輻射方法中的紅外光譜法是目前研究比較多的一種方法,但這種技術(shù)單純使用紅外技術(shù)來檢測血糖值,由于數(shù)據(jù)單一,同時測量環(huán)境對紅外線的干擾較大,導(dǎo)致單純的紅外光無法準(zhǔn)確穩(wěn)定地檢測血糖值。

從國內(nèi)外的公開研究資料和報道來看,基于多波長的無創(chuàng)血糖檢測具有良好的應(yīng)用前景,且在重復(fù)性和穩(wěn)定性等方面還有很多值得我們?nèi)パ芯?本文對一種多信息融合無創(chuàng)血糖檢測數(shù)據(jù)處理方法進(jìn)行研究。

1 無創(chuàng)血糖檢測參數(shù)選取

如何采用無創(chuàng)血糖測試方法準(zhǔn)確測量血糖值一直受到人們的廣泛關(guān)注(包括測量那些參數(shù)來標(biāo)定血糖值,以及如何準(zhǔn)確測量各種參數(shù)的方法)[2-4]。葡萄糖在被人體消耗的過程中,由于會消耗血液中氧氣的同時還生成相應(yīng)的熱量,因此可以采用能量守恒的辦法間接測量血糖值。其中代謝產(chǎn)生能量的主要去向是與外界發(fā)生熱交換[4-5],通過可見光和近紅外雙波段與人體血液中相互作用,得出一個參數(shù)值,用于血糖指標(biāo)檢測。目前,無創(chuàng)血糖檢測參數(shù)選取主要有血氧飽和度、心率(血液流速)、被測部位體表溫度、被測部位熱輻射溫度及環(huán)境溫度。血氧飽和度(SpO2)表示血液中血氧的濃度,它是呼吸循環(huán)的重要生理參數(shù),如果血糖異常影響呼吸功能的話,會出現(xiàn)血氧飽和度偏低的情況;心跳快慢與耗氧速度直接有關(guān);血糖代謝會釋放能量,產(chǎn)生熱量,使體表有溫度;體溫會受到環(huán)境的影響,測量環(huán)境溫度可降低干擾;光電檢測中的紅光和紅外光有輻射,因此還需考慮環(huán)境輻射溫度,被測部位輻射溫度。因此選取這些參數(shù)作為血糖濃度的數(shù)據(jù)來源,并對這些數(shù)據(jù)進(jìn)行多信息融合算法探討及處理。

2 數(shù)據(jù)預(yù)處理

多信息融和算法主要涉及數(shù)據(jù)預(yù)處理和模型建立。對采集的數(shù)據(jù)通過MATLAB進(jìn)行處理,以剔除原始數(shù)據(jù)中異常點造成結(jié)果不精確的影響。

由文獻(xiàn)[6-7]可知,選定的血氧飽和度指標(biāo)與紅外光(ir)和紅光(red)透射率有關(guān),其中利用氧合血紅蛋白(HbO2)和血紅蛋白(Hb)對紅光和紅外光吸收的差異,可以得出兩束光的吸收比率,即特征值

R=ACred·DCir/(ACir·DCred)

(1)

式中:ACred和DCred分別為red透射光的交流成分和直流成分;ACir和DCir分別為ir透射光的交流成分和直流成分。SpO2血氧飽和度作為應(yīng)變量與自變量R滿足一元二次方程關(guān)系,因此我們通過采集紅外光和紅光經(jīng)血樣后的透射能量以獲得R值。

在數(shù)據(jù)采集后,發(fā)現(xiàn)每個樣本的紅外光透射能量在29 000個數(shù)據(jù)點中前5 000個和后5 000個數(shù)據(jù)都存在較大的異常波動。下面對MATLAB導(dǎo)出的一個樣本的紅外光透射圖進(jìn)行分析。

圖1為某一樣本紅外透射能量圖,由圖可以看到29 000個數(shù)據(jù)點中,中段(5 500～24 000)數(shù)據(jù)較為平穩(wěn),前后5 000個數(shù)據(jù)點都出現(xiàn)較大范圍的波動,由于同一個樣本數(shù)據(jù)采集的是同一人,因此數(shù)據(jù)點偏差不應(yīng)該過大。產(chǎn)生波動的原因是測量開始及結(jié)束時,手指在傳感器中在一個相對運(yùn)動的狀態(tài),此時記錄數(shù)據(jù)為非正常工作狀態(tài)。所以,我們將輸入紅外光透射的前后5 000個數(shù)據(jù)點進(jìn)行刪除,只截取中段數(shù)據(jù),最后得到有效紅外光透射能量圖,如圖2所示。

圖1 指尖紅外光透射能量圖(全部數(shù)據(jù)點)Fig.1 The infrared transmission energy diagram of the fingers(all data points)

圖2 指尖紅外光透射能量(有效數(shù)據(jù)點)Fig.2 The infrared transmission energy diagram of the fingers(valid data points)

由圖2可以看出,除了個別異常點外,相對圖1,此時的數(shù)據(jù)點波動顯然較小。

我們同時也將同一個樣本的紅光透射譜用MATLAB導(dǎo)出,得到了圖3所示的手指透射red能量的全部數(shù)據(jù)。在圖3的29 000個數(shù)據(jù)點中,中段(5 500～24 000)數(shù)據(jù)較為平穩(wěn),前后5 000的數(shù)據(jù)點都出現(xiàn)較大范圍的波動,由于同一個樣本數(shù)據(jù)采集的是同一人,因此數(shù)據(jù)點偏差不應(yīng)該過大。因此我們將前后5 000個數(shù)據(jù)點進(jìn)行刪除,只截取中段數(shù)據(jù),最后得到圖4所示的有效紅光透射能量圖。由圖4可以看出,除了個別異常點外,相對圖3,此時的透射能量波動顯然較小,僅在(1±0.025)×106間波動。

圖3 指尖紅光透射能量(全部數(shù)據(jù)點)Fig.3 The red light transmission energy diagram of the fingers(all data points)

圖4 指尖紅光透射能量(有效數(shù)據(jù)點)Fig.4 The red light transmission energy diagram of the fingers(valid data points)

現(xiàn)在我們將有效數(shù)據(jù)點的紅外光和紅光透射能量繪制在同一張圖中,如圖5所示。圖中透射能量仍然存在一定的波動,以下對其進(jìn)一步處理。

首先標(biāo)定紅外光與紅光透射能量的峰值點,圖6為紅光和紅外光峰值點標(biāo)定圖。

圖5 紅外光和紅光透射能量譜圖(有效數(shù)據(jù)點)Fig.5 The infrared and red light transmission energy diagram of the fingers(valid data points)

圖6 紅光、紅外光透射能量峰值點標(biāo)定(有效數(shù)據(jù)點)Fig.6 The peaks of the red light and infrared transmission energy calibration(valid data points)

由圖6可以看出,前后兩個數(shù)據(jù)點不如中段數(shù)據(jù)穩(wěn)定,因此,我們同樣剔除前后兩個峰值點。接著用差分函數(shù)求相鄰兩個峰值點能量差求平均,以及計算最大峰值與所有數(shù)據(jù)點透射能量的偏差后也用差分函數(shù)求相鄰能量差求平均。最終的波動偏差由上述兩種方法得出的結(jié)果取平均,此時能量誤差達(dá)到了4.5%左右。于是在對整體數(shù)據(jù)取平均時,不會存在較大誤差。

我們將測量的指尖溫度作為體表溫度。由于體溫會受到環(huán)境的影響,所以還需測量環(huán)境溫度,圖7是所測的體表溫度與環(huán)境溫度曲線。

體表溫度應(yīng)為此時人體指尖溫度的一個測量值,通過圖7我們可以發(fā)現(xiàn),在最后10個點時溫度趨向穩(wěn)定,因此我們體表溫度用最后10個點的均值表示。環(huán)境溫度波動情況較小,我們對其整段取均值。

與體表溫度、環(huán)境溫度分析同理,測得的體表輻射能量和環(huán)境輻射能量如圖8所示。

圖7 體表溫度與環(huán)境溫度Fig.7 Surface temperature and environmental temperature

圖8 體表輻射能量與環(huán)境輻射能量Fig.8 The radiation energy of body and environment

體表輻射能量是人體指尖輻射能量的一個測量值,通過圖8我們可以發(fā)現(xiàn),在最后10個點時輻射能量趨向穩(wěn)定,因此我們體表輻射能量用最后10個點的均值表示。環(huán)境輻射能量波動情況較小,我們對其整段取均值。

3 多信息融合無創(chuàng)血糖檢測數(shù)據(jù)處理方法

針對預(yù)處理后的血氧飽和度、心率、被測部位體表溫度、被測部位熱輻射溫度、環(huán)境溫度建立數(shù)學(xué)模型。由于變量之間存在相關(guān)性,因此我們將對變量進(jìn)行相關(guān)性分析,之后進(jìn)行主成分分析提取得分向量將結(jié)構(gòu)簡化[8],最終建模。

首先將這些參數(shù)和傳統(tǒng)有創(chuàng)血糖儀檢測數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化:

(1)

式中:m為自變量;n為數(shù)據(jù)點點數(shù);E0為標(biāo)準(zhǔn)化后的自變量(選取和血糖有關(guān)的6個參數(shù))矩陣;F0為標(biāo)準(zhǔn)化后的因變量(血糖)矩陣;xj1(j=1,2,…,n)為血氧飽和度;xj2(j=1,2,…,n)為心率;xj3(j=1,2,…,n)為被測部位體表溫度;xj4(j=1,2,…,n)為被測部位熱輻射溫度;xj5(j=1,2,…,n)為環(huán)境溫度;xj6(j=1,2,…,n)為環(huán)境輻射溫度。然后進(jìn)行以下5個步驟。

(1) 計算E0第一成分的得分向量,記為t1:

(2)

式中w1為權(quán)重,w1=[w11…w1m]T。

t1和F0的協(xié)方差Cov(t1,F0)可用t1和F0的內(nèi)積計算,因而求解可化為Lagrange乘數(shù)法求條件極值問題,于是有

(3)

(2) 假定回歸模型為

(4)

式中:E1和F1為殘差矩陣;α1、β1為模型效應(yīng)負(fù)荷量,它們的估計值為

(5)

(3) 用E1和F1代替E0和F0重復(fù)(2)步驟。

(4) 設(shè)r(E0)

(6)

即得多信息融合算法表達(dá)式

y=a1x1+…+amxm

(7)

(5) 進(jìn)行交叉有效性檢驗。

多信息融合算法可像主成分分析一樣,僅通過前l(fā)個成分(l≤r),即可得到較好預(yù)測能力的回歸模型?？赏ㄟ^交叉有效性檢驗來確定提取的主成分個數(shù)l。

(8)

(9)

4 測試結(jié)果與分析

為驗證多信息無創(chuàng)血糖檢測數(shù)據(jù)處理方法的正確性,在附近醫(yī)院對20～30歲之間的18名健康志愿者分別進(jìn)行了測試。被測者在測量前一天晚上保證良好的睡眠,在早上8:30左右接受測試。測量時,要求被測者手指潔凈干燥,保持靜坐姿勢,心情平靜,被測手指無大幅運(yùn)動,通過醫(yī)用指夾將傳感器固定在手指上,并保持手掌平直。用無創(chuàng)血糖儀檢測各項特征信息并經(jīng)數(shù)據(jù)處理得到血糖濃度,同時還采用醫(yī)院的血液分析儀測量18名志愿者的血糖濃度并將其作為標(biāo)準(zhǔn)值,測試數(shù)據(jù)如表1所示。

表1 有創(chuàng)測量和無創(chuàng)測量測試數(shù)據(jù)Tab.1 Measured data of invasive and noninvasive measurement

圖9為無創(chuàng)測量與有創(chuàng)測量結(jié)果對比曲線。

通過對健康人體的測試發(fā)現(xiàn),利用本文所建數(shù)據(jù)處理模型得到無創(chuàng)血糖測量結(jié)果與傳統(tǒng)的有創(chuàng)測量結(jié)果之間的相關(guān)系數(shù)達(dá)到了85.4%,同時驗證了血糖濃度與所選的幾項人體特征信息具有相關(guān)性,無創(chuàng)血糖檢測系統(tǒng)可有效反映血糖濃度的動態(tài)變化趨勢。

圖9 無創(chuàng)測量與有創(chuàng)測量結(jié)果對比曲線Fig.9 Contrast of the non-invasive measurement and invasive measurements

5 結(jié) 論

本文以生物特征檢測技術(shù)為基礎(chǔ),研究了一種無創(chuàng)血糖檢測數(shù)據(jù)處理方法,經(jīng)無創(chuàng)血糖檢測與傳統(tǒng)有創(chuàng)血糖儀檢測對比,同一測試者兩種檢測結(jié)果的相關(guān)系數(shù)達(dá)到了85.4%。目前無創(chuàng)血糖技術(shù)還不成熟還需要較長的研究,后續(xù)將會繼續(xù)對算法進(jìn)行研究,使測量準(zhǔn)確性進(jìn)一步提高。

[1] 焦繪春.無創(chuàng)血糖檢測方法研究[D].哈爾濱:哈爾濱工程大學(xué),2004.

[2] 孫斐,孔德義,梅濤,等.無創(chuàng)傷法測量血糖濃度的技術(shù)及其發(fā)展前景[J].生物醫(yī)學(xué)工程學(xué)雜志,2005,22(1):171-174.

[3] 陳文亮,徐可欣,杜振輝,等.人體無創(chuàng)血糖檢測技術(shù)[J].儀器儀表學(xué)報,2003,24(S1):258-261,265.

[4] 劉志宏.人體無創(chuàng)血糖檢測方法的研究[J].中國科技信息,2006(2):91.

[5] 張紅艷.近紅外光譜技術(shù)在人體血糖無創(chuàng)檢測中的應(yīng)用研究[D].長春:中國科學(xué)院長春光學(xué)與物理研究所,2005.

[6] 劉創(chuàng).近紅外無創(chuàng)血糖檢測校正模型研究[D].西安:西北工業(yè)大學(xué)生物醫(yī)學(xué)工程,2006.

[7] 劉慶珍,徐可欣,蔣誠志.近紅外光譜技術(shù)無創(chuàng)測量人體血糖口[J].激光生物學(xué)報,2004,13(2):129-135.

[8] 姜啟源.數(shù)學(xué)模型[M].第3版.北京:高等教育出版社,2008:294-308.

Aresearchondataprocessingmethodfornon-invasivebloodglucosebasedondetectioninformationfusion

BAO Ying, SHAN Xinzhi, WANG Guanxue, HONG Ruijin, ZHANG Dawei, GAO Xiumin

(School of Optical-Electrical and Computer Engineering, University of Shanghai for Science and Technology, Shanghai 200093, China)

This paper takes the biometric detection technology as a foundation to analyze the data of blood sugar examination basing on the constructed no blood sugar detection system. Taking the photoelectric detection technology as the primary position,combined with other features of human body information detection,a large amount of data is obtained by the noninvasive testing to the human body. The mathematics model is set up to calculate the blood sugar concentration. Compared with a common information of the blood glucose meter traumatic inspection,the correlation coefficient of the test results for one tester achieves 0.854,which provides a kind of prospective data-processing method for the systems of non-invasive blood glucose detection with multi-information fusion.

multi-information fusion; non-invasive; blood glucose detection

1005－5630（2017）05－0001－07

2017-05-28

國家自然科學(xué)基金項目(61378035)

包 穎(1992—),女,碩士研究生,主要從事光電檢測系統(tǒng)設(shè)計方面的研究。E-mail:827215334@qq.com

高秀敏(1978—),男,研究員,主要從事光學(xué)性質(zhì)可調(diào)控性、光譜檢測技術(shù)、光學(xué)成像等方面的研究。E-mail:897188391@qq.com

TN 216

A

10.3969/j.issn.1005-5630.2017.05.001

(編輯:劉鐵英)