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

    基于外周彈性約束的礦坑密集點(diǎn)云濾波方法

    2023-10-19 08:35:44熊何喜孫久運(yùn)閆志剛張新耐
    金屬礦山 2023年9期

    熊何喜 孫久運(yùn) 閆志剛 張新耐

    (中國礦業(yè)大學(xué)環(huán)境與測繪學(xué)院,江蘇 徐州 221116)

    高分辨率數(shù)字高程模型可為礦山安全開采、礦區(qū)沉降監(jiān)測、土地修復(fù)、生態(tài)恢復(fù)等提供基礎(chǔ)的數(shù)據(jù)支撐[1-2]。因此,近年來礦區(qū)DEM 提取獲得了業(yè)內(nèi)學(xué)者的廣泛的關(guān)注[3]。礦坑作為開采活動(dòng)的一種特殊遺留形式,由于其地形復(fù)雜、坑壁陡峭、灌叢密布,從DEM 構(gòu)建數(shù)據(jù)采集來看,傳統(tǒng)航空航天采集方式無法有效采集坑壁數(shù)據(jù),人工RTK 采集困難、危險(xiǎn);從數(shù)據(jù)處理來看,目前點(diǎn)云濾波方法效果不甚理想,因而礦坑精細(xì)DEM 構(gòu)建成為一個(gè)難點(diǎn)。

    近年來,快速發(fā)展的傾斜攝影測量由于其寬基線、密集匹配的特點(diǎn)[4-5],密集點(diǎn)云精度、密度能與Li-DAR 點(diǎn)云媲美,數(shù)據(jù)采集更加靈活,為礦坑精細(xì)DEM構(gòu)建提供了新的方向。LiDAR 點(diǎn)云和傾斜攝影測量獲取的密集點(diǎn)云均包含地面點(diǎn)與地物點(diǎn),將兩者分離的濾波工作尤為關(guān)鍵,國內(nèi)外學(xué)者提出了多種濾波算法提取點(diǎn)云中的地面點(diǎn)[6-7]。當(dāng)前濾波算法主要依據(jù)鄰近點(diǎn)云間的高程突變,即局部不連續(xù)進(jìn)行的,處理方法大致分為基于坡度的方法、基于數(shù)學(xué)形態(tài)學(xué)的方法和基于地表的方法。其中,基于坡度的濾波方法,通常假設(shè)地面坡度變化是漸變的,非地面點(diǎn)與地面點(diǎn)之間的變化較大,由此VOSSELMAN 提出通過計(jì)算最大坡度和垂直測量來確定高差閾值進(jìn)行點(diǎn)云濾波,目前基于坡度濾波的研究主要圍繞坡度閾值[8]、結(jié)合地表特征[9]、區(qū)域生長[10]等方面展開?;跀?shù)學(xué)形態(tài)學(xué)的濾波方法,由斯圖加特大學(xué)LINDENBERGER于1993年首次提出,基于膨脹和腐蝕兩項(xiàng)幾何學(xué)基本操作,并將兩個(gè)基礎(chǔ)算子根據(jù)濾波需要進(jìn)行組合,形成開運(yùn)算和閉運(yùn)算,用于不同場景的點(diǎn)云數(shù)據(jù)處理,不少學(xué)者主要圍繞濾波窗口[11]、閾值設(shè)置[12]、結(jié)合語義分割[13]等方面開展研究。前面兩種濾波方法利用坡度或形態(tài)學(xué)算子,可以直接從原始點(diǎn)云中分離地面點(diǎn)和非地面點(diǎn),而基于地表模擬的濾波方法,需要先對(duì)地面進(jìn)行模擬得到一個(gè)近似裸露的地表,再根據(jù)模擬地表分離點(diǎn)云。目前,基于地表模擬濾波方法主要包括三角網(wǎng)濾波法[14]、移動(dòng)曲面濾波法[15]、布料模擬法[16]等。其中,2016年ZHANG 等首次提出布料模擬濾波法,就被廣泛應(yīng)用于各種點(diǎn)云場景,目前主要圍繞結(jié)合地形因子[17]、加密布料[18]等對(duì)布料模擬濾波算法進(jìn)行改進(jìn)。

    1 方法

    1.1 基本思路

    布料模擬濾波使用一塊剛性布料覆蓋在倒置的點(diǎn)云上,通過模擬重力作用下的布料下落過程,考慮布料網(wǎng)格節(jié)點(diǎn)與對(duì)應(yīng)點(diǎn)云之間的相互作用,得到近似地表,從而提取出地面點(diǎn)云,如圖1所示。然而,該方法使用高程最低點(diǎn)作為初始的布料,在地形起伏較大區(qū)域,布料初始化效果差,無法對(duì)地面進(jìn)行有效模擬。

    圖1 布料模擬示意Fig.1 Schematic of cloth simulation filtering

    借鑒布料模擬思路,在下凹地形區(qū)域使用外周彈性約束,用皮筋形變約束過程代替布料重力下落過程,實(shí)現(xiàn)對(duì)礦坑地面的有效模擬。面向下凹式地形,按照一定的高差將點(diǎn)云切割為多個(gè)類似圓形的計(jì)算單元,對(duì)每個(gè)計(jì)算單元使用皮筋進(jìn)行約束,假設(shè)皮筋足夠柔軟,即將點(diǎn)云直接加到皮筋中,得到的皮筋網(wǎng)為礦坑數(shù)字表面模型(Digital Surface Model,DSM);如果將周長略小于類圓計(jì)算單元的彈性皮筋套在外側(cè)進(jìn)行約束,得到的皮筋網(wǎng)為數(shù)字高程模型(Digital Elevation Model,DEM)。為此,本研究提出了下凹地形的外周彈性約束濾波算法,從礦坑密集點(diǎn)云中提取地面點(diǎn),如圖2所示。首先利用密集點(diǎn)云顏色信息和坡度信息對(duì)點(diǎn)云進(jìn)行快速初分類,分類結(jié)果用于整個(gè)皮筋網(wǎng)基本形狀的初始化,使用計(jì)算單元最外側(cè)點(diǎn)云初始化皮筋網(wǎng)格,通過分析每個(gè)皮筋網(wǎng)格單元與對(duì)應(yīng)點(diǎn)云之間的相互作用,加密優(yōu)化皮筋,得到與礦坑地面近似的皮筋網(wǎng),從而提取出原始密集點(diǎn)云中的地面點(diǎn)。

    圖2 彈性約束濾波示意Fig.2 Schematic of the elastic constraint filtering

    1.2 基本原理

    1.2.1 皮筋構(gòu)造

    面向下凹式地形,首先按照一定的高差將點(diǎn)云切割為多個(gè)類似圓形的計(jì)算單元,切割數(shù)目N可根據(jù)下式計(jì)算:

    式中,zmax為礦坑最大高度,m;zmin為最小高度,m;Δh為指定高差,m;rand(·)為取整函數(shù)。

    每個(gè)類圓單元Ci的重心Gi(x,y,z)和平均半徑Ri的計(jì)算公式分別為

    式中,Ti為Ci中的點(diǎn)云總數(shù);Pk為Ci中的三維點(diǎn);d(Gi,Pk)為兩點(diǎn)之間的距離。

    根據(jù)GARCH(1,1)模型的回歸結(jié)果,通過Eviews8.0軟件可以迭代計(jì)算出第t天條件方差的預(yù)測值,然后將開方得到的標(biāo)準(zhǔn)差σt代入VaRt= Pt-1Zασt,這樣根據(jù)第(t-1)天的股票價(jià)格Pt-1可以預(yù)測得到第t天的VaR值。在以上公式中,Zα為置信水平c下的分位數(shù),c=1-α,相應(yīng)的顯著性水平α。為了驗(yàn)證結(jié)果的準(zhǔn)確性,本文計(jì)算了每一天的VaR值(見表5),由于樣本數(shù)據(jù)太多,這里只列出一部分預(yù)測計(jì)算結(jié)果。

    根據(jù)平均半徑Ri,將Ci切割成長度為Δw的單元,得到大小為Δw ×Δh的立面網(wǎng)格,每個(gè)類圓單元中的網(wǎng)格數(shù)量Mi的計(jì)算公式為

    式中,2π 為圓的弧度值;Δw/Ri為每個(gè)立面網(wǎng)格對(duì)應(yīng)的弧度值。此時(shí),可以用立面網(wǎng)格Sij來表示礦坑f(x,y,z),公式為

    立面網(wǎng)格切割過程如圖3所示。考慮到下凹地形特征的立面網(wǎng)格劃分方法,與常規(guī)平面網(wǎng)格劃分相比,更有利于濾波單元的有效計(jì)算。

    圖3 立面網(wǎng)格切割過程Fig.3 Process of cutting vertical grids

    1.2.2 皮筋彈性約束

    基于點(diǎn)云的立面網(wǎng)格劃分,進(jìn)行皮筋彈性約束過程模擬,可以得到與下凹式地形近似的皮筋網(wǎng)。將周長略小于類圓計(jì)算單元的皮筋套在外側(cè),皮筋會(huì)被計(jì)算單元點(diǎn)云撐開產(chǎn)生形變。為了模擬皮筋在每個(gè)立面網(wǎng)格中的形狀,計(jì)算皮筋節(jié)點(diǎn)位置,對(duì)形變過程進(jìn)行簡化,即將每個(gè)網(wǎng)格中的皮筋節(jié)點(diǎn)視為一個(gè)整體分析,并且單個(gè)皮筋網(wǎng)格只受鄰近網(wǎng)格影響,如圖4所示。

    圖4 單個(gè)皮筋單元形變過程Fig.4 Deformation process of a single elastic cell

    在皮筋約束過程中,首先計(jì)算每個(gè)網(wǎng)格中離皮筋中心點(diǎn)最遠(yuǎn)的點(diǎn)云作為初始化皮筋網(wǎng)格,使用初分類點(diǎn)云固定基本皮筋網(wǎng)格,然后考慮鄰近網(wǎng)格影響,更新的皮筋網(wǎng)格位置如圖5所示。

    1.3 彈性約束濾波的實(shí)現(xiàn)

    彈性約束濾波的主要步驟為:首先將點(diǎn)云初進(jìn)行分類,使用分類結(jié)果中高置信度的地面點(diǎn)云初始化皮筋形狀,如果皮筋網(wǎng)格內(nèi)點(diǎn)云數(shù)量滿足要求,則固定該網(wǎng)格;若不滿足要求,則使用其他點(diǎn)云加密優(yōu)化該網(wǎng)格,從而得到與礦坑地面近似的皮筋網(wǎng)。根據(jù)皮筋網(wǎng),計(jì)算原始密集點(diǎn)云到皮筋節(jié)點(diǎn)的距離,如果距離小于指定的閾值,則將其劃分為地面點(diǎn);否則,歸類為非地面點(diǎn)。

    濾波算法實(shí)現(xiàn)過程為:① 預(yù)處理原始點(diǎn)云,去掉重復(fù)點(diǎn)、剔除異常值;② 按照濾波網(wǎng)格分辨率,劃分出立面網(wǎng)格;③ 使用密集點(diǎn)云顏色信息或地形信息提取高地面置信度的點(diǎn)云,用于皮筋網(wǎng)形狀初始化;④ 加密皮筋節(jié)點(diǎn),主要使用遠(yuǎn)離皮筋中心的點(diǎn)加密皮筋節(jié)點(diǎn);⑤ 約束皮筋網(wǎng)格,目標(biāo)網(wǎng)格受到鄰近網(wǎng)格影響,朝著皮筋中心方向發(fā)生位移;⑥ 計(jì)算距離差值,分離點(diǎn)云,根據(jù)得到的皮筋網(wǎng),計(jì)算原始點(diǎn)云到皮筋節(jié)點(diǎn)之間的距離,若小于指定閾值,則將其劃分為地面點(diǎn);否則,劃分為非地面點(diǎn)。

    2 試驗(yàn)與討論

    2.1 試驗(yàn)區(qū)概況及數(shù)據(jù)采集

    本研究試驗(yàn)區(qū)為江蘇省徐州市某廢棄露天礦坑,研究區(qū)域東西長約340 m,南北寬約300 m,礦坑深度約130 m。區(qū)內(nèi)主要有裸露巖體、土地、植被以及道路。數(shù)據(jù)采集設(shè)備為大疆Phantom4 RTK 無人機(jī),搭載數(shù)碼相機(jī)DJI FC6310R,相機(jī)傳感器尺寸為12.833 3 mm,焦距 8.8 mm,拍攝影像大小54 723 648。共采集了礦坑影像877 幅,影像地面分辨率為3.6 cm/pixel,使用PhotoScan 軟件對(duì)齊照片、進(jìn)行密集匹配,生成礦坑密集點(diǎn)云,如圖6所示。

    圖6 礦坑點(diǎn)云圖Fig.6 Point cloud diagram of the pit

    對(duì)原始點(diǎn)云進(jìn)行預(yù)處理,即去掉重復(fù)點(diǎn)和剔除異常值,得到礦坑密集點(diǎn)云共14 751 528 個(gè)點(diǎn),點(diǎn)云密度約155 點(diǎn)/m2。通過目視解譯得到礦坑真實(shí)的地面點(diǎn)云和非地面點(diǎn)云,具體信息見表1。

    表1 參考點(diǎn)云信息Table 1 Information of reference point cloud

    2.2 彈性約束濾波

    礦坑場景中雖然地面起伏較大,但地物類型比較單一,可以利用密集點(diǎn)云的顏色信息和坡度信息快速地對(duì)點(diǎn)云進(jìn)行初分類。通過選取少量典型地面樣本為參考,設(shè)置顏色參數(shù)和坡度參數(shù),能快速得到地面置信度較高的點(diǎn)云,用于皮筋形狀的初始化,再使用立面網(wǎng)格中其他點(diǎn)云加密約束皮筋網(wǎng)格,最終得到與礦坑地面近似的皮筋網(wǎng),從而過濾出地面點(diǎn)云。礦坑皮筋約束濾波細(xì)節(jié)和主要成果如圖7、圖8所示。

    圖7 彈性約束濾波局部細(xì)節(jié)Fig.7 Local details of elastic constraint filtering

    2.3 濾波精度

    為了驗(yàn)證濾波的有效性,選取礦坑中具有不同特征的樣本來測試本研究算法的性能,樣本基本信息見表2。

    表2 樣本描述Table 2 The description of samples

    本研究主要采用準(zhǔn)確率(Accuracy)和ISPRS(International Society for Photogrammetry and Remote Sensing)提出的交叉表格[19]來量化點(diǎn)云濾波結(jié)果和評(píng)價(jià)點(diǎn)云濾波精度。其中,準(zhǔn)確率即為被正確分類的點(diǎn)云數(shù)目與點(diǎn)云總數(shù)的比值,準(zhǔn)確率越高,濾波效果越好。ISPRS 交叉表格包括:點(diǎn)云濾波精度Ⅰ類誤差(地面點(diǎn)云被誤分類為非地面點(diǎn)的比例)、Ⅱ類誤差(非地面點(diǎn)云被誤分類為地面點(diǎn)云的比例)、總誤差(錯(cuò)誤分類點(diǎn)云占比)。經(jīng)過皮筋約束濾波后,計(jì)算樣本的Ⅰ類誤差、Ⅱ類誤差、總誤差和準(zhǔn)確率見表3。由表3 可知:在Ⅰ類誤差方面,皮筋約束濾波在樣本21 中有較好的濾波結(jié)果,能夠有效地提取地面點(diǎn)云;在Ⅱ類誤差方面,地面起伏較小的樣本12 中,皮筋約束能很好地將非地面點(diǎn)云準(zhǔn)確分類;在總誤差和準(zhǔn)確率方面,皮筋約束濾波在樣本21 中表現(xiàn)最佳,其余樣本中表現(xiàn)相當(dāng),說明皮筋約束對(duì)地面起伏較大區(qū)域有著獨(dú)特的優(yōu)勢??傮w來說,綜合考慮露天礦坑密集點(diǎn)云數(shù)據(jù)特點(diǎn)和下凹地形特征,使用彈性約束濾波的方式,濾波精度受地形劇烈起伏干擾較少,但皮筋網(wǎng)格的初始化會(huì)影響濾波效果。在低矮植物覆蓋的地面區(qū)域,地面點(diǎn)云和非地面點(diǎn)云較難用特定的指標(biāo)進(jìn)行劃分,這同樣會(huì)帶來較大誤差。

    表3 樣本精度評(píng)價(jià)結(jié)果Table 3 Evaluation results of sample accuracy %

    本研究選取露天礦坑中常見的地形(樣本12、樣本21),進(jìn)行皮筋約束濾波并計(jì)算樣本的Ⅰ類誤差和Ⅱ類誤差,濾波結(jié)果如圖9 和圖10所示。

    圖9 樣本12 濾波結(jié)果Fig.9 Filtering results of sample 12

    圖10 樣本21 濾波結(jié)果Fig.10 Filtering results of sample 21

    由圖9 和圖10 可知:無論是樣本區(qū)域起伏較小的樣本12,還是地表起伏大的樣本21,皮筋約束都能有效地提取出地面點(diǎn)云,并且能夠成功保留地表細(xì)節(jié)。同時(shí),樣本12、樣本21 的Ⅰ類誤差主要分布在地面與低矮植物相連區(qū)域,一些高程低的植被被誤分類為地面點(diǎn)云;Ⅱ類誤差主要分布地面顏色不明顯,并且有少量植被覆蓋區(qū)域,這主要是因?yàn)檫M(jìn)行皮筋形狀初始化時(shí),該范圍選取的皮筋節(jié)點(diǎn)稀疏,后續(xù)使用皮筋網(wǎng)分離點(diǎn)云時(shí),會(huì)被誤分類為非地面點(diǎn),但這對(duì)構(gòu)建礦坑精細(xì)DEM 影響較小。

    為了定量分析皮筋約束濾波算法的準(zhǔn)確性,使用點(diǎn)云魔方V2.0 軟件中的坡度濾波、數(shù)學(xué)形態(tài)學(xué)濾波與基于地表的濾波(布料模擬、三角網(wǎng)、移動(dòng)曲面)進(jìn)行比較分析,濾波結(jié)果如表4 和表5所示。

    表4 濾波方法在不同地形條件下的總誤差Table 4 Total errors of filtering methods under different terrain conditions%

    表5 濾波方法在礦坑場景中的總誤差和準(zhǔn)確度統(tǒng)計(jì)Table 5 Statistics of total error and accuracy of filtering methods in the pit%

    由表4 和表5 可知:在有植被覆蓋的露天礦坑區(qū)域,使用三角網(wǎng)模擬地面的濾波方法效果不理想,后續(xù)討論略過該方法。在地面起伏較小的區(qū)域(樣本11、樣本12),每種濾波方法處理效果均較理想。在地面起伏較大的區(qū)域(樣本21、樣本22),坡度濾波方法很難實(shí)現(xiàn)坡度閾值的自適應(yīng)調(diào)整,濾波誤差較大;形態(tài)學(xué)濾波方法在植被較少、地面起伏大的區(qū)域樣本21 中,能夠提取到較多的地面種子點(diǎn)進(jìn)行形態(tài)學(xué)操作,誤差較小,而在植被覆蓋較高、地面起伏較大的區(qū)域樣本22 中,誤差較大;通過將剛性布料放置在倒置的點(diǎn)云上,分析布料下落過程和布料節(jié)點(diǎn)之間的關(guān)系來模擬地面的布料模擬濾波算法,在地面起伏較小的區(qū)域效果較理想,但面對(duì)高差達(dá)到百米的露天礦坑區(qū)域,地面模擬效果差,誤差較大;通過擬合函數(shù)進(jìn)行地面模擬的移動(dòng)曲面濾波方法,濾波精度相比于平坦區(qū)域稍有下降,但使用地面擬合函數(shù)很難區(qū)分礦坑中植被與突起巖體之間的差別,導(dǎo)致礦坑中較多地面細(xì)節(jié)丟失,總體濾波精度一般??傮w來看,本研究提出的礦坑外周彈性約束濾波方法,無論是礦坑平坦區(qū)域還是地形起伏較大的坑壁區(qū)域都有不錯(cuò)的表現(xiàn),點(diǎn)云濾波準(zhǔn)確率提高了約10.15%,在起伏較大區(qū)域表現(xiàn)良好。該方法在有植被覆蓋的礦坑區(qū)域,可以有效提取出礦坑地面點(diǎn)云,從而為建立礦坑精細(xì)DEM 提供數(shù)據(jù)基礎(chǔ)。

    2.4 DEM 生成

    使用不規(guī)則三角網(wǎng)對(duì)彈性約束濾波得到的地面點(diǎn)進(jìn)行插值,從鄰近點(diǎn)組成的多個(gè)三角形共同形成的表面提取網(wǎng)格單元高程值,生成精細(xì)礦坑DEM,如圖11所示。生成的礦坑DEM 可以有效展示礦坑整體高程分布情況,并且在地形起伏較大區(qū)域能保留大量地面細(xì)節(jié),在植被覆蓋率高、缺少地面點(diǎn)的區(qū)域,使用距離皮筋節(jié)點(diǎn)近的數(shù)據(jù)進(jìn)行計(jì)算,避免了大面積無地面數(shù)據(jù)區(qū)域的插值,有利于礦坑精細(xì)DEM 數(shù)據(jù)的生成。

    圖11 生成的礦坑DEMFig.11 Generate DEM of the pit

    隨機(jī)抽取地面點(diǎn),計(jì)算插值高程點(diǎn)和真實(shí)高程點(diǎn)之間擬合優(yōu)度R2、均方根誤差RMSE值,從而對(duì)生成的DEM 數(shù)據(jù)進(jìn)行定量評(píng)價(jià),結(jié)果見表6。

    表6 檢查點(diǎn)高程統(tǒng)計(jì)Table 6 Elevation statistics of checkpoints m

    根據(jù)表6 可以求得RMSE為0.048,擬合優(yōu)度R2為0.830。由此可知:礦坑復(fù)雜場景中使用彈性約束濾波可以有效提取地面點(diǎn),通過不規(guī)則三角網(wǎng)插值生成的DEM 與真實(shí)地面擬合較好,能夠?yàn)榈V山提供可靠的高精度DEM 數(shù)據(jù)。

    3 結(jié)語

    面向復(fù)雜的礦坑場景,借鑒布料模擬濾波思路,提出了用水平向內(nèi)的皮筋彈性約束代替布料沿重力方向的自然下落,在豎直方向上將礦坑水平剖分為若干形似“皮筋”的網(wǎng)格,進(jìn)行彈性約束模擬,擬合出礦坑的外輪廓,并剔除內(nèi)部植被點(diǎn)云的濾波方法。利用該方法對(duì)徐州某露天礦坑進(jìn)行了濾波試驗(yàn),結(jié)果表明:該方法與常規(guī)方法相比,準(zhǔn)確度有明顯提高,得到了均方根誤差為0.048 的廢棄礦坑DEM,能夠滿足精細(xì)DEM 構(gòu)建的需求,為礦山高精度礦坑DEM 提取提供了一種思路。

    成人国产麻豆网| 欧美zozozo另类| 69av精品久久久久久| 亚洲五月天丁香| 男插女下体视频免费在线播放| 国产片特级美女逼逼视频| 日韩 亚洲 欧美在线| 三级毛片av免费| 欧美另类亚洲清纯唯美| 又爽又黄a免费视频| 波多野结衣巨乳人妻| 国产成人freesex在线| 日产精品乱码卡一卡2卡三| 99热只有精品国产| 中文精品一卡2卡3卡4更新| 国产精品久久久久久亚洲av鲁大| 欧美成人精品欧美一级黄| 高清日韩中文字幕在线| .国产精品久久| 少妇熟女欧美另类| 波多野结衣高清作品| 亚洲最大成人av| 国产免费一级a男人的天堂| 日本-黄色视频高清免费观看| 欧美+亚洲+日韩+国产| 色噜噜av男人的天堂激情| 国产成年人精品一区二区| 欧美一区二区国产精品久久精品| 天美传媒精品一区二区| 欧美成人a在线观看| 欧美色视频一区免费| 亚洲一区高清亚洲精品| 99久久成人亚洲精品观看| 精品国内亚洲2022精品成人| 国产精品1区2区在线观看.| 看非洲黑人一级黄片| 成人特级av手机在线观看| 久久6这里有精品| 免费观看人在逋| 男的添女的下面高潮视频| 欧美3d第一页| 日韩精品有码人妻一区| 真实男女啪啪啪动态图| 精品人妻熟女av久视频| 青青草视频在线视频观看| 久久久精品欧美日韩精品| 观看美女的网站| 69人妻影院| 性色avwww在线观看| 麻豆国产av国片精品| 秋霞在线观看毛片| 国产精品伦人一区二区| 亚洲精品乱码久久久v下载方式| 一本精品99久久精品77| 久久久久久久亚洲中文字幕| 亚洲精华国产精华液的使用体验 | 日本免费一区二区三区高清不卡| 色吧在线观看| 一区二区三区四区激情视频 | 日韩成人伦理影院| 免费观看精品视频网站| 国模一区二区三区四区视频| 国产亚洲5aaaaa淫片| 国产成人精品一,二区 | 国产色婷婷99| 三级毛片av免费| 最近的中文字幕免费完整| 麻豆av噜噜一区二区三区| 精品99又大又爽又粗少妇毛片| 99热这里只有是精品50| 少妇的逼水好多| 国产探花极品一区二区| 亚洲经典国产精华液单| 午夜福利在线观看免费完整高清在 | 一区二区三区高清视频在线| 亚洲国产精品sss在线观看| 国产蜜桃级精品一区二区三区| 国产精品女同一区二区软件| 校园春色视频在线观看| 午夜福利在线观看免费完整高清在 | 午夜福利视频1000在线观看| 毛片一级片免费看久久久久| 两个人的视频大全免费| 国产日本99.免费观看| 中国国产av一级| 亚洲自拍偷在线| 麻豆久久精品国产亚洲av| 免费av观看视频| 三级经典国产精品| 久久这里只有精品中国| 亚洲欧美日韩高清专用| 亚洲,欧美,日韩| 99国产极品粉嫩在线观看| 久久久精品欧美日韩精品| 成人鲁丝片一二三区免费| 亚洲欧美日韩东京热| 黄片wwwwww| 嫩草影院入口| 欧美潮喷喷水| 深夜a级毛片| 偷拍熟女少妇极品色| 久久久色成人| 哪里可以看免费的av片| 亚洲最大成人手机在线| 最后的刺客免费高清国语| 亚洲成人中文字幕在线播放| av黄色大香蕉| 亚洲五月天丁香| 国产成人91sexporn| 成人三级黄色视频| 国产精品无大码| 久久中文看片网| 亚洲av一区综合| 99riav亚洲国产免费| 男人狂女人下面高潮的视频| 欧美日韩综合久久久久久| 99热只有精品国产| 成年免费大片在线观看| 女同久久另类99精品国产91| 国产成人福利小说| 91精品国产九色| 人人妻人人澡人人爽人人夜夜 | 97超碰精品成人国产| 美女大奶头视频| 欧美成人a在线观看| 久久久a久久爽久久v久久| 精品一区二区三区人妻视频| 亚洲av一区综合| 内地一区二区视频在线| 成人毛片a级毛片在线播放| 亚洲精品国产成人久久av| 国产单亲对白刺激| 亚洲国产高清在线一区二区三| 床上黄色一级片| 亚洲成a人片在线一区二区| 亚洲最大成人中文| 国产一区亚洲一区在线观看| a级毛色黄片| 人妻制服诱惑在线中文字幕| 深夜精品福利| 在线免费观看不下载黄p国产| 精品久久久噜噜| 日日摸夜夜添夜夜添av毛片| 日本熟妇午夜| 欧美+日韩+精品| 精品人妻一区二区三区麻豆| 午夜福利视频1000在线观看| 久久久精品欧美日韩精品| 最近视频中文字幕2019在线8| а√天堂www在线а√下载| 美女内射精品一级片tv| 特大巨黑吊av在线直播| 老司机福利观看| 亚洲国产精品成人久久小说 | 国产亚洲精品久久久久久毛片| 久久国产乱子免费精品| 我要看日韩黄色一级片| 99热这里只有是精品在线观看| 亚洲,欧美,日韩| 少妇人妻一区二区三区视频| 国产男人的电影天堂91| 一个人观看的视频www高清免费观看| 久久国内精品自在自线图片| 99久久人妻综合| .国产精品久久| 久久久久久九九精品二区国产| 国产精品一及| 亚洲va在线va天堂va国产| 尤物成人国产欧美一区二区三区| 国产精品国产高清国产av| 2022亚洲国产成人精品| 麻豆国产97在线/欧美| 日本欧美国产在线视频| 乱码一卡2卡4卡精品| 亚洲色图av天堂| av国产免费在线观看| 亚洲七黄色美女视频| 成人鲁丝片一二三区免费| 赤兔流量卡办理| 亚洲精品亚洲一区二区| 久久婷婷人人爽人人干人人爱| 一个人免费在线观看电影| 一夜夜www| 一进一出抽搐动态| 青春草视频在线免费观看| 欧美日韩一区二区视频在线观看视频在线 | 在线观看美女被高潮喷水网站| 99久久成人亚洲精品观看| 国产乱人视频| 男人的好看免费观看在线视频| 亚洲人成网站在线观看播放| 国产高清激情床上av| 亚洲欧美中文字幕日韩二区| 一级黄色大片毛片| 国产精品精品国产色婷婷| 久久久久久久久久久丰满| 成人毛片a级毛片在线播放| 精品久久久久久久久av| 高清毛片免费观看视频网站| 日本黄色片子视频| 婷婷亚洲欧美| 免费av毛片视频| 国产高清有码在线观看视频| 少妇人妻精品综合一区二区 | 久久99蜜桃精品久久| 丰满乱子伦码专区| 日韩亚洲欧美综合| 一级毛片aaaaaa免费看小| av.在线天堂| 久久人人爽人人爽人人片va| 1024手机看黄色片| 精品久久久久久久久久免费视频| 天堂中文最新版在线下载 | 97超视频在线观看视频| 综合色丁香网| 国产一级毛片在线| 久久久久国产网址| 国产蜜桃级精品一区二区三区| 亚洲五月天丁香| 99久久人妻综合| 女同久久另类99精品国产91| 日日撸夜夜添| 男人的好看免费观看在线视频| 免费人成在线观看视频色| 国产精品久久久久久亚洲av鲁大| 观看美女的网站| 校园春色视频在线观看| h日本视频在线播放| 久久这里只有精品中国| 亚洲婷婷狠狠爱综合网| 色尼玛亚洲综合影院| 全区人妻精品视频| 99在线人妻在线中文字幕| 亚洲人与动物交配视频| 日韩成人伦理影院| 99国产极品粉嫩在线观看| 国产精品女同一区二区软件| 精品日产1卡2卡| 婷婷亚洲欧美| 丰满的人妻完整版| 最近手机中文字幕大全| 亚洲在线观看片| 永久网站在线| 精品久久国产蜜桃| 国产亚洲欧美98| 欧美激情久久久久久爽电影| 国产成人91sexporn| 一本一本综合久久| 亚洲人成网站在线播| 人妻制服诱惑在线中文字幕| 亚洲美女视频黄频| 禁无遮挡网站| 国产激情偷乱视频一区二区| 日本与韩国留学比较| 69人妻影院| 精品久久久久久久久久久久久| 不卡一级毛片| 我要看日韩黄色一级片| 欧美成人免费av一区二区三区| 日本成人三级电影网站| 国产国拍精品亚洲av在线观看| 秋霞在线观看毛片| 久久久久国产网址| 色吧在线观看| 国产黄色视频一区二区在线观看 | 欧美区成人在线视频| 长腿黑丝高跟| 青春草国产在线视频 | 免费电影在线观看免费观看| 日本欧美国产在线视频| 欧美日本视频| 黄片无遮挡物在线观看| 国产精品爽爽va在线观看网站| 国产成人aa在线观看| 亚洲在久久综合| 日日撸夜夜添| 国产中年淑女户外野战色| 黄色一级大片看看| 国产精品人妻久久久影院| 日本一二三区视频观看| 黄色欧美视频在线观看| 日韩人妻高清精品专区| 国产精品女同一区二区软件| 一级av片app| 日本熟妇午夜| 99久久精品热视频| 男人的好看免费观看在线视频| 国产成人a区在线观看| 寂寞人妻少妇视频99o| 亚洲精品色激情综合| 国产亚洲av嫩草精品影院| 国产精品久久电影中文字幕| 国产白丝娇喘喷水9色精品| 在线免费观看的www视频| 最近手机中文字幕大全| 国产欧美日韩精品一区二区| 日产精品乱码卡一卡2卡三| 亚洲婷婷狠狠爱综合网| 大又大粗又爽又黄少妇毛片口| 亚洲真实伦在线观看| 少妇猛男粗大的猛烈进出视频 | 少妇高潮的动态图| 黑人高潮一二区| 熟女人妻精品中文字幕| 特级一级黄色大片| avwww免费| 国产精品嫩草影院av在线观看| 久久亚洲国产成人精品v| 热99在线观看视频| 久99久视频精品免费| 免费看光身美女| 亚洲精品色激情综合| 成年免费大片在线观看| 在线观看免费视频日本深夜| 欧美一区二区国产精品久久精品| 欧美极品一区二区三区四区| 成人特级黄色片久久久久久久| 欧美一区二区国产精品久久精品| 婷婷色av中文字幕| 少妇高潮的动态图| 丰满人妻一区二区三区视频av| 欧美色视频一区免费| 尾随美女入室| 国产成人91sexporn| 国产黄色视频一区二区在线观看 | 亚洲国产高清在线一区二区三| 毛片一级片免费看久久久久| 亚洲av二区三区四区| 亚洲欧美日韩卡通动漫| 亚洲,欧美,日韩| 中文字幕熟女人妻在线| 真实男女啪啪啪动态图| 波多野结衣高清无吗| 免费看美女性在线毛片视频| 校园春色视频在线观看| 欧美一区二区精品小视频在线| 寂寞人妻少妇视频99o| 狠狠狠狠99中文字幕| 天堂网av新在线| 国产三级在线视频| 国产成人影院久久av| 在线播放无遮挡| 久久久久久伊人网av| 欧美色欧美亚洲另类二区| 在线观看66精品国产| 亚洲精品国产成人久久av| 亚洲美女搞黄在线观看| 日韩一本色道免费dvd| 性插视频无遮挡在线免费观看| 亚洲欧美日韩高清在线视频| 国产精品无大码| 国产亚洲av嫩草精品影院| 麻豆成人av视频| 亚洲自偷自拍三级| 中文字幕制服av| 国产爱豆传媒在线观看| 亚洲最大成人av| 成人亚洲精品av一区二区| 欧美精品一区二区大全| 日韩av在线大香蕉| 欧美一级a爱片免费观看看| 欧美日韩在线观看h| 午夜福利视频1000在线观看| 国产成人一区二区在线| 日韩欧美三级三区| 干丝袜人妻中文字幕| 国产一级毛片七仙女欲春2| 天美传媒精品一区二区| 天堂网av新在线| av国产免费在线观看| 成人av在线播放网站| 国产黄片美女视频| 国产视频首页在线观看| 老师上课跳d突然被开到最大视频| 免费看a级黄色片| 99久国产av精品国产电影| 日本熟妇午夜| 色哟哟·www| 又爽又黄a免费视频| 黄片无遮挡物在线观看| 99在线视频只有这里精品首页| 看十八女毛片水多多多| 综合色av麻豆| 国产高清三级在线| 在线观看一区二区三区| 国产人妻一区二区三区在| 天堂中文最新版在线下载 | 亚洲成人久久性| 麻豆国产av国片精品| 99视频精品全部免费 在线| 欧美性感艳星| 国产 一区 欧美 日韩| 日产精品乱码卡一卡2卡三| 国产精华一区二区三区| 少妇熟女aⅴ在线视频| 国产成人精品婷婷| 欧美高清性xxxxhd video| 99国产极品粉嫩在线观看| 99久国产av精品国产电影| 色综合色国产| 国产精品三级大全| 又爽又黄a免费视频| 亚洲真实伦在线观看| 亚洲四区av| 成年版毛片免费区| 亚洲精品国产成人久久av| 真实男女啪啪啪动态图| 黄片wwwwww| 桃色一区二区三区在线观看| 亚洲精品成人久久久久久| 蜜桃久久精品国产亚洲av| 国产成人91sexporn| 又黄又爽又刺激的免费视频.| 国产亚洲91精品色在线| 在线播放无遮挡| 嫩草影院新地址| 99久久九九国产精品国产免费| 久久午夜福利片| 全区人妻精品视频| 国产精品美女特级片免费视频播放器| 久久人人精品亚洲av| 人妻久久中文字幕网| 热99在线观看视频| 国产日韩欧美在线精品| 日韩强制内射视频| 国产精品乱码一区二三区的特点| 亚洲婷婷狠狠爱综合网| 三级国产精品欧美在线观看| 韩国av在线不卡| 91av网一区二区| 欧洲精品卡2卡3卡4卡5卡区| 免费电影在线观看免费观看| 国产伦精品一区二区三区四那| 嫩草影院精品99| 不卡视频在线观看欧美| 午夜老司机福利剧场| 日本黄大片高清| 久久精品夜色国产| 婷婷色综合大香蕉| 久久人妻av系列| 有码 亚洲区| 亚洲无线观看免费| 嫩草影院入口| 99久久久亚洲精品蜜臀av| 久久人人爽人人片av| 亚洲国产精品sss在线观看| 亚洲美女搞黄在线观看| 国产亚洲91精品色在线| 国产伦理片在线播放av一区 | 国产成人午夜福利电影在线观看| 两个人视频免费观看高清| 一级黄色大片毛片| 人妻少妇偷人精品九色| 91av网一区二区| 国产乱人偷精品视频| 亚洲18禁久久av| 国产成人午夜福利电影在线观看| 99国产精品一区二区蜜桃av| 精品熟女少妇av免费看| 一进一出抽搐动态| 3wmmmm亚洲av在线观看| 99热这里只有是精品50| 中国国产av一级| 老熟妇乱子伦视频在线观看| 久久久色成人| 身体一侧抽搐| 一本一本综合久久| 国产极品精品免费视频能看的| 亚洲va在线va天堂va国产| 不卡视频在线观看欧美| 噜噜噜噜噜久久久久久91| 国产欧美日韩精品一区二区| 久久久成人免费电影| 欧洲精品卡2卡3卡4卡5卡区| 色综合色国产| 哪里可以看免费的av片| 国产高清三级在线| 国产私拍福利视频在线观看| 日本三级黄在线观看| 少妇丰满av| 如何舔出高潮| 深夜a级毛片| 99热这里只有是精品在线观看| 亚洲精品日韩在线中文字幕 | 国产黄片美女视频| 日日摸夜夜添夜夜爱| 黄色配什么色好看| 一级毛片久久久久久久久女| 亚洲,欧美,日韩| 舔av片在线| 国产精品人妻久久久久久| 国产精品三级大全| 国产欧美日韩精品一区二区| 国产精品三级大全| 少妇的逼水好多| or卡值多少钱| 久久久久免费精品人妻一区二区| 嫩草影院精品99| 欧美日韩国产亚洲二区| 观看美女的网站| 国产精品乱码一区二三区的特点| 国产黄片视频在线免费观看| 少妇熟女aⅴ在线视频| 国产高清有码在线观看视频| 久久久久国产网址| 欧美在线一区亚洲| 岛国毛片在线播放| 欧美在线一区亚洲| 亚洲av二区三区四区| 一区二区三区免费毛片| 欧美日韩综合久久久久久| 插阴视频在线观看视频| 免费观看在线日韩| 亚洲国产色片| 亚洲精品乱码久久久久久按摩| 久久久久性生活片| 自拍偷自拍亚洲精品老妇| 美女高潮的动态| 18禁黄网站禁片免费观看直播| 亚洲av不卡在线观看| 成人国产麻豆网| 精品人妻一区二区三区麻豆| 国产精品国产三级国产av玫瑰| 久久久精品大字幕| 国产精品三级大全| 国产亚洲精品av在线| 午夜亚洲福利在线播放| 白带黄色成豆腐渣| 最后的刺客免费高清国语| 天堂中文最新版在线下载 | 国产91av在线免费观看| 久久亚洲精品不卡| 日本在线视频免费播放| 亚洲av男天堂| 十八禁国产超污无遮挡网站| 黑人高潮一二区| 国产成人a区在线观看| 国产av不卡久久| 一级毛片aaaaaa免费看小| 一边亲一边摸免费视频| 小蜜桃在线观看免费完整版高清| 舔av片在线| 免费看光身美女| 少妇猛男粗大的猛烈进出视频 | 天美传媒精品一区二区| 精品一区二区三区人妻视频| 午夜福利成人在线免费观看| 国产一区二区激情短视频| 伦理电影大哥的女人| 国产美女午夜福利| 色综合亚洲欧美另类图片| 97超碰精品成人国产| 亚洲va在线va天堂va国产| 国产午夜精品久久久久久一区二区三区| 十八禁国产超污无遮挡网站| 日本爱情动作片www.在线观看| 国产免费一级a男人的天堂| 成熟少妇高潮喷水视频| 久久久久久国产a免费观看| 波野结衣二区三区在线| 久久久久性生活片| 久久精品夜夜夜夜夜久久蜜豆| 中文字幕人妻熟人妻熟丝袜美| 精品午夜福利在线看| 久久这里有精品视频免费| 最近视频中文字幕2019在线8| 色5月婷婷丁香| 欧美在线一区亚洲| 国产熟女欧美一区二区| 午夜老司机福利剧场| 一级黄色大片毛片| 日韩欧美精品免费久久| 免费av不卡在线播放| 亚洲婷婷狠狠爱综合网| av福利片在线观看| 成人无遮挡网站| 久久久成人免费电影| 乱码一卡2卡4卡精品| 国产国拍精品亚洲av在线观看| 人人妻人人看人人澡| 最近中文字幕高清免费大全6| 夫妻性生交免费视频一级片| 国产亚洲精品久久久com| 99久久久亚洲精品蜜臀av| 99热全是精品| 欧美bdsm另类| 日本一本二区三区精品| 三级男女做爰猛烈吃奶摸视频| а√天堂www在线а√下载| 国语自产精品视频在线第100页| 一本久久精品| 1000部很黄的大片| 青春草国产在线视频 | 成年女人永久免费观看视频| 91狼人影院| 男女视频在线观看网站免费| 午夜福利在线观看免费完整高清在 | av黄色大香蕉| 九色成人免费人妻av| 欧美日韩国产亚洲二区| 日日摸夜夜添夜夜爱| 午夜福利在线观看免费完整高清在 | 免费人成在线观看视频色| 免费大片18禁| 一级黄色大片毛片| 国产片特级美女逼逼视频| 国产极品天堂在线| 女的被弄到高潮叫床怎么办| av免费在线看不卡| 日韩欧美在线乱码| 国产成人a∨麻豆精品| 欧美一区二区精品小视频在线| 免费看日本二区| 国产不卡一卡二| 欧美极品一区二区三区四区|