• <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一区二区三区| 亚洲欧美日韩高清在线视频| 精品少妇一区二区三区视频日本电影| 亚洲一区二区三区色噜噜 | 亚洲精品在线观看二区| 美女高潮喷水抽搐中文字幕| www.www免费av| 国产亚洲欧美精品永久| 国产免费男女视频| 色在线成人网| 久久99一区二区三区| 久99久视频精品免费| 国产在线精品亚洲第一网站| 亚洲va日本ⅴa欧美va伊人久久| 韩国精品一区二区三区| 日韩视频一区二区在线观看| 夜夜夜夜夜久久久久| 久久国产精品人妻蜜桃| 成人永久免费在线观看视频| 一级毛片高清免费大全| 一个人免费在线观看的高清视频| 天堂√8在线中文| 男人舔女人下体高潮全视频| 婷婷六月久久综合丁香| 久久精品aⅴ一区二区三区四区| 国产免费av片在线观看野外av| 亚洲自拍偷在线| 久久狼人影院| 国产av一区二区精品久久| a级片在线免费高清观看视频| 黄片小视频在线播放| 久久精品影院6| 亚洲性夜色夜夜综合| 国产成人精品久久二区二区91| 国产精品 国内视频| 国产成人影院久久av| 亚洲人成电影免费在线| xxxhd国产人妻xxx| 午夜免费成人在线视频| 国产伦人伦偷精品视频| 欧美成人免费av一区二区三区| 欧美人与性动交α欧美软件| 日韩精品免费视频一区二区三区| 视频在线观看一区二区三区| 亚洲精品久久午夜乱码| 国产精品影院久久| 51午夜福利影视在线观看| 97超级碰碰碰精品色视频在线观看| 一级毛片高清免费大全| 丁香欧美五月| 一级,二级,三级黄色视频| 精品久久久久久成人av| 欧美成人午夜精品| 男女床上黄色一级片免费看| 免费在线观看黄色视频的| 久久国产亚洲av麻豆专区| 91字幕亚洲| 国产亚洲欧美在线一区二区| 成人国语在线视频| 久久精品亚洲熟妇少妇任你| 亚洲视频免费观看视频| 国产熟女午夜一区二区三区| 麻豆一二三区av精品| 18禁观看日本| 免费观看人在逋| 女人被狂操c到高潮| 亚洲一区二区三区不卡视频| 女人高潮潮喷娇喘18禁视频| 亚洲七黄色美女视频| 大码成人一级视频| 亚洲五月色婷婷综合| 91麻豆精品激情在线观看国产 | 在线十欧美十亚洲十日本专区| 精品久久久久久,| 免费女性裸体啪啪无遮挡网站| 少妇的丰满在线观看| 国产免费男女视频| 亚洲成人国产一区在线观看| 别揉我奶头~嗯~啊~动态视频| 亚洲精品av麻豆狂野| 亚洲av成人不卡在线观看播放网| av在线天堂中文字幕 | 夜夜看夜夜爽夜夜摸 | 久久伊人香网站| 久热爱精品视频在线9| 精品国产美女av久久久久小说| 久久精品成人免费网站| 久久国产精品人妻蜜桃| 99国产极品粉嫩在线观看| 嫩草影院精品99| av片东京热男人的天堂| 777久久人妻少妇嫩草av网站| 丰满的人妻完整版| 成人黄色视频免费在线看| 人人澡人人妻人| 亚洲国产欧美日韩在线播放| 男女之事视频高清在线观看| 国产高清国产精品国产三级| 欧美一区二区精品小视频在线| 国产一区二区激情短视频| aaaaa片日本免费| 好看av亚洲va欧美ⅴa在| 国产三级在线视频| 久久精品91无色码中文字幕| 成人精品一区二区免费| 亚洲国产看品久久| 日韩人妻精品一区2区三区| 伦理电影免费视频| 亚洲 国产 在线| 大陆偷拍与自拍| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品99久久99久久久不卡| 国产一卡二卡三卡精品| 人人妻人人澡人人看| 男男h啪啪无遮挡| 日本五十路高清| 黄频高清免费视频| 国产精品久久电影中文字幕| 亚洲专区中文字幕在线| 亚洲aⅴ乱码一区二区在线播放 | 免费高清在线观看日韩| 中文字幕色久视频| 国产99白浆流出| 91成年电影在线观看| 精品久久久久久电影网| 免费在线观看日本一区| 岛国在线观看网站| 999久久久精品免费观看国产| 日韩 欧美 亚洲 中文字幕| 宅男免费午夜| 亚洲七黄色美女视频| 天天添夜夜摸| 中文欧美无线码| 黄色丝袜av网址大全| 99久久99久久久精品蜜桃| 亚洲精品成人av观看孕妇| 国产男靠女视频免费网站| 日本欧美视频一区| 欧美老熟妇乱子伦牲交| 国产精品日韩av在线免费观看 | 夜夜爽天天搞| 日本一区二区免费在线视频| 91大片在线观看| 国产免费男女视频| 亚洲全国av大片| 国产一卡二卡三卡精品| 久久99一区二区三区| 国产成人av教育| 我的亚洲天堂| 国产精品一区二区精品视频观看| 亚洲成人精品中文字幕电影 | 欧美乱妇无乱码| 熟女少妇亚洲综合色aaa.| 免费观看精品视频网站| 亚洲精品成人av观看孕妇| 天天影视国产精品| 国产精品爽爽va在线观看网站 | 国产精品亚洲一级av第二区| 亚洲第一av免费看| 琪琪午夜伦伦电影理论片6080| 久久久国产成人免费| 午夜91福利影院| 国产不卡一卡二| 桃红色精品国产亚洲av| 不卡av一区二区三区| 在线免费观看的www视频| 国产精品美女特级片免费视频播放器 | 高清av免费在线| 午夜两性在线视频| 国产午夜精品久久久久久| 99在线视频只有这里精品首页| 欧美中文综合在线视频| 精品国产国语对白av| 亚洲成人免费电影在线观看| 亚洲第一av免费看| 亚洲中文av在线| 国产精品成人在线| 狂野欧美激情性xxxx| 国产成人免费无遮挡视频| 热re99久久精品国产66热6| 一区福利在线观看| 日韩av在线大香蕉| 美女 人体艺术 gogo| 99国产精品99久久久久| 黄色 视频免费看| 宅男免费午夜| 男女之事视频高清在线观看| 国产国语露脸激情在线看| 国产av一区在线观看免费| 又黄又爽又免费观看的视频| 在线视频色国产色| videosex国产| 欧美人与性动交α欧美精品济南到| 成人影院久久| 老熟妇乱子伦视频在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 久久久久久久午夜电影 | 亚洲成人久久性| 亚洲av成人av| √禁漫天堂资源中文www| 亚洲成人久久性| 欧美成人性av电影在线观看| 亚洲 欧美一区二区三区| 成人三级黄色视频| 波多野结衣一区麻豆| 亚洲欧美一区二区三区久久| 波多野结衣av一区二区av| 欧美激情极品国产一区二区三区| 国产精品久久视频播放| 99久久精品国产亚洲精品| 亚洲欧美一区二区三区久久| 97碰自拍视频| 欧美黑人欧美精品刺激| 欧美在线一区亚洲| √禁漫天堂资源中文www| 不卡av一区二区三区| 国产亚洲精品第一综合不卡| 淫秽高清视频在线观看| 黄色片一级片一级黄色片| 亚洲 欧美一区二区三区| 国产一区二区三区视频了| 国产一区二区激情短视频| 美女扒开内裤让男人捅视频| 色在线成人网| 日本wwww免费看| 欧美成人午夜精品| 亚洲一码二码三码区别大吗| 丝袜美足系列| 亚洲精品美女久久久久99蜜臀| 久久久久国产一级毛片高清牌| 欧美日韩黄片免| 久久99一区二区三区| 国产精品野战在线观看 | 日韩一卡2卡3卡4卡2021年| 久久久久久久久免费视频了| 精品久久蜜臀av无| 欧美不卡视频在线免费观看 | 男女之事视频高清在线观看| 欧美久久黑人一区二区| 天天躁夜夜躁狠狠躁躁| 中文字幕人妻丝袜制服| 高清在线国产一区| 国产高清视频在线播放一区| 国产成人精品无人区| 两个人看的免费小视频| 日韩免费av在线播放| 一个人免费在线观看的高清视频| 黄色片一级片一级黄色片| 色老头精品视频在线观看| 精品乱码久久久久久99久播| 久久久精品国产亚洲av高清涩受| 美女福利国产在线| 久久人人97超碰香蕉20202| 国产熟女xx| 韩国精品一区二区三区| 国产一区二区三区在线臀色熟女 | 天天影视国产精品| 母亲3免费完整高清在线观看| 国产日韩一区二区三区精品不卡| 老汉色∧v一级毛片| 成在线人永久免费视频| 国产精品 欧美亚洲| 国产成人欧美| 天天添夜夜摸| 桃色一区二区三区在线观看| 日韩国内少妇激情av| 精品人妻在线不人妻| 一边摸一边抽搐一进一小说| 国内久久婷婷六月综合欲色啪| 一边摸一边做爽爽视频免费| 日本精品一区二区三区蜜桃| 亚洲aⅴ乱码一区二区在线播放 | 脱女人内裤的视频| 欧美一区二区精品小视频在线| 国产一区二区在线av高清观看| 无人区码免费观看不卡| 在线国产一区二区在线| 精品人妻1区二区| 欧美激情高清一区二区三区| 91麻豆av在线| 国产精品国产高清国产av| 国产主播在线观看一区二区| 中文欧美无线码| 看黄色毛片网站| 国产欧美日韩一区二区三| 国产成人欧美在线观看| 新久久久久国产一级毛片| tocl精华| 亚洲中文字幕日韩| 一级片免费观看大全| 咕卡用的链子| 嫩草影院精品99| 啦啦啦在线免费观看视频4| 国产精品国产av在线观看| 国产成人欧美在线观看| 精品久久久久久,| 99久久国产精品久久久| 久久精品91蜜桃| 嫩草影院精品99| 人人妻人人澡人人看| 精品久久久久久久毛片微露脸| 成人国语在线视频| 欧美 亚洲 国产 日韩一| 在线免费观看的www视频| 欧美中文日本在线观看视频| 国产精品一区二区精品视频观看| 巨乳人妻的诱惑在线观看| 一级,二级,三级黄色视频| 欧美久久黑人一区二区| 日韩精品青青久久久久久| 久久人妻福利社区极品人妻图片| 欧美中文日本在线观看视频| 久热爱精品视频在线9| 视频区欧美日本亚洲| bbb黄色大片| 精品国产亚洲在线| 一个人观看的视频www高清免费观看 | 女人爽到高潮嗷嗷叫在线视频| 欧美日本亚洲视频在线播放| 久久人妻熟女aⅴ| 亚洲 欧美一区二区三区| 男女做爰动态图高潮gif福利片 | 99久久国产精品久久久| 又紧又爽又黄一区二区| 国产精品av久久久久免费| 我的亚洲天堂| 波多野结衣av一区二区av| 后天国语完整版免费观看| 午夜91福利影院| 欧美在线黄色| 欧美成狂野欧美在线观看| 男男h啪啪无遮挡| 巨乳人妻的诱惑在线观看| 如日韩欧美国产精品一区二区三区| 又黄又爽又免费观看的视频| 人人妻人人爽人人添夜夜欢视频| 一进一出抽搐动态| 亚洲精品成人av观看孕妇| 亚洲国产精品一区二区三区在线| 可以在线观看毛片的网站| 亚洲成人免费av在线播放| 成人手机av| 手机成人av网站| 国产精品av久久久久免费| 免费一级毛片在线播放高清视频 | 日本a在线网址| 真人做人爱边吃奶动态| 国产精品香港三级国产av潘金莲| 久9热在线精品视频| aaaaa片日本免费| 欧美一级毛片孕妇| 国产高清视频在线播放一区| 18美女黄网站色大片免费观看| 亚洲av成人一区二区三| 久久精品国产亚洲av香蕉五月| 在线观看日韩欧美| 国产精品99久久99久久久不卡| 日本免费a在线| 亚洲 欧美一区二区三区| 亚洲精品一二三| 大码成人一级视频| 午夜两性在线视频| 叶爱在线成人免费视频播放| 我的亚洲天堂| 久热这里只有精品99| 两人在一起打扑克的视频| 国产乱人伦免费视频| 国产精品自产拍在线观看55亚洲| 高清av免费在线| 亚洲欧洲精品一区二区精品久久久| 我的亚洲天堂| 亚洲美女黄片视频| 视频在线观看一区二区三区| 久久国产亚洲av麻豆专区| 80岁老熟妇乱子伦牲交| 欧美日韩视频精品一区| 精品国产国语对白av| 男女床上黄色一级片免费看| 伦理电影免费视频| 日韩精品免费视频一区二区三区| 狂野欧美激情性xxxx| 欧美+亚洲+日韩+国产| 中文欧美无线码| 精品国产亚洲在线| 欧美日韩视频精品一区| www.自偷自拍.com| av片东京热男人的天堂| 国产又爽黄色视频| 欧美成人免费av一区二区三区| 亚洲国产看品久久| 色尼玛亚洲综合影院| 久久天堂一区二区三区四区| 黄色怎么调成土黄色| 欧美激情极品国产一区二区三区| 午夜91福利影院| 中文字幕人妻丝袜制服| 9热在线视频观看99| a级毛片在线看网站| 欧美成人午夜精品| 精品一区二区三区视频在线观看免费 | 久久精品亚洲av国产电影网| 免费在线观看视频国产中文字幕亚洲| 日韩 欧美 亚洲 中文字幕| 国产亚洲精品第一综合不卡| 久久精品91无色码中文字幕| 99久久99久久久精品蜜桃| 男女做爰动态图高潮gif福利片 | 亚洲国产看品久久| 国产欧美日韩综合在线一区二区| 国产精品国产高清国产av| 在线天堂中文资源库| 少妇的丰满在线观看| 999精品在线视频| 亚洲精品中文字幕一二三四区| 精品国产乱子伦一区二区三区| 巨乳人妻的诱惑在线观看| 真人一进一出gif抽搐免费| 亚洲国产精品合色在线| 国产一区二区激情短视频| 夫妻午夜视频| 日韩欧美国产一区二区入口| 男人的好看免费观看在线视频 | 亚洲成人免费电影在线观看| 亚洲七黄色美女视频| 99精品在免费线老司机午夜| 宅男免费午夜| 亚洲avbb在线观看| 天堂俺去俺来也www色官网| 女人被躁到高潮嗷嗷叫费观| 久久99一区二区三区| 国产激情欧美一区二区| 一区福利在线观看| 亚洲第一欧美日韩一区二区三区| 亚洲成人久久性| 久久精品成人免费网站| 成人精品一区二区免费| 亚洲,欧美精品.| 一级片'在线观看视频| 一级作爱视频免费观看| 丝袜美足系列| 十八禁网站免费在线| 成人18禁高潮啪啪吃奶动态图| 精品人妻在线不人妻| xxx96com| 色综合站精品国产| 国产亚洲精品综合一区在线观看 | 国产伦人伦偷精品视频| 久久久久久久久中文| 成年人免费黄色播放视频| 9色porny在线观看| 免费av毛片视频| 婷婷丁香在线五月| 丁香欧美五月| 性欧美人与动物交配| 欧美一区二区精品小视频在线| 这个男人来自地球电影免费观看| 欧美 亚洲 国产 日韩一| 国产三级黄色录像| 国产日韩一区二区三区精品不卡| 久久精品aⅴ一区二区三区四区| 久久久久国产一级毛片高清牌| 最好的美女福利视频网| 国产成人一区二区三区免费视频网站| 成人国产一区最新在线观看| 精品日产1卡2卡| 亚洲成人国产一区在线观看| 热re99久久国产66热| 丝袜美腿诱惑在线| netflix在线观看网站| 90打野战视频偷拍视频| 久久 成人 亚洲| 亚洲精品美女久久久久99蜜臀| 最新美女视频免费是黄的| 日韩欧美一区二区三区在线观看| 欧美成人免费av一区二区三区| 精品久久久久久,| 女人被狂操c到高潮| 午夜免费成人在线视频| 老司机午夜十八禁免费视频| 一区二区日韩欧美中文字幕| 窝窝影院91人妻| 精品免费久久久久久久清纯| 激情在线观看视频在线高清| 人成视频在线观看免费观看| 99久久精品国产亚洲精品| 久久精品亚洲精品国产色婷小说| 免费高清视频大片| 欧美久久黑人一区二区| 80岁老熟妇乱子伦牲交| 欧洲精品卡2卡3卡4卡5卡区| 国产精品久久久av美女十八| 极品人妻少妇av视频| 三上悠亚av全集在线观看| 波多野结衣一区麻豆| 热re99久久精品国产66热6| 黄片大片在线免费观看| 一级a爱视频在线免费观看| 欧美在线一区亚洲| 亚洲色图av天堂| av超薄肉色丝袜交足视频| 激情在线观看视频在线高清| 少妇裸体淫交视频免费看高清 | 国产一卡二卡三卡精品| 黑人操中国人逼视频| 久久天躁狠狠躁夜夜2o2o| 18禁观看日本| 多毛熟女@视频| 国产亚洲欧美在线一区二区| 老司机在亚洲福利影院| 中文字幕人妻熟女乱码| 一a级毛片在线观看| 精品一区二区三区av网在线观看| 久久中文字幕人妻熟女| 亚洲中文字幕日韩| 成人国语在线视频| 在线免费观看的www视频| 搡老熟女国产l中国老女人| 一边摸一边抽搐一进一小说| 国产av在哪里看| 成人18禁在线播放| 亚洲人成77777在线视频| 亚洲午夜精品一区,二区,三区| 久久久国产成人免费| 啦啦啦免费观看视频1| 日韩大尺度精品在线看网址 | 熟女少妇亚洲综合色aaa.| 欧美在线黄色| 啦啦啦在线免费观看视频4| 欧美在线一区亚洲| 亚洲全国av大片| 欧美中文综合在线视频| 久久久久久大精品| 黑人巨大精品欧美一区二区mp4| 人人妻人人澡人人看| 午夜福利一区二区在线看| 国产成人精品久久二区二区免费| 亚洲一区二区三区欧美精品| 精品人妻1区二区| 久久影院123| 少妇裸体淫交视频免费看高清 | 一边摸一边抽搐一进一小说| 日本五十路高清| 国产aⅴ精品一区二区三区波| 极品教师在线免费播放| 久久久久久免费高清国产稀缺| 久久草成人影院| 日本黄色日本黄色录像| 亚洲国产精品合色在线| 一进一出抽搐动态| 亚洲视频免费观看视频| 一本综合久久免费| 久久久久久久精品吃奶| 中文亚洲av片在线观看爽| 国产熟女午夜一区二区三区| 亚洲男人的天堂狠狠| 欧美日本亚洲视频在线播放| 国产单亲对白刺激| 天天影视国产精品| 欧洲精品卡2卡3卡4卡5卡区| 国产精品免费一区二区三区在线| 最新美女视频免费是黄的| 国产精品二区激情视频| 日日摸夜夜添夜夜添小说| 欧美激情极品国产一区二区三区| 看片在线看免费视频| 妹子高潮喷水视频| 男人操女人黄网站| 伊人久久大香线蕉亚洲五| 在线av久久热| 日本欧美视频一区| 丰满人妻熟妇乱又伦精品不卡| 在线免费观看的www视频| 中文字幕av电影在线播放| 91成年电影在线观看| 日韩国内少妇激情av| 亚洲精华国产精华精| 999久久久精品免费观看国产| 后天国语完整版免费观看| av在线播放免费不卡| 在线国产一区二区在线| 国产精华一区二区三区| 国产精品久久久av美女十八| 看片在线看免费视频| 日韩大码丰满熟妇| 国产又爽黄色视频| 精品高清国产在线一区| 国产精华一区二区三区| 一a级毛片在线观看| 国产亚洲精品综合一区在线观看 | 一级,二级,三级黄色视频| 久久午夜综合久久蜜桃| 男女床上黄色一级片免费看| xxxhd国产人妻xxx| 国产免费现黄频在线看| 无限看片的www在线观看| 69av精品久久久久久| 欧美不卡视频在线免费观看 | 国产成人一区二区三区免费视频网站| 久久久国产欧美日韩av| 制服人妻中文乱码| 国产成人精品无人区| 久久久久亚洲av毛片大全| 精品第一国产精品| 日韩av在线大香蕉| 欧美日本亚洲视频在线播放| 日日夜夜操网爽| 欧美日韩乱码在线| 亚洲精品成人av观看孕妇| 一a级毛片在线观看| 亚洲国产中文字幕在线视频| 久久精品国产清高在天天线|