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

    基于相位光流法的結(jié)構(gòu)振動微小位移測量①

    2022-05-10 02:29:48王浩林
    計算機系統(tǒng)應(yīng)用 2022年2期
    關(guān)鍵詞:光流法光流像素

    王浩林,劉 源,王 蕓

    (長安大學(xué) 信息工程學(xué)院,西安 710064)

    隨著我國基礎(chǔ)設(shè)施建設(shè)的需求和數(shù)量的不斷擴大,巨大數(shù)量的基礎(chǔ)設(shè)施需要定期檢查評估以確保其安全性和可靠性.近年來利用結(jié)構(gòu)健康監(jiān)測(SHM)技術(shù)管理基礎(chǔ)設(shè)施成為研究人員關(guān)注的熱點.以模態(tài)分析為中心的基于結(jié)構(gòu)動力學(xué)的振動方法是民用、機械和航空結(jié)構(gòu)健康監(jiān)測、模型表征和更新最廣泛使用的技術(shù)之一[1].

    橋梁結(jié)構(gòu)的模態(tài)分析需要通過分析傳感器的振動數(shù)據(jù)來確定模態(tài)特性,傳統(tǒng)的傳感器多為接觸式傳感器,如加速度計,通過直接固定在結(jié)構(gòu)表面上進行測量.這種方法通常非常精確并具有高動態(tài)范圍,人們也可以通過使用應(yīng)變儀來測量動態(tài)應(yīng)變[2].但接觸式傳感器在許多人類難以到達的結(jié)構(gòu)處安裝困難、費時且昂貴,并且對于長期監(jiān)測而言,傳感器的傳輸電纜或電池維護也是一個挑戰(zhàn).最近非接觸式傳感器系統(tǒng)用于結(jié)構(gòu)振動測量,如微波遙感[3],數(shù)字圖像處理[4],激光多普勒技術(shù)[5,6]和基于視頻序列的光流分析[7].這些方法的優(yōu)點在于傳感器不需固定在結(jié)構(gòu)上,并且隨著消費級相機的普及和算法的改進,基于計算機視覺的位移測量技術(shù)越來越受到研究人員關(guān)注.

    目前在結(jié)構(gòu)位移檢測領(lǐng)域內(nèi)應(yīng)用較多的方法有模板匹配、數(shù)字圖像相關(guān)技術(shù)(DIC)、SIFT 特征點匹配、傳統(tǒng)Lucas-Kanade 光流神經(jīng)、基于深度學(xué)習(xí)的光流等方法.例如Ye 等人[8]通過在橋梁上使用專用的視覺標(biāo)記如二維碼或LED 燈進行模板匹配監(jiān)測大跨度橋梁的結(jié)構(gòu)位移,這種方法的限制是仍需在結(jié)構(gòu)表面添加人工目標(biāo).Mazzoleni 等人[9]通過DIC 第一次在土木結(jié)構(gòu)上估算人群引起的垂直動荷載.Khuc 等人[10]研究使用尺度不變特征變換SIFT 方法進行結(jié)構(gòu)位移測量,和傳統(tǒng)傳感器之間的最大差異小于5%.Zhao 等人[11]結(jié)合支持相關(guān)濾波算法與Kanade-Lucas-Tomasi (KLT)算法測量斜拉橋的位移,但僅在KLT 算法上較為魯棒.Dong 等人[7]提出一種基于神經(jīng)網(wǎng)絡(luò)Flownet2的全場光流法位移計算方法,比傳統(tǒng)的光流算法具有更高的精度.

    人眼感知運動的范圍是有限度的,對于視頻中許多幅度微小的運動,人眼無法感知.最近,一種用于可視化肉眼看不見的微小運動的運動放大算法受到關(guān)注,基于歐拉的運動放大和基于相位的運動放大算法類似于流體流動分析的歐拉框架[12–14],即跟蹤圖像中每一像素在不同幀中亮度的變化.Wadhwa 等人[13]通過復(fù)數(shù)可操縱金字塔將圖像分解為多尺度的相位和幅值圖,再通過放大系數(shù)放大不同幀之間的相位差,將放大后的相位疊加對圖像進行重構(gòu),放大了微小位移.2015年,Chen等人[15]結(jié)合基于局部相位的光流技術(shù)[16],將基于相位的運動放大算法拓展為一種可以量化微小運動的光流法,并用于懸臂梁的模態(tài)識別.2017年,Chen 等人[17]測量了校園建筑頂部距離超過175 m的天線塔,識別出0.21 mm的諧振振幅,小于1/170 像素.Cha 等人[18]通過卡爾曼濾波降低了基于相位的光流法的噪聲.Shang 等人[19]通過基于運動放大的光流法進行人行天橋的多點模態(tài)頻率和振型測量,誤差均小于4%.基于相位的光流法在微小位移測量方面顯示出巨大的應(yīng)用潛力.

    目前基于相位的光流法在橋梁結(jié)構(gòu)振動檢測方面應(yīng)用較少,本文受以上研究啟發(fā),并根據(jù)Diamond 等人[2]對基于相位運動放大的光流法精確度的驗證研究,提出了一種基于相位的改進光流法計算結(jié)構(gòu)微小位移的方法.本文選擇基于相位的運動放大方法中的構(gòu)造的復(fù)數(shù)可操縱金字塔中的第一層y方向濾波器來計算局部相位,利用局部相位輪廓不變性求得視頻中運動結(jié)構(gòu)的像素位移時間歷程,再通過比例因子將像素位移轉(zhuǎn)換為實際位移.并在實驗室通過搭建簡支梁模型進行振動測量實驗并與傳感器進行對比,驗證了本文方法的性能.

    1 結(jié)構(gòu)微小振動位移計算方法設(shè)計

    1.1 位移測量系統(tǒng)流程

    本文提出的基于相位的光流法計算結(jié)構(gòu)微小位移實驗流程如圖1所示,大致可分為4 個步驟:數(shù)據(jù)采集、數(shù)據(jù)預(yù)處理、基于相位的光流法計算位移和數(shù)據(jù)后處理.首先進行相機標(biāo)定,建立二維圖像和三維真實世界之間的對應(yīng)關(guān)系,并通過相機采集結(jié)構(gòu)振動視頻數(shù)據(jù).接著進行數(shù)據(jù)預(yù)處理,采集的視頻需要進行裁剪選擇出感興趣(region of interest,ROI)區(qū)域,目的是降低運算量.并將視頻拆解為幀,因為位移計算是計算第一幀和第n幀之間的位移.第3 步就是計算像素位移,采用改進的基于相位的光流法計算振動的像素位移.最后是數(shù)據(jù)的后處理,算法要對像素位移數(shù)據(jù)進行降噪處理,并將像素單位的位移轉(zhuǎn)換為物理單位的位移并進行數(shù)據(jù)結(jié)果分析.

    1.2 相機標(biāo)定

    相機標(biāo)定是為了估計投影過程.在基于視覺的結(jié)構(gòu)位移測量研究中,常用的方法主要有3 種:比例因子法、全投影矩陣法和平面單應(yīng)矩陣法[7].當(dāng)相機和鏡頭的軸垂直于運動平面時,比例因子計算為真實世界物體尺寸與圖像尺寸的比率.如果在基于視覺的位移測量系統(tǒng)中應(yīng)用比例因子,則徑向失真可以忽略不計的假設(shè)必須成立.Dong 等人[7]建議使用徑向畸變可忽略不計的變焦鏡頭,以避免現(xiàn)場應(yīng)用中繁瑣的攝像機固有矩陣校準(zhǔn)工作.本文中所使用的數(shù)碼相機具有失真校正功能,且本文實驗中相機鏡頭光軸垂直于物體運動平面,因此選擇比例因子法將運動結(jié)構(gòu)的像素運動轉(zhuǎn)換為毫米的物理單位.

    1.3 濾波器構(gòu)造

    可操縱金字塔(steerable pyramid)是一種過完備變換,根據(jù)空間比例、方向和位置分解圖像.變換的基函數(shù)類似于Gabor 小波,由高斯包絡(luò)加窗的正弦曲線,并且是可操縱的[20].可操縱金字塔的定向頻帶 ψl,θ中的頻域轉(zhuǎn)換函數(shù)是基本濾波器的不同尺度和不同方向的副本,由尺度l和方向θ 決定.如圖2所示,每一塊區(qū)域均代表一組濾波器副本.頻率軸范圍從?π 到 π,K=4,分為4 個方向和多層尺度[21].基本函數(shù)通過尺度變換和方向變換組成多組濾波器副本(初始高通帶和最終低通子帶除外).

    可控金字塔是通過遞歸地將圖像分割成一組定向子帶和一個低通殘留帶來實現(xiàn)的.轉(zhuǎn)換的系統(tǒng)圖如圖3所示,H表示高通濾波,L表示低通濾波,B表示帶通濾波,低通部分通過遞歸的方式進行子項分解.

    在傅里葉域中,濾波器構(gòu)造公式[22]如下:

    帶有徑向和角度部分:

    其中,r,θ是極坐標(biāo)下的極徑和極角,K表示總共分解的幾個方向.

    復(fù)數(shù)可操縱金字塔(complex steerable pyramid,CSP)的轉(zhuǎn)換函數(shù)只包含相應(yīng)的實部可控金字塔濾波器的正頻率,即 2 cos(ωx)=eiωx+e?iωx中eiωx的響應(yīng),因此存在幅度和相位的概念[13].

    本文通過時域卷積定理進行濾波,即時域內(nèi)的卷積對應(yīng)頻域內(nèi)的乘積,如式(3)所示:

    在基于相位的視頻運動放大中,局部相位需要進行放大疊加以形成重構(gòu)圖像.但在本文中,計算的是局部相位輪廓隨時間偏移形成的位移,實驗設(shè)計為y方向垂直振動的簡支梁振動模型,因此選擇復(fù)數(shù)可操縱金字塔第一層y方向的濾波器,即 θ=π/2.濾波器的可操縱性基本上沒有被使用,因為本文中主要計算的是縱向振動位移.

    1.4 基于相位的光流法

    基于相位的光流法是由Chen 等人[15]新提出的一種結(jié)構(gòu)位移計算方法,適用于無背景雜波的微小像素位移計算.Wadhwa 等人[13]首次提出了一種基于相位的運動放大算法,通過構(gòu)造的復(fù)數(shù)可操縱金字塔濾波器對圖像亮度進行濾波,計算出局部空間相位和幅值,局部空間相位信號在時間上被傅立葉變換分解成一系列代表諧波運動的正弦曲線.相位信號隨后在時間上被帶通濾波、放大并重新組合以形成運動放大視頻.Wadhwa 等人[14]又通過里斯變換進行改進加快了運算速度.Chen 等人[15]受基于相位的運動放大啟發(fā),利用局部相位在時間上的恒定輪廓的運動[16,23]來計算相應(yīng)的位移信號.

    空間局部相位和局部振幅是類似于傅立葉級數(shù)系數(shù)的相位和振幅的量.相位控制基函數(shù)的位置,而振幅控制其強度.在傅里葉變換的情況下,相位對應(yīng)于全局運動.局部相位提供了一種計算局部運動的方法[16].對于視頻數(shù)據(jù),如式(4)所示,由I(x,y,t0)表示在空間位置(x,y) 處和t0時刻的圖像亮度,通過使用復(fù)數(shù)濾波器Gθ2+iH2θ對視頻幀進行卷積運算,獲得時間t0的幀的方位θ的局部相位?θ(x,y,t0)和局部振幅Aθ(x,y,t0).

    位移信號是從局部相位在時間上的恒定輪廓的運動中獲得的[16,23],可以表示為:

    其中,m表示常數(shù),意味著隨著時間的推移相位的局部相位恒定輪廓不變,某個方向的位移信號就是局部相位輪廓在第一幀和當(dāng)前幀之間的移動距離,在實際計算中僅考慮和坐標(biāo)軸x,y平行的兩個方向的分量,式(6)是通過關(guān)于時間的微分獲得的方程式.

    其中,u和v表示圖像x,y軸方向上的速度.通過求解式(6)可得:

    其中,? ?0(x,y,t)/?y=0和? ?π/2(x,y,t)/?x=0.由于計算相鄰幀的速度進行積分得到的位移會累計誤差,因此選擇初始幀圖像和當(dāng)前幀圖像之間進行計算,本文對原始公式進行優(yōu)化變形,計算結(jié)構(gòu)振動在y方向的位移?y,得到式(8):

    其中,η表示一個極小的正數(shù),目的是防止分母為0.α表示位移系數(shù),用于將速度轉(zhuǎn)換為位移,β為下采樣倍率,根據(jù)Diamond 等人[2]的實驗,在應(yīng)用復(fù)數(shù)濾波器之前,對輸入圖像執(zhí)行下采樣.這樣做是為了降低圖像中的噪聲.Chen 等人建議選擇視頻中結(jié)構(gòu)邊緣或者有紋理區(qū)域中的點進行運算[15],這是因為無紋理區(qū)域的運動是局部模糊的,在本文中不做討論.

    1.5 數(shù)據(jù)后處理

    通過多次實驗發(fā)現(xiàn),在算法計算過程中,由于圖像噪聲的影響,位移時程數(shù)據(jù)中會出現(xiàn)某些超出正常范圍值得尖峰噪聲,因此選擇中值濾波進行降噪處理,濾波器尺寸為1×3.最后通過比例因子將像素單位的位移轉(zhuǎn)換為以毫米為單位的真實物理位移,并于其他方法進行比較.在接下來的章節(jié)中,將通過實驗室實驗來驗證所提出的改進算法的可行性和性能.

    2 實驗與結(jié)果分析

    2.1 實驗設(shè)置

    本文的實驗環(huán)境為實驗室搭建的簡支梁模型,如圖4所示,圖4(a)為簡支梁模型,左側(cè)為固定端,右側(cè)為滑動端;圖4(b)為傳感器固定位置,為簡支梁左側(cè)1/4 位置處,傳感器下方箭頭指向的像素點為選擇進行計算的點,位于簡支梁的邊緣,將傳感器的數(shù)據(jù)作為參考標(biāo)準(zhǔn);圖4(c)為激勵源,使用橡皮錘在簡支梁右側(cè)10 cm 處輕輕敲擊,產(chǎn)生振動源.兩個底座固定在同一木板上,簡支梁是一組總長150 cm,寬30 mm,厚度為2 mm的鋼條組成,在簡支梁中部可通過改變螺絲松緊度來模擬橋梁損傷,但本文實驗?zāi)繕?biāo)是驗證振動位移測量方法的性能,因此將螺絲旋轉(zhuǎn)至固定.

    本文中進行視頻采集的設(shè)備為Canon EOS-1D X Mark II 相機,攝像模式參數(shù)為1920×1080的分辨率和120 fps的幀率,光照條件為穩(wěn)定的室內(nèi)光源.相機通過三腳架和云臺固定,與簡支梁和傳感器處于同一水平面上,且鏡頭正對簡支梁和傳感器,距離簡支梁模型1.3 m,以免因夾角產(chǎn)生計算誤差.本文中使用的為891-4 型超低頻測振儀垂直傳感器,傳感器的采樣頻率為200 Hz,用來測量豎直方向的速度作為參考數(shù)據(jù),并在最后的實驗比對階段進行積分處理,振動數(shù)據(jù)采集設(shè)備如圖5所示,圖5(a)為視頻采集設(shè)備,圖5(b)為振動傳感器數(shù)據(jù)采集儀和采集軟件.

    本文使用的位移真值近似值由用于超低頻或低頻振動測量的891-4 型拾振器獲得,通過傳感器測量的速度積分計算位移,該傳感器的速度分辨率為1.6×10–6m/s,速度靈敏度為1.738 mV/mm/s,位移分辨率為1.6×10–6m,傳感器采樣率為200 Hz.本文認為由891-4型拾振器采集的原始位移數(shù)據(jù)屬于更高一級的位移測量方法,因此選擇此數(shù)據(jù)作為振動位移的真值近似值.

    2.2 實驗結(jié)果與分析

    為了驗證本文方法的性能,選取了兩種基于神經(jīng)網(wǎng)絡(luò)的光流方法進行對比,文獻[24]中的Flownet2.0網(wǎng)絡(luò)是一種經(jīng)典的基于神經(jīng)網(wǎng)絡(luò)的光流方法,性能追平了傳統(tǒng)光流方法,文獻[25]中的基于循環(huán)全對場變換(recurrent all-pairs field transforms,RAFT)的神經(jīng)網(wǎng)絡(luò)是截止2020年最為先進的深度學(xué)習(xí)光流方法,兩種網(wǎng)絡(luò)具體結(jié)構(gòu)見文獻[24,25],幾種方法均選擇同一像素進行位移計算.

    在簡支梁上選取點并放置振動傳感器作為參照,使用本文方法及深度學(xué)習(xí)光流方法計算簡支梁振動位移,本文方法與其他光流方法的數(shù)據(jù)可視化對比如圖6所示,橫坐標(biāo)為圖像幀數(shù),縱坐標(biāo)為真實位移,單位為mm.視頻中簡支梁結(jié)構(gòu)振動振幅不超過1 個像素.將傳感器數(shù)據(jù)降采樣為120 Hz,手動進行傳感器和相機的時間對齊,將兩者數(shù)據(jù)進行同步,選擇了1 300 個數(shù)據(jù)點進行對比.圖6(a)為總的1 300 個數(shù)據(jù),為了更直觀的展示方法差異,選擇部分?jǐn)?shù)據(jù)進行展示,圖6(b)為200–400幀數(shù)據(jù)可視化對比,圖6(c)為800–1 000 幀的可視化對比,從中可以看出本文方法位移時程序列明顯更加平滑.

    本文采用兩種評價參數(shù)作為性能指標(biāo),一是皮爾森相關(guān)系數(shù)(Pearson correlation coefficient),是兩組位移時程數(shù)據(jù)的協(xié)方差和標(biāo)準(zhǔn)差的商,用于衡量算法數(shù)據(jù)與傳感器數(shù)據(jù)的相關(guān)程度,相關(guān)系數(shù)越大,表示相關(guān)程度越高:

    二是平均絕對誤差(mean absolute error,MAE),用來度量兩條位移時程序列間各點的平均偏差程度,MAE越小,表示偏差程度越小:

    圖7展示了3 種方法與傳感器數(shù)據(jù)的誤差分析可視化圖,左側(cè)為絕對誤差的時程曲線,右側(cè)為絕對誤差的分布直方圖.方法的誤差服從正態(tài)分布,選擇誤差的均值、方差、和標(biāo)準(zhǔn)差進行分析.均值公式、方差、標(biāo)準(zhǔn)差公式分別為:

    Chen 等人建議在使用此方法時量程不要超過2 個像素單位[15],在本文實驗中,經(jīng)過驗證,本方法在位移為3 個像素單位內(nèi)測量較為準(zhǔn)確,因此本文方法量程為3 個像素單位,即單側(cè)位移不超過1.5 個像素單位.測量方法的準(zhǔn)確度等級為(最大絕對誤差/量程)×100,本文方法測量的最大絕對誤差為0.267 8 個像素單位,最大量程為3 個像素單位,因此本文方法的精確度等級為8.9.量程較小屬于方法的固有缺陷,這限制了本文方法的應(yīng)用場景,但本文方法的優(yōu)勢在于對于微小位移的測量的相關(guān)度和誤差分析結(jié)果均優(yōu)于現(xiàn)有光流方法.

    本文實驗通過橡皮錘敲擊簡支梁產(chǎn)生單次振動,并重復(fù)10 次實驗,每次振動的峰-峰值均接近最大量程,本文方法的10 次實驗結(jié)果如表1.每次均選擇從振動開始的1 300 個數(shù)據(jù)點(每秒120 個數(shù)據(jù)點)進行分析.

    表1 基于相位光流法的10 次振動實驗數(shù)據(jù)分析和誤差分析

    接著使用其他兩種方法對10 次實驗進行數(shù)據(jù)分析和誤差分析對比,結(jié)果如表2所示,本文的方法在穩(wěn)定光照條件下測得的簡支梁振動數(shù)據(jù)各項指標(biāo)均優(yōu)于其他兩種光流方法.

    表2 不同方法的實驗數(shù)據(jù)和誤差分析10 組數(shù)據(jù)均值

    綜合以上實驗結(jié)果,本文基于相位運動放大的結(jié)構(gòu)振動微小位移測量方法在微小振動位移的精確度和穩(wěn)定性上具有明顯優(yōu)勢,并且通過實驗室簡支梁振動實驗驗證了方法的有效性.

    3 結(jié)論與展望

    結(jié)構(gòu)健康監(jiān)測技術(shù)對于維護土木建筑結(jié)構(gòu)的安全和可靠性具有重要意義,基于計算機視覺的技術(shù)可以實現(xiàn)在較少用戶參與的情況下,實現(xiàn)土木結(jié)構(gòu)的非接觸位移監(jiān)測.本文提出了一種改進的基于相位的光流法,與傳統(tǒng)的位移和振動監(jiān)測測量技術(shù)相比,該方法需要更簡單、不太復(fù)雜且更具成本效益的硬件.本文將方法應(yīng)用于一個實驗室搭建的簡支梁模型,它很好地代表了許多類似的相對細長的結(jié)構(gòu).本文方法相比于其他光流方法在微小位移檢測上具有更好的精確度,但當(dāng)視頻中的結(jié)構(gòu)出現(xiàn)大的運動(大于3 個像素)或者光照劇烈變化時,該方法并不適用,這限制了方法的適用范圍.在未來的研究中,要加強方法對于外界干擾的魯棒性研究,擴展方法的適用范圍,提高方法的實用性.

    猜你喜歡
    光流法光流像素
    趙運哲作品
    藝術(shù)家(2023年8期)2023-11-02 02:05:28
    利用掩膜和單應(yīng)矩陣提高LK光流追蹤效果
    像素前線之“幻影”2000
    “像素”仙人掌
    基于物理學(xué)的改善粒子圖像測速穩(wěn)健光流方法研究
    基于背景分類的監(jiān)控視頻中的運動目標(biāo)檢測算法綜述
    Matlab下視頻處理系統(tǒng)設(shè)計與實現(xiàn)
    權(quán)重系數(shù)自適應(yīng)光流法運動目標(biāo)檢測
    高像素不是全部
    CHIP新電腦(2016年3期)2016-03-10 14:22:03
    初始段安控圖像目標(biāo)識別方法研究
    国产高清不卡午夜福利| 啦啦啦韩国在线观看视频| eeuss影院久久| .国产精品久久| 久久久欧美国产精品| 亚洲久久久久久中文字幕| 色吧在线观看| 三级国产精品欧美在线观看| 国产午夜福利久久久久久| 午夜福利网站1000一区二区三区| 久久99精品国语久久久| 亚洲av免费在线观看| 青青草视频在线视频观看| 婷婷色综合www| 简卡轻食公司| 狂野欧美激情性xxxx在线观看| 少妇人妻一区二区三区视频| 免费大片黄手机在线观看| 国产色婷婷99| 观看免费一级毛片| 午夜爱爱视频在线播放| 国产成人a∨麻豆精品| 美女xxoo啪啪120秒动态图| 街头女战士在线观看网站| 免费看不卡的av| 亚洲精品一区蜜桃| 91精品国产九色| 国产精品国产三级国产专区5o| 好男人在线观看高清免费视频| 亚洲自拍偷在线| 99视频精品全部免费 在线| 亚洲欧美中文字幕日韩二区| 丰满人妻一区二区三区视频av| 国产白丝娇喘喷水9色精品| 国产日韩欧美在线精品| 国产亚洲av嫩草精品影院| 国产av不卡久久| 亚洲av成人精品一二三区| 少妇人妻一区二区三区视频| 亚洲av国产av综合av卡| 国产极品天堂在线| 成人一区二区视频在线观看| 欧美xxxx性猛交bbbb| 欧美成人a在线观看| 一夜夜www| 亚洲18禁久久av| av播播在线观看一区| 国产精品久久久久久久电影| 久久久久久久久久人人人人人人| 国产精品一区二区三区四区免费观看| 成人午夜精彩视频在线观看| 高清av免费在线| 91久久精品电影网| 国产精品一区二区在线观看99 | 夜夜爽夜夜爽视频| 99久国产av精品| 国产精品爽爽va在线观看网站| 免费高清在线观看视频在线观看| 18+在线观看网站| 精品人妻一区二区三区麻豆| 青春草国产在线视频| 18禁裸乳无遮挡免费网站照片| 国产精品爽爽va在线观看网站| 亚洲伊人久久精品综合| 国产乱人视频| 黄色日韩在线| 在线天堂最新版资源| 免费观看无遮挡的男女| 久久这里只有精品中国| av在线观看视频网站免费| 成人毛片60女人毛片免费| 国产免费福利视频在线观看| 日日啪夜夜撸| 国产精品.久久久| 可以在线观看毛片的网站| 精品一区二区三卡| 亚洲av中文字字幕乱码综合| 久久精品国产鲁丝片午夜精品| 人妻少妇偷人精品九色| 国产精品久久视频播放| 国产精品女同一区二区软件| 国产精品女同一区二区软件| 免费观看无遮挡的男女| 欧美成人一区二区免费高清观看| 日韩欧美精品免费久久| 69av精品久久久久久| 五月玫瑰六月丁香| 少妇丰满av| 国产极品天堂在线| 免费看a级黄色片| 亚洲自拍偷在线| 亚洲成人中文字幕在线播放| 在线观看av片永久免费下载| .国产精品久久| 高清视频免费观看一区二区 | 精品酒店卫生间| 国产亚洲精品av在线| 免费黄频网站在线观看国产| 午夜福利成人在线免费观看| a级毛色黄片| 亚洲国产精品成人久久小说| 亚洲国产精品成人久久小说| 99热全是精品| 亚洲av不卡在线观看| 99久久人妻综合| 免费在线观看成人毛片| 插阴视频在线观看视频| 国模一区二区三区四区视频| 国产在视频线精品| 精品欧美国产一区二区三| 91精品伊人久久大香线蕉| 免费大片18禁| 搡老乐熟女国产| 日本免费a在线| 久久人人爽人人爽人人片va| 搡老乐熟女国产| 国产高清有码在线观看视频| 26uuu在线亚洲综合色| 汤姆久久久久久久影院中文字幕 | 伦理电影大哥的女人| av在线亚洲专区| 国产一级毛片七仙女欲春2| 男女下面进入的视频免费午夜| 亚洲综合色惰| 一级片'在线观看视频| 国产精品日韩av在线免费观看| 久久久久久久久久黄片| 国产真实伦视频高清在线观看| 亚洲av中文av极速乱| 中文天堂在线官网| 国产成人精品久久久久久| 亚洲成人av在线免费| 看免费成人av毛片| 熟妇人妻久久中文字幕3abv| 久久亚洲国产成人精品v| 国产激情偷乱视频一区二区| 又黄又爽又刺激的免费视频.| 欧美高清性xxxxhd video| 精品熟女少妇av免费看| 亚洲av成人精品一区久久| 久久亚洲国产成人精品v| 伦理电影大哥的女人| 日韩欧美 国产精品| 国产一级毛片七仙女欲春2| 精品欧美国产一区二区三| 男女下面进入的视频免费午夜| 少妇人妻一区二区三区视频| 国产精品一区二区三区四区久久| 国产av国产精品国产| 国产精品久久久久久av不卡| 一级毛片我不卡| 91在线精品国自产拍蜜月| 久久99热这里只频精品6学生| 校园人妻丝袜中文字幕| 日本免费a在线| 91在线精品国自产拍蜜月| 亚洲精品久久午夜乱码| 久久97久久精品| 久久久色成人| 哪个播放器可以免费观看大片| 欧美+日韩+精品| 国产伦在线观看视频一区| 少妇的逼水好多| 亚洲欧美日韩东京热| 国产精品一区二区在线观看99 | 国产乱人偷精品视频| 亚洲不卡免费看| 精品久久久久久成人av| 蜜桃久久精品国产亚洲av| 欧美成人精品欧美一级黄| 免费在线观看成人毛片| 欧美xxxx性猛交bbbb| 春色校园在线视频观看| 夫妻午夜视频| 久久精品久久久久久噜噜老黄| 国产乱人视频| 亚洲av中文av极速乱| 亚洲真实伦在线观看| 自拍偷自拍亚洲精品老妇| 91aial.com中文字幕在线观看| 国产精品无大码| 亚洲成人中文字幕在线播放| 亚洲精品成人久久久久久| eeuss影院久久| 婷婷六月久久综合丁香| 精品久久久久久电影网| 我的女老师完整版在线观看| 爱豆传媒免费全集在线观看| 色吧在线观看| 午夜激情久久久久久久| 国产精品三级大全| 亚洲婷婷狠狠爱综合网| 啦啦啦中文免费视频观看日本| 18禁在线无遮挡免费观看视频| 欧美xxxx黑人xx丫x性爽| 人妻系列 视频| 国产成人精品一,二区| 黄色配什么色好看| 日韩av不卡免费在线播放| 日日撸夜夜添| 久久久精品欧美日韩精品| 青青草视频在线视频观看| 午夜福利视频精品| 成人av在线播放网站| 精品久久久久久成人av| 久久亚洲国产成人精品v| 欧美日韩综合久久久久久| 免费无遮挡裸体视频| 大陆偷拍与自拍| 91午夜精品亚洲一区二区三区| 日本wwww免费看| 国产午夜精品久久久久久一区二区三区| 神马国产精品三级电影在线观看| 免费电影在线观看免费观看| 狂野欧美激情性xxxx在线观看| 91久久精品国产一区二区成人| 18禁裸乳无遮挡免费网站照片| 亚洲国产色片| 99久国产av精品| 又黄又爽又刺激的免费视频.| 亚洲精华国产精华液的使用体验| 国产白丝娇喘喷水9色精品| 国内精品一区二区在线观看| 一个人看的www免费观看视频| 亚洲第一区二区三区不卡| 国产乱人视频| av免费观看日本| 水蜜桃什么品种好| 国产精品久久视频播放| 日韩欧美精品免费久久| 国产黄色视频一区二区在线观看| 亚洲一区高清亚洲精品| 男女那种视频在线观看| 可以在线观看毛片的网站| 欧美最新免费一区二区三区| 26uuu在线亚洲综合色| 欧美激情在线99| 久久人人爽人人片av| 伊人久久国产一区二区| 成年女人在线观看亚洲视频 | 日本三级黄在线观看| 国产精品人妻久久久久久| 国产成人精品福利久久| 两个人的视频大全免费| 久久国产乱子免费精品| a级毛片免费高清观看在线播放| 毛片一级片免费看久久久久| 韩国高清视频一区二区三区| 日韩一区二区视频免费看| 国产av码专区亚洲av| 2021天堂中文幕一二区在线观| 日韩欧美三级三区| 啦啦啦中文免费视频观看日本| 熟女人妻精品中文字幕| 成年av动漫网址| 欧美变态另类bdsm刘玥| 插逼视频在线观看| 在线天堂最新版资源| 亚洲精品影视一区二区三区av| 国产成人精品婷婷| av一本久久久久| 视频中文字幕在线观看| 欧美日韩国产mv在线观看视频 | 久久人人爽人人片av| 校园人妻丝袜中文字幕| 欧美成人一区二区免费高清观看| 成人无遮挡网站| 麻豆精品久久久久久蜜桃| 久久久久久久大尺度免费视频| 亚洲国产成人一精品久久久| 久久久久久久亚洲中文字幕| 尾随美女入室| 欧美激情久久久久久爽电影| 日韩欧美 国产精品| 国产三级在线视频| 亚洲国产精品专区欧美| 我要看日韩黄色一级片| 一级毛片电影观看| a级毛片免费高清观看在线播放| 亚洲av福利一区| 嫩草影院精品99| 男插女下体视频免费在线播放| 亚洲精品中文字幕在线视频 | 淫秽高清视频在线观看| 在线 av 中文字幕| 又黄又爽又刺激的免费视频.| 白带黄色成豆腐渣| 日韩中字成人| 亚洲精品久久午夜乱码| 国产人妻一区二区三区在| 青春草国产在线视频| 欧美一级a爱片免费观看看| 日韩,欧美,国产一区二区三区| 天堂av国产一区二区熟女人妻| 99久久九九国产精品国产免费| 亚洲欧美精品专区久久| 成人特级av手机在线观看| 亚洲综合精品二区| 免费在线观看成人毛片| 大又大粗又爽又黄少妇毛片口| 99久久九九国产精品国产免费| 亚洲精品日韩av片在线观看| 国产精品女同一区二区软件| 在现免费观看毛片| 国产一级毛片在线| 国产高清三级在线| 亚洲精品456在线播放app| 秋霞伦理黄片| 亚洲av日韩在线播放| 成人毛片a级毛片在线播放| 国产一区亚洲一区在线观看| 亚洲欧洲日产国产| 高清欧美精品videossex| 久久精品夜色国产| 亚洲国产精品国产精品| 最近最新中文字幕大全电影3| 亚洲精华国产精华液的使用体验| 网址你懂的国产日韩在线| 亚洲电影在线观看av| 亚洲国产精品sss在线观看| 在线免费观看的www视频| 少妇猛男粗大的猛烈进出视频 | 国产成人aa在线观看| 成人欧美大片| 日韩av在线免费看完整版不卡| 亚洲欧洲日产国产| 免费在线观看成人毛片| 久久久久久久久久成人| 3wmmmm亚洲av在线观看| 观看美女的网站| 蜜桃亚洲精品一区二区三区| av专区在线播放| 久久99蜜桃精品久久| www.色视频.com| 精品国产三级普通话版| 最近中文字幕高清免费大全6| 一区二区三区免费毛片| 精品久久久久久成人av| 日韩一区二区视频免费看| 99热这里只有是精品50| 大又大粗又爽又黄少妇毛片口| 伦理电影大哥的女人| 国产 一区精品| 国产亚洲最大av| 成人性生交大片免费视频hd| 亚洲天堂国产精品一区在线| 一个人看的www免费观看视频| 国产 一区精品| 亚洲av一区综合| 色视频www国产| 日日摸夜夜添夜夜添av毛片| 国产高潮美女av| 99久久精品热视频| 最后的刺客免费高清国语| 极品少妇高潮喷水抽搐| 少妇高潮的动态图| 国产黄频视频在线观看| 久久精品综合一区二区三区| 午夜激情欧美在线| 亚洲欧美成人精品一区二区| 国产 一区精品| videossex国产| 中文字幕久久专区| 亚洲激情五月婷婷啪啪| 丰满乱子伦码专区| 三级毛片av免费| 国产一区二区三区综合在线观看 | 免费av不卡在线播放| 亚洲怡红院男人天堂| 插阴视频在线观看视频| 中文字幕人妻熟人妻熟丝袜美| 日本熟妇午夜| 国产精品久久久久久精品电影| 久久国内精品自在自线图片| 一级毛片aaaaaa免费看小| 天天躁夜夜躁狠狠久久av| 青春草亚洲视频在线观看| 国产精品久久久久久久久免| 亚洲成色77777| 国产精品人妻久久久影院| a级毛色黄片| 久久鲁丝午夜福利片| 久久这里只有精品中国| 成年人午夜在线观看视频 | 日本-黄色视频高清免费观看| 国产成人91sexporn| 国产黄片视频在线免费观看| 人妻系列 视频| 午夜免费观看性视频| 寂寞人妻少妇视频99o| 草草在线视频免费看| 久久97久久精品| 美女高潮的动态| 久久久欧美国产精品| 日韩欧美 国产精品| 一边亲一边摸免费视频| 禁无遮挡网站| 亚洲精品一区蜜桃| 你懂的网址亚洲精品在线观看| av一本久久久久| 日本三级黄在线观看| 三级国产精品欧美在线观看| 激情五月婷婷亚洲| 蜜桃亚洲精品一区二区三区| 一级毛片黄色毛片免费观看视频| 91久久精品国产一区二区三区| 成人二区视频| 精品国产露脸久久av麻豆 | 国内精品宾馆在线| 在线播放无遮挡| 看非洲黑人一级黄片| 日韩亚洲欧美综合| 亚洲av中文字字幕乱码综合| 两个人的视频大全免费| 菩萨蛮人人尽说江南好唐韦庄| 少妇裸体淫交视频免费看高清| 午夜激情欧美在线| 午夜免费男女啪啪视频观看| 亚洲自拍偷在线| 嫩草影院入口| 日韩视频在线欧美| 免费高清在线观看视频在线观看| 视频中文字幕在线观看| av线在线观看网站| 寂寞人妻少妇视频99o| 亚洲av电影在线观看一区二区三区 | 嘟嘟电影网在线观看| 看免费成人av毛片| 日韩不卡一区二区三区视频在线| 日韩欧美精品免费久久| 黄色欧美视频在线观看| 国产一区二区三区av在线| 久久99热这里只频精品6学生| 直男gayav资源| 人人妻人人澡人人爽人人夜夜 | 色综合站精品国产| 亚洲图色成人| 亚洲精品亚洲一区二区| 一区二区三区乱码不卡18| 亚洲,欧美,日韩| 亚洲av福利一区| 亚洲精品国产av蜜桃| 天堂中文最新版在线下载 | 91午夜精品亚洲一区二区三区| 日日干狠狠操夜夜爽| 国产在线一区二区三区精| 日本免费在线观看一区| 一级二级三级毛片免费看| 国产单亲对白刺激| 亚洲熟妇中文字幕五十中出| 欧美日韩国产mv在线观看视频 | 国产成人福利小说| 免费黄网站久久成人精品| 男女啪啪激烈高潮av片| 一个人免费在线观看电影| 国产一区有黄有色的免费视频 | 亚洲国产精品专区欧美| 高清欧美精品videossex| 成人综合一区亚洲| av播播在线观看一区| 成人午夜高清在线视频| 亚洲av成人精品一区久久| 人妻夜夜爽99麻豆av| 91久久精品国产一区二区三区| 国产精品国产三级国产专区5o| 伊人久久精品亚洲午夜| 久久精品夜夜夜夜夜久久蜜豆| 大香蕉97超碰在线| 日韩,欧美,国产一区二区三区| 亚洲熟妇中文字幕五十中出| 欧美高清性xxxxhd video| 久久久午夜欧美精品| 男插女下体视频免费在线播放| 成人高潮视频无遮挡免费网站| 精品久久国产蜜桃| 精品人妻一区二区三区麻豆| 欧美日韩在线观看h| 国产人妻一区二区三区在| 欧美日韩综合久久久久久| 观看美女的网站| 国产在线男女| 美女黄网站色视频| 亚洲精品456在线播放app| 日韩一本色道免费dvd| 一区二区三区免费毛片| 97超碰精品成人国产| 国产精品久久久久久精品电影| 国产精品久久久久久精品电影小说 | 最近中文字幕2019免费版| a级毛片免费高清观看在线播放| 亚洲久久久久久中文字幕| 久久久午夜欧美精品| 日韩欧美精品免费久久| 搞女人的毛片| 黑人高潮一二区| 国产精品蜜桃在线观看| 亚洲av中文av极速乱| 国产亚洲精品久久久com| www.av在线官网国产| 丰满乱子伦码专区| 男女那种视频在线观看| 26uuu在线亚洲综合色| 亚洲成人久久爱视频| 夜夜看夜夜爽夜夜摸| 亚洲精品久久午夜乱码| 可以在线观看毛片的网站| 久久这里只有精品中国| 日韩一区二区三区影片| 永久免费av网站大全| 欧美日韩国产mv在线观看视频 | 国产黄a三级三级三级人| 国产精品蜜桃在线观看| 2021天堂中文幕一二区在线观| 人妻少妇偷人精品九色| 成年女人看的毛片在线观看| 97热精品久久久久久| 国产成人精品福利久久| 一个人免费在线观看电影| 亚洲欧美成人精品一区二区| 少妇的逼好多水| 精品久久久噜噜| 日韩欧美三级三区| 精品亚洲乱码少妇综合久久| 91av网一区二区| 非洲黑人性xxxx精品又粗又长| 国产成人精品久久久久久| 天天躁日日操中文字幕| 亚洲人成网站在线观看播放| 亚洲一级一片aⅴ在线观看| 熟妇人妻不卡中文字幕| 深夜a级毛片| 国产精品一二三区在线看| 男人舔奶头视频| 亚洲国产色片| 国产成人精品久久久久久| av国产久精品久网站免费入址| 亚洲欧美日韩东京热| 一级二级三级毛片免费看| 日韩在线高清观看一区二区三区| 亚洲国产欧美人成| 精品久久久噜噜| 国产欧美另类精品又又久久亚洲欧美| 精品久久久久久久久久久久久| 国产精品av视频在线免费观看| 国产精品久久久久久久久免| 日本与韩国留学比较| 青春草国产在线视频| 免费高清在线观看视频在线观看| 国产黄色视频一区二区在线观看| 亚洲欧美一区二区三区黑人 | 国产精品综合久久久久久久免费| 午夜老司机福利剧场| 精品少妇黑人巨大在线播放| 一级a做视频免费观看| 欧美一级a爱片免费观看看| 最新中文字幕久久久久| 久久精品夜色国产| 特大巨黑吊av在线直播| 精品人妻偷拍中文字幕| 插阴视频在线观看视频| 亚洲经典国产精华液单| 最近视频中文字幕2019在线8| 免费观看精品视频网站| 久久久久久伊人网av| 午夜福利成人在线免费观看| 五月伊人婷婷丁香| 久久草成人影院| 国产真实伦视频高清在线观看| 最近的中文字幕免费完整| 91久久精品电影网| 久久久色成人| 身体一侧抽搐| 国产精品久久久久久精品电影小说 | 国产成人免费观看mmmm| 亚洲图色成人| 日日撸夜夜添| 国产高清三级在线| 国产男人的电影天堂91| 国产探花在线观看一区二区| 亚洲内射少妇av| 国产午夜精品一二区理论片| 狂野欧美激情性xxxx在线观看| 18+在线观看网站| 欧美性猛交╳xxx乱大交人| 尾随美女入室| 国产黄色小视频在线观看| 精品少妇黑人巨大在线播放| 搡老妇女老女人老熟妇| kizo精华| 国产黄色小视频在线观看| 欧美激情国产日韩精品一区| 99九九线精品视频在线观看视频| 看黄色毛片网站| 乱人视频在线观看| 波多野结衣巨乳人妻| a级毛片免费高清观看在线播放| 亚洲人成网站在线观看播放| 舔av片在线| 国产一区二区三区综合在线观看 | 成年人午夜在线观看视频 | 亚洲自偷自拍三级| 国产淫语在线视频| 男女国产视频网站| 91狼人影院| 一区二区三区四区激情视频| 成人亚洲精品一区在线观看 | 亚洲天堂国产精品一区在线| 亚洲综合色惰| 久久久a久久爽久久v久久| 五月玫瑰六月丁香| 日韩电影二区| 色综合站精品国产| 亚洲性久久影院| 最近中文字幕高清免费大全6| 久久久久国产网址|