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

    基于地面激光掃描的典型海岸帶鹽沼潮灘地形反演

    2018-12-06 07:48:26周云軒錢偉偉湛玉劍黃蓋先
    關(guān)鍵詞:白茅潮灘鹽沼

    魏 偉,周云軒,田 波,錢偉偉,湛玉劍,黃蓋先

    華東師范大學(xué)河口海岸學(xué)國家重點實驗室,上海 200062

    0 引言

    潮灘是指海岸帶地區(qū)高潮時淹沒、低潮時露出的區(qū)域[1],該區(qū)域常伴隨著耐鹽且喜水的陸生植被即鹽沼植被的出現(xiàn)。鹽沼植被的發(fā)達(dá)根系有助于保持潮灘泥沙,其對潮水的緩流作用和對懸浮細(xì)顆粒泥沙的黏附、捕捉作用也促進(jìn)了泥沙沉降和灘面淤高[2];同時,潮灘高程的增加加快了鹽沼植被群落的演替和向海延伸,植株密度和生物量也隨之增加[2-3]。高精度潮灘地形數(shù)據(jù)可為鹽沼植被對潮灘動力、地貌、沉積過程的影響機(jī)理研究以及鹽沼生物量的獲取提供定量支持,探測和量測潮灘地形變化對于制定海岸帶綜合管理計劃具有重要的現(xiàn)實意義和參考價值。然而潮灘灘面寬廣、潮溝縱橫、灘面質(zhì)地松軟以及鹽沼植被的覆蓋等干擾情況,給鹽沼潮灘地形的精準(zhǔn)快速測量帶來較大難度。

    目前,常用的地形測量方式主要有InSAR(interferometric synthetic aperture radar)、LiDAR(light detection and ranging)以及地面觀測如全站儀、GPS-RTK(global positioning system - real time kinematic)技術(shù)等。InSAR、LiDAR主要用于較大區(qū)域范圍的測量,全站儀、GPS-RTK主要用于小范圍的測量。這些方式方法和手段在應(yīng)用于潮灘時,分別存在著重訪周期長、時效性差、成本高、單點測量效率低下等問題。近年來,TLS(terrestrial laser scanner)在潮灘區(qū)域的地形測量中得到良好的應(yīng)用,它能獲取目標(biāo)物的紋理信息和顏色信息。相對于其他地面測量手段,TLS獲取地形數(shù)據(jù)的效率更高,覆蓋范圍可達(dá)2~5 km,能夠保證野外連續(xù)實時測量。近年來,國內(nèi)外學(xué)者使用TLS在潮灘地形測量和反演方面取得了不少研究成果。常直楊等[4]使用TLS在江蘇潮灘獲取灘涂地形數(shù)據(jù),評估了在該條件下使用的水平和高程偏差,結(jié)果表明該技術(shù)可以高精度獲取潮灘微地貌特征數(shù)據(jù)。Hannam等[5]利用TLS對帕迪拉灣區(qū)域展開連續(xù)3年監(jiān)測,以了解沉積地貌的時空變化模式。Guarnieri等[6]使用TLS獲取意大利威尼斯瀉湖鹽沼區(qū)域的地形數(shù)據(jù),構(gòu)建植被濾除算法,分別得出TLS在較密和較稀疏植被區(qū)域使用的精度。

    TLS在潮灘光灘地區(qū)獲取地形數(shù)據(jù)具有良好的效果,然而在鹽沼潮灘區(qū)域,植被對激光發(fā)射信號有遮蔽效應(yīng),TLS獲取的激光反射信號包含了植被、地面信息等數(shù)據(jù),因此,要獲取真實的鹽沼地形數(shù)據(jù),必須去除植被數(shù)據(jù)的干擾。本項研究以長江口4種典型鹽沼植被(蘆葦、白茅、互花米草、海三棱藨草)的樣方為研究目標(biāo),采用先稀釋植被蓋度后TLS掃描地形數(shù)據(jù)的步驟,構(gòu)建分離植被與地面激光點云的植被濾除算法,研究TLS在鹽沼區(qū)域使用的適用性及典型鹽沼植被與反演地形精度的關(guān)系,為使用TLS反演鹽沼潮灘地形演變量化研究提供參考。

    1 研究區(qū)域概況

    長江口是長江在我國東部的入???,長江上游豐富的泥沙資源在河口的大量淤積使得沿岸潮灘處于動態(tài)的淤漲變化過程中。作為典型的淤泥質(zhì)河口,長江口的灘坡十分平緩,這在一定程度上有利于植被的生長和擴(kuò)散,使得長江口鹽沼植被資源非常豐富。長江口鹽沼植被主要有蘆葦群落、藨草群落、海三棱藨草群落、糙葉苔草群落以及外來物種互花米草群落等[7]。

    本研究區(qū)域選取崇明東灘、南匯邊灘兩塊長江口潮灘鹽沼區(qū)域(圖1a)。崇明東灘位于上海市崇明島東部,南臨長江入??冢瑬|瀕東海,潮溝發(fā)育豐富,走向大致與灘地延伸方向一致,在平面上呈樹狀向外散射[8]。目前該區(qū)域鹽沼植被分布的一般趨勢為:低潮帶的鹽漬藻類帶;中潮帶及高潮帶下部的海三棱藨草和藨草群落帶;高潮帶及潮上帶的蘆葦群落帶。此外還散生一些白茅、糙葉苔草等斑塊狀群落。崇明東灘濕地互花米草生態(tài)控制與鳥類棲息地優(yōu)化工程的順利開展使得大片互花米草種群基本消失。由于不同類型植被的平均株高和密度不同,本文分別在崇明東灘南部地區(qū)選取3塊典型類型的蘆葦區(qū)域、白茅與低矮蘆葦?shù)幕焐鷧^(qū)域(后續(xù)稱白茅區(qū)域)、海三棱藨草區(qū)域建立樣方,分別分布于圖1b中的S1、S2、S3位置,平均株高分別為150、50、30 cm。南匯邊灘地處長江口和杭州灣的交匯處,該區(qū)域潮灘鹽沼植被主要由低潮區(qū)的海三棱藨草(包括糙葉苔草和藨草群落帶),中、高潮混生的互花米草群落及高潮帶內(nèi)側(cè)的蘆葦群落三大類組成[9]。本文在該地區(qū)選取互花米草分布帶建立樣方,分布于圖1c中的S4位置,平均植株高度為150 cm。

    圖1 研究區(qū)及樣方位置示意圖Fig.1 The study area and quadrat location

    2 數(shù)據(jù)獲取與預(yù)處理

    2.1 點云數(shù)據(jù)獲取及植被蓋度稀釋

    本研究使用的地面激光掃描儀型號為RIEGL VZ 4000。該掃描儀的有效掃描距離為4 000 m,單次測量精度為15 mm,重復(fù)測量精度為10 mm,激光發(fā)射波長為單色近紅外波段。4種植被的實驗樣方大小均為2 m×2 m,儀器架設(shè)距離樣方5 m,掃描儀距離地面高度2.5 m,該條件下掃描實際覆蓋掃描距離為950 m左右。完成TLS掃描后,修剪稀疏植被植株,以減小樣方內(nèi)的植被蓋度,再使用TLS對樣方實施稀疏后的測站掃描;重復(fù)該步驟,至修剪完植被的樣方掃描數(shù)據(jù)為真實地形數(shù)據(jù)。同時,TLS針對每個樣方進(jìn)行重復(fù)掃描測量實驗時,不移動測站位置,確保同一樣方的激光掃描測量擁有相同的自身坐標(biāo)系統(tǒng)SOCS(scanner’s own coordinate system)——以掃描儀測站位置為基礎(chǔ)的空間三維坐標(biāo)系。

    2.2 植被蓋度獲取及估算

    各樣方植被蓋度信息從數(shù)碼相片上獲取,植被蓋度定義為觀測區(qū)域內(nèi)植被垂直投影面積占地表面積的百分比[10]。數(shù)碼相機(jī)架設(shè)在距地面2.5 m高的位置,垂直對樣方拍攝。由于樣方的尺寸較大,人為修剪一次后,用數(shù)碼相機(jī)分別拍攝樣方的4個拐角,最后取蓋度均值作為樣方植被蓋度。數(shù)碼相片的像元值代表植被和土壤像素。一副相片的像素可分為植被和非植被兩種像素的二值圖。數(shù)碼相片上的植物蓋度通過計算植被像素占總像素的比例算得。本文使用Envi軟件監(jiān)督分類模塊進(jìn)行相片二值化處理,從原始到多次稀釋的結(jié)果為:蘆葦蓋度估值分別為50%、45%、40%、35%、30%、15%、0;白茅蓋度估值分別為70%、55%、40%、30%、20%、0;互花米草蓋度估值分別為65%、50%、40%、30%、15%、0;海三棱藨草蓋度估值分別為65%、60%、50%、40%、30%、20%、0。

    2.3 點云噪聲點去除

    使用TLS在野外獲取數(shù)據(jù)時,由于物體的遮擋以及其他外部條件的影響,所獲取的點云數(shù)據(jù)不可避免地會帶有噪聲點[11],這些噪聲點對構(gòu)建DEM(digital elevation model)不可避免地產(chǎn)生較大的干擾。因此,在獲取原始點云數(shù)據(jù)后,必須先進(jìn)行點云數(shù)據(jù)去噪處理,刪除噪聲點的工作通過點云可視化處理軟件Riscan pro完成。

    3 植被濾除算法構(gòu)建與反演精度評估

    3.1 植被濾除算法構(gòu)建

    點云數(shù)據(jù)包含三維坐標(biāo)、回光振幅、回光反射率以及顏色等屬性信息。本文所使用的植被濾除算法在移動窗口法[12]的基礎(chǔ)上,結(jié)合點云屬性信息對植被噪聲點進(jìn)行濾除。移動窗口法已經(jīng)被許多研究者用來去除植被影響[13-15],其思想是提取窗口內(nèi)的最低高程點為地面點,獲取到地面點的可能性與窗口的尺寸大小密切相關(guān)。當(dāng)窗口很小時,TLS返回的激光信號高程值可能高于地面實際高程,隨著窗口尺度的增加,該值越來越有可能是激光穿透植被冠層到達(dá)地面的返還點。然而窗口過大時,又有可能平滑或去掉一些小的地形不連續(xù)部分[16]。不同的植被樣方因植株高度與密度的不同,其最佳窗口尺度也不同。本文分別對4塊樣方點云數(shù)據(jù)使用移動窗口法過濾,窗口尺寸從1 cm到20 cm遞增,增量為1 cm。同時,激光掃描儀返回信號里包含的信息可以被用于地表不同材質(zhì)物體的分類,例如振幅屬性[17];采用非監(jiān)督分類K-means的方法去除分類結(jié)果為植被的點云。先移動窗口后去除植被點云的方法相比單一使用移動窗口法精度有進(jìn)一步的提升[6]。K-means算法是一種以平均值作為聚類中心的分割聚類算法。RIEGL VZ 4000 TLS激光反射信號為近紅外波段,對水體沒有回波反射信號,因此在研究區(qū)域內(nèi)激光反射信號只有植被和土壤兩類。因此,算法根據(jù)屬性值(本例中指點云振幅屬性)將簇群的n個目標(biāo)(本例中指TLS掃描點云數(shù)據(jù))分成k類(本例中指土壤和植被)。

    該算法的具體步驟如下:

    1)從數(shù)據(jù)集中隨機(jī)選擇k個對象C1,C2,…,Ck作為初始的聚類中心。

    2)把每個對象分配到與之距離最小的聚合,每個聚合用其中所有對象的均值來代表。對于每個點Vi(i=1,2, …,n),找出一個質(zhì)心Cj(j=1,2, …,k),使它們之間的距離d(Vi,Cj)最小,并把Vi分配到第j組。

    3)把所有的點都分配到相應(yīng)的組之后重新計算每個組的質(zhì)心Cj。

    4)循環(huán)執(zhí)行第2)步和第3)步,直到數(shù)據(jù)的劃分不再發(fā)生變化。

    選取裸土和覆蓋植被的土壤表面(以白茅植被為例)兩類點云數(shù)據(jù),兩種點云振幅屬性頻率分布如圖2所示。K-means算法對白茅植被分類結(jié)果為12,振幅小于12時為土壤點云數(shù)據(jù)、大于12時為鹽沼植被點云數(shù)據(jù)。該分類結(jié)果與裸土振幅頻率分布圖(圖2a)有較好的一致性,故去除振幅值大于12的點云。由于本文實驗樣方較小,同時部分植被存在激光穿不透現(xiàn)象,為確保樣方有足夠多點云構(gòu)建DEM,本方法僅用于去除個別高程異常高點,使得構(gòu)建樣方DEM點云不少于100個。

    3.2 植被樣方DEM構(gòu)建與精度評估

    以不同蓋度構(gòu)建的樣方DEM定量分析反演4種不同鹽沼樣方區(qū)域的地形精度。在完成上述去除低質(zhì)量點以及植被濾除后,為了便于后續(xù)精度評估,對點云數(shù)據(jù)進(jìn)行插值建模。本文基于Arcgis10.2軟件,使用克里金插值法構(gòu)建掃描樣方DEM。

    通過生成的DEM與完全去除植被后的DEM對比不同蓋度情況下的地形反演精度。不同情境下的實驗完全獨立進(jìn)行,使用均方根誤差(root-mean-square error,RMSE)來衡量算法生成DEM的精度。

    4 結(jié)果與分析

    4.1 測量精度與植被蓋度的關(guān)系

    根據(jù)實驗結(jié)果建立TLS反演4種不同植被樣方地表地形精度與蓋度、窗口尺度的關(guān)系,如圖3所示。從圖3可以看出,在窗口尺度為0~20 cm的范圍內(nèi),隨著窗口尺度的增加,TLS反演不同蓋度的4種植被樣方地表地形精度不斷提高。

    1)從圖3a可以看出:當(dāng)窗口尺度為20 cm、蘆葦蓋度為50%時,TLS反演蘆葦樣方地表地形精度最低,RMSE為4.5 cm;在窗口尺度為20 cm、蘆葦蓋度小于50%時,RMSE均在4.5 cm以內(nèi)。因此,激光基本上可以穿透蓋度為50%的蘆葦冠層到達(dá)地面。

    2)從圖3b可以看出:在窗口尺度大于15 cm時,RMSE雖有降低,但降低的幅度較小,幾乎平行于坐標(biāo)軸。在窗口尺度為20 cm、白茅蓋度為70%時,TLS反演白茅樣方地表地形精度最低,RMSE為22 cm,激光無法穿透植被冠層到達(dá)地表。在該情況下,以10 cm窗口濾除算法構(gòu)建的DEM及原始數(shù)據(jù)地表DEM見圖4(以掃描儀的SOCS空間坐標(biāo)系為參考)。從圖4中可以看出:沿入射方向,激光可以部分穿透植被,但有效距離很短;沿著該方向,隨著與TLS距離的增加,激光穿透的能力越來越弱,呈現(xiàn)一個緩坡。樣方內(nèi)的白茅株高為50 cm左右,在樣方后部,頂層DEM距離地面DEM約為40 cm,激光穿透深度為10 cm。當(dāng)窗口尺寸為20 cm、RMSE為5 cm時,白茅蓋度為30%左右。

    3)從圖3c可以看出:在窗口尺度為20 cm、互花米草蓋度為65%時,TLS反演互花米草樣方地表地形精度最低,RMSE為22 cm,激光無法穿透植被冠層到達(dá)地表。在窗口尺度為20 cm、互花米草蓋度小于等于50%時,互花米草曲線關(guān)系圖與蘆葦曲線關(guān)系圖有一定的相似性。在窗口尺度為20 cm,互花米草、蘆葦蓋度為50%時,互花米草RMSE為5.7 cm,稍高于蘆葦(4.5 cm);當(dāng)兩者蓋度小于50%時精度相近。當(dāng)窗口尺寸為20 cm、RMSE為5.7 cm時,互花米草蓋度為50%左右。

    4)從圖3d可以看出:在窗口尺度為20 cm、海三棱藨草蓋度為65%時,TLS反演海三棱藨草樣方地表地形精度最低,RMSE為8.6 cm,激光無法穿透植被冠層到達(dá)地表。當(dāng)窗口尺寸為20 cm、RMSE為5 cm時,海三棱藨草蓋度為40%左右。

    a.裸土;b. 覆蓋白茅的土壤。圖2 點云振幅屬性頻率分布圖Fig.2 Frequency distribution of amplitude attributes of point cloud

    a. 蘆葦;b. 白茅;c. 互花米草;d. 海三棱藨草。圖3 植被過濾算法窗口尺度與RMSE關(guān)系圖Fig.3 Window sizes of vegetation filtering algorithm vs. RMSE

    a. 10 cm窗口尺度過濾算法構(gòu)建的DEM;b. 原始數(shù)據(jù)地表DEM。圖4 白茅樣方DEM示意圖Fig.4 DEM of Imperata cylindrica quadrat

    4.2 激光穿透植被分析

    當(dāng)蓋度大于50%時,白茅、互花米草、海三棱藨草樣方原始蓋度分別為70%、65%、65%,RMSE分別為22.0、22.0、8.6 cm,TLS激光均不能穿透植被冠層到達(dá)地面。以白茅樣方為例,圖5a為白茅樣方原始蓋度掃描結(jié)果,圖5b為經(jīng)過3次人為稀釋植被蓋度后的掃描結(jié)果:隨著植被蓋度的稀釋,激光穿透能力增強(qiáng)。當(dāng)激光穿不透時,任何算法都不能在原始數(shù)據(jù)缺失的情況下提取地表、地形、地貌特征。

    黃色點云為掃描儀掃描結(jié)果,藍(lán)色點云為修剪完植被的真實地形數(shù)據(jù)。圖5 白茅植被稀釋前掃描數(shù)據(jù)(a)與3次稀釋后掃描數(shù)據(jù)(b)Fig.5 Imperata cylindrica of TLS scanning before coverage reduced (a) vs. three coverage reduced (b)

    當(dāng)蓋度為50%時,TLS反演4種鹽沼植被樣方地形的誤差為蘆葦<互花米草<海三棱藨草<白茅,RMSE分別為4.5、5.7、6.6、16.0 cm。白茅樣方

    反演精度較差,地形反演精度與蘆葦相差3.56倍,海三棱藨草地形反演精度與蘆葦相差1.47倍,互花米草地形反演精度與蘆葦?shù)匦畏囱菥认嘟?。造成低矮植被比較高植被反演精度差的原因是植物構(gòu)型的影響。由圖6可見:蘆葦和互花米草(圖6a、c)植株較高,莖稈光桿長度較長,葉片主要在上部,下部的莖稈對激光穿透影響很小,因此能夠穿過蘆葦、互花米草莖稈縫隙到達(dá)地面;白茅(圖6b)樣方內(nèi)為白茅與小蘆葦?shù)幕焐鷧^(qū),該區(qū)域植被高度較低,莖稈光桿長度很短,上部分布了多層葉片,在同樣的TLS視角下,激光不但被表層植被葉子遮擋且很難穿過被葉片遮擋的較短光桿部分,致使誤差較高即圖4a現(xiàn)象;海三棱藨草(圖6d)株高同樣較低,表層根狀莖對激光穿透作用較蘆葦、互花米草莖稈影響更大,出現(xiàn)同樣蓋度情況下反演精度比蘆葦和互花米草差的現(xiàn)象。海三棱藨草比白茅誤差小的原因可能有兩種:一是白茅的平均株高比海三棱藨草高;二是海三棱藨草的根須莖比白茅葉片對激光穿透的影響更小。

    當(dāng)蓋度小于50%時,由于激光能夠穿透蘆葦、互花米草,TLS反演兩種鹽沼植被地形精度量級變化不明顯;TLS反演海三棱藨草精度小幅度提升;由于激光難穿透白茅,隨著蓋度的稀釋,TLS反演精度提升較為明顯。在同種蓋度情況下,海三棱藨草誤差小于白茅誤差。

    綜上所述,鹽沼植被蓋度越高,地形反演精度越低,兩者呈現(xiàn)負(fù)相關(guān)關(guān)系。值得注意的是,當(dāng)起始窗口尺度為1 cm時,蘆葦誤差≈互花米草誤差>白茅誤差>海三棱藨草誤差。這是算法所致,較小窗口獲取點云高程最小值時易獲取其值相對較大的點,導(dǎo)致整體樣方誤差偏高。

    a. 蘆葦;b. 白茅;c. 互花米草;d. 海三棱藨草。圖6 TLS掃描樣方示意圖Fig.6 Quadrat of TLS scanning

    5 結(jié)論與建議

    本文以長江口鹽沼潮灘區(qū)域典型4種植被:蘆葦、白茅、互花米草、海三棱藨草樣方的TLS掃描點云為數(shù)據(jù)源,構(gòu)建植被點云過濾算法,使用控制變量的方法固定儀器的高度、與樣方的距離和掃描分辨率;通過人為改變鹽沼植被蓋度的方法,研究植被蓋度與點云數(shù)據(jù)反演鹽沼潮灘區(qū)域地形精度之間的關(guān)系,探究TLS在反演鹽沼潮灘區(qū)域地形的適用性。通過初步分析有以下認(rèn)識:鹽沼植被蓋度越高,地形反演精度越低,兩者呈現(xiàn)負(fù)相關(guān)關(guān)系;反演4種鹽沼植被地形的精度達(dá)到5 cm左右的蓋度臨界值分別為:蘆葦、互花米草等于50%左右,白茅等于30%左右,海三棱藨草等于40%左右。當(dāng)蓋度等于50%時,TLS反演蘆葦、互花米草、白茅、海三棱藨草地形的精度分別為4.5、5.7、16.0、6.6 cm。當(dāng)蓋度小于50%時,TLS激光更容易穿透株高較高的禾本科鹽沼植被如蘆葦、互花米草以及株高較低的莎草科植被海三棱藨草,激光對株高較低的禾本科植被白茅的反演效果隨著蓋度的降低精度提高較多。

    本文通過TLS反演小樣方鹽沼植被的實驗,探究了TLS在反演鹽沼潮灘區(qū)域地形的適用性。然而鹽沼潮灘地區(qū)植被大面積分布,大片鹽沼植被蓋度均大于50%,距離TLS越遠(yuǎn),激光點云密度越小,衰減作用越明顯,激光穿透植被能力越弱,遠(yuǎn)處激光點云數(shù)據(jù)可能只在植被表層。在鹽沼潮灘地區(qū)使用TLS反演地形時,增加TLS架設(shè)高度、對同一區(qū)域多方位反復(fù)掃描可能有助于提高地形反演效果。

    猜你喜歡
    白茅潮灘鹽沼
    白茅:與愛有關(guān),與美有關(guān)
    杭州灣北岸臨港新片區(qū)岸段潮灘沉積特征
    海洋通報(2022年5期)2022-11-30 12:03:48
    崇明東北潮灘泥沙輸運過程的研究
    海洋通報(2021年1期)2021-07-23 01:55:40
    潮灘干濕轉(zhuǎn)換的地貌發(fā)育物理模型及動力機(jī)制
    海洋通報(2020年3期)2020-10-20 09:26:46
    天空之鏡
    玻利維亞——烏尤尼鹽沼
    白茅純束
    江南榿木與蘆葦在潮灘濕地的適應(yīng)性對比研究
    国产精品国产三级国产av玫瑰| 午夜免费激情av| 国产男靠女视频免费网站| av天堂在线播放| 亚洲中文字幕日韩| 国产三级中文精品| 日日摸夜夜添夜夜添av毛片| 亚洲第一区二区三区不卡| 日韩av不卡免费在线播放| 22中文网久久字幕| 欧美日韩一区二区视频在线观看视频在线 | 国产精品久久视频播放| 国产蜜桃级精品一区二区三区| 黄色日韩在线| 小蜜桃在线观看免费完整版高清| 久久久久性生活片| 国产麻豆成人av免费视频| 国产片特级美女逼逼视频| 国产成人a∨麻豆精品| 日韩欧美免费精品| 国产亚洲欧美98| 好男人在线观看高清免费视频| 亚洲av不卡在线观看| 直男gayav资源| 亚洲av二区三区四区| av国产免费在线观看| 深爱激情五月婷婷| 中文字幕精品亚洲无线码一区| 国产乱人偷精品视频| 三级经典国产精品| 美女黄网站色视频| 亚洲高清免费不卡视频| 美女 人体艺术 gogo| 久久久色成人| 欧美成人免费av一区二区三区| 日韩av在线大香蕉| 听说在线观看完整版免费高清| 18禁黄网站禁片免费观看直播| 老熟妇乱子伦视频在线观看| 在线国产一区二区在线| 免费看a级黄色片| 亚洲无线在线观看| 国产69精品久久久久777片| 在线国产一区二区在线| 亚洲国产精品sss在线观看| 乱人视频在线观看| 51国产日韩欧美| 又黄又爽又刺激的免费视频.| 亚洲av五月六月丁香网| 亚洲中文字幕日韩| 亚洲av二区三区四区| 亚洲熟妇中文字幕五十中出| 免费av不卡在线播放| 精品不卡国产一区二区三区| h日本视频在线播放| 少妇的逼水好多| 淫妇啪啪啪对白视频| 在线a可以看的网站| 美女免费视频网站| 午夜老司机福利剧场| 亚洲天堂国产精品一区在线| 精品一区二区三区视频在线观看免费| 国产精品福利在线免费观看| 中文字幕av成人在线电影| 久久久久九九精品影院| 亚洲国产欧洲综合997久久,| 少妇丰满av| 免费av毛片视频| 黑人高潮一二区| 国产亚洲精品久久久com| 精品久久久久久久久久久久久| 一级毛片aaaaaa免费看小| 国产激情偷乱视频一区二区| 人妻夜夜爽99麻豆av| 美女黄网站色视频| 欧美xxxx性猛交bbbb| 国产精品野战在线观看| 特级一级黄色大片| 亚洲国产精品国产精品| 深夜精品福利| 精品人妻偷拍中文字幕| 欧美激情久久久久久爽电影| 日本三级黄在线观看| 别揉我奶头 嗯啊视频| 欧美不卡视频在线免费观看| 精品午夜福利在线看| 国产精品嫩草影院av在线观看| 国产精品久久久久久久久免| 特级一级黄色大片| 婷婷六月久久综合丁香| 亚洲四区av| 中文字幕久久专区| 最好的美女福利视频网| 男女视频在线观看网站免费| 99热这里只有是精品50| 此物有八面人人有两片| 久久6这里有精品| 欧洲精品卡2卡3卡4卡5卡区| 身体一侧抽搐| 深爱激情五月婷婷| 老师上课跳d突然被开到最大视频| 黄色一级大片看看| 亚州av有码| 亚洲欧美中文字幕日韩二区| 欧美绝顶高潮抽搐喷水| 亚洲一级一片aⅴ在线观看| 给我免费播放毛片高清在线观看| 丝袜喷水一区| 男人舔奶头视频| 午夜亚洲福利在线播放| 亚洲美女搞黄在线观看 | 亚洲一级一片aⅴ在线观看| 国产欧美日韩一区二区精品| 国产精品亚洲美女久久久| 97人妻精品一区二区三区麻豆| 最新中文字幕久久久久| 长腿黑丝高跟| 精品久久国产蜜桃| 听说在线观看完整版免费高清| 九九热线精品视视频播放| 欧美zozozo另类| 亚洲五月天丁香| 人妻少妇偷人精品九色| 如何舔出高潮| 亚洲av熟女| 久久久国产成人免费| 一级av片app| 国产午夜精品论理片| 国产高清视频在线播放一区| 亚洲精品色激情综合| 国产美女午夜福利| 一个人看的www免费观看视频| 久久婷婷人人爽人人干人人爱| 中国美女看黄片| 久久久精品大字幕| 97人妻精品一区二区三区麻豆| 亚洲无线在线观看| 嫩草影院新地址| 国产欧美日韩精品一区二区| 欧美丝袜亚洲另类| 最近在线观看免费完整版| 99热这里只有是精品在线观看| 俺也久久电影网| 国产v大片淫在线免费观看| 欧美日本视频| 婷婷精品国产亚洲av在线| 亚洲一级一片aⅴ在线观看| 国产精品日韩av在线免费观看| 亚洲精品久久国产高清桃花| 精品99又大又爽又粗少妇毛片| 国产v大片淫在线免费观看| 国产精品电影一区二区三区| 99久久九九国产精品国产免费| 国产极品精品免费视频能看的| 国产精品久久久久久久电影| 日本 av在线| 校园人妻丝袜中文字幕| 欧美xxxx性猛交bbbb| 美女被艹到高潮喷水动态| 高清午夜精品一区二区三区 | 国产精品不卡视频一区二区| 大香蕉久久网| 欧美+亚洲+日韩+国产| 免费电影在线观看免费观看| 最近中文字幕高清免费大全6| 性插视频无遮挡在线免费观看| 一进一出抽搐动态| 国产大屁股一区二区在线视频| 中文字幕av在线有码专区| h日本视频在线播放| 午夜福利18| 久久精品影院6| 岛国在线免费视频观看| 久久午夜福利片| 97超级碰碰碰精品色视频在线观看| 久久精品国产亚洲av天美| 看免费成人av毛片| 中文亚洲av片在线观看爽| 不卡视频在线观看欧美| 九九在线视频观看精品| 国产男靠女视频免费网站| 亚洲精品色激情综合| 免费人成视频x8x8入口观看| 麻豆av噜噜一区二区三区| av专区在线播放| 日韩精品中文字幕看吧| 欧美三级亚洲精品| 午夜福利在线在线| 国产高清视频在线观看网站| 国产伦在线观看视频一区| 三级男女做爰猛烈吃奶摸视频| 国产女主播在线喷水免费视频网站 | 久久国内精品自在自线图片| 亚洲精品456在线播放app| 国产aⅴ精品一区二区三区波| 老师上课跳d突然被开到最大视频| 伦精品一区二区三区| 波多野结衣高清无吗| 国产伦精品一区二区三区视频9| 国产午夜精品久久久久久一区二区三区 | 日本与韩国留学比较| 91麻豆精品激情在线观看国产| 天美传媒精品一区二区| 麻豆av噜噜一区二区三区| 午夜日韩欧美国产| 精品国产三级普通话版| 中文字幕精品亚洲无线码一区| 亚洲欧美日韩东京热| 免费电影在线观看免费观看| 欧美极品一区二区三区四区| 亚洲激情五月婷婷啪啪| 亚洲最大成人手机在线| 国产精品亚洲美女久久久| 美女黄网站色视频| 亚洲久久久久久中文字幕| 两个人的视频大全免费| 日韩欧美精品免费久久| 亚洲成av人片在线播放无| 久久婷婷人人爽人人干人人爱| 国产淫片久久久久久久久| 少妇熟女aⅴ在线视频| 长腿黑丝高跟| 成人亚洲精品av一区二区| 12—13女人毛片做爰片一| 日韩成人av中文字幕在线观看 | 长腿黑丝高跟| 久久久久精品国产欧美久久久| 国产高潮美女av| 丝袜美腿在线中文| 欧美+日韩+精品| 丰满人妻一区二区三区视频av| 最近的中文字幕免费完整| 亚洲色图av天堂| 亚洲无线在线观看| 中文字幕人妻熟人妻熟丝袜美| 黄色配什么色好看| 狂野欧美激情性xxxx在线观看| 亚洲av不卡在线观看| 22中文网久久字幕| 久久久色成人| 网址你懂的国产日韩在线| 色哟哟·www| 亚洲国产精品国产精品| 久久精品国产自在天天线| 好男人在线观看高清免费视频| 赤兔流量卡办理| 97人妻精品一区二区三区麻豆| 美女黄网站色视频| 狂野欧美白嫩少妇大欣赏| 国产探花极品一区二区| 日韩欧美一区二区三区在线观看| 18禁在线无遮挡免费观看视频 | 久久精品国产亚洲av香蕉五月| 嫩草影院入口| 99热6这里只有精品| 欧美+日韩+精品| 91av网一区二区| 亚洲高清免费不卡视频| 亚洲美女搞黄在线观看 | 免费不卡的大黄色大毛片视频在线观看 | 日韩一本色道免费dvd| 国产一区二区在线观看日韩| 在现免费观看毛片| 欧美高清成人免费视频www| avwww免费| 亚洲国产精品成人久久小说 | 精品久久久久久久人妻蜜臀av| 两个人的视频大全免费| 欧美日韩乱码在线| 国产精品久久电影中文字幕| 国产精品1区2区在线观看.| 久久亚洲精品不卡| 欧美xxxx黑人xx丫x性爽| 中文字幕久久专区| 国产熟女欧美一区二区| av国产免费在线观看| 日本a在线网址| 精品日产1卡2卡| 大香蕉久久网| 91麻豆精品激情在线观看国产| 校园人妻丝袜中文字幕| 亚洲美女搞黄在线观看 | 国产日本99.免费观看| 在线天堂最新版资源| 91在线精品国自产拍蜜月| 国产精品野战在线观看| 国产亚洲91精品色在线| 中文在线观看免费www的网站| 国产精品久久久久久亚洲av鲁大| 丰满乱子伦码专区| 99久国产av精品| 亚洲中文字幕一区二区三区有码在线看| 国产精品一区二区三区四区久久| 最近2019中文字幕mv第一页| 亚洲丝袜综合中文字幕| 久久精品夜夜夜夜夜久久蜜豆| 亚洲欧美清纯卡通| 人妻久久中文字幕网| 一级黄色大片毛片| 国产一区二区激情短视频| 一级av片app| 在线观看66精品国产| 国产精品三级大全| 国产精品久久久久久久久免| 国产成人freesex在线 | 十八禁国产超污无遮挡网站| 日本 av在线| 国产单亲对白刺激| av天堂在线播放| 简卡轻食公司| 日韩人妻高清精品专区| 精品一区二区三区视频在线| 国产女主播在线喷水免费视频网站 | 全区人妻精品视频| 午夜亚洲福利在线播放| 日日摸夜夜添夜夜添av毛片| 深夜精品福利| 精品一区二区三区av网在线观看| 一区福利在线观看| 欧美成人精品欧美一级黄| 最好的美女福利视频网| 最近2019中文字幕mv第一页| 日韩一本色道免费dvd| 99在线视频只有这里精品首页| av天堂中文字幕网| 99国产精品一区二区蜜桃av| 两性午夜刺激爽爽歪歪视频在线观看| 久久久久久国产a免费观看| 久久久精品大字幕| 长腿黑丝高跟| 啦啦啦韩国在线观看视频| 日韩欧美三级三区| 国产在线精品亚洲第一网站| 国产高清视频在线观看网站| 亚洲三级黄色毛片| 色哟哟·www| 国产色婷婷99| 床上黄色一级片| 99热这里只有是精品50| 国产亚洲91精品色在线| 成人一区二区视频在线观看| 国产黄a三级三级三级人| 丝袜喷水一区| 亚洲国产欧洲综合997久久,| 啦啦啦韩国在线观看视频| 国产精品一区二区三区四区久久| 国产在视频线在精品| 日韩制服骚丝袜av| 国产在视频线在精品| 精品一区二区免费观看| 成年av动漫网址| 一级a爱片免费观看的视频| 麻豆av噜噜一区二区三区| av中文乱码字幕在线| 久99久视频精品免费| 亚洲第一区二区三区不卡| 男人的好看免费观看在线视频| 国产黄片美女视频| 最后的刺客免费高清国语| 亚洲精品一卡2卡三卡4卡5卡| 日本爱情动作片www.在线观看 | 丝袜美腿在线中文| 黄色日韩在线| 欧美xxxx黑人xx丫x性爽| 国产精品1区2区在线观看.| 少妇猛男粗大的猛烈进出视频 | 国产不卡一卡二| 成年女人毛片免费观看观看9| 国产乱人视频| 麻豆乱淫一区二区| 亚洲国产精品sss在线观看| 女人被狂操c到高潮| 最好的美女福利视频网| 99久国产av精品| 亚洲精品成人久久久久久| 能在线免费观看的黄片| 日本五十路高清| a级毛色黄片| 午夜福利高清视频| 日本撒尿小便嘘嘘汇集6| 午夜福利在线在线| 亚洲精品456在线播放app| 一级a爱片免费观看的视频| 如何舔出高潮| 日韩制服骚丝袜av| 亚洲人成网站高清观看| 久久久久久久久大av| 九九爱精品视频在线观看| 在线观看av片永久免费下载| 日日摸夜夜添夜夜爱| 国产激情偷乱视频一区二区| 男女做爰动态图高潮gif福利片| 国产精品三级大全| 国内久久婷婷六月综合欲色啪| 亚洲人成网站在线播放欧美日韩| 国产精品爽爽va在线观看网站| 亚洲美女搞黄在线观看 | 乱系列少妇在线播放| 一级黄色大片毛片| 午夜精品在线福利| 国产精品综合久久久久久久免费| 一进一出抽搐动态| 2021天堂中文幕一二区在线观| 亚洲人与动物交配视频| 日产精品乱码卡一卡2卡三| 久久精品国产亚洲网站| 婷婷色综合大香蕉| 亚洲第一电影网av| 国产精品人妻久久久影院| 欧美潮喷喷水| 嫩草影院精品99| 欧美+日韩+精品| 欧美激情在线99| 在线观看66精品国产| 嫩草影院新地址| 麻豆一二三区av精品| 国产一区二区亚洲精品在线观看| 啦啦啦观看免费观看视频高清| 国产真实伦视频高清在线观看| 午夜精品在线福利| 国产三级中文精品| 欧美成人精品欧美一级黄| 男女做爰动态图高潮gif福利片| 亚洲图色成人| 国产成年人精品一区二区| 天堂网av新在线| 国产一区二区在线av高清观看| 联通29元200g的流量卡| or卡值多少钱| 少妇高潮的动态图| 最近手机中文字幕大全| 日韩人妻高清精品专区| 亚洲精品456在线播放app| 精品国产三级普通话版| 天堂网av新在线| 欧美性猛交黑人性爽| 国产伦一二天堂av在线观看| 毛片一级片免费看久久久久| 欧美日韩一区二区视频在线观看视频在线 | 亚洲av熟女| 亚洲精品一卡2卡三卡4卡5卡| 亚洲aⅴ乱码一区二区在线播放| 大又大粗又爽又黄少妇毛片口| 99久久久亚洲精品蜜臀av| 日韩av不卡免费在线播放| 卡戴珊不雅视频在线播放| 国产欧美日韩精品一区二区| 久久久久久久久大av| 美女 人体艺术 gogo| 少妇熟女aⅴ在线视频| 久久精品影院6| 久久人妻av系列| 成人三级黄色视频| 欧美成人a在线观看| 国产成人freesex在线 | 国产精品日韩av在线免费观看| 成年女人毛片免费观看观看9| 色哟哟哟哟哟哟| 一本久久中文字幕| 亚洲av二区三区四区| 最近中文字幕高清免费大全6| 日日摸夜夜添夜夜添小说| 久久久久久久久久黄片| 久久久久久久久久久丰满| 亚洲第一电影网av| 国产一区二区在线av高清观看| 国产欧美日韩精品一区二区| 午夜亚洲福利在线播放| 欧美日韩综合久久久久久| 国产真实乱freesex| 亚洲美女黄片视频| 级片在线观看| 有码 亚洲区| 欧美激情国产日韩精品一区| 最后的刺客免费高清国语| 国产爱豆传媒在线观看| 久久精品国产99精品国产亚洲性色| av国产免费在线观看| 老师上课跳d突然被开到最大视频| 国产一区亚洲一区在线观看| 校园春色视频在线观看| 国产不卡一卡二| 国产精品一区二区免费欧美| 国产欧美日韩精品亚洲av| 国产成人一区二区在线| 色播亚洲综合网| 在线天堂最新版资源| 欧美xxxx性猛交bbbb| 欧美日韩精品成人综合77777| 日本三级黄在线观看| 国产欧美日韩精品亚洲av| av在线亚洲专区| 嫩草影院入口| 大型黄色视频在线免费观看| 中文字幕免费在线视频6| 99久久中文字幕三级久久日本| 欧美一级a爱片免费观看看| av在线天堂中文字幕| 在线观看66精品国产| 别揉我奶头~嗯~啊~动态视频| 国产男靠女视频免费网站| 高清毛片免费看| 久久精品国产清高在天天线| 亚洲中文日韩欧美视频| 我要看日韩黄色一级片| 亚洲成人av在线免费| 51国产日韩欧美| 国内揄拍国产精品人妻在线| 国产欧美日韩精品一区二区| 99国产精品一区二区蜜桃av| 日韩欧美 国产精品| 热99re8久久精品国产| 日韩精品有码人妻一区| 人妻丰满熟妇av一区二区三区| 久久久久久大精品| 亚洲国产精品合色在线| 免费观看的影片在线观看| 亚洲精品一区av在线观看| 校园春色视频在线观看| 国产一区二区在线av高清观看| 国产午夜精品论理片| 欧美xxxx黑人xx丫x性爽| 麻豆久久精品国产亚洲av| 狂野欧美白嫩少妇大欣赏| 国产av在哪里看| 蜜桃亚洲精品一区二区三区| 国产爱豆传媒在线观看| 亚洲精品成人久久久久久| 亚洲精品456在线播放app| 国产高清三级在线| 久久中文看片网| 久久久久久久亚洲中文字幕| 天堂动漫精品| 性欧美人与动物交配| 午夜福利18| 99热这里只有精品一区| 国产精品av视频在线免费观看| 国产日本99.免费观看| 亚洲欧美日韩无卡精品| 两个人视频免费观看高清| av天堂中文字幕网| 好男人在线观看高清免费视频| 国内精品美女久久久久久| 美女内射精品一级片tv| 变态另类成人亚洲欧美熟女| 99riav亚洲国产免费| 国产午夜福利久久久久久| 好男人在线观看高清免费视频| 国产高潮美女av| 久久精品夜色国产| 91午夜精品亚洲一区二区三区| 国产淫片久久久久久久久| 此物有八面人人有两片| 成人鲁丝片一二三区免费| 色噜噜av男人的天堂激情| 午夜久久久久精精品| 丰满人妻一区二区三区视频av| 国产aⅴ精品一区二区三区波| 精品一区二区免费观看| 成人美女网站在线观看视频| 久久久久久国产a免费观看| 欧美人与善性xxx| 亚洲av熟女| 亚洲电影在线观看av| 女人被狂操c到高潮| 久久人人爽人人爽人人片va| 亚洲欧美日韩东京热| 久99久视频精品免费| 日韩欧美一区二区三区在线观看| 少妇的逼水好多| 高清毛片免费观看视频网站| 久久国产乱子免费精品| 国产av麻豆久久久久久久| 美女内射精品一级片tv| 国产精品福利在线免费观看| 日韩大尺度精品在线看网址| 欧美zozozo另类| 男人舔奶头视频| 国产色爽女视频免费观看| 嫩草影院新地址| 日韩人妻高清精品专区| 一本精品99久久精品77| 日韩成人av中文字幕在线观看 | 久久这里只有精品中国| 欧美日本视频| 国产精品嫩草影院av在线观看| 深夜a级毛片| 国产欧美日韩精品一区二区| 性插视频无遮挡在线免费观看| 伊人久久精品亚洲午夜| 最近2019中文字幕mv第一页| 国产精品,欧美在线| 看免费成人av毛片| 中文字幕精品亚洲无线码一区| 九九久久精品国产亚洲av麻豆| 欧美潮喷喷水| 日韩高清综合在线| 国产午夜精品久久久久久一区二区三区 | 久久久欧美国产精品| 一级毛片我不卡| 亚洲精品在线观看二区| 亚洲人成网站高清观看| 99久久无色码亚洲精品果冻| 自拍偷自拍亚洲精品老妇| 精品国内亚洲2022精品成人| 欧美成人一区二区免费高清观看| 大型黄色视频在线免费观看| 国产 一区 欧美 日韩| 日韩成人av中文字幕在线观看 | 国产精品伦人一区二区| 亚洲内射少妇av| 在线天堂最新版资源| 国产91av在线免费观看| 18禁在线播放成人免费| 午夜福利高清视频|