高偉翔 李 擎* 楊玉竹
1(北京信息科技大學自動化學院 北京 100192)
2(北京科技大學自動化學院 北京 100083)
從近些年的情況來看,通過一些技術(shù),人們可以自行檢測身體的心率等生理參數(shù),在發(fā)生異常時及時就醫(yī)便可大概率抑制病發(fā)。關(guān)于衡量心臟健康程度,心率無疑是非常重要的生理指標參數(shù),它是否穩(wěn)定標志著被試者心臟的功能的優(yōu)良程度,同時它也是診斷和治療疾病最重要的依據(jù)之一。
醫(yī)學行業(yè)里往往通過相關(guān)測量儀器去測量心率,比如指夾式脈搏血氧儀等,但是大部分都要和人的皮膚接觸,可能會因為被試者的心理波動造成數(shù)據(jù)不準確,也可能會因為被試者的皮膚問題導致皮膚感染等病變現(xiàn)象。此外,在一些特殊時期,比如傳染性較高的疫情期間,人為接觸式測量可能會涉及對人力資源的要求,以及為確保不傳染,對接觸式測量儀器的無菌性要求,對于測量與被測量人員的安全性要求。
與傳統(tǒng)心率檢測設(shè)備相比,借助攝像機進行非接觸式和網(wǎng)絡(luò)在線的心率檢測,既能適應(yīng)疫情期間快速檢測以及無法或不想進入醫(yī)院等情況的人體參數(shù)檢測,又能適應(yīng)大數(shù)據(jù)分析。
近些年來人們提出用非接觸的心率測量方式去檢測這一指標。這種心率測量方式的核心在于遠程光容積描記技術(shù)(rPPG)。這是一種無接觸測量心率的方法,經(jīng)過相機設(shè)備提取到面部PPG信號的周期性變化,再檢測血液容積的周期性變化,最終得以測得心率數(shù)據(jù),此方法常用在遠程醫(yī)療等領(lǐng)域。
Haan等[4]提出了基于色度模型的心率檢測方法,大大改善了光照對檢測到的心率的結(jié)果的影響。
Wang等[5]提出了一種新的心率檢測框架來改善rPPG的運動魯棒性。其基本思想來源于觀察到攝像機可以同時對多個皮膚區(qū)域進行并行采樣,并且每個皮膚區(qū)域都可以作為獨立的脈沖測量傳感器進行處理。因此可以利用圖像傳感器的空間冗余來區(qū)分脈沖信號和運動引起的噪聲。為此,該論文構(gòu)造基于像素的rPPG傳感器以使用運動補償?shù)南袼氐较袼孛}沖提取、空間修剪和時間濾波來估計魯棒的脈沖信號,該系統(tǒng)大大改善了對于運動引起的膚色變化而導致的心率不穩(wěn)的情況。
Wang等[6]提出了一個數(shù)學模型,該模型結(jié)合了皮膚反射的相關(guān)光學和生理特性,目的是增加對遠程光電容積描記法(rPPG)背后的算法原理的理解。該模型用于解釋現(xiàn)有rPPG脈沖提取方法的不同選擇。來自該模型的理解可用于設(shè)計穩(wěn)健或特定應(yīng)用的rPPG解決方案。他們通過設(shè)計一種替代的rPPG方法來說明這一點,其中與膚色正交的投影平面用于脈沖提取。各種rPPG方法的一個大基準表明,它們的相對優(yōu)點確實可以從提出的模型中理解。
Macwan等[7]提出了一種新的半盲源提取方法,用于使用自相關(guān)作為周期性度量的多目標優(yōu)化方法來測量Rppg。
Yu等[8]提出了一種新穎的不顯眼的多模態(tài)傳感器,用于監(jiān)測生理參數(shù),包括電容式心電圖(cECG)、反射光電容積描記圖(rPPG)和單個傳感器中的磁感應(yīng)監(jiān)測(MI)。傳感器系統(tǒng)包括設(shè)計和優(yōu)化的傳感器節(jié)點,用于集成到床中的多個傳感器的網(wǎng)格狀陣列和用于數(shù)據(jù)收集和處理的中央控制器盒中。因此,它在應(yīng)用中具有高度通用性,適用于在專業(yè)環(huán)境和家庭護理環(huán)境中對生命體征進行不顯眼的監(jiān)測。所提出的硬件設(shè)計考慮了cECG和MI之間的模間干擾以及由于緊鄰的兩個MI傳感器之間的串擾引起的模內(nèi)干擾。對于心跳間隔,該系統(tǒng)使用cECG實現(xiàn)了2.3 ms的RMSE和0.99的相關(guān)性。類似的,使用rPPG,該系統(tǒng)實現(xiàn)了18.9 ms的RMSE,其相關(guān)性為0.99。關(guān)于源自MI的呼吸間隔,該系統(tǒng)實現(xiàn)了1.12 s的RMSE和0.90的相關(guān)性。
張家成等[12]提出了一種改進色度模型的非接觸式心率檢測與估計方法,他們考慮到在基于色度模型的視頻圖像心率估計算法未考慮到實驗中光照變化產(chǎn)生的信號干擾。針對該問題,他們提出了一種基于改進色度模型的非接觸心率估計算法,通過定義感興趣(ROI)區(qū)域和背景區(qū)域的反射光模型,利用兩種區(qū)域顏色通道信息的比值來消除光照干擾,從而得到消除光照干擾的色度模型,并進行脈搏波的提取與心率估計。對10名在校大學生采集到的人臉圖像進行分析,發(fā)現(xiàn)該算法估計心率的相對誤差為3.52%、2.53%,原始色度模型算法的相對誤差為7.23%、2.82%,兩者對比有顯著性差異(P<0.01)。該算法能有效抑制光照變化所引起的干擾,使所提取脈搏波信號的信噪比平均提高1.52 dB,這對更精確地進行非接觸式心率估計與脈搏波信號提取具有重要意義。
當前,非接觸式心率檢測大致有三類:小波多普勒雷達方法、熱成像方法和基于計算機視覺的圖像方法。但是第一種辦法還有問題未解決,比如魯棒性有待提高、信噪比不滿足實驗所需。熱紅外成像方法的心率檢測成本較為昂貴。
所以本文采用基于視頻的圖像方法,首先,圖像的特征信息比較豐富,這種豐富的信息有利于全面準確地獲取心率信息,但是其噪聲相對來說也比較多,如何對圖像得到的心率進行濾波也是一個需要解決的問題,另外,基于視頻的心率檢測也比較容易受到光照的影響,但是其成本相對來說比較廉價,可以大規(guī)模地在人群中使用,考慮到目前通過普通攝像機測心率準確性低的問題,本文提出用PTZ攝像機的云臺控制組成心率檢測系統(tǒng),以測得心率數(shù)據(jù),即將基于攝像頭的非接觸式的心率檢測應(yīng)用于此??蛇@樣理解,這種非接觸式的心率檢測方式能做到讓人在不察覺的狀態(tài)下被檢測,由此可以規(guī)避由接觸式檢測帶來的心理波動以及不舒服的感覺,此外,這種方式依賴于攝像頭,只要有攝像頭能安裝的地點即可,攝像頭可以安裝的地點條件基本沒有要求,所以可認為這種方法基本不受場地限制。綜上所述,用攝像頭檢測心率的適用范圍極其廣。
首要在于人臉信息檢測的準確性,本文采用的PTZ相機用以實現(xiàn)智能機器人檢測心率的目的。該相機系統(tǒng)可自動調(diào)焦,結(jié)合云臺控制系統(tǒng)得到分辨率適宜的視頻圖像(捕捉目標至視頻中心處)。用海康威視PTZ相機進行人臉的捕捉放大和追蹤,從而得以進行對于長距離非約束的人臉捕捉。之后由系統(tǒng)自動進行調(diào)節(jié)焦距,捕捉到目標再將其放大,之后檢測心率。
為了進行實驗系統(tǒng)的準確性校驗,分別在被試者靜止狀態(tài)與上半身運動狀態(tài)進行多組檢測,且在用相機測心率的同時,用血壓儀檢測得到真實實時心率。通過計算求出兩個狀態(tài)下實驗系統(tǒng)與血壓儀所測心率結(jié)果的平均絕對誤差A、誤差的標準差B、均方根誤差C,且用medcalc軟件繪制實驗系統(tǒng)與血壓儀這兩種方法的B-A圖去分析兩種辦法的一致性,進而檢驗基于PTZ相機的非接觸式心率檢測這一方法的準確程度。
為了闡述基于PTZ相機非接觸式心率檢測實驗系統(tǒng)的檢測原理,本文從基礎(chǔ)開始定義了一個rPPG模型,表示皮膚反射的光學和生理特性,如圖1所示,采取一個光源去照明一部分皮膚組織(該組織內(nèi)要包含人體動脈血液),同時由一個遠程攝像機記錄這幅圖像。
圖1 皮膚反射模型示意圖
在實驗中,隨著人的心臟跳動(心室的收縮和舒張),人體血管內(nèi)的血液容積會產(chǎn)生周期性的變化,則由光源照射到皮膚,經(jīng)皮膚反射到PTZ相機傳感器的光照強度也產(chǎn)生周期性的變化。采用面部作為研究點,是因為面部皮膚的毛細血管較為豐富,會很好地體現(xiàn)出血液容積的周期性變化,利用PTZ攝像機的云臺控制、PTZ自動巡航檢測功能,使其捕捉人臉目標并跟隨,將人臉區(qū)域放大,也就是使面部皮膚的毛細血管內(nèi)的血液容積的周期性的變化更明顯,經(jīng)過算法去分析提取到面部PPG信號的周期性變化,最終得以測得心率數(shù)據(jù)。
1.2.1實驗儀器
因為人的狀態(tài)不一定恒為靜止,為了檢驗實驗的科學性,既要做靜止組實驗也要做運動組實驗。針對于運動組實驗,考慮到實際生活中人的運動往往具有多角度隨機性,則選用PTZ相機,控制其云臺全方位(左右/上下)移動及鏡頭變倍、變焦,實現(xiàn)多角度清晰跟蹤式測量,PTZ相機如圖2所示。在電腦上編寫程序進行調(diào)用。為了后續(xù)檢驗實驗系統(tǒng)的準確度,選用血壓儀提取同步的心率數(shù)據(jù),血壓儀如圖3所示。
圖2 實驗所選PTZ相機 圖3 實驗所選血壓儀
1.2.2實驗操作
為了檢驗基于PTZ相機的非接觸式心率檢測系統(tǒng),實驗環(huán)節(jié)安排被試者在不同的相機和人臉距離(0.5 m、1.0 m、1.5 m、2.0 m、2.5 m、3.0 m、3.5 m、4.0 m)下測量心率,以分析使用距離對于本實驗系統(tǒng)的影響。之后由控制變量法設(shè)定相機的焦距為2倍,光照強度和距離不變時,讓被試者分別以靜止和上半身運動的狀態(tài)進行30 s時長的測量。例如相機和人臉距離為2.0 m時,現(xiàn)場布置如圖4所示。
圖4 2.0 m距離時實驗布置示意圖
為了讓相機測出無察覺狀態(tài)下被試者的心率,靜止狀態(tài)只規(guī)定被試者頭部保持靜止,且不能以后腦勺朝向鏡頭,其余的面部及眼睛部位信息不作要求。
上半身運動狀態(tài)規(guī)定:由上半身運動帶動頭部位置變化,導致相機鏡頭和眼睛連線由居中到傾斜至左側(cè)60°,再到傾斜至右側(cè)60°,最后回到居中。規(guī)定整個運動周期為10 s,每個距離的設(shè)定下連續(xù)運動三個周期。
特別說明,為了確保實驗的科學性,所有實驗均讓被試者在被相機檢測的同時,用血壓儀檢測出真實同步心率數(shù)據(jù)。例如讓被試者提前佩戴血壓儀,在靜止組的相機與人臉距離為1.5 m實驗過程中,運行程序。交互界面的結(jié)果如圖5所示,血壓儀檢測結(jié)果如圖6所示。
由圖5可知,此時實驗系統(tǒng)所測心率為64.1次/min。由圖6可知,此時血壓儀所測心率為66次/min。最后分別將靜止組和運動組所得的各距離下兩種心率測量方法的實驗數(shù)據(jù)記錄整合,再進行后續(xù)的數(shù)據(jù)處理與分析。
1.3.1人臉檢測
利用PTZ相機,為了實現(xiàn)能夠捕捉定位并放大動態(tài)的人臉,就要在視頻中精準檢測到人臉形成定位捕捉框,并檢測出盡可能多的關(guān)鍵點以便后續(xù)進行人臉關(guān)鍵點的預處理,去篩選出真正的面部皮膚像素點以便于后面的膚色檢測。
參考于仕琪老師的人臉檢測開源庫,對于視頻中的被試者進行人臉檢測,得到68個人臉特征點。并將被試者的面部進行捕捉定位,得到人臉定位矩形框。
1.3.2皮膚檢測
(1) 簡介。經(jīng)過人臉檢測,得到68個人臉特征點以及人臉定位框。但是由于人臉定位框所選中的區(qū)域并不完全是想要的皮膚像素點,即包括眼鏡和眼睛以及眉毛等無關(guān)信息。因此需要做一定的圖像處理,初步篩選出明顯不屬于人臉皮膚的區(qū)域,然后再進行膚色檢測。這樣可以使得到的效果更加準確,也大大提高了算法的效率。
(2) 算法的比較與選用。對于所要用的皮膚檢測的方法,經(jīng)前期調(diào)研,可認為基于橢圓皮膚模型的皮膚檢測在非機器學習方法中是比較準確的,另外,本文還參考了機器學習方法OCSVM(即單類支持向量機)。使用這一模型將獲得的樣本集合投影到高維空間里,基于此,再通過計算使得目標與坐標圖的原點能大限度剝離開。
在不同的光照環(huán)境下測試OCSVM與橢圓皮膚模型兩種皮膚檢測算法,分析得出橢圓皮膚模型(即閾值法)的分割更為連續(xù),且分割效果相較于OCSVM而言更好,則本文選用橢圓膚色模型進行膚色檢測的研究。
(3) 橢圓膚色模型。首先,使用68個人臉特征點求凸包,即將最外側(cè)人臉特征點連接成一個封閉區(qū)域。其次,以鼻子處特征點為中心做滿水填充,生成所分割區(qū)域的圖像掩碼。最后得到預處理之后的區(qū)域部分,完成對圖像中人面部的預處理。
關(guān)于皮膚檢測,本文參考了傳統(tǒng)的橢圓膚色模型,此算法針對像素點的分布情況,判斷待檢測像素點的位置是否不在限定區(qū)域(橢圓)之外去判斷像素是否是皮膚像素??梢哉J為,本質(zhì)上來說該算法屬于一種閾值法。
通過此模型進行人臉皮膚關(guān)鍵點的篩選,去除掉眼鏡和眼睛以及眉毛等無關(guān)信息,只保留毛細血管豐富的面部皮膚像素點區(qū)域,結(jié)果如圖7所示。
圖7 皮膚檢測結(jié)果
1.3.3心率信號的提取
(1) 光流法提取rPPG信號。針對于前面皮膚檢測步驟得到的被試者的面部皮膚像素點區(qū)域,本文用光流法去看其圖里像素點活動的瞬時速度,檢查時域內(nèi)圖像像素的轉(zhuǎn)變以及相鄰各圖像間的一些聯(lián)系。也可用這種辦法去搜羅上一幀圖像和目前幀的圖像之間的那種關(guān)系,進而再由計算推廣得出各相鄰圖像之間的目標各移動狀態(tài)伴隨的數(shù)據(jù)信息。通過光流法,能夠捕捉上一幀圖像中每個皮膚的像素在下一幀圖像中的估計位置,從而獲得兩幀之間圖像的相對位移量。
通過調(diào)研,選出兩種光流算法,即以待選用。即Farneback光流算法和DIS光流算法。
Farneback光流算法會借由公式:f(x)=xTAx+bTx+c去對圖像進行逼近,其中:x為像素矩陣;A、b、c為參數(shù)。對于每幅圖像中的每個像素,要先定義一個有效的鄰域,緊接著使用已定義的該鄰域里(2n+1)×(2n+1)個像素點當成最小二乘法的樣本點集合,就可以由此來調(diào)整所獲得的圖像的中心像素點的系數(shù)數(shù)據(jù)。對于某個像素點,假設(shè)原始灰度值是f(x),并且全局位移之后的值增加至f(x+d)。假設(shè)f(x)=f(x+d)。然后利用多項式的對應(yīng)系數(shù),能夠求解出這個像素的全局位移d,最后獲得圖像的密集光流場。
DIS光流算法主要包含有3個方面:逆搜索,得到稠密的光流場,精細化搜索。
DIS算法最大的優(yōu)點在于運算速度的提升,在單核CPU上,1 024×436分辨率可以達到300 Hz~600 Hz運行速率;包含圖像預處理。OpenCV4中提供了DIS光流法的接口函數(shù),其使用方式與Farneback算法在OpenCV中的接口函數(shù)類似[17]。
(2) 算法的比較與選用。針對于所選的Farneback光流算法和DIS光流算法,進行二者的比較和最終選用。
通過由DIS光流法進行計算后,會得到一個雙通道圖像(即光流flow),其中各通道分別包括x及y兩個方向的位移,通過munsell顏色系統(tǒng)將其實現(xiàn)可視化。
DIS光流算法和Farneback光流算法的演示結(jié)果分別如圖8和圖9所示??梢钥吹?由DIS光流法得到的人臉區(qū)域光流信息較Farneback光流法更全面,且該算法運行的速度相較于Farneback光流法更快,則本文選用DIS光流算法。
圖8 DIS光流法 圖9 Farneback光流法
1.3.4PPG信號提取的精確化
通過使用DIS光流法得到相鄰幀之間的像素對,然后對它們進行級聯(lián),在一個長的時間序列上將會得到若干個原始心率信號。但是,這種方法還是有弊端,頭部的某些運動(比如轉(zhuǎn)頭)有可能導致某些像素在下一幀中丟失。
為了實現(xiàn)PPG信號的精確化提取,完善DIS光流算法的問題,參考論文《Exploiting Spatial-redundancy of Image Sensor for Motion Robust rPPG》。該文提到,在恒定的光照環(huán)境下,皮膚不同位置的像素呈現(xiàn)相同的相對PPG幅值。這意味著如果每個像素對的脈沖誘導的顏色變化被時域歸一化,它們可以按照任意的順序連在一起,從而得到一個長期的信號。相鄰兩幀之間的時域歸一化的像素量可理解為:
(1)
借用Chrom法將時域歸一化的RGB差值投影至色度面上(把脈沖信號定義成RGB通道的線性組合),如下:
(2)
(3)
(4)
式中:l是由滑動窗口定義的色度跟蹤的間隔長度。本系統(tǒng)中,把l定為128幀。之后求滑動窗中的ppG信號為:
(5)
式中:α代表標準差。
(6)
注意,在上述過程中,局部像素的空間平均可以減少時間顏色正?;陂g的量化誤差。因此從人臉檢測步驟開始對人臉ROI進行下采樣(在OpenCV中,調(diào)用resize接口函數(shù)即可),這不僅降低了基于像素的rPPG傳感器的噪聲靈敏度,而且提高了DIS光流法的處理速度。選擇下采樣時,要斟酌關(guān)于尺寸的研究。考慮到精度和效率的大小,將人臉ROI區(qū)域均勻下采樣至寬度為40像素的比例,下采樣的比例由第一幀圖像的人臉尺寸確定。
1.3.5利用PCA降維
在數(shù)據(jù)處理時,特征維度很可能比樣本數(shù)量多很多,那樣的話效果不一定利于實操。有兩方面的原因:(1) 噪聲影響機器運算的效果;(2) 因為無法避免地造成計算量增加,會導致時間復雜度和空間復雜度不必要地增加。為了解決這些問題,就要進行數(shù)據(jù)轉(zhuǎn)換:關(guān)于特征維度高的問題就要進行降低維度的處理;為了減少噪聲并且降低計算量,就需要消除一些特征之間的相關(guān)性。則可以選用PCA技術(shù)。將原始數(shù)據(jù)的特征空間投影至較小的子空間,以盡可能喪失少的信息去表達。為了更好地去表達數(shù)據(jù),就進行降維并從子空間內(nèi)得到數(shù)據(jù)信息。
之前提取出來的PPG心率信號是k個面部皮膚像素點組成的多維k×54矩陣,不利于實操分析。則本實驗系統(tǒng)利用PCA技術(shù)將這些多維心率信號降成一維的心率信號。即將k×54矩陣進行數(shù)據(jù)的降維,降維到矩陣X(一個1×54矩陣)。
利用PCA技術(shù)得到的在一段時間內(nèi)的心率信號如圖10所示。
圖10 一段時間內(nèi)得到的心率信號
1.3.6用傅里葉變換計算心率信號
由上一步的PCA技術(shù)得到一段時間內(nèi)的心率信號是處于時域狀態(tài)的結(jié)果,因為要求得被試者心率數(shù)據(jù),要得到頻率值,則要將得到的時域心率信號轉(zhuǎn)換到頻域進行分析,尋找最高峰值處的數(shù)值即為被試者的心率數(shù)值。
由傅里葉變換去達到分析信號分量的目的,也能夠操縱這些分量合成某信號,并使用正弦波作為信號分量。關(guān)于f(t)的傅里葉變換如式(7)所示。
(7)
式中:f(t)是F(ω)的原像。
在OpenCV中dft函數(shù)可以實現(xiàn)關(guān)于離散信號的傅里葉變換過程,將心率信號轉(zhuǎn)換到頻率域去算心率信息。首先找到頻域中的峰值pmax,然后利用公式bpm=pmax×fps/total×SEC_PER_MIN,其中:fps為視頻幀率;total為處理一組視頻的幀數(shù)(在本系統(tǒng)中固定為54);SEC_PER_MIN為常數(shù)60。
經(jīng)過傅里葉變換后得到頻域的心率信號如圖11所示。
圖11 心率信號的頻域結(jié)果
由圖11可見,頻域圖中的最高點即為心率數(shù)據(jù)采集點,在該峰值處的心率數(shù)值傳至含有人臉的視頻頁面,顯示此時被試者的心率數(shù)值。
2.1.1實驗測量
根據(jù)前面介紹的方案,在不同的相機與人臉距離下,做靜止組實驗并將心率數(shù)據(jù)統(tǒng)計整合,結(jié)果如表1所示。
表1 靜止狀態(tài)下的心率測量結(jié)果
2.1.2實驗系統(tǒng)的精確度分析
由表11中數(shù)據(jù)計算出各距離下,兩種心率測量方式測量結(jié)果的平均絕對誤差為1.30次/min,誤差的標準差是2.073 5次/min以及均方根誤差為1.554 8次/min。
結(jié)果表明,在不同的機臉距離下,兩種測量方法的平均絕對誤差小于2.00次/min,不超過中華人民共和國醫(yī)藥業(yè)的誤差標準5次/min,誤差的標準差小于2.1次/min,均方根誤差小于1.7次/min,這表明被試者靜止狀態(tài)下,實驗系統(tǒng)和血壓儀測量的結(jié)果非常接近,具有一定的準確性。
2.1.3一致性分析
在醫(yī)學領(lǐng)域里,要畫出B-A散點圖去評估實驗系統(tǒng)和血壓儀兩種心率測量方法的一致性,即檢驗實驗系統(tǒng)的科學性。
借助MedCalc軟件繪制兩種心率測量方法的B-A圖,結(jié)果如圖12所示。
圖12 靜止狀態(tài)的Bland-Altman數(shù)據(jù)分析圖
關(guān)于圖12,Mean線表示誤差平均值線,SD代表標準差,±1.96 SD表示95%的置信區(qū)間。一般而言,從指標顯示,二者的差值平均數(shù)越接近0,95%置信區(qū)間的范圍程度越小時,可認為本實驗系統(tǒng)和血壓儀這兩種心率測量方法的一致性越高,即本實驗系統(tǒng)的準確度越高,科學性越強。
從圖12可知,在被試者靜止狀態(tài)的情況下,實驗系統(tǒng)和血壓儀這兩種心率測量方法的心率結(jié)果的差值平均數(shù)為-1.3次/min,取絕對值后為1.3次/min,95%置信區(qū)間為-3.1~+0.5 次/min,則可認為被試者靜止狀態(tài)下,本實驗系統(tǒng)和血壓儀這兩種心率測量方法的一致性較高、即本實驗系統(tǒng)的準確度較高,科學性較強。
2.2.1實驗測量
根據(jù)前面介紹的方案,在不同的機臉距離下,做運動組實驗并將心率數(shù)據(jù)統(tǒng)計整合,結(jié)果如表2所示。
表2 上半身運動狀態(tài)下的心率測量結(jié)果
2.2.2實驗系統(tǒng)的精確度分析
由表13中數(shù)據(jù)計算出各距離下,兩種心率測量方式測量結(jié)果的平均絕對誤差為0.987 5次/min,誤差的標準差是2.379 3次/min以及均方根誤差為1.184 8次/min。
結(jié)果表明,在不同的相機與被試者面部距離下,兩種測量方法的平均絕對誤差小于1.00次/min,不超過中華人民共和國醫(yī)藥行業(yè)的規(guī)定誤差標準(即誤差控制在不超過5次/min),誤差的標準差小于2.4次/min,均方根誤差小于1.2次/min,這表明被試者在運動狀態(tài)下,實驗系統(tǒng)具有較高的準確性。
2.2.3一致性分析
借助MedCalc軟件繪制運動狀態(tài)下,兩種心率測量方法的B-A圖去分析二者的一致性,結(jié)果如圖13所示。
圖13 上半身運動狀態(tài)的Bland-Altman數(shù)據(jù)分析圖
從圖13可知,在被試者上半身狀態(tài)的情況下,實驗系統(tǒng)和血壓儀這兩種心率測量方法的心率結(jié)果的差值平均數(shù)為-1.0次/min,取絕對值后約為1.0次/min,95%置信區(qū)間為-2.4~+0.4次/min,則可認為被試者上半身運動狀態(tài)下,本實驗系統(tǒng)和血壓儀這兩種心率測量方法的一致性較高,實驗系統(tǒng)科學性較強。
通過前面統(tǒng)計的心率數(shù)據(jù),比較靜止狀態(tài)與上半身運動狀態(tài)這兩種被試者實驗狀態(tài)組的平均絕對誤差、誤差的標準差、均方根誤差,發(fā)現(xiàn)均符合測量標準,且變動不大。
通過Bland-Altman圖,比較靜止狀態(tài)與上半身運動狀態(tài)這兩種被試者實驗狀態(tài)組的差值平均值的絕對值以及95%置信區(qū)間,發(fā)現(xiàn)變動不大,且符合檢測標準。
則本實驗系統(tǒng)對于靜止以及上半身運動的人來說,均可進行精確度較高的心率測量,對非接觸式測心率的研究有參考意義。
醫(yī)學行業(yè)測量心率,大部分都要和人的皮膚接觸,可能會因為被試者的心理波動造成數(shù)據(jù)不準確,還可能會引發(fā)被測人員的皮膚問題。在一些特殊時期,比如傳染性較高的疫情期間,人為接觸式測量可能會涉及對人力資源的要求。以及為確保不傳染,對接觸式測量儀器的無菌性要求,對于測量與被測量人員的安全性要求。則非接觸式的心率檢測具有很大的應(yīng)用前景。
本文中采用PTZ相機作為非接觸式檢測心率的媒介。該相機系統(tǒng)可自動調(diào)焦,結(jié)合云臺控制系統(tǒng)得到分辨率適宜的視頻圖像(捕捉目標至視頻中心處)。在操作相機系統(tǒng)的自動巡航功能捕捉到人臉目標后,進行視頻處理。借助經(jīng)比較后所選的橢圓膚色模型將得到的人臉框和特征點區(qū)域進行皮膚檢測。針對于提取到的皮膚像素點,結(jié)合經(jīng)過預先挑選比較后的采用的DIS光流法精確提取其PPG信號。再由PCA技術(shù)降維整合得到時域的心率信號,最后借助傅里葉變化得到頻域的心率信號,讀取得到心率數(shù)據(jù)。分別在靜止與上半身運動這兩種被試者狀態(tài)下,各進行八組不同距離的對比實驗,求出各組的實驗值與血壓儀測量值的平均絕對誤差、誤差的標準差、均方根誤差,并繪制了B-A一致性分析圖去判斷95%置信區(qū)間的范圍,由數(shù)據(jù)分析證明出本實驗系統(tǒng)與血壓儀測量這兩種測量方法的一致性較強,即證明實驗系統(tǒng)具有一定的準確性。
目前關(guān)于非接觸式心率檢測,一些學者采用自適應(yīng)去噪進行心率檢測。本文的準確性較之更為優(yōu)良,且不亞于同類的多數(shù)其他文章里的實驗系統(tǒng),存在科學應(yīng)用性。
同類大多數(shù)文章所采用的攝像頭均為固定式攝像頭或者某個固定攝像頭平臺,只能捕捉一定角度范圍內(nèi)的人臉目標。本文選用PTZ相機作為視頻處理的媒介就具有很大的應(yīng)用優(yōu)勢。因為應(yīng)用于現(xiàn)實中,所測之人當時不一定保持靜止,則基于PTZ相機,編寫程序調(diào)用其360°全方位自動巡航的功能就可以全方位捕捉在附近運動任意角度的人的面部,具有更大的應(yīng)用型優(yōu)勢。
同類一些文章的實驗中,被試者直視攝像頭進行信息采集,一定程度上會對被試者的心理造成波動,得到的心率數(shù)據(jù)不一定準確。本文在用相機采集人臉時,被試者多數(shù)以側(cè)面朝向鏡頭,且注意力不在相機鏡頭上,一定程度實現(xiàn)了無察覺的心率測量,且得到的數(shù)據(jù)具有一定的準確度。應(yīng)用到現(xiàn)實中就可以將PTZ相機進行隱蔽處理,從而使人在未知的情況下被采集信息,解決了強制直視鏡頭進行信息采集導致數(shù)據(jù)說服力可能不足的問題,從而提升了最終數(shù)據(jù)的生態(tài)效度。
對于應(yīng)用距離的遠近程度,本文的實驗系統(tǒng)的檢測距離由0.5 m到4.0 m。通過數(shù)據(jù)可知,具備一定的準確性。通過猜想和分析,當應(yīng)用距離即相機與人臉間距大于4.0 m時依然存在可觀的準確性,但是還未經(jīng)證實,目前仍在研究,已獲得的成果已達預期的目的,實現(xiàn)了顯著的效果,且具有一定的科學性。基于PTZ相機的非接觸式心率檢測有了很大優(yōu)勢。