• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    自適應(yīng)人臉多區(qū)域分析的視頻心率檢測(cè)①

    2022-11-07 09:07:38楊學(xué)志李龍偉王定良
    關(guān)鍵詞:區(qū)域信號(hào)檢測(cè)

    陳 宇,楊學(xué)志,方 帥,李龍偉,王定良

    1(合肥工業(yè)大學(xué) 計(jì)算機(jī)與信息學(xué)院,合肥 230009)

    2(工業(yè)安全與應(yīng)急技術(shù)安徽省重點(diǎn)實(shí)驗(yàn)室,合肥 230009)

    3(合肥工業(yè)大學(xué) 軟件學(xué)院,合肥 230009)

    4(中國(guó)科學(xué)技術(shù)大學(xué)附屬第一醫(yī)院,合肥 230036)

    心率作為人的基本生命體征,不僅可以反應(yīng)人體的健康狀態(tài),也是評(píng)估人體心腦血管疾病最重要和有效的參數(shù)之一.其在心腦血管疾病、高血壓、心理和精神疾病等領(lǐng)域都被證實(shí)是一種有效的診斷手段[1].

    臨床上傳統(tǒng)的心率監(jiān)測(cè)技術(shù)主要是利用心電圖機(jī)(electrocardiogram,ECG),其通過設(shè)備上的十二導(dǎo)聯(lián)綁定在體表不同的裸露部位獲取精確的心電數(shù)據(jù)進(jìn)行分析處理[2].這種方式可以得到較為精準(zhǔn)的心率,但是其存在成本高昂、操作繁雜、專業(yè)要求高,普通家庭無(wú)法獨(dú)立使用.除此以外,家用的脈搏血氧儀也被廣泛用來(lái)測(cè)量受試者心率,其是利用光電傳感器夾在受試者指尖、耳垂等體表末梢通過接收反射或透射的光信號(hào)進(jìn)行心率測(cè)量.該方法需要測(cè)試部位和設(shè)備緊密貼合且不能夠有相對(duì)運(yùn)動(dòng),使用不便,舒適度較差不利于長(zhǎng)時(shí)間實(shí)時(shí)檢測(cè).針對(duì)傳統(tǒng)方法的在實(shí)際應(yīng)用中的局限性,研究一種使用方便的非接觸式的心率檢測(cè)技術(shù)就顯得十分重要.

    基于視頻的光電容積描記術(shù)(video photoplethysmography,VPPG),是一種非接觸式心率檢測(cè)技術(shù).其可以利用普通彩色相機(jī)采集面部或指尖等區(qū)域表皮視頻,通過檢測(cè)皮下毛細(xì)血管由于心臟周期泵血使得血液血紅蛋白濃度變化產(chǎn)生的微弱顏色變化,提取血液容積脈沖信號(hào)(blood volume pulse,BVP).通過分析該信號(hào)的頻率即可進(jìn)行心率估計(jì).Verkruysse 等人[3]最早使用普通彩色攝像頭,在可見光作為光源的條件下證明了使用VPPG 技術(shù)進(jìn)行非接觸式心率檢測(cè)的可行性.隨后,Poh 等人[4]提出對(duì)視頻幀進(jìn)行人臉感興趣區(qū)域檢測(cè)并計(jì)算RGB 三通道像素均值,然后使用獨(dú)立成分分析(ICA)進(jìn)行BVP 信號(hào)提取和心率估計(jì).盡管這些算法都可以比較好的進(jìn)行心率估計(jì),但需要測(cè)試者全程保持靜止或人臉和環(huán)境背景的穩(wěn)定才能保證良好的實(shí)驗(yàn)結(jié)果,這不僅會(huì)給測(cè)試者帶來(lái)不便,而且會(huì)限制算法的使用場(chǎng)景.由于人臉運(yùn)動(dòng)(如說(shuō)話、表情變化)、頭部晃動(dòng)和光照不均勻等干擾因素的存在,會(huì)影響VPPG技術(shù)BVP 信號(hào)的提取以及心率估計(jì)的準(zhǔn)確性.因此,如何保證在真實(shí)的生活場(chǎng)景中提升心率檢測(cè)算法的穩(wěn)定性和準(zhǔn)確性,已經(jīng)成為國(guó)內(nèi)外學(xué)者研究的重點(diǎn).

    研究者們主要做出的改進(jìn),有如下幾個(gè)方面: 其一,是根據(jù)皮膚光照反射特點(diǎn)采用不同的顏色空間進(jìn)行BVP 信號(hào)提取,通過不同顏色通道信號(hào)的差異性和互補(bǔ)性來(lái)增強(qiáng)原始BVP 信號(hào).如de Haan 等人[5]利用不同顏色通道信號(hào)之間的互補(bǔ)性,提出了基于色度的顏色空間投影分解算法(chrominance-based rPPG,CHROM)以減輕頭部晃動(dòng)產(chǎn)生的影響.其二,是根據(jù)BVP 信號(hào)特性以及先驗(yàn)知識(shí)等進(jìn)行噪聲濾除和信號(hào)增強(qiáng)等操作.如Wang 等人[6]通過分析BVP 信號(hào)頻域上心動(dòng)信號(hào)和噪聲信號(hào)的區(qū)別,采取頻域修正的方式來(lái)去除非心率范圍內(nèi)的干擾.Tasli 等人[7]通過估計(jì)頭部運(yùn)動(dòng)頻率并進(jìn)行剔除,取得了更魯棒的心率估計(jì).

    除此以外,近年來(lái)研究表明通過同步分析多個(gè)不同的面部感興趣區(qū)域的信號(hào)可以用來(lái)提升原始BVP信號(hào)的提取效果.如Kumar 等人[8]使用信噪比來(lái)對(duì)不同的感興趣區(qū)域信號(hào)進(jìn)行加權(quán).Feng 等人[9]將ROI 區(qū)域分塊后設(shè)計(jì)并提取了兩種特征,最后使用K均值算法自適應(yīng)選擇有效的感興趣區(qū)域.Li 等人[10]使用人臉特征點(diǎn)裁剪獲得固定的人臉皮膚區(qū)域.Favilla 等人[11]使用ICA 算法分析面部3 個(gè)區(qū)域的多通道信息進(jìn)行心率檢測(cè).周雙等人[12]采用矩陣低秩分解來(lái)對(duì)不同感興趣區(qū)域的信號(hào)進(jìn)行降噪處理.然而,由于以上的多區(qū)域分析的劃分均是在人臉檢測(cè)框的基礎(chǔ)上進(jìn)行簡(jiǎn)單的矩形分割,此方式雖然操作簡(jiǎn)單,但實(shí)際上并沒有考慮到分塊內(nèi)人臉輪廓、陰影和皮膚紋理等信息,以及視頻幀之間分塊區(qū)域位于人臉位置的一致性,ROI 分塊劃分存在一定的不合理之處.

    為了解決上述ROI 分塊劃分和挑選中存在的問題,本文受傳統(tǒng)的超像素分割算法的啟發(fā),提出一種基于人臉視頻的自適應(yīng)超像素分割和ROI 優(yōu)選方法.使得分塊區(qū)域擁有相同或近似的特征集,一方面提升子塊劃分質(zhì)量,從源頭提升了BVP 信號(hào)提取的質(zhì)量,另一方面通過優(yōu)選算法尋找面部信號(hào)最佳區(qū)域.增強(qiáng)了面部運(yùn)動(dòng)干擾以及環(huán)境光分布不均時(shí)算法的抗干擾能力.

    1 心率檢測(cè)方法

    本文心率檢測(cè)算法是以ROI 超像素分割和脈搏矩陣信號(hào)優(yōu)選為基礎(chǔ).通過裁切視頻幀內(nèi)面部ROI,隨后根據(jù)輪廓、皮膚紋理和陰影等特征將面部ROI 通過超像素聚類分割成互不重疊的子塊,再對(duì)分塊結(jié)果進(jìn)行色度BVP 矩陣提取,最后進(jìn)行多指標(biāo)綜合分析優(yōu)選出最佳面部子區(qū)域進(jìn)行信號(hào)處理與心率估計(jì).

    本方法的核心步驟包括以下幾個(gè)部分: 1)人臉與特征點(diǎn)檢測(cè); 2)人臉自適應(yīng)超像素分割; 3)色度BVP矩陣提取; 4)脈搏信號(hào)評(píng)估優(yōu)選與心率估計(jì).具體算法框圖如圖1 所示.

    圖1 算法流程及示意圖

    1.1 人臉與特征點(diǎn)檢測(cè)

    在真實(shí)生活場(chǎng)景中,由于檢測(cè)過程中頭部晃動(dòng)以及表情變化等會(huì)影響攝像頭對(duì)面部顏色周期性變化的提取效果,進(jìn)而導(dǎo)致BVP 信號(hào)提取質(zhì)量的大幅降低.因此對(duì)原始視頻進(jìn)行人臉與特征點(diǎn)檢測(cè)算法,來(lái)獲取穩(wěn)定的人臉ROI 區(qū)域(圖像幀).首先根據(jù)視頻第1 幀圖像,使用Viola 等人[13]的人臉檢測(cè)器輸出視頻中人臉檢測(cè)矩形區(qū)域的頂點(diǎn)坐標(biāo); 然后采用魯棒判別響應(yīng)圖擬合 (DRMF)[14]方法對(duì)人臉坐標(biāo)區(qū)域內(nèi)主要器官進(jìn)行特征點(diǎn)定位,記錄人臉66 個(gè)特征點(diǎn)位置C1,包含眼睛、鼻子、眉毛、臉頰和嘴巴等區(qū)域.之后對(duì)后續(xù)的視頻幀進(jìn)行KLT 跟蹤特征點(diǎn)的位置,獲取連續(xù)幀面部特征點(diǎn)矩陣C,記作:

    其中,i為視頻幀數(shù),m為特征點(diǎn)個(gè)數(shù),范圍為1-66.

    1.2 人臉自適應(yīng)超像素分割

    傳統(tǒng)的ROI 分割算法均是將人臉檢測(cè)框內(nèi)區(qū)域簡(jiǎn)單的分割成m行n列的矩形.如圖2 所示.這樣分割,首先并沒有考慮到頭部運(yùn)動(dòng)導(dǎo)致分塊對(duì)應(yīng)人臉位置的變化,其次也沒有考慮到分塊內(nèi)像素間的特征聯(lián)系,ROI 分割存在不合理之處.例如: 平面內(nèi)頭部偏轉(zhuǎn)會(huì)造成檢測(cè)框改變,分塊對(duì)應(yīng)面部位置變化,導(dǎo)致連續(xù)幀之間分塊不連續(xù); 光照分布不均時(shí)陰影區(qū)域和正常亮度區(qū)域若被隨機(jī)劃分為一個(gè)子塊,會(huì)導(dǎo)致子塊內(nèi)像素均值增大,顏色變化幅度減小進(jìn)一步降低BVP 信號(hào)提取質(zhì)量.

    圖2 傳統(tǒng)ROI 分割效果演示

    超像素分割算法是通過像素間共享特征(如亮度、紋理、陰影、輪廓等)將圖像分割為包含特定圖像結(jié)構(gòu)的子區(qū)域[15].該算法可用于對(duì)面部區(qū)域進(jìn)行分割,但是傳統(tǒng)的超像素分割僅對(duì)圖像進(jìn)行處理,未考慮連續(xù)幀之間的分塊個(gè)數(shù)及位置的連續(xù)性.故針對(duì)心率檢測(cè)場(chǎng)景對(duì)傳統(tǒng)方法進(jìn)行改進(jìn)提出人臉自適應(yīng)超像素分割算法.

    方法首先針對(duì)視頻圖像中人臉位置運(yùn)動(dòng)做出位置自適應(yīng)處理,再結(jié)合簡(jiǎn)單的線性迭代聚類(simple linear iterative clustering,SLIC)[16]實(shí)現(xiàn)連續(xù)幀間人臉區(qū)域分塊的一致性.能夠很好地應(yīng)對(duì)視頻中人臉位置移動(dòng)、遠(yuǎn)近距離變化和平面內(nèi)旋轉(zhuǎn)等情況下人臉分塊的穩(wěn)定性,保證子塊BVP 信號(hào)提取區(qū)域?qū)?yīng)人臉位置的穩(wěn)定性和連續(xù)性.具體步驟為:

    (1)視頻人臉獲取.根據(jù)第1.1 節(jié)獲取的視頻人臉位置對(duì)原始圖形進(jìn)行裁切獲取原始人臉區(qū)域圖像幀序列F.

    (2)自適應(yīng)處理.傾斜校正: 如圖3 所示根據(jù)獲取的圖像幀序列F以及第1.1 節(jié)獲取的人臉特征點(diǎn)位置C使用戚剛[17]的傾斜校正方法校正視頻中旋轉(zhuǎn)的人臉獲得校正人臉序列F1.

    圖3 人臉傾斜校正

    人臉尺寸一致處理: 由于頭部運(yùn)動(dòng)會(huì)影響圖像中面部像素區(qū)域大小[18],若采用固定的分割個(gè)數(shù)可能會(huì)使得分割子塊像素區(qū)域過小,導(dǎo)致BVP 信號(hào)提取質(zhì)量的急速下降[19].因此通過設(shè)置初始化超像素分塊尺寸,根據(jù)視頻初始幀人臉像素尺寸自適應(yīng)決定分塊個(gè)數(shù).由于人臉五官位置在不同尺寸圖片中比例保持相對(duì)穩(wěn)定,因此根據(jù)面部特征點(diǎn)C對(duì)傾斜校正序列F1進(jìn)行對(duì)比縮放,獲得尺寸一致的人臉圖像序列F2,以確保分塊個(gè)數(shù)一致.

    (3)初始化種子點(diǎn).設(shè)序列F2中圖像像素個(gè)數(shù)為K,為了確保分塊BVP 提取質(zhì)量,設(shè)置初始超像素大小為S=50×50的矩形區(qū)域,則超像素個(gè)數(shù)M=K/S,超像素種子點(diǎn)間的距離L=sqrt(K/S).

    (4)相似度衡量.針對(duì)序列F2圖像幀中的每個(gè)像素,分別計(jì)算與距離最近的種子點(diǎn)間的相似程度,為該像素點(diǎn)賦予最相似種子點(diǎn)標(biāo)簽.相似度衡量公式如下所示:

    其中,dlab為像素點(diǎn)之間的顏色差異,dxy是像素點(diǎn)間的空間距離.Di用來(lái)衡量?jī)上袼亻g的相似程度;L為種子點(diǎn)的間距,e為平衡參數(shù),用來(lái)反應(yīng)顏色差異與空間距離在相似度衡量中的比重.Di越大,表明兩像素之間相似度越高.

    自適應(yīng)超像素分割迭代過程如圖4,通過不斷迭代上述過程直至達(dá)到收斂,分塊之間充分考慮到面部像素之間距離和顏色特征等聯(lián)系,并且經(jīng)過了面部?jī)A斜校正和一致性處理,可以確保連續(xù)幀之間分塊位置的一致性,獲得最終自適應(yīng)分割結(jié)果標(biāo)簽Label,范圍1-M.

    圖4 自適應(yīng)超像素分割迭代過程

    1.3 分塊色度脈搏矩陣提取

    對(duì)獲取的連續(xù)幀人臉自適應(yīng)分塊結(jié)果,根據(jù)Label標(biāo)記區(qū)域進(jìn)行BVP 信號(hào)提取.研究表明[20],HSV 顏色空間H通道色度特征對(duì)面部運(yùn)動(dòng)具有良好的魯棒性,可以抵抗面部運(yùn)動(dòng)光線反射角度變化造成的皮膚亮度和顏色變化.因此采用H通道進(jìn)行色度脈搏提取.首先根據(jù)分塊Label 結(jié)果,分別計(jì)算M個(gè)ROI 分割子區(qū)域像素均值,記為bm(N)其中m表示第m個(gè)分塊區(qū)域,N表示的視頻總幀數(shù).隨后將計(jì)算得到的M個(gè)區(qū)域的像素均值矩陣B通過顏色空間轉(zhuǎn)換到H通道獲得H通道色度BVP 信號(hào)矩陣,記為:

    其中,hm(N)表示第m個(gè)分塊區(qū)域色度脈搏提取信號(hào).

    H通道顏色轉(zhuǎn)換方法在以下公式中進(jìn)行描述.使用以下公式對(duì)RGB 值進(jìn)行標(biāo)準(zhǔn)化與顏色空間轉(zhuǎn)換:

    其中,R′,G′,B′為RGB三通道的歸一化值,范圍0-1;Cmax和Cmin分別為歸一化RGB通道的最大值和最小值,以及歸一化RGB通道的差值Δ;H即為HSV 顏色空間的色度分量.

    1.4 脈搏信號(hào)篩選與心率估計(jì)

    獲取到面部M個(gè)子塊色度BVP 信號(hào)后,對(duì)各子塊的原始脈搏信號(hào)質(zhì)量進(jìn)行評(píng)估與優(yōu)選.在進(jìn)行脈搏波質(zhì)量評(píng)估的特征選取時(shí),必須要符合以下準(zhǔn)則.首先,因?yàn)檎鎸?shí)場(chǎng)景下無(wú)心率參考真值,質(zhì)量評(píng)估指標(biāo)不依賴真值對(duì)比; 其次,特征數(shù)值大小需體現(xiàn)信號(hào)質(zhì)量?jī)?yōu)劣; 最后,特征值應(yīng)方便從簡(jiǎn)短的波形獲得,確保特征及時(shí)更新實(shí)現(xiàn)ROI 優(yōu)選的時(shí)效性.考慮到上述條件,因此選取幅度、頻域和非線性3 個(gè)方面指標(biāo)綜合評(píng)判脈搏波質(zhì)量.

    (1)幅度特征: 分塊區(qū)域亮度變化標(biāo)準(zhǔn)差(light std),表示分塊內(nèi)像素亮度與均值之間的偏離程度,可以體現(xiàn)分塊內(nèi)亮度變化程度.公式如式(8):

    其中,light(m)表示的是第m個(gè)分塊的亮度變化程度指標(biāo),lvaluemi代表第m分塊中第i幀的亮度值,μ代表分塊中所有幀的亮度均值,N為視頻幀總數(shù).

    (2)頻域特征: 分塊脈搏波信噪比(signal noise ratio,SNR),表示脈搏信號(hào)平均功率與噪聲平均功率的比值.可以反應(yīng)心率信號(hào)能量占正常心率區(qū)間內(nèi)的總能量大小,公式如式(9):

    其中,SNR(m)表示的是第m個(gè)分塊的信噪比指標(biāo),fHR表示信號(hào)頻譜峰值頻率,fs表示信號(hào)采樣率,函數(shù)S表示頻譜密度函數(shù).

    (3)非線性特征: 分塊內(nèi)樣本熵.熵值作為一種非線性特征,可以體現(xiàn)隨機(jī)變量不確定性的程度,反應(yīng)信號(hào)的規(guī)則與不規(guī)則性.近來(lái)已經(jīng)被越來(lái)越多人應(yīng)用于對(duì)生理信號(hào)分析上[21],樣本熵(sample entropy,SampEn)對(duì)心率信號(hào)的高斯噪聲尤其敏感,故選取作為波形質(zhì)量評(píng)估特征.

    以上3 種指標(biāo)由于評(píng)價(jià)尺度不一無(wú)法直接對(duì)比數(shù)值大小使用,現(xiàn)分別對(duì)指標(biāo)進(jìn)行有效性歸一化處理,公式如下:

    其中,min 和max 分別為各分塊區(qū)域中對(duì)應(yīng)指標(biāo)的最小值和最大值.Lightnorm(m)、SNRnorm(m)、SampEnnorm(m)分別為第m塊有效性歸一化后的指標(biāo),數(shù)值越大表明分塊內(nèi)BVP 信號(hào)質(zhì)量越好.

    綜合有效性歸一化后的指標(biāo),得到分?jǐn)?shù)Score值,公式如下:

    其中,Score(m)為第m塊分區(qū)BVP 質(zhì)量最終分?jǐn)?shù),分?jǐn)?shù)越高說(shuō)明綜合質(zhì)量越好.

    如圖5 所示,對(duì)Score(m)數(shù)值排序后選取最佳分塊區(qū)域?qū)?yīng)BVP 信號(hào),對(duì)BVP 信號(hào)進(jìn)行帶通濾波,選取在0.8-3.5 Hz區(qū)間最大功率對(duì)應(yīng)的信號(hào)頻率fHR得到受試者心率.其中對(duì)應(yīng)關(guān)系如下:

    圖5 優(yōu)選BVP 信號(hào)及功率譜密度

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

    2.1 實(shí)驗(yàn)設(shè)備及數(shù)據(jù)

    本文實(shí)驗(yàn)和數(shù)據(jù)采集是在白色LED 照明環(huán)境中進(jìn)行.使用普通彩色攝像頭(羅技C920)連接微軟Surface 筆記本(Windows 10 ,16 GB 內(nèi)存,512 GB 硬盤)進(jìn)行面部視頻采集,同步使用指夾式脈搏傳感器(紫金電子 PWS-20D )采集左手食指指尖ppg 信號(hào).攝像頭設(shè)置采集分辨率為1080 P、幀率30 fps,每次采集受試者30 s 的人臉視頻,并以RGB 顏色空間保存成mp4 格式.脈搏傳感器設(shè)置采樣率為200 Hz,采樣精度為10 位作為脈搏真值保存.共有20 名受試者參與本次實(shí)驗(yàn),其中男性12 名、女性8 名,年齡分別為20-43 歲、19-35 歲.數(shù)據(jù)采集方式和采集界面如圖6 所示,受試者需正對(duì)攝像頭并且保持距離40-60 cm 左右而坐,同時(shí)采集面部視頻與指尖ppg 信號(hào).

    圖6 采集場(chǎng)景演示

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

    本文針對(duì)算法的效果復(fù)現(xiàn)了文獻(xiàn)[10-12] 中的3 種心率檢測(cè)算法作為對(duì)比,主要從實(shí)際場(chǎng)景下對(duì)算法進(jìn)行評(píng)估.其中采用了5 種誤差指標(biāo)對(duì)準(zhǔn)確性進(jìn)行比較.第1 個(gè)指標(biāo)為平均絕對(duì)誤差,代表心率檢測(cè)值與真實(shí)值之間誤差的絕對(duì)平均值,記為Me; 第2 個(gè)指標(biāo)為誤差標(biāo)準(zhǔn)差SDe,用來(lái)反映誤差之間相互離散的程度;第3 個(gè)指標(biāo)為誤差均方根記作RMSE,用來(lái)評(píng)估面部測(cè)量值和指尖ppg 參考值之間的偏離程度; 第4 個(gè)指標(biāo)為心率平均準(zhǔn)確率HRac,主要是反映心率檢測(cè)結(jié)果的平均準(zhǔn)確性,具體公式如式(15):

    其中,T為實(shí)驗(yàn)樣本總數(shù),為第t組面部心率檢測(cè)的心率測(cè)量值,為第t組脈搏傳感器心率測(cè)量的參考值; 第5 個(gè)參數(shù)是皮爾森相關(guān)系數(shù)記作r,是用來(lái)反應(yīng)心率測(cè)量值和真實(shí)值之間線性相關(guān)性的統(tǒng)計(jì)量.r值通常范圍在-1 至1 之間,且絕對(duì)值越大說(shuō)明兩變量之間的相關(guān)程度越高.

    為了全面驗(yàn)證本文方法在實(shí)際場(chǎng)景下的性能,設(shè)計(jì)了3 個(gè)典型場(chǎng)景驗(yàn)證算法的準(zhǔn)確性.首先分析受試者保持面部靜止場(chǎng)景,其次是受試者在鏡頭內(nèi)存在頭部運(yùn)動(dòng)和表情變化等動(dòng)態(tài)干擾場(chǎng)景.隨后使用不同光照?qǐng)鼍跋乱曨l驗(yàn)證算法的穩(wěn)定性,最后對(duì)受試者進(jìn)行10 min 的連續(xù)心率檢測(cè).

    2.2.1 靜止?fàn)顟B(tài)下心率檢測(cè)

    在本次實(shí)驗(yàn)中,為檢驗(yàn)靜止?fàn)顟B(tài)下本文算法心率檢測(cè)的準(zhǔn)確性,受試者需端坐在采集攝像頭前50 cm左右位置,保持環(huán)境光照穩(wěn)定及受試者頭部靜止采集2 組人臉視頻.每組視頻時(shí)長(zhǎng)為30 s,實(shí)驗(yàn)樣本總數(shù)為40 組.實(shí)驗(yàn)結(jié)果如表1 所示,其中單位為bpm (beat per minute)即心臟每分鐘跳動(dòng)次數(shù).

    表1 靜態(tài)場(chǎng)景下算法準(zhǔn)確性分析

    如表1 所示,展示了4 種方法在受試者保持靜止時(shí)的準(zhǔn)確性情況.可以發(fā)現(xiàn)4 種方法都可以很好的應(yīng)對(duì)靜態(tài)場(chǎng)景下的心率檢測(cè),其中,Me均低于1.6 bpm,SDe和RMSE低于2 bpm,HRac可以達(dá)到98% ,r接近1.通過觀察可以發(fā)現(xiàn),本文方法中的Me最低、r最高,分析原因可能是: 在視頻采集中,受試者雖然身體保持靜止但是由于呼吸動(dòng)作的存在會(huì)引入輕微的起伏干擾,其次眨眼等動(dòng)作的存在也會(huì)導(dǎo)致面部某些區(qū)域存在輕微的干擾.而本文采用的人臉自適應(yīng)超像素分割和綜合指標(biāo)ROI 優(yōu)選方法可以從面部區(qū)域中自動(dòng)尋找最佳區(qū)域進(jìn)行心率檢測(cè),表現(xiàn)出來(lái)更好的準(zhǔn)確性.

    2.2.2 動(dòng)態(tài)干擾場(chǎng)景下心率檢測(cè)

    由于檢測(cè)過程中,頭部的晃動(dòng)和表情變化等,會(huì)使得鏡頭內(nèi)人臉位置和背景發(fā)生變化造成皮膚與光源反射角度的變化,最終導(dǎo)致BVP 信號(hào)的失真進(jìn)一步影響心率檢測(cè)的準(zhǔn)確性.因此本文特地設(shè)計(jì)一組頭部移動(dòng)和表情變化干擾場(chǎng)景.在這次的實(shí)驗(yàn)中,運(yùn)動(dòng)干擾包括在鏡頭內(nèi)以一定的速率水平或豎直移動(dòng)頭部以及進(jìn)行平面內(nèi)旋轉(zhuǎn)等操作; 表情變化為視頻采集時(shí)受試者隨機(jī)朗讀1-9 號(hào)數(shù)字.在本組實(shí)驗(yàn)中,每名受試者采集3 組數(shù)據(jù),樣本總數(shù)為60 組.具體的實(shí)驗(yàn)結(jié)果如表2 所示.為了直觀展示實(shí)驗(yàn)結(jié)果,實(shí)驗(yàn)數(shù)據(jù)通過散點(diǎn)圖和Bland-Altman 圖進(jìn)行展示,如圖7 所示.

    圖7 運(yùn)動(dòng)結(jié)果散點(diǎn)和Bland-Altman

    由表2 結(jié)果可以得知,由于運(yùn)動(dòng)干擾和表情變化的引入導(dǎo)致4 種方法的準(zhǔn)確性均表現(xiàn)出不同程度的降低.經(jīng)過對(duì)比發(fā)現(xiàn),文獻(xiàn)[10-11]方法的性能明顯偏低,其中,Me、SDe、RMSE均大于 7 bpm,HRac下降至90% 以下,相比較靜止?fàn)顟B(tài)下最差下降幅度高達(dá)11.2%,與心率的相關(guān)性明顯降低; 文獻(xiàn)[12]的效果雖然有稍許提升,但是Me依然大于5 bpm 超過實(shí)際場(chǎng)景下心率檢測(cè)的誤差上限.而本文的方法,在4 種方法中受到運(yùn)動(dòng)干擾和表情變化的干擾幅度最小,尤其是心率準(zhǔn)確性指標(biāo)HRac達(dá)到了95.34%相關(guān)系數(shù)r高達(dá)0.956.上述數(shù)據(jù)表明本文的人臉自適應(yīng)超像素分割方法可以有效地從運(yùn)動(dòng)的人臉中提取信號(hào)最佳的區(qū)域進(jìn)行心率估計(jì),具有良好的性能表現(xiàn).

    表2 動(dòng)態(tài)干擾場(chǎng)景下算法準(zhǔn)確性分析

    通過對(duì)比圖7 第1 行中4 種方法的散點(diǎn)圖可以看出文獻(xiàn)[10-12]的方法存在較多的數(shù)據(jù)誤差較大明顯偏離線性回歸線,而本文方法的數(shù)據(jù)集中在線性回歸線附近,且線性回歸線的斜率接近于1,表明本文方法的測(cè)量值與心率參考值的相關(guān)最強(qiáng).

    圖7 第2 行是運(yùn)動(dòng)結(jié)果的Bland-Altman 圖,其在醫(yī)學(xué)領(lǐng)域經(jīng)常被用來(lái)分析新方法與標(biāo)準(zhǔn)方法測(cè)量結(jié)果的一致程度.虛線內(nèi)的區(qū)域是95%置信區(qū)間,中心實(shí)線代表心率測(cè)量值與參考值之間的相對(duì)平均誤差.從圖中可知,本文方法的虛區(qū)域遠(yuǎn)遠(yuǎn)窄于其他3 種方法且中心實(shí)線最接近于0,表明本文方法較其他3 種方法的測(cè)量值與標(biāo)準(zhǔn)參考值有著更好的一致性.進(jìn)一步分析原因可能如下:

    1)文獻(xiàn)[10]采用的自適應(yīng)濾波算法,能夠較好地處理高斯白噪聲干擾.但由于使用的是整張人臉作為ROI,而運(yùn)動(dòng)時(shí)人臉引起的尖銳干擾并不能被自適應(yīng)濾波器很好地濾除.因此該方法無(wú)法提取到純凈的BVP 信號(hào),進(jìn)而無(wú)法在運(yùn)動(dòng)干擾條件下進(jìn)行準(zhǔn)確的心率估計(jì).

    2)文獻(xiàn)[11]采用的多區(qū)域分析方法,主要選取額頭和兩頰3 個(gè)獨(dú)立區(qū)域作為心率檢測(cè)的ROI.該方法從一定程度上可以從運(yùn)動(dòng)干擾中提取優(yōu)質(zhì)BVP 信號(hào)進(jìn)行心率估計(jì).但是由于選取的是固定區(qū)域,固定區(qū)域無(wú)法確保不受到偶然干擾影響,并且額頭部位在女性中容易受到劉海覆蓋干擾,影響實(shí)際應(yīng)用場(chǎng)景下的表現(xiàn).

    3)文獻(xiàn)[12]也是通過將面部進(jìn)行分割提取不同區(qū)塊的色度特征信號(hào)進(jìn)行心率估計(jì).但是,首先其對(duì)人臉區(qū)域僅進(jìn)行簡(jiǎn)單的矩形分割沒有考慮區(qū)域內(nèi)像素間關(guān)系; 其次由于固定了分割矩形塊大小為20×20 的像素區(qū)域,并沒有對(duì)不同幀內(nèi)人臉圖像尺寸進(jìn)行一致性處理,無(wú)法確保運(yùn)動(dòng)中分塊個(gè)數(shù)的穩(wěn)定; 最后由于沒有對(duì)人臉旋轉(zhuǎn)進(jìn)行校正,難以保證傾斜狀態(tài)下分塊位置的連續(xù)性.

    4)本文方法與文獻(xiàn)[12]采取的方法類似,都是通過將臉部區(qū)域分割成不同的子塊通過對(duì)子塊信號(hào)的分析進(jìn)行心率估計(jì).不同的是,本文在簡(jiǎn)單矩形分割的基礎(chǔ)上引入了自適應(yīng)超像素聚類算法,確保區(qū)域內(nèi)像素間的特征相關(guān)性.同時(shí)本文算法針對(duì)人臉視頻分割后分塊位置的一致性和連續(xù)性進(jìn)行了人臉校正操作.一方面提升原始BVP 信號(hào)的提取質(zhì)量,另一方面通過綜合指標(biāo)優(yōu)選ROI 方法進(jìn)一步提升了算法在實(shí)際場(chǎng)景中的抗干擾能力.

    2.2.3 不同光照?qǐng)鼍跋滦穆蕶z測(cè)

    考慮到真實(shí)場(chǎng)景中視頻采集環(huán)境的復(fù)雜性,位于不同光照條件下心率檢測(cè)的準(zhǔn)確性會(huì)受到影響.因此本次實(shí)驗(yàn)分別采集了受試者靜止?fàn)顟B(tài)4 種光照?qǐng)鼍跋碌囊曨l進(jìn)行對(duì)比實(shí)驗(yàn),具體為面向光源(正)、側(cè)面對(duì)光源(側(cè))、逆向光源(逆)、以及暗光條件(暗).在本次實(shí)驗(yàn)中,每名受試者均在特定場(chǎng)景下采集2 組視頻,視頻時(shí)長(zhǎng)30 s,樣本總數(shù)為160 組.實(shí)驗(yàn)的結(jié)果如圖8 所示.

    圖8 不同光照?qǐng)鼍跋碌男穆蕼?zhǔn)確性變化

    由圖8 可知,在正面向光源的場(chǎng)景中4 種方法的準(zhǔn)確率均可以達(dá)到97%以上,表明在光照條件良好的狀態(tài)下4 種方法均可以準(zhǔn)確地進(jìn)行心率檢測(cè).但是針對(duì)側(cè)光和逆光這類光照分布不均勻的場(chǎng)景,文獻(xiàn)[10-12]方法的準(zhǔn)確性均有較大幅度的下降,其中文獻(xiàn)[10]更是分別下降了7.8%、12.4%.本文方法展示出較高的準(zhǔn)確性優(yōu)勢(shì),其中側(cè)光準(zhǔn)確率達(dá)到94.3%,逆光準(zhǔn)確率達(dá)到92.1%,對(duì)比最差方法提升了5.2%與8.2%,但是針對(duì)暗光環(huán)境4 種方法均表現(xiàn)不佳.分析原因?yàn)閭?cè)光、逆光時(shí)光照分布不均勻,本文方法可以根據(jù)光照分布的情況從面部選取像素特征相近區(qū)塊進(jìn)行分析并挑選質(zhì)量最佳子塊進(jìn)行心率估計(jì).面對(duì)暗光這樣的光照?qǐng)鼍?由于射入皮膚的光照強(qiáng)度變低導(dǎo)致面部反射的光強(qiáng)變低,無(wú)法從ROI 挑選上對(duì)光照進(jìn)行增強(qiáng),故以上方法在暗光場(chǎng)景中均表現(xiàn)欠佳.

    2.2.4 連續(xù)時(shí)間心率檢測(cè)

    長(zhǎng)時(shí)間連續(xù)的心率檢測(cè),可以用來(lái)評(píng)估用戶的心率變化趨勢(shì)以及對(duì)一些心血管疾病的產(chǎn)生進(jìn)行早期預(yù)警.為了最終評(píng)價(jià)本算法的準(zhǔn)確性,設(shè)計(jì)一組連續(xù)10 min的心率檢測(cè)實(shí)驗(yàn).受試者正坐在攝像機(jī)前50 cm 的位置,保持10 min 內(nèi)面部一直處于視頻采集畫面內(nèi),期間不限制受試者頭部移動(dòng)和表情等變化.并同時(shí)使用脈搏血傳感器每10 s 記錄一次實(shí)時(shí)心率值做對(duì)比分析,實(shí)驗(yàn)結(jié)果如圖9 所示.

    圖9 心率變化折線圖

    在圖9 中的紅線和藍(lán)線分別代表指尖PPG 傳感器采集的心率參考值和面部心率檢測(cè)的測(cè)量值.從整體上看,藍(lán)線的變化趨勢(shì)隨著紅線有著相似的起伏波動(dòng),說(shuō)明測(cè)量值和參考值之間有著良好的一致性(r=0.9558).綜上所述,本文的算法可以有效地進(jìn)行心率長(zhǎng)期檢測(cè),提升真實(shí)場(chǎng)景下算法的魯棒性.

    3 結(jié)論與展望

    近年來(lái)由于視頻圖像處理技術(shù)的快速發(fā)展,基于視頻的非接觸心率檢測(cè)技術(shù)發(fā)展迅速,但在實(shí)際場(chǎng)景應(yīng)用仍有諸多挑戰(zhàn).本文針對(duì)面部運(yùn)動(dòng)和光照變化干擾問題,提出一種人臉自適應(yīng)分割和多區(qū)域多指標(biāo)綜合評(píng)價(jià)算法,提升脈搏波質(zhì)量,提高心率檢測(cè)可靠性,增強(qiáng)實(shí)際場(chǎng)景魯棒性.主要由以下幾部分組成: 面部區(qū)域自適應(yīng)超像素分割為互不重疊的子塊; 提取色度特征構(gòu)建脈搏信號(hào)矩陣; 根據(jù)綜合指標(biāo)評(píng)估各分塊質(zhì)量,最后對(duì)最優(yōu)子塊進(jìn)行心率估計(jì).實(shí)驗(yàn)結(jié)果表明,相比傳統(tǒng)ROI 分割方法,本文算法分割后脈搏矩陣信號(hào)質(zhì)量有顯著提升.在面部運(yùn)動(dòng)和光照不均勻場(chǎng)景下心率檢測(cè)準(zhǔn)確性和穩(wěn)定性均有明顯提升.然而,在暗光和環(huán)境光照變化劇烈等場(chǎng)景下,本文心率檢測(cè)算法性能出現(xiàn)比較明顯的下降.在以后的工作中,將對(duì)環(huán)境光照變化造成的干擾問題進(jìn)行研究,進(jìn)一步提升復(fù)雜光照?qǐng)鼍跋碌男穆蕶z測(cè)魯棒性.

    猜你喜歡
    區(qū)域信號(hào)檢測(cè)
    “不等式”檢測(cè)題
    “一元一次不等式”檢測(cè)題
    “一元一次不等式組”檢測(cè)題
    信號(hào)
    鴨綠江(2021年35期)2021-04-19 12:24:18
    完形填空二則
    基于FPGA的多功能信號(hào)發(fā)生器的設(shè)計(jì)
    電子制作(2018年11期)2018-08-04 03:25:42
    小波變換在PCB缺陷檢測(cè)中的應(yīng)用
    關(guān)于四色猜想
    分區(qū)域
    基于LabVIEW的力加載信號(hào)采集與PID控制
    国产一区二区三区av在线 | 中文字幕av成人在线电影| 国产精品无大码| 国产三级在线视频| 午夜福利高清视频| 国产免费一级a男人的天堂| 日本-黄色视频高清免费观看| 亚洲中文字幕日韩| 亚洲国产高清在线一区二区三| 看免费成人av毛片| 91精品国产九色| 真人一进一出gif抽搐免费| 亚洲性夜色夜夜综合| 免费电影在线观看免费观看| 日本五十路高清| 可以在线观看毛片的网站| 真人一进一出gif抽搐免费| 精华霜和精华液先用哪个| 精品一区二区三区人妻视频| 少妇熟女aⅴ在线视频| 波多野结衣巨乳人妻| 久久精品国产鲁丝片午夜精品 | 亚洲va日本ⅴa欧美va伊人久久| 精品久久久噜噜| 中出人妻视频一区二区| 亚洲国产色片| 在线观看美女被高潮喷水网站| 香蕉av资源在线| 国产三级在线视频| 三级国产精品欧美在线观看| 老熟妇乱子伦视频在线观看| 国产色婷婷99| 精品久久久久久久久久久久久| 国产中年淑女户外野战色| 国模一区二区三区四区视频| 少妇裸体淫交视频免费看高清| 国产精品人妻久久久影院| 精品久久久久久久久亚洲 | 久久久久久久午夜电影| 18禁黄网站禁片免费观看直播| 国产免费一级a男人的天堂| 国产白丝娇喘喷水9色精品| 国产在线男女| 看片在线看免费视频| 99热这里只有精品一区| 韩国av在线不卡| 午夜福利视频1000在线观看| 少妇裸体淫交视频免费看高清| 国产成人福利小说| www日本黄色视频网| 日本免费一区二区三区高清不卡| 国产黄片美女视频| 韩国av一区二区三区四区| 久久精品国产亚洲av天美| 免费在线观看成人毛片| 亚洲精华国产精华精| 美女xxoo啪啪120秒动态图| 日本精品一区二区三区蜜桃| 精品久久久久久久久亚洲 | 国产成人影院久久av| 欧美日韩乱码在线| 亚洲精品国产成人久久av| 99九九线精品视频在线观看视频| 欧美成人a在线观看| 最近在线观看免费完整版| 成人鲁丝片一二三区免费| 亚洲图色成人| 免费一级毛片在线播放高清视频| 两性午夜刺激爽爽歪歪视频在线观看| 精品不卡国产一区二区三区| 亚洲 国产 在线| 欧美bdsm另类| 亚洲国产欧洲综合997久久,| 琪琪午夜伦伦电影理论片6080| 国产男人的电影天堂91| 亚州av有码| 长腿黑丝高跟| 国产一区二区激情短视频| 男女边吃奶边做爰视频| 国产一区二区激情短视频| 美女高潮喷水抽搐中文字幕| 国产伦精品一区二区三区视频9| 日本精品一区二区三区蜜桃| 精品一区二区免费观看| 亚洲18禁久久av| 国产精品三级大全| 国产一区二区三区在线臀色熟女| 噜噜噜噜噜久久久久久91| 最近中文字幕高清免费大全6 | 女人被狂操c到高潮| 亚洲成人久久爱视频| 简卡轻食公司| av女优亚洲男人天堂| 亚洲四区av| 桃红色精品国产亚洲av| 国产精品久久久久久av不卡| 少妇猛男粗大的猛烈进出视频 | 成人无遮挡网站| 男女做爰动态图高潮gif福利片| 亚洲成人久久性| 特级一级黄色大片| 特级一级黄色大片| 熟妇人妻久久中文字幕3abv| 国产伦在线观看视频一区| 国产精品人妻久久久久久| 18+在线观看网站| 亚洲国产日韩欧美精品在线观看| 在线国产一区二区在线| 国产白丝娇喘喷水9色精品| 永久网站在线| 国产白丝娇喘喷水9色精品| 亚洲专区国产一区二区| 伦精品一区二区三区| 啦啦啦观看免费观看视频高清| 亚洲图色成人| 老司机午夜福利在线观看视频| 国产人妻一区二区三区在| 欧美xxxx黑人xx丫x性爽| 国产精品女同一区二区软件 | 欧美极品一区二区三区四区| 国内精品久久久久久久电影| 春色校园在线视频观看| 69av精品久久久久久| av在线天堂中文字幕| 国产成年人精品一区二区| 免费av毛片视频| 精华霜和精华液先用哪个| 2021天堂中文幕一二区在线观| 欧美一区二区精品小视频在线| 久久精品91蜜桃| 一个人看的www免费观看视频| 一个人看视频在线观看www免费| 国产精品不卡视频一区二区| 看片在线看免费视频| a级一级毛片免费在线观看| 看片在线看免费视频| 在线观看舔阴道视频| 亚洲va日本ⅴa欧美va伊人久久| 国内精品一区二区在线观看| 不卡视频在线观看欧美| 老司机深夜福利视频在线观看| 白带黄色成豆腐渣| 九九在线视频观看精品| 国产精品一区二区性色av| 在线观看一区二区三区| 欧美性猛交╳xxx乱大交人| 麻豆成人午夜福利视频| 国产久久久一区二区三区| 91在线观看av| 女同久久另类99精品国产91| 日韩欧美 国产精品| 免费看a级黄色片| 色哟哟·www| 国产亚洲精品久久久com| 亚洲国产精品久久男人天堂| 在线播放无遮挡| 不卡一级毛片| 91av网一区二区| 亚洲国产欧洲综合997久久,| 两个人视频免费观看高清| 国产欧美日韩一区二区精品| 成年女人毛片免费观看观看9| 日本一本二区三区精品| 色综合色国产| 亚洲精品日韩av片在线观看| 欧美一区二区国产精品久久精品| videossex国产| 欧美色欧美亚洲另类二区| 久久国产精品人妻蜜桃| 色5月婷婷丁香| 天天一区二区日本电影三级| 亚洲第一区二区三区不卡| av中文乱码字幕在线| 美女高潮喷水抽搐中文字幕| 亚洲精华国产精华精| 日本欧美国产在线视频| 啦啦啦啦在线视频资源| 一卡2卡三卡四卡精品乱码亚洲| 联通29元200g的流量卡| 亚洲国产高清在线一区二区三| 国产精品久久久久久精品电影| 狂野欧美激情性xxxx在线观看| 精品欧美国产一区二区三| 精品久久久久久成人av| 欧美日本视频| 国产免费一级a男人的天堂| 国产高清激情床上av| 男女边吃奶边做爰视频| 欧美成人一区二区免费高清观看| 亚洲av不卡在线观看| 午夜视频国产福利| 国产高潮美女av| 久久久久久久久大av| 久久久久久久精品吃奶| 午夜激情福利司机影院| 99久久精品一区二区三区| 小说图片视频综合网站| 婷婷色综合大香蕉| 久久久精品欧美日韩精品| 亚洲精品久久国产高清桃花| 91麻豆av在线| 国产精品不卡视频一区二区| 97碰自拍视频| 深夜精品福利| 国产精华一区二区三区| 欧美+日韩+精品| 国产精品免费一区二区三区在线| 亚洲欧美精品综合久久99| 国产极品精品免费视频能看的| 国产精品女同一区二区软件 | 搡女人真爽免费视频火全软件 | 国产单亲对白刺激| 91久久精品电影网| 亚洲欧美日韩卡通动漫| 成人一区二区视频在线观看| 午夜a级毛片| 亚洲成人久久爱视频| 少妇高潮的动态图| 国产 一区精品| 日韩欧美三级三区| 中文亚洲av片在线观看爽| 搡女人真爽免费视频火全软件 | 午夜福利欧美成人| 在线观看美女被高潮喷水网站| 精品日产1卡2卡| 亚洲色图av天堂| 色吧在线观看| 欧美激情在线99| 中文字幕人妻熟人妻熟丝袜美| 国产熟女欧美一区二区| 精品一区二区免费观看| 日本a在线网址| 久99久视频精品免费| 国产 一区 欧美 日韩| 国产视频一区二区在线看| 香蕉av资源在线| 丰满人妻一区二区三区视频av| 日本色播在线视频| 俄罗斯特黄特色一大片| 成熟少妇高潮喷水视频| 免费观看在线日韩| 日韩强制内射视频| 国产人妻一区二区三区在| 亚洲人成网站高清观看| 91久久精品电影网| 免费观看人在逋| 久久草成人影院| 一个人看视频在线观看www免费| 亚洲精品日韩av片在线观看| 最近最新免费中文字幕在线| 最近最新中文字幕大全电影3| 成人特级黄色片久久久久久久| 国产精品国产三级国产av玫瑰| 久久久久久久精品吃奶| 国产麻豆成人av免费视频| 国内精品久久久久精免费| 欧美日韩综合久久久久久 | 久久精品国产清高在天天线| 欧美3d第一页| 又粗又爽又猛毛片免费看| 国产精品久久视频播放| 黄色一级大片看看| 俄罗斯特黄特色一大片| 99热这里只有是精品50| 精品一区二区免费观看| 中文字幕免费在线视频6| 神马国产精品三级电影在线观看| 毛片女人毛片| 在线天堂最新版资源| 日本爱情动作片www.在线观看 | 午夜福利高清视频| 在现免费观看毛片| 夜夜看夜夜爽夜夜摸| а√天堂www在线а√下载| 免费一级毛片在线播放高清视频| 免费av毛片视频| 男人狂女人下面高潮的视频| 午夜福利高清视频| 国产欧美日韩一区二区精品| 日本成人三级电影网站| 久久天躁狠狠躁夜夜2o2o| 搞女人的毛片| 国产精品女同一区二区软件 | 狠狠狠狠99中文字幕| 亚洲av中文字字幕乱码综合| 变态另类丝袜制服| 黄色丝袜av网址大全| 久久国产乱子免费精品| 欧美3d第一页| 午夜福利18| 精品人妻熟女av久视频| 99久久久亚洲精品蜜臀av| 最后的刺客免费高清国语| 熟妇人妻久久中文字幕3abv| 舔av片在线| 精品久久久噜噜| 日韩欧美一区二区三区在线观看| 亚洲国产精品久久男人天堂| 婷婷色综合大香蕉| 成年免费大片在线观看| av在线观看视频网站免费| 男人和女人高潮做爰伦理| 亚洲av二区三区四区| 制服丝袜大香蕉在线| 国产欧美日韩精品一区二区| 亚洲av.av天堂| 精品午夜福利视频在线观看一区| 精品久久久久久,| 国产成人aa在线观看| 成年女人看的毛片在线观看| 中国美女看黄片| 欧美激情国产日韩精品一区| 欧美一区二区亚洲| 亚洲无线观看免费| 久久久精品欧美日韩精品| 日日夜夜操网爽| 日韩欧美三级三区| 亚洲av成人精品一区久久| 国产精品久久久久久精品电影| 人妻久久中文字幕网| 人人妻人人澡欧美一区二区| 97超级碰碰碰精品色视频在线观看| 老熟妇仑乱视频hdxx| 毛片女人毛片| 亚洲一区二区三区色噜噜| 国产成人影院久久av| 国产精品98久久久久久宅男小说| 啪啪无遮挡十八禁网站| 亚洲经典国产精华液单| 精品久久久久久久末码| av福利片在线观看| 成人无遮挡网站| 亚洲经典国产精华液单| 国产精品久久视频播放| 国产老妇女一区| 欧美成人a在线观看| 国内精品一区二区在线观看| 69av精品久久久久久| 搞女人的毛片| 国产高清有码在线观看视频| 成人国产综合亚洲| 亚洲av一区综合| 九色国产91popny在线| 亚洲综合色惰| 久久国产精品人妻蜜桃| 亚洲自偷自拍三级| 亚洲三级黄色毛片| 精品久久国产蜜桃| 中文字幕高清在线视频| 久久久精品欧美日韩精品| 成人精品一区二区免费| 国产一区二区亚洲精品在线观看| 日日摸夜夜添夜夜添小说| 免费不卡的大黄色大毛片视频在线观看 | 亚洲va在线va天堂va国产| h日本视频在线播放| 一夜夜www| 欧美日韩国产亚洲二区| 舔av片在线| 人人妻人人看人人澡| 亚洲av日韩精品久久久久久密| 联通29元200g的流量卡| 国产精品日韩av在线免费观看| 国产探花在线观看一区二区| 久久精品影院6| 精品国内亚洲2022精品成人| 久久精品影院6| 精品久久久噜噜| 亚洲自偷自拍三级| 国内精品久久久久精免费| 88av欧美| 欧美最黄视频在线播放免费| 国产精品三级大全| 亚洲精华国产精华液的使用体验 | 99久久中文字幕三级久久日本| 亚洲第一电影网av| 国产精品久久久久久精品电影| 久久精品国产亚洲av天美| 国产高清激情床上av| 亚洲久久久久久中文字幕| 国产一区二区在线观看日韩| 亚洲精品乱码久久久v下载方式| 日韩中字成人| 国产精品久久电影中文字幕| 色综合亚洲欧美另类图片| 蜜桃久久精品国产亚洲av| 日韩精品有码人妻一区| 欧洲精品卡2卡3卡4卡5卡区| 日韩欧美一区二区三区在线观看| 日本一本二区三区精品| 少妇的逼好多水| av在线老鸭窝| 日韩欧美三级三区| 欧美最黄视频在线播放免费| 综合色av麻豆| 欧美日韩国产亚洲二区| 校园人妻丝袜中文字幕| 亚洲最大成人中文| 性插视频无遮挡在线免费观看| 国产毛片a区久久久久| 国产伦人伦偷精品视频| 亚洲经典国产精华液单| 亚洲四区av| 成人毛片a级毛片在线播放| 国产成人aa在线观看| 国产精品久久久久久亚洲av鲁大| 成人二区视频| 国产精品av视频在线免费观看| 中文字幕免费在线视频6| 成人国产麻豆网| 久久久久久久久久久丰满 | 日本黄色片子视频| 日韩中字成人| 久久久久精品国产欧美久久久| 日日撸夜夜添| 黄色欧美视频在线观看| 狂野欧美激情性xxxx在线观看| 日韩大尺度精品在线看网址| 香蕉av资源在线| 91精品国产九色| 日韩精品青青久久久久久| а√天堂www在线а√下载| 天天躁日日操中文字幕| 熟妇人妻久久中文字幕3abv| 国产一区二区三区在线臀色熟女| 亚洲三级黄色毛片| 99热只有精品国产| 色尼玛亚洲综合影院| 欧美日韩精品成人综合77777| 久久99热6这里只有精品| 亚洲精品亚洲一区二区| 好男人在线观看高清免费视频| 免费av不卡在线播放| 免费观看人在逋| netflix在线观看网站| 有码 亚洲区| 日韩一区二区视频免费看| 欧美精品国产亚洲| 亚洲欧美日韩高清在线视频| 长腿黑丝高跟| 国产精品久久电影中文字幕| 免费看a级黄色片| 国产人妻一区二区三区在| 波野结衣二区三区在线| 国产高清有码在线观看视频| 日韩在线高清观看一区二区三区 | 俄罗斯特黄特色一大片| 精品人妻1区二区| 欧美成人a在线观看| 久久精品国产亚洲av涩爱 | av福利片在线观看| 免费大片18禁| 久久久久久久久久成人| 亚洲av熟女| 欧美激情国产日韩精品一区| 久久久久久大精品| 尾随美女入室| 嫩草影院新地址| 免费av不卡在线播放| 给我免费播放毛片高清在线观看| 国产爱豆传媒在线观看| 久久久国产成人精品二区| eeuss影院久久| 亚洲第一电影网av| 九色成人免费人妻av| 天堂av国产一区二区熟女人妻| 九九热线精品视视频播放| 欧美+日韩+精品| 久久久久久久精品吃奶| 精品一区二区三区视频在线观看免费| 熟女人妻精品中文字幕| 亚洲欧美日韩无卡精品| 日韩欧美在线二视频| 久久久久久久久久黄片| 国产伦精品一区二区三区四那| 国产精品女同一区二区软件 | 成年人黄色毛片网站| 亚洲专区中文字幕在线| 性欧美人与动物交配| 色噜噜av男人的天堂激情| 日韩欧美国产在线观看| 变态另类丝袜制服| 少妇熟女aⅴ在线视频| 亚洲精品久久国产高清桃花| 欧美+亚洲+日韩+国产| 午夜福利成人在线免费观看| 亚洲自偷自拍三级| 欧美+日韩+精品| 最后的刺客免费高清国语| 国产乱人视频| 人妻久久中文字幕网| 五月玫瑰六月丁香| 欧美+日韩+精品| www.www免费av| 国产精品综合久久久久久久免费| 老司机午夜福利在线观看视频| 99久久中文字幕三级久久日本| 亚洲中文日韩欧美视频| 免费人成视频x8x8入口观看| 国产高清视频在线播放一区| 老司机福利观看| 亚州av有码| 在线观看舔阴道视频| 欧美xxxx黑人xx丫x性爽| 中文字幕av成人在线电影| 欧美另类亚洲清纯唯美| 有码 亚洲区| 欧美日韩瑟瑟在线播放| 麻豆成人午夜福利视频| 97人妻精品一区二区三区麻豆| 中文在线观看免费www的网站| 亚洲无线在线观看| 天天一区二区日本电影三级| 亚洲欧美日韩高清专用| 国产极品精品免费视频能看的| 亚洲国产色片| 久久精品国产自在天天线| 国产在线精品亚洲第一网站| 我要搜黄色片| 国产极品精品免费视频能看的| 麻豆av噜噜一区二区三区| 久久精品国产清高在天天线| 亚洲成av人片在线播放无| 久久亚洲真实| 国内少妇人妻偷人精品xxx网站| 成年女人毛片免费观看观看9| 99热只有精品国产| 韩国av一区二区三区四区| 日韩高清综合在线| av在线观看视频网站免费| 乱码一卡2卡4卡精品| 欧美性猛交╳xxx乱大交人| 三级男女做爰猛烈吃奶摸视频| 无人区码免费观看不卡| 色综合亚洲欧美另类图片| 国产精品人妻久久久久久| videossex国产| 人人妻,人人澡人人爽秒播| 国产久久久一区二区三区| 亚洲av日韩精品久久久久久密| 人妻制服诱惑在线中文字幕| 成人无遮挡网站| 嫩草影视91久久| 别揉我奶头 嗯啊视频| 日韩强制内射视频| 在线观看66精品国产| 人妻少妇偷人精品九色| 国产一区二区在线观看日韩| 免费在线观看日本一区| 校园人妻丝袜中文字幕| 亚洲成人久久性| 亚洲人成网站在线播| 欧美一区二区国产精品久久精品| av福利片在线观看| 夜夜夜夜夜久久久久| 亚洲av免费在线观看| 久久精品国产亚洲av涩爱 | 日本 欧美在线| 日日啪夜夜撸| 久久国产精品人妻蜜桃| 日本a在线网址| 精品一区二区三区人妻视频| 干丝袜人妻中文字幕| 最新中文字幕久久久久| 午夜精品一区二区三区免费看| 黄色视频,在线免费观看| 午夜视频国产福利| 日韩av在线大香蕉| 日日啪夜夜撸| 久久久久久久亚洲中文字幕| 九九热线精品视视频播放| 内射极品少妇av片p| 我要搜黄色片| 久久精品91蜜桃| 国产成人av教育| 亚洲精品日韩av片在线观看| 日本撒尿小便嘘嘘汇集6| 久久亚洲精品不卡| 亚洲欧美日韩高清在线视频| 夜夜爽天天搞| 欧美一区二区精品小视频在线| 直男gayav资源| 国产高清视频在线观看网站| 伦理电影大哥的女人| videossex国产| av在线观看视频网站免费| 精品久久久噜噜| 我要搜黄色片| 国产精品综合久久久久久久免费| 老司机午夜福利在线观看视频| 免费不卡的大黄色大毛片视频在线观看 | 嫩草影视91久久| 中文字幕免费在线视频6| 国产一区二区三区在线臀色熟女| 波野结衣二区三区在线| 久久久色成人| av.在线天堂| 精品无人区乱码1区二区| 国产午夜精品久久久久久一区二区三区 | www日本黄色视频网| 欧美成人免费av一区二区三区| 在线观看一区二区三区| 性欧美人与动物交配| 久久国产乱子免费精品| 日本在线视频免费播放| 精品国产三级普通话版| 久久国产精品人妻蜜桃| 亚洲自偷自拍三级| 亚洲精华国产精华液的使用体验 | av专区在线播放| 欧美性猛交╳xxx乱大交人| 91狼人影院| 在线观看午夜福利视频| 老熟妇乱子伦视频在线观看| 在线观看午夜福利视频|