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

    基于機載LiDAR的點云數(shù)據(jù)濾波方法

    2022-04-29 06:17:36黃代軍
    測繪標準化 2022年1期
    關(guān)鍵詞:偏度曲面高程

    黃代軍

    (浙江省測繪科學技術(shù)研究院 浙江杭州 310030)

    機載LiDAR技術(shù)經(jīng)過20多年的發(fā)展,其硬件設(shè)備已經(jīng)取得了長足的進步。作為主動式數(shù)據(jù)采集手段,機載LiDAR技術(shù)已經(jīng)廣泛應(yīng)用在場景數(shù)據(jù)采集中,如建筑物數(shù)據(jù)采集、電線巡檢數(shù)據(jù)采集和中小型區(qū)域地形植被數(shù)據(jù)采集等。相比于傳統(tǒng)的航空攝影測量,機載LiDAR不僅抗干擾能力強,而且采集的點云數(shù)據(jù)信息豐富,可以全天候作業(yè)。點云數(shù)據(jù)的優(yōu)點在于可以真實反映場景中地物的細節(jié)信息,如強度信息和位置信息等,但也存在一些缺點,如采樣粒度不均勻、冗余高和數(shù)據(jù)量大等。

    點云數(shù)據(jù)的后處理問題一直是測繪行業(yè)研究的難題之一,相比于點云采集的硬件設(shè)備,點云數(shù)據(jù)的處理技術(shù)仍然相對滯后,其中點云數(shù)據(jù)的濾波更是點云數(shù)據(jù)處理中的關(guān)鍵點和難點之一[1]。經(jīng)過幾十年的發(fā)展,許多專家學者提出并驗證了多種點云濾波算法。目前,常用的點云數(shù)據(jù)濾波算法主要包括基于不規(guī)則三角網(wǎng)濾波算法、移動窗口濾波算法、數(shù)學形態(tài)學濾波算法和移動曲面擬合濾波算法等。本文在傳統(tǒng)的移動曲面擬合濾波算法的基礎(chǔ)上,引入偏度概念,對原始算法進行了改進,并采用國際攝影測量與遙感學會提供的數(shù)據(jù),通過試驗驗證改進算法的可行性,同時與原始算法進行對比,判斷兩種算法的效果。

    1 常用的機載LiDAR點云數(shù)據(jù)濾波方法

    1.1 基于不規(guī)則三角網(wǎng)濾波算法

    Terrasolid軟件使用基于不規(guī)則三角網(wǎng)濾波算法?;诓灰?guī)則三角網(wǎng)濾波算法的機理:首先,對所有數(shù)據(jù)進行計算得到初始參數(shù),再通過選取地面種子點構(gòu)成不規(guī)則三角網(wǎng);其次,設(shè)置條件,通過向不規(guī)則三角網(wǎng)中添加滿足條件的地面點來加密不規(guī)則三角網(wǎng);最后,當?shù)孛鏇]有滿足加入三角網(wǎng)的點集時,結(jié)束運算。

    基于不規(guī)則三角網(wǎng)濾波算法僅能計算一些連續(xù)表面環(huán)境的點云,具有一定的局限性,對于城市中的一些不是連續(xù)表面中的點會因為不能滿足算法設(shè)置的條件而被遺漏。為了解決這一問題,李卉等[2]使用鏡像的方法將遺漏的點加入到不規(guī)則三角網(wǎng)中。此方法雖有一定的可行性,但是理論依據(jù)不足。為此,他們在算法中融入了一種融合區(qū)域增長思想,大大減少了漏選點。

    1.2 移動窗口濾波算法

    移動窗口濾波算法的機理:首先,在某個范圍內(nèi)取一個比較大的窗口,找到最低點;其次,設(shè)置一定的閾值,如果最低點與窗口中最高點的高差不大于此閾值,那么可將此點視為地面點。通過不斷減少窗口大小和閾值,尋找到更為精確的地面點。同時,對利用不同窗口大小得到的地面點進行配權(quán),以加權(quán)平均的方式得到DEM成果。

    移動窗口濾液算法在使用時,如果設(shè)置的窗口太大,可能導(dǎo)致表面不連續(xù)的點被濾掉,窗口太小有可能會將建筑物等識別為地面點。設(shè)置的閾值過小,可能會將表面不連續(xù)的點識別為非地面點;設(shè)置的閾值過大,可能會將植被等點識別為地面點。所以使用移動窗口濾波算法的關(guān)鍵在于參數(shù)的選取是否合適。

    1.3 數(shù)學形態(tài)學濾波算法

    數(shù)學形態(tài)學濾波算法在數(shù)字圖像處理中的應(yīng)用相當廣泛,其機理是通過具有一定形態(tài)的結(jié)構(gòu)實現(xiàn)圖像中對應(yīng)形狀的提取,達到圖像識別和圖像分析的目的。

    數(shù)字形態(tài)學濾波算法包括兩個不同的運算流程:開運算與閉運算,通過膨脹與腐蝕兩個基本運算實現(xiàn)。先進行腐蝕操作,再進行膨脹操作的是開運算;先進行膨脹操作,再進行腐蝕操作的是閉運算。開運算可以實現(xiàn)對機載LiDAR點云的數(shù)據(jù)濾波。

    數(shù)學形態(tài)學濾波算法受窗口尺寸選取影響較大,在實際使用前,需要對點云數(shù)據(jù)所在環(huán)境具有一定的先驗知識,才能設(shè)置合適的窗口尺寸,從而實現(xiàn)對機載LiDAR點云數(shù)據(jù)的濾波[3]。

    1.4 移動曲面擬合濾波算法

    移動曲面擬合濾波算法可以實現(xiàn)對數(shù)據(jù)的濾波,計算速度快且簡單明了。

    移動曲面擬合濾波算法的機理:首先,在種子區(qū)域內(nèi)找到最低的3個角點作為初始地面點,對其進行二維排序,通過求解平面方程的方式擬合初始平面,在平面方程中加入相鄰角點的平面坐標,得到擬合高程值;其次,比較擬合高程與觀測高程的大小,如果二者之間的差值在給定的閾值范圍內(nèi),則此相鄰角點為地面點,否則為非地面點;再次,將判斷為地面點的相鄰角點與組成初始平面的3個角點擬合為1個新的曲面;最后,按照上述方法繼續(xù)判斷相鄰角點是否為地面點,當判斷得到的地面點數(shù)量達到6個時,在增加一個新地面點的同時去掉一個最遠的角點。

    重復(fù)上述操作,通過簡單的地形曲面擬合方式完成對整個測區(qū)的過濾[4]。

    2 基于偏度平衡的移動曲面擬合濾波算法

    本文在移動曲面擬合濾波算法的基礎(chǔ)上,引入統(tǒng)計學中的偏度概念。偏度(skewness),也稱偏度系數(shù),是對變量的概率分布密度函數(shù)曲線偏離對稱度的統(tǒng)計,用以表現(xiàn)概率密度曲線尾部的相對長度。

    假設(shè)有一隨機變量X,同時擁有三階和四階矩陣,那么偏度Sk可表示為[5]

    (1)

    標準正態(tài)分布的偏度為0,當Sk<0時稱為負偏態(tài),Sk>0時稱為正偏態(tài);Sk的絕對值越大,表明其偏離程度越大。本文引入偏度概念的目的在于剔除干擾點。

    試驗前,假設(shè)數(shù)據(jù)中的地面點概率密度分布函數(shù)服從標準正態(tài)分布,非地面點呈偏態(tài)分布,對數(shù)據(jù)進行干擾,即通過不斷地從點云中剔除非地面點的方式使得偏態(tài)分布更接近于正態(tài)分布。同時,為了使樣本的分布逐漸趨于平衡,還需要不斷地對偏度值進行校正。

    移動曲面擬合濾波算法實現(xiàn)原理較為簡單,但對濾波數(shù)據(jù)的要求較高,只有在點云密度足夠高的前提下才能得到較好的效果。該算法通過判斷相鄰角點是否為地面點的方式不斷擬合新的曲面,同時不斷加入新點,擬合地形不斷變化,誤差不斷累積與傳遞就會導(dǎo)致濾波精度降低。

    另外,該算法的閾值不可以隨著地形起伏的變化而自動調(diào)整,導(dǎo)致對坡度較大的地形特征點不能進行有效分離;其次,該算法對場景中的所有點都是單獨進行比較,會大大降低濾波效率。本文主要是對移動曲面擬合濾波算法中的誤差累積進行改進,改進后的算法流程為[6]:

    1)對機載LiDAR點云數(shù)據(jù)進行預(yù)處理后,對點云區(qū)域進行虛擬網(wǎng)格劃分,并為不同網(wǎng)格中的激光點建立高程序列索引。

    2)確定初始窗口尺寸,在相鄰3×3窗口中選取6個高程最低點,構(gòu)成點集合;然后采用偏度校正的方式剔除虛擬網(wǎng)格中高程最低的非地面點,剩余的點組成種子點集。

    3)使用最小二乘法通過種子點擬合點云數(shù)據(jù)區(qū)塊。

    4)遍歷場景中所有點云,計算遍歷的激光角點的擬合高程與實際高程的差值Δhi,設(shè)置一定的閾值Δh。如果某個點的Δhi在閾值外,那么可認定此點為非地面點,反之則為地面點,同時初步標記地面點與非地面點。設(shè)置合適的閾值對點云濾波非常重要,直接影響濾波精度。本文閾值設(shè)定的方式為

    Threshold=(Hmax-Hmin)×Scale

    (2)

    式中:Threshold為高差閾值,Hmax和Hmin分別為區(qū)域內(nèi)點云的最大高程與最小高程,Scale為窗口尺寸。

    5)擴大窗口尺寸,重復(fù)步驟3)和步驟4),得到分類更加精細的激光角點。當窗口尺寸達到閾值時,結(jié)束算法運行,實現(xiàn)濾波。

    基于偏度平衡的移動曲面擬合濾波算法的具體流程如圖1所示。

    圖1 基于偏度平衡的移動曲面擬合濾波算法流程Fig.1 Flow Chart of Moving Surface Fitting Filtering Algorithm Based on Skewness Balance

    3 試驗驗證與分析

    3.1 數(shù)據(jù)來源

    為了驗證基于偏度平衡的移動曲面擬合濾波算法的可行性,本文采用國際攝影測量與遙感學會公布的CSite2首次回波數(shù)據(jù)及利用不同濾波方法得到的參考濾波數(shù)據(jù)作為試驗數(shù)據(jù)對改進后算法進行驗證與分析。

    國際攝影測量與遙感學會提供的數(shù)據(jù)包含激光脈沖首次、末次的反射強度和反射位置,其中,Csite2的點密度可達到0.67點/平方米。

    3.2 誤差指標

    濾波后的數(shù)據(jù)除了需要定性分析外,還需要使用誤差指標對濾波結(jié)果進行定量分析。目前常用的誤差指標有兩種:

    1)激光點的誤分率。誤分率可分為總誤差、一類誤差和二類誤差。其中,一類誤差主要是在進行地面點分類時剔除部分地面點產(chǎn)生的誤差;二類誤差是在不能有效剔除地物點時產(chǎn)生的誤差;總誤差則是一類誤差與二類誤差通過加權(quán)求和得到的誤差。

    改進后算法的適用性可利用一類誤差和二類誤差進行反映,通常是選擇有效的方法降低二類誤差,以得到更加全面的地形數(shù)據(jù)[7]。本文主要使用激光點的誤分率(表1)對改進算法的效果進行驗證。

    表1 濾波誤分率Tab.1 Filtering Error Rates參考結(jié)果濾波結(jié)果地面點個數(shù)地物點個數(shù)參考數(shù)據(jù)點個數(shù)誤分率/%地面點個數(shù)abe=a+bp= b / e地物點個數(shù)cdf=c+dq= c / f濾波后點數(shù)g=a+ch=b+dm=a+b+c+dr=(b+c)/m

    表1中,a和d分別為在濾波過程中正確分類的地面點和地物點的個數(shù),b為將地面點錯分為地物點的個數(shù),c為將非地面點錯分為地面點的個數(shù),e為正確地面點的個數(shù),f為正確地物點個數(shù),m為場景中所有的LiDAR點個數(shù),p為一類誤差,q為二類誤差,r為總誤差。

    2)在不同區(qū)域選取離散樣本點,計算濾波后地面模型的高程值與原始數(shù)據(jù)高程值的均方根誤差。將濾波前后高程值的均方根誤差和高程差值作為評價指標對改進算法的濾波效果進行評價,均方根誤差越小,高程差值越小,表明濾波效果越好,精度越高[8]。

    3.3 試驗過程

    為了更加直觀地反映基于偏度平衡的移動曲面擬合濾波算法的濾波效果,使用CSite2中的樣本數(shù)據(jù)Samp31和Samp23對改進算法進行驗證。由于樣本數(shù)據(jù)Samp31和Samp23中既有原始數(shù)據(jù)又有濾波參考信息,因此可同時進行定性和定量比對。

    通過MATLAB軟件來運行和實現(xiàn)改進算法,并對傳統(tǒng)算法與改進算法得到的試驗結(jié)果進行比對。在試驗過程中,濾波效果受網(wǎng)格大小的影響較為嚴重,經(jīng)過試驗比較,采用2 m×2 m的網(wǎng)格得到的效果較好。

    在每個網(wǎng)格中選取3個高程最低點組成種子點集,并通過偏度平衡過濾掉網(wǎng)格中高程較高的點。由于網(wǎng)格面積較小,可能會覆蓋到建筑物表面,從而將建筑物點加入到地面點中,在進行平面擬合時由于含有非地面點,導(dǎo)致誤差累積。剔除高程較大值后,得到所有最低點高程,對不同最低點密度進行統(tǒng)計得到種子點直方圖,如圖2所示。

    圖2 種子點高頻直方圖Fig.2 High Frequency Histogram of Seed Points

    通過種子點進行曲面擬合,以Samp31數(shù)據(jù)為例得到接近于地面點的擬合曲面如圖3所示。

    圖3 擬合曲面Fig.3 Fitting Surface

    使用改進濾波算法與基于不規(guī)則三角網(wǎng)濾波算法對Samp31和Samp23試驗數(shù)據(jù)進行濾波,結(jié)果如圖4所示。

    圖4 不同濾波方法試驗結(jié)果Fig.4 Experimental Results of Different Filtering Methods

    圖4(a)是采用TerraScan軟件的算法對Samp31數(shù)據(jù)進行濾波處理得到的地面點集,可以看出其算法基本分離出了地面點,視覺上可以明顯看到地物點和地面點的區(qū)別,但實際是地面點中仍然含有少量地物點。

    圖4(b)是采用基于偏度平衡的移動曲面擬合濾波算法對Samp31數(shù)據(jù)進行濾波處理得到的地面點集,可以明顯看到空白區(qū)域為建筑物所在區(qū)域,過濾掉存在的小型植被區(qū)域點云。圖4(c)是采用Terrasolid軟件的算法對Samp23數(shù)據(jù)進行濾波處理得到的地面點集,圖4(d)是采用基于偏度平衡的移動曲面擬合濾波算法對Samp23數(shù)據(jù)進行濾波處理得到的地面點集,可以看到提取效果明顯好于對Samp31數(shù)據(jù)進行濾波處理的效果,但是仍存在一類誤差和二類誤差。其中,圖4(d)中的點云更能反映出地面區(qū)域特征,更明確地表達出建筑物及植被的信息。

    3.4 結(jié)果分析

    為了定量分析基于偏度平衡的移動曲面擬合濾波算法的濾波效果,本文使用誤分率這一評價指標對濾波效果進行評價。

    由于點云數(shù)據(jù)量大,采用人工選取數(shù)據(jù)進行濾波效果對比的方式不實際,因此,本文選擇與Csite樣本數(shù)據(jù)中標記的參考濾波數(shù)據(jù)進行對比的方式,分別使用移動曲面擬合濾波法和基于偏度平衡的移動曲面擬合濾波算法對Samp31數(shù)據(jù)和Samp23數(shù)據(jù)進行濾波處理,并定量評價處理效果,結(jié)果如表2所示。

    表2 濾波效果Tab.2 Filtering Effects 參考結(jié)果Samp31數(shù)據(jù)濾波結(jié)果地面點/個地物點/個參考數(shù)據(jù)點/個誤分率/%Samp23數(shù)據(jù)濾波結(jié)果地面點/個地物點/個參考數(shù)據(jù)點/個誤分率/%地面點/個15 23674513 1905.6516 34365713 1904.98地物點/個1 19912 18816 1057.441 05311 28416 1056.54濾波后點/個16 43512 93329 2956.6417 39611 94129 2955.84

    根據(jù)表2可知,移動曲面擬合濾波法和基于偏度平衡的移動曲面擬合濾波算法的濾波效果都在10%以內(nèi),一類誤差小于二類誤差,二類誤差和總誤差稍偏大,但都在可接受范圍內(nèi),并且基于偏度平衡的移動曲面擬合濾波算法的一類誤差、二類誤差與總誤差均小于移動曲面擬合濾波算法。試驗結(jié)果表明,基于偏度平衡的移動曲面擬合濾波算法適用于城市地區(qū)點云濾波,結(jié)果較為可靠,精度更高。

    4 結(jié) 語

    本文首先介紹了幾種常用的點云濾波方法,在分析移動曲面擬合濾波算法的基礎(chǔ)上,提出了一種基于偏度的移動曲面擬合濾波算法。該算法在確定種子點的過程中加入偏度概念,以降低種子點誤差,達到在上一次濾波完成之后繼續(xù)對非地面點和地面點進行濾波,有效分離非地面點與地面點,從而提高濾波精度的目的。同時,使用國際攝影測量與遙感學會提供的數(shù)據(jù)對改進算法進行驗證,結(jié)果表明,基于偏度的移動曲面擬合濾波算法的一類誤差和二類誤差都在10%以內(nèi),且濾波效果優(yōu)于移動曲面擬合濾波算法,可以減少一類誤差、二類誤差和總誤差?;谄鹊囊苿忧鏀M合濾波算法適用于城市地區(qū)點云濾波,可取得較為可靠的濾波效果和更高的濾波精度。

    猜你喜歡
    偏度曲面高程
    對稱分布的矩刻畫
    8848.86m珠峰新高程
    當代陜西(2020年23期)2021-01-07 09:24:44
    相交移動超曲面的亞純映射的唯一性
    圓環(huán)上的覆蓋曲面不等式及其應(yīng)用
    GPS控制網(wǎng)的高程異常擬合與應(yīng)用
    基于偏度的滾動軸承聲信號故障分析方法
    基于曲面展開的自由曲面網(wǎng)格劃分
    考慮偏度特征的動態(tài)多響應(yīng)穩(wěn)健參數(shù)設(shè)計與優(yōu)化
    基于偏度、峰度特征的BPSK信號盲處理結(jié)果可信性評估
    電子器件(2015年5期)2015-12-29 08:42:56
    SDCORS高程代替等級水準測量的研究
    国产亚洲欧美在线一区二区| 久久av网站| 色播在线永久视频| 亚洲天堂av无毛| 他把我摸到了高潮在线观看 | 91九色精品人成在线观看| 国产三级黄色录像| 国产一区二区 视频在线| a级片在线免费高清观看视频| 午夜久久久在线观看| 日本欧美视频一区| 一级毛片电影观看| 天堂中文最新版在线下载| 亚洲中文字幕日韩| 嫁个100分男人电影在线观看| 欧美精品亚洲一区二区| 国内毛片毛片毛片毛片毛片| 国产精品免费视频内射| 亚洲天堂av无毛| 国产xxxxx性猛交| 亚洲一区中文字幕在线| 久久久精品94久久精品| 岛国毛片在线播放| 国产精品99久久99久久久不卡| 男女边摸边吃奶| 国产av一区二区精品久久| 美国免费a级毛片| 成人手机av| 久久久精品免费免费高清| 久久精品亚洲精品国产色婷小说| 国产极品粉嫩免费观看在线| kizo精华| 久久久久网色| 一本色道久久久久久精品综合| 黑丝袜美女国产一区| 精品国产超薄肉色丝袜足j| 大片电影免费在线观看免费| 国产一区二区激情短视频| 老司机亚洲免费影院| 老司机午夜十八禁免费视频| 99精品欧美一区二区三区四区| 免费日韩欧美在线观看| 欧美亚洲 丝袜 人妻 在线| 一本一本久久a久久精品综合妖精| 国产精品99久久99久久久不卡| 欧美日韩一级在线毛片| 亚洲精品粉嫩美女一区| 久久久久久人人人人人| 50天的宝宝边吃奶边哭怎么回事| 欧美在线一区亚洲| 伊人久久大香线蕉亚洲五| 丰满迷人的少妇在线观看| 亚洲人成电影免费在线| 在线永久观看黄色视频| 国产精品98久久久久久宅男小说| 亚洲伊人久久精品综合| 精品少妇黑人巨大在线播放| 欧美精品亚洲一区二区| 国产精品一区二区在线不卡| 久久久久久久久久久久大奶| 国产野战对白在线观看| 老熟妇乱子伦视频在线观看| 亚洲精品国产色婷婷电影| 中文欧美无线码| 亚洲熟女毛片儿| 99热网站在线观看| 免费观看人在逋| 宅男免费午夜| 日韩欧美一区视频在线观看| 久久久久久久大尺度免费视频| 老司机深夜福利视频在线观看| 欧美另类亚洲清纯唯美| a级毛片在线看网站| 一个人免费在线观看的高清视频| 精品国产一区二区三区四区第35| 两性午夜刺激爽爽歪歪视频在线观看 | 成年人黄色毛片网站| 日本欧美视频一区| 丝袜在线中文字幕| 国产精品久久久久成人av| 久久99热这里只频精品6学生| 国产人伦9x9x在线观看| 久久久久久亚洲精品国产蜜桃av| 亚洲av欧美aⅴ国产| 久久久水蜜桃国产精品网| 亚洲全国av大片| 变态另类成人亚洲欧美熟女 | 国产精品自产拍在线观看55亚洲 | 美女视频免费永久观看网站| 日日摸夜夜添夜夜添小说| 老司机福利观看| 如日韩欧美国产精品一区二区三区| 国产在线视频一区二区| 国产在线观看jvid| 午夜精品国产一区二区电影| 老熟女久久久| 在线观看舔阴道视频| 熟女少妇亚洲综合色aaa.| 亚洲黑人精品在线| 99在线人妻在线中文字幕 | 黄色 视频免费看| 最新美女视频免费是黄的| 国产熟女午夜一区二区三区| 啪啪无遮挡十八禁网站| 国产在线视频一区二区| 国产主播在线观看一区二区| 成人特级黄色片久久久久久久 | 9热在线视频观看99| 亚洲人成电影免费在线| 18禁美女被吸乳视频| 国产欧美亚洲国产| 十八禁网站免费在线| 999久久久国产精品视频| 真人做人爱边吃奶动态| 精品欧美一区二区三区在线| 日韩大码丰满熟妇| www.999成人在线观看| 男人舔女人的私密视频| 久久精品亚洲熟妇少妇任你| 免费在线观看日本一区| 午夜激情久久久久久久| 波多野结衣一区麻豆| av有码第一页| 日韩成人在线观看一区二区三区| 黄色成人免费大全| 日日摸夜夜添夜夜添小说| 欧美另类亚洲清纯唯美| 免费不卡黄色视频| 成年人黄色毛片网站| 亚洲欧美一区二区三区黑人| 欧美黄色片欧美黄色片| 每晚都被弄得嗷嗷叫到高潮| 免费在线观看视频国产中文字幕亚洲| 国产成人av激情在线播放| 中文字幕色久视频| 免费人妻精品一区二区三区视频| 亚洲中文av在线| 国产免费av片在线观看野外av| 亚洲五月色婷婷综合| 中文字幕精品免费在线观看视频| 久久久久久久久免费视频了| 欧美一级毛片孕妇| 日本vs欧美在线观看视频| 日本黄色视频三级网站网址 | 久久久久国产一级毛片高清牌| 中文欧美无线码| 欧美在线一区亚洲| 欧美日韩中文字幕国产精品一区二区三区 | 欧美日本中文国产一区发布| 国产日韩欧美亚洲二区| 侵犯人妻中文字幕一二三四区| 久久精品亚洲精品国产色婷小说| 在线观看免费高清a一片| 久久人人爽av亚洲精品天堂| 视频区图区小说| 一区二区三区精品91| 一边摸一边抽搐一进一小说 | 久久精品国产a三级三级三级| 欧美日韩黄片免| 亚洲av片天天在线观看| 欧美在线黄色| 99国产极品粉嫩在线观看| 欧美乱码精品一区二区三区| 乱人伦中国视频| 精品高清国产在线一区| www.熟女人妻精品国产| 亚洲国产欧美在线一区| 久久香蕉激情| 一级,二级,三级黄色视频| 男女边摸边吃奶| 亚洲精品中文字幕一二三四区 | 午夜久久久在线观看| 男女之事视频高清在线观看| 国产成人精品久久二区二区91| cao死你这个sao货| 欧美激情 高清一区二区三区| 999久久久国产精品视频| 免费一级毛片在线播放高清视频 | 在线观看人妻少妇| 久久久久久免费高清国产稀缺| 在线十欧美十亚洲十日本专区| 99久久99久久久精品蜜桃| 丁香六月欧美| 黑人欧美特级aaaaaa片| 国产高清videossex| 免费久久久久久久精品成人欧美视频| 亚洲成a人片在线一区二区| 老司机午夜十八禁免费视频| 国产黄频视频在线观看| 在线观看免费午夜福利视频| 99香蕉大伊视频| 亚洲天堂av无毛| 在线观看免费日韩欧美大片| 日韩大码丰满熟妇| 天天操日日干夜夜撸| 少妇裸体淫交视频免费看高清 | 一本大道久久a久久精品| 国产精品电影一区二区三区 | 女人精品久久久久毛片| 国产精品美女特级片免费视频播放器 | 国产成人精品无人区| 亚洲av第一区精品v没综合| 国产有黄有色有爽视频| 精品国产乱码久久久久久男人| 亚洲人成77777在线视频| 色94色欧美一区二区| 男女下面插进去视频免费观看| 亚洲av日韩精品久久久久久密| 人人妻人人澡人人爽人人夜夜| 露出奶头的视频| 国产一区二区激情短视频| 亚洲熟妇熟女久久| 一个人免费在线观看的高清视频| 视频区欧美日本亚洲| 国产欧美亚洲国产| 深夜精品福利| 国产xxxxx性猛交| 少妇被粗大的猛进出69影院| 十八禁网站网址无遮挡| 高清视频免费观看一区二区| 一二三四社区在线视频社区8| 五月开心婷婷网| 国产午夜精品久久久久久| 啦啦啦 在线观看视频| 视频区欧美日本亚洲| 免费观看av网站的网址| 欧美在线一区亚洲| 69av精品久久久久久 | 国产成人av激情在线播放| 精品第一国产精品| 久久人人爽av亚洲精品天堂| 久久久久精品国产欧美久久久| 在线播放国产精品三级| 多毛熟女@视频| 日日摸夜夜添夜夜添小说| 精品卡一卡二卡四卡免费| 成人永久免费在线观看视频 | 91精品三级在线观看| 久久精品人人爽人人爽视色| 69av精品久久久久久 | 一本综合久久免费| 欧美人与性动交α欧美软件| 中文亚洲av片在线观看爽 | 亚洲专区字幕在线| 国产精品久久久久久精品电影小说| 一级毛片精品| 色婷婷久久久亚洲欧美| 国产精品99久久99久久久不卡| 欧美 日韩 精品 国产| 亚洲美女黄片视频| 日韩欧美一区二区三区在线观看 | 丰满少妇做爰视频| 精品福利观看| 美女视频免费永久观看网站| 老汉色∧v一级毛片| 久久久久久久国产电影| 啦啦啦免费观看视频1| 亚洲av美国av| 天天操日日干夜夜撸| 一本—道久久a久久精品蜜桃钙片| 欧美日韩亚洲国产一区二区在线观看 | 两性午夜刺激爽爽歪歪视频在线观看 | 久久久国产精品麻豆| 国产精品九九99| 成在线人永久免费视频| 日本a在线网址| 欧美亚洲 丝袜 人妻 在线| 最黄视频免费看| 日韩中文字幕视频在线看片| 国产精品欧美亚洲77777| 久久久精品94久久精品| 成人精品一区二区免费| 香蕉久久夜色| 午夜两性在线视频| 日本五十路高清| 一进一出好大好爽视频| 久久久欧美国产精品| 亚洲中文av在线| 啪啪无遮挡十八禁网站| 国产成人一区二区三区免费视频网站| 色综合婷婷激情| 制服诱惑二区| 无人区码免费观看不卡 | 亚洲精品在线美女| 精品人妻在线不人妻| 99在线人妻在线中文字幕 | 日韩欧美一区视频在线观看| 国产91精品成人一区二区三区 | 午夜成年电影在线免费观看| 亚洲中文日韩欧美视频| 精品久久久精品久久久| 不卡一级毛片| av在线播放免费不卡| 久久久国产精品麻豆| 最近最新中文字幕大全免费视频| 欧美日韩视频精品一区| 国产精品免费大片| 91精品国产国语对白视频| 亚洲午夜精品一区,二区,三区| 在线亚洲精品国产二区图片欧美| 久久狼人影院| 91av网站免费观看| 美女主播在线视频| 国产高清videossex| 欧美人与性动交α欧美软件| 另类精品久久| 99re在线观看精品视频| videosex国产| 亚洲五月婷婷丁香| 自拍欧美九色日韩亚洲蝌蚪91| 久久国产亚洲av麻豆专区| 一个人免费看片子| 国产在线观看jvid| 深夜精品福利| 母亲3免费完整高清在线观看| 男人操女人黄网站| 国产激情久久老熟女| 国产精品久久久久成人av| 久久久久久久国产电影| 久久毛片免费看一区二区三区| 视频区欧美日本亚洲| 亚洲一卡2卡3卡4卡5卡精品中文| 在线播放国产精品三级| 精品少妇内射三级| 黄色 视频免费看| 国产主播在线观看一区二区| 亚洲第一青青草原| cao死你这个sao货| 国产成人欧美在线观看 | 国产97色在线日韩免费| 精品第一国产精品| 99国产精品一区二区三区| 成人国产一区最新在线观看| 99精品欧美一区二区三区四区| 亚洲一码二码三码区别大吗| 69精品国产乱码久久久| 高清av免费在线| 亚洲av成人一区二区三| 亚洲国产欧美日韩在线播放| 18禁国产床啪视频网站| 亚洲国产精品一区二区三区在线| 欧美成狂野欧美在线观看| av网站在线播放免费| 在线观看免费高清a一片| 麻豆乱淫一区二区| 丝袜喷水一区| 国产精品九九99| 亚洲色图 男人天堂 中文字幕| 成人影院久久| 亚洲精品一卡2卡三卡4卡5卡| 99国产极品粉嫩在线观看| 高清视频免费观看一区二区| 久久亚洲真实| 高清av免费在线| 又黄又粗又硬又大视频| 男人舔女人的私密视频| 精品乱码久久久久久99久播| 亚洲七黄色美女视频| 亚洲免费av在线视频| 久久99热这里只频精品6学生| 黄色成人免费大全| 日韩欧美一区二区三区在线观看 | 91九色精品人成在线观看| 另类精品久久| 精品人妻在线不人妻| 亚洲国产精品一区二区三区在线| 精品熟女少妇八av免费久了| 制服人妻中文乱码| www日本在线高清视频| av天堂在线播放| 99在线人妻在线中文字幕 | 国产成人精品在线电影| 久久香蕉激情| 精品国内亚洲2022精品成人 | 一区二区av电影网| 黄片播放在线免费| 精品国产超薄肉色丝袜足j| 日本精品一区二区三区蜜桃| 国产熟女午夜一区二区三区| 精品少妇久久久久久888优播| 欧美日韩黄片免| 天堂动漫精品| 麻豆国产av国片精品| 2018国产大陆天天弄谢| 一进一出好大好爽视频| 一级片'在线观看视频| 精品人妻在线不人妻| www.精华液| 巨乳人妻的诱惑在线观看| 水蜜桃什么品种好| 久久久久久人人人人人| 美女高潮喷水抽搐中文字幕| 一区二区三区精品91| 久久影院123| 亚洲全国av大片| 美女高潮到喷水免费观看| 女警被强在线播放| 精品国产国语对白av| 一个人免费看片子| 国产福利在线免费观看视频| 老司机福利观看| 午夜激情av网站| 国产精品国产av在线观看| 又大又爽又粗| 精品国产乱子伦一区二区三区| 午夜精品国产一区二区电影| 亚洲av日韩精品久久久久久密| 人人妻人人澡人人看| 搡老熟女国产l中国老女人| 少妇精品久久久久久久| 国产福利在线免费观看视频| 国产精品欧美亚洲77777| 精品国产亚洲在线| 亚洲精品美女久久久久99蜜臀| 国产精品美女特级片免费视频播放器 | 黄色成人免费大全| 国产欧美日韩一区二区三| 如日韩欧美国产精品一区二区三区| 真人做人爱边吃奶动态| 丝瓜视频免费看黄片| 黄色视频,在线免费观看| 国产精品98久久久久久宅男小说| 久久中文看片网| 777米奇影视久久| 韩国精品一区二区三区| 亚洲av欧美aⅴ国产| 国产精品免费大片| 久热爱精品视频在线9| 宅男免费午夜| 久9热在线精品视频| 亚洲自偷自拍图片 自拍| 国产高清激情床上av| 久久久国产成人免费| 咕卡用的链子| 日韩欧美一区视频在线观看| 99精品久久久久人妻精品| 国产成+人综合+亚洲专区| 99re在线观看精品视频| 成人黄色视频免费在线看| 中文字幕高清在线视频| 国产不卡一卡二| 国产精品 国内视频| 91老司机精品| 色在线成人网| 超色免费av| 一本综合久久免费| 两个人免费观看高清视频| 日本一区二区免费在线视频| 老熟妇乱子伦视频在线观看| 19禁男女啪啪无遮挡网站| 久久久久久亚洲精品国产蜜桃av| 亚洲欧美日韩另类电影网站| 青草久久国产| 正在播放国产对白刺激| 免费在线观看影片大全网站| 男女之事视频高清在线观看| 淫妇啪啪啪对白视频| 色在线成人网| 国产精品成人在线| 久久精品亚洲精品国产色婷小说| 欧美激情 高清一区二区三区| 国产成人免费观看mmmm| 中文亚洲av片在线观看爽 | 天堂动漫精品| 熟女少妇亚洲综合色aaa.| 欧美日韩黄片免| 一夜夜www| 制服人妻中文乱码| 精品亚洲乱码少妇综合久久| 男女边摸边吃奶| 91老司机精品| 一区福利在线观看| 久久香蕉激情| 又黄又粗又硬又大视频| a级片在线免费高清观看视频| 国产成人免费观看mmmm| 精品国产亚洲在线| 少妇裸体淫交视频免费看高清 | 国产在线免费精品| 色综合欧美亚洲国产小说| 每晚都被弄得嗷嗷叫到高潮| 国产免费福利视频在线观看| 亚洲精品久久午夜乱码| 日韩成人在线观看一区二区三区| 最新的欧美精品一区二区| 国产成人精品久久二区二区免费| 久久毛片免费看一区二区三区| e午夜精品久久久久久久| 丰满饥渴人妻一区二区三| 亚洲av日韩在线播放| 国产成人免费无遮挡视频| 欧美激情 高清一区二区三区| 婷婷丁香在线五月| 一级黄色大片毛片| 91av网站免费观看| 啦啦啦中文免费视频观看日本| 精品视频人人做人人爽| 久久久精品国产亚洲av高清涩受| 亚洲中文字幕日韩| 久久久久久久精品吃奶| 亚洲av成人一区二区三| 十八禁网站免费在线| av一本久久久久| 啦啦啦免费观看视频1| 国产麻豆69| 少妇粗大呻吟视频| a在线观看视频网站| 在线 av 中文字幕| 女人久久www免费人成看片| 91大片在线观看| 啦啦啦免费观看视频1| 免费不卡黄色视频| 每晚都被弄得嗷嗷叫到高潮| 一级,二级,三级黄色视频| 一本久久精品| 精品人妻熟女毛片av久久网站| svipshipincom国产片| 久久久国产欧美日韩av| 超色免费av| 下体分泌物呈黄色| 99精国产麻豆久久婷婷| 黄色毛片三级朝国网站| 久久av网站| 高清在线国产一区| 国产不卡一卡二| 亚洲伊人久久精品综合| 三上悠亚av全集在线观看| 欧美精品一区二区免费开放| 成年版毛片免费区| 久久人妻福利社区极品人妻图片| 91精品国产国语对白视频| 亚洲成av片中文字幕在线观看| 国产精品欧美亚洲77777| 啦啦啦免费观看视频1| 色婷婷av一区二区三区视频| 欧美日韩亚洲高清精品| 十八禁人妻一区二区| 亚洲专区国产一区二区| 国产在线免费精品| 成人国语在线视频| 日韩免费av在线播放| 超色免费av| 丝袜美足系列| 久久av网站| 国产精品一区二区在线不卡| 欧美黑人精品巨大| av又黄又爽大尺度在线免费看| 国产精品影院久久| 国产精品美女特级片免费视频播放器 | 纯流量卡能插随身wifi吗| 狠狠婷婷综合久久久久久88av| 99久久99久久久精品蜜桃| 黄色视频不卡| 脱女人内裤的视频| 99re6热这里在线精品视频| 黄色丝袜av网址大全| 在线观看一区二区三区激情| 国产精品秋霞免费鲁丝片| 亚洲精品在线美女| 亚洲人成电影观看| 亚洲精品一卡2卡三卡4卡5卡| 久久久久久久精品吃奶| 美国免费a级毛片| 老司机影院毛片| 久久国产亚洲av麻豆专区| 捣出白浆h1v1| 亚洲色图综合在线观看| 18禁国产床啪视频网站| 亚洲欧洲精品一区二区精品久久久| 国产精品香港三级国产av潘金莲| 少妇的丰满在线观看| 国产日韩欧美亚洲二区| 亚洲人成伊人成综合网2020| 婷婷成人精品国产| 亚洲五月婷婷丁香| 9热在线视频观看99| 啪啪无遮挡十八禁网站| 曰老女人黄片| 黄频高清免费视频| 美女福利国产在线| 久久久水蜜桃国产精品网| 正在播放国产对白刺激| 嫩草影视91久久| 中文字幕人妻熟女乱码| 欧美 日韩 精品 国产| 欧美成狂野欧美在线观看| 超碰成人久久| 亚洲av第一区精品v没综合| 久久精品亚洲熟妇少妇任你| 免费女性裸体啪啪无遮挡网站| 久久精品亚洲精品国产色婷小说| 香蕉丝袜av| 老司机亚洲免费影院| 亚洲av第一区精品v没综合| 757午夜福利合集在线观看| 纵有疾风起免费观看全集完整版| 国产精品国产高清国产av | 捣出白浆h1v1| 一级毛片女人18水好多| 国产欧美日韩精品亚洲av| 人人妻人人添人人爽欧美一区卜| 精品人妻熟女毛片av久久网站| 久久久久久久久免费视频了| 亚洲男人天堂网一区| 看免费av毛片| 亚洲中文av在线| 欧美日韩国产mv在线观看视频| 叶爱在线成人免费视频播放| 老司机亚洲免费影院| 天天操日日干夜夜撸| 黄色视频在线播放观看不卡| 日韩人妻精品一区2区三区| 中文亚洲av片在线观看爽 | 18禁美女被吸乳视频| 色老头精品视频在线观看|