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

    基于Landsat-8數據的洞庭湖區(qū)地表水體提取方法評價

    2023-08-22 07:53:20葉雨陽呂獻林羅鍇澍操華雙
    水利水電快報 2023年8期
    關鍵詞:豐水期水體精度

    葉雨陽,呂獻林,羅鍇澍,操華雙,陳 剛

    (1.中國地質大學 海洋學院,湖北 武漢 430074; 2.中國電建集團河南省電力勘測設計院有限公司,河南 鄭州 450007; 3.湖北大學 資源環(huán)境學院,湖北 武漢 430062; 4.竹山縣水利和湖泊局,湖北 十堰 442200)

    0 引 言

    湖泊作為陸地水圈的重要組成部分,參與自然界的水分循環(huán),使得地球表層系統(tǒng)各圈層相互聯(lián)結。同時湖泊也是環(huán)境變化重要的指示器,近年來由于氣候變化以及人類活動的影響[1],洞庭湖趨于萎縮,洪澇干旱災害頻繁發(fā)生,水環(huán)境形勢不容樂觀,進而影響洞庭湖區(qū)及周邊生態(tài)環(huán)境。因此,精準、快速識別湖泊水體信息并掌握其水體時空分布特征對于水資源監(jiān)測與應用具有重要意義[2],且可為有關部門制定保護政策從而改善湖泊周邊生態(tài)環(huán)境提供可靠依據。

    衛(wèi)星遙感技術因具有宏觀、綜合、動態(tài)、快速等無可比擬的優(yōu)點[3],已經成為分析和傳遞不同地球資源特別是地表水變化數據的重要信息來源。目前,Landsat系列衛(wèi)星數字產品因其較高的空間分辨率、覆蓋范圍全球化及易獲取性被廣泛應用于水體信息提取。利用遙感數據提取水體信息的方法一般大致分為兩類:一類為光譜分析方法,即考慮水體在多光譜波段的輻射特性提取水體信息。例如,Xu[4]基于各地物在不同波段光譜特性的差異提出改進的歸一化差異水體指數(Modified Normalized Difference Water Index,MNDWI),可以有效分割建筑區(qū)水體;Feyisa等[5]提出的自動水體提取指數(Automated Water Extraction Index,AWEIsh)能有效抑制陰影以正確識別水體特征。另一類為影像分類方法,即采用特征提取和分類技術進行水域檢測和定量分析。例如,Paul等[6]的研究表明,即使僅用少量的標記樣本進行訓練,支持向量機(Support Vector Machine,SVM)分類方法也能較好地從遙感影像中提取水體信息;王雪等[7]構建的全卷積神經網絡(Fully Convolution Neural Networks,FCN)模型可以從抽象特征中恢復出每個像素所屬類別來提取水體;王新宇[8]利用隨機森林算法(Random Forest,RF)挖掘特征組合信息,提升遙感影像水體范圍的提取效果。

    本文對比分析了MNDWI,AWEIsh,SVM,ANNs和RF等5種典型方法在洞庭湖水體提取中的表現(xiàn)。以Landsat-8枯水期影像數據為基礎,豐水期影像數據為輔,將通過目視解譯的水體信息作為真值,據此建立混淆矩陣并評價各方法提取精度。通過對比枯水期5種方法提取結果細節(jié)差異及枯、豐水期水體提取精度等,為各提取方法在不同應用及制圖需求中的靈活運用及基于中等空間分辨率遙感數據的湖泊動態(tài)監(jiān)測精準化研究等方面提供一定參考[9]。

    1 研究區(qū)概況及數據來源

    1.1 研究區(qū)概況

    洞庭湖區(qū)位于長江中游(27°39′N~29°51′N,111°19′E~113°34′E;如圖1所示),是中國第二大淡水湖,大型通江湖泊,具有防旱蓄洪、調節(jié)湖區(qū)周邊濕地生態(tài)環(huán)境等重要作用[10]。東洞庭湖向南連接湘江,向西連通資水、沅江和澧水,經調蓄后由城陵磯匯入長江[1],形成錯綜復雜的吞吐調蓄性湖泊。其東、南、西三面環(huán)山,呈河網平原地貌景觀,屬于典型的亞熱帶季風濕潤氣候。1~3月為枯水期,7~9月為豐水期,枯、豐水期水域分布變化極為強烈,素有“洪水一大片,枯水幾條線”之景。鑒于研究區(qū)基礎影像是洞庭湖區(qū)枯水期影像且湖區(qū)內水系面域較大[11],本文選取背景地物復雜、水體類型多樣的東洞庭湖區(qū)及水系分布復雜的湖區(qū),共2塊具有代表性的區(qū)域,見圖1(b)、(c),作精度評定。

    注:黑色方框為精度驗證區(qū)所在位置,水體邊界以黃色線條描繪。圖1 研究區(qū)示意

    1.2 數據來源及預處理

    本文以Landsat-8為數據源,枯水期影像數據為基礎、豐水期影像數據為輔,其中枯水期選取成像時間為2020年3月19日和2020年4月29日各一景影像,軌道號為124/40和123/40;豐水期選取成像時間為2020年10月22日的兩景影像,軌道號為123/39和123/40。所選影像皆為云量低于5%的高質量影像。

    因獲取影像存在地物間對比度低,受大氣散射、吸收、反射影響等問題,本文進行了如下預處理操作。首先使用ENVI5.3軟件對影像數據進行線性拉伸,突出顯示水體信息[12];其次進行輻射定標及FLAASH大氣校正;然后使用三次卷積法進行無縫鑲嵌;最后使用研究區(qū)矢量圖裁剪影像,以獲得研究區(qū)范圍枯、豐水期數據。

    2 研究方法

    為全面對比分析湖泊水體信息提取方法在洞庭湖的應用效果,本文從水體指數法和機器學習法中分別選取多種代表性分類方法,通過總體精度、Kappa系數、漏提率及過提率等4種評價指標進行評價分析。

    2.1 水體指數法

    水體指數法因其提取速度快、可擴展性強等優(yōu)點被廣泛用于水體信息提取。本文選取改進的歸一化差異水體指數(Modified Normalized Difference Water Index,MNDWI)和自動水體提取指數(Automated Water Extraction Index,AWEIsh)兩種方法作為典型代表用于后文對比分析。通過研究發(fā)現(xiàn),閾值的大小會受亞像素土地覆蓋組成成分占比影響,當土地覆蓋由水、土壤和植被組成時,水分占比越高,選取的最優(yōu)分割閾值越接近理論分割閾值0[13]。因此,本文通過對比分析水體、植被、農用地等地物灰度值,將理論分割閾值0進行調整,得到最優(yōu)分割閾值。且通過試驗發(fā)現(xiàn),水體指數法提取時間皆小于30 s。

    2.1.1 MNDWI

    Xu[4]改進了歸一化差異水體指數(Normalized Difference Water Index,NDWI)[14],將近紅外波段(NIR)替換為短波紅外波段(SWIR1),提出MNDWI,即

    (1)

    式中:Green(0.525~0.600 μm)為綠光波段;SWIR1(1.560~1.660 μm)為短波紅外1波段。由于水體與建設用地等地物反射特性在NIR波段相似但在SWIR1波段差異較大,MNDWI能更有效地抑制甚至去除建成區(qū)、植被及土壤噪聲。根據單波段影像灰度直方圖(圖2)選取閾值發(fā)現(xiàn),當閾值設定為0.464 2和0.459 4時,水體信息提取效果良好且受背景地物影響較小。因此,本文將大于0.464 2的像元劃分為枯水期水體,將大于0.459 4的像元劃分為豐水期水體。

    圖2 MNDWI法枯、豐水期水體灰度直方圖

    2.1.2 AWEIsh

    Feyisa等[5]利用Landsat-5 TM數據提出了AWEIsh指數,表達式為

    AWEIsh=Blue+2.5×Green-

    1.5×(NIR+SWIR1)-0.25×SWIR2

    (2)

    式中:Blue(0.450~0.515 μm)為藍光波段;SWIR2(2.100~2.300 μm)為短波紅外2波段;下標“sh”表明旨在有效消除陰影像元,提高存在山體陰影和其他暗面區(qū)域的水體提取精度。本文根據單波段影像灰度直方圖(圖3)選取閾值,最終選定0.064 8和0.093 0作為最優(yōu)分割閾值,將大于0.064 8的像元劃分為枯水期水體,大于0.093 0的像元劃分為豐水期水體。

    圖3 AWEIsh法枯、豐水期水體灰度直方圖

    2.2 機器學習法

    機器學習分為傳統(tǒng)機器學習和深度學習,近年來在水體信息提取方面的應用受到廣泛關注,該方法可以提取出更精細的水體細部信息。本文選取傳統(tǒng)機器學習支持向量機(Support Vector Machine,SVM)、隨機森林(Random Forest,RF)以及深度學習人工神經網絡(Artificial Neural Networks,ANNs)共3種方法用于后文對比分析,且均使用5,4,3波段(NIR,Red,Green)作為輸入波段。

    2.2.1 SVM

    SVM是以統(tǒng)計學為基礎,實現(xiàn)結構風險最小化原則的智能算法[15]。在水體提取中,由于SVM能有效解決影像混合像元等問題而被廣泛應用。其性能主要取決于核函數類型選擇和參數設置[16]。徑向基(RBF)核函數相較于線性核函數、多項式核函數及Sigmoid核函數識別率更高,性能更好,且訓練集減少時分類性能最穩(wěn)定[11],因此本文選取RBF作核函數。設置不同的Gamma值相當于調整模型復雜度,當Gamma值過大時,模型過擬合;Gamma值過小時,模型欠擬合。本文測試了Gamma值為[1,20]時的提取效果,最終將Gamma值設置為10,懲罰參數C為100。訓練時間為10 min。

    2.2.2 RF

    RF由Breiman[17]于2001年提出,是一種基于分類樹的機器學習算法[18],其思路是改進早前提出的Bagging算法[19]在節(jié)點特征選擇部分引入了隨機過程,它提高了簡單決策樹的預測精度。本文使用基尼不純度確定節(jié)點最優(yōu)條件,測試樹數量為[100,150]時的提取效果,通過對比發(fā)現(xiàn)當子樹數量為130時提取效果最佳,訓練時間為3 min。

    2.2.3 ANNs

    人工神經網絡簡稱為“神經網絡”(Neural Networks,NNs),本文采用的神經網絡由3層組成:輸入層、隱藏層和輸出層。輸入層負責數據接收和類型轉化等操作,隱含層將抽象層處理后的數據進行數據屬性抽象描述[20],輸出層負責將特征映射到特定維度并輸出預測結果。本文使用209個訓練樣本,146個測試樣本,Logistic函數作為激活函數[21],測試隱含層數量為[1,2],學習率為[0.2,0.5]時的提取效果,最終將隱含層設置為1層,學習率為0.2,迭代次數為1 000次,訓練時間為7 min。

    2.3 精度評價指標

    本文基于混淆矩陣[22]計算總體分類精度、Kappa系數、漏提率和過提率等4個指標進行提取結果精度評價。

    2.3.1 總體分類精度(Overall Accuracy,OA)

    OA表示被正確分為水體的像元數占像元總數的百分比,其表達式為

    (3)

    式中:TP表示實際為水體像元且檢測結果也為水體像元的像元數;TN表示實際為非水體像元且檢測結果也為非水體像元的像元數;FP表示實際為非水體像元而檢測結果為水體像元的像元數;FN表示實際為水體像元而檢測結果為非水體像元的像元數。

    2.3.2 Kappa系數

    Kappa系數是衡量一致性的常用統(tǒng)計方法,其表達式為

    (4)

    (5)

    式中:P0為總體分類精度;Pe為期望一致率,即兩次提取結果由于偶然機會所造成的一致率;n為總像元數。

    2.3.3 漏提率及誤提率

    本文采用“漏提率”和“過提率”對水體提取效果差異進行量化描述,公式為

    (6)

    (7)

    2.4 驗證數據

    為驗證各提取方法精度,本文通過如下方法提取驗證數據集,制作流程如圖4所示。

    圖4 驗證數據集制作流程

    (1) 分別選取枯、豐水期影像為底圖,在ArcGIS軟件平臺通過人工目視解譯獲得精度驗證區(qū)水體真值;

    (2) 在人工目視解譯結果的基礎上,以枯、豐水期水體分布作為約束條件,創(chuàng)建隨機水體樣本點[23];

    (3) 同理,生成非水體樣本點。

    最后,本文共收集了6 512個水體樣本點,6 512個非水體樣本點。

    3 結果與分析

    本研究通過分析水體提取精度及影響因素揭示各提取方法魯棒性并針對提取精度、提取時長及適用范圍對各方法作出評價。

    3.1 精度評價

    本文使用MNDWI,AWEIsh,SVM,ANNs和RF法對研究區(qū)水體進行提取,5種方法提取結果總體精度及漏提、過提等情況如表1所示;精度驗證區(qū)如圖1(b),(c)所示??梢钥闯?區(qū)域b中,SVM法漏提率最低,總體精度最高。區(qū)域c中,雖然AWEIsh法過提率最高但總體精度也最高,這是由于它對湖體的提取較其他方法更加完整,即漏提率最低;SVM法過提率最低但漏提率高于AWEIsh法和RF法,因此總體精度略低于這兩種方法。2個區(qū)域中,RF法過提率雖相對AWEIsh法更低,漏提率卻更高,因此總體精度略低于AWEIsh法。由此可見,漏提率指標在總體精度評定中起到更重要的參考作用。

    3.2 水體提取精度影響因素分析

    3.2.1 水體類型及背景地物影響

    本文將選取精度驗證區(qū)b,c中的差異水體類型及背景地物共4個重點區(qū)域,通過視覺比較和定量指標進一步評估5種方法的水體提取性能并分析影響提取精度的原因。其中,因本研究提取對象為洞庭湖區(qū)內由湖泊、水庫及細小河流組成的水系,魚場雖為水體但會對提取結果造成一定影響,因此將其視為背景地物。具體分類結果如圖5~6所示。

    圖5 5種算法在水庫、細小河流的提取效果差異

    通過與原始影像數據目視比較,提取效果差異主要表現(xiàn)為水庫、細小河流的漏提現(xiàn)象,農用地、魚場的過提現(xiàn)象。其中,MNDWI法、AWEIsh法對水庫漏提較嚴重,SVM法、AWEIsh法對農用地及魚場過提較嚴重,RF法最輕且對湖體間淺水處水體提取效果最好。為了定量評價水體分類的準確性及差異性,本文計算了水庫、細小河流、農用地及魚場的漏提率或過提率,結果如表2所示。

    表2 5種方法在4個區(qū)域的提取效果評價

    從表2中可知,MNDWI法對水庫漏提最嚴重,其原因可能是水庫與背景地物通過波段計算后DN值相似。AWEIsh法對細小河流的提取效果最好,但對背景地物過提嚴重,通過分析發(fā)現(xiàn)農用地、魚場等背景地物的光譜特性與部分水體相近容易造成過提。SVM法提取的水庫最貼近真實形態(tài),且通過對比效果圖發(fā)現(xiàn),該方法提取水體細部信息的效果較ANNs法更好。RF法對背景環(huán)境敏感性最低,即過提率最低??傮w來說,水體指數法對細小水體提取效果較好,機器學習法對于水體輪廓復雜度適應性良好,對背景環(huán)境敏感度基本低于水體指數方法,魯棒性更強。且綜合各水體類型及背景地物的漏提率、過提率排序來看,可與表1中結果相互印證。

    3.2.2 水域季節(jié)性變化影響

    由于洞庭湖區(qū)水系水體分布隨季節(jié)變化巨大,分為枯水期與豐水期(圖7中淡藍色部分為豐水期水體分布,深藍色虛線為枯水期水體分布范圍),本文通過計算4種定量評價指標分析水域變化對各方法提取精度的影響程度,結果見表3。

    通過對比發(fā)現(xiàn),水域面積季節(jié)性變化對各方法漏提率及過提率影響較大但對總體精度及Kappa系數影響不大。除ANNs法豐水期提取精度外,各方法枯、豐水期總體精度高的漏提率較低,且過提率相近時漏提率越高,總體精度越低。分析豐水期過提率高于枯水期的原因,發(fā)現(xiàn)豐水期水體分布更加復雜且裸露河床及淺水區(qū)較多易過提。各方法中,SVM法枯、豐水期總體提取精度均最高(95.02%,95.76%),表明其在應對湖泊水域面積變化時具有更好的魯棒性。精度表現(xiàn)最差的是MNDWI法(94.61%,93.06%),該方法枯、豐水期平均總體精度最低,且兩者差異最大,這表明指數型方法在區(qū)域化應用時泛化能力較差。因此綜合來看,表現(xiàn)更好的是機器學習法,其中SVM法效果最好,而水體指數法表現(xiàn)略遜一籌。

    3.3 提取方法評價

    通過前文對5種水體提取方法的綜合分析得知,MNDWI法和AWEIsh法提取結果受光譜特性影響較大,因此,在有大面積農用地等背景地物時不建議選擇水體指數法,但該種方法操作簡單、提取速度快,可用于自動化應急監(jiān)測。RF與ANNs相比:從提取水體類型來看,RF法適于提取細小河流,ANNs法適于提取面域較大的水體;從提取時長、背景環(huán)境影響來看,RF法速度更快且對背景環(huán)境敏感性更低,適用于背景地物復雜地區(qū)水域動態(tài)監(jiān)測。SVM法能較好地適應水域面積變化并能保證輪廓復雜水體的提取精度要求,適用于提取精度要求較高、水體類型多樣時的水體提取并制作地表水體分布圖。但該方法提取速度受樣本數量影響較大,因此,選取有效訓練樣本至關重要,這可以提高正確識別率及分類速度。

    4 結論與展望

    本文基于Landsat-8影像數據采用5種方法提取洞庭湖區(qū)枯、豐水期水體信息,得到以下結論:

    (1) 漏提率指標在精度評價中具有更高參考價值。

    (2) 基于Landsat-8影像,洞庭湖區(qū)水體信息提取精度及魯棒性排序為SVM>RF>AWEIsh>ANNs>MNDWI;提取時長排序為SVM>ANNs>RF>水體指數法。

    (3) 機器學習方法中,SVM法提取精度最高、魯棒性強,適用于多尺度、多類型水體精細化提取,能夠基本滿足高精度地表水體制圖需求;RF法對背景地物敏感度最低,適用于背景地物復雜地區(qū)的水體提取;ANNs法較適用于大面積水域提取。水體指數法受光譜特征影響較大,適用于城區(qū)內小型湖泊提取。

    本文研究僅選取同一衛(wèi)星兩個不同時期的洞庭湖水系進行分析,下一步將選取Landsat系列衛(wèi)星,使用多時相、長時間序列遙感影像數據分析洞庭湖水系面積季節(jié)變化及年際變化,為洞庭湖資源合理利用規(guī)劃提供科學依據。

    猜你喜歡
    豐水期水體精度
    豐水期巢湖水體與沉積物重金屬污染及分配特征研究
    陽谷縣農村地區(qū)枯水期和豐水期生活飲用水微生物檢測結果分析
    膠東沿海某地區(qū)生活飲用水總α、總β放射性水平調查分析
    農村黑臭水體治理和污水處理淺探
    生態(tài)修復理念在河道水體治理中的應用
    基于DSPIC33F微處理器的采集精度的提高
    電子制作(2018年11期)2018-08-04 03:25:38
    廣元:治理黑臭水體 再還水清岸美
    GPS/GLONASS/BDS組合PPP精度分析
    滇池流域入湖河流水文形貌特征對豐水期大型底棲動物群落結構的影響
    改進的Goldschmidt雙精度浮點除法器
    国产精品精品国产色婷婷| 特级一级黄色大片| 观看免费一级毛片| 色5月婷婷丁香| 国产精品乱码一区二三区的特点| 国产高清不卡午夜福利| 国产三级中文精品| 热99在线观看视频| 亚洲精品日韩在线中文字幕 | 又粗又硬又长又爽又黄的视频 | 亚洲精品国产成人久久av| 美女xxoo啪啪120秒动态图| 69人妻影院| 欧美xxxx黑人xx丫x性爽| 国产精品久久视频播放| 国产黄色视频一区二区在线观看 | 中文字幕熟女人妻在线| 99riav亚洲国产免费| 一级毛片电影观看 | 男女下面进入的视频免费午夜| 99久久精品热视频| 麻豆国产97在线/欧美| 亚洲av第一区精品v没综合| 中国美女看黄片| 能在线免费观看的黄片| 免费黄网站久久成人精品| 亚洲人成网站高清观看| 老师上课跳d突然被开到最大视频| 精品久久久久久久人妻蜜臀av| 日本熟妇午夜| 国产精品一区二区性色av| 国产精品久久久久久亚洲av鲁大| 精华霜和精华液先用哪个| 久久久久久久久大av| 日韩三级伦理在线观看| 一级毛片电影观看 | 高清午夜精品一区二区三区 | 边亲边吃奶的免费视频| 人妻久久中文字幕网| 我的女老师完整版在线观看| 亚洲国产欧洲综合997久久,| 亚洲激情五月婷婷啪啪| 91av网一区二区| 精品一区二区免费观看| 国产一区二区在线av高清观看| av又黄又爽大尺度在线免费看 | 在线播放国产精品三级| 成人午夜精彩视频在线观看| 国产成人精品久久久久久| 一区福利在线观看| 精品久久久久久久末码| 日本撒尿小便嘘嘘汇集6| 久久人妻av系列| 夫妻性生交免费视频一级片| 国产私拍福利视频在线观看| 亚洲欧美清纯卡通| 国产 一区 欧美 日韩| 波多野结衣巨乳人妻| 亚洲欧美成人精品一区二区| 日韩 亚洲 欧美在线| 国产在视频线在精品| 婷婷精品国产亚洲av| 亚洲精品日韩av片在线观看| 色播亚洲综合网| 嘟嘟电影网在线观看| 国产成人午夜福利电影在线观看| 国产三级在线视频| 国产黄片视频在线免费观看| 成年av动漫网址| 欧美激情久久久久久爽电影| 婷婷六月久久综合丁香| 午夜免费激情av| 国产av一区在线观看免费| 精品久久久久久久人妻蜜臀av| 久久鲁丝午夜福利片| 国产精品一区二区在线观看99 | 国内精品美女久久久久久| 寂寞人妻少妇视频99o| 人人妻人人看人人澡| 国产 一区 欧美 日韩| 成人午夜精彩视频在线观看| 白带黄色成豆腐渣| 欧美一区二区精品小视频在线| 成人三级黄色视频| 国产精品爽爽va在线观看网站| 99久久精品国产国产毛片| 嘟嘟电影网在线观看| 久久久久久久亚洲中文字幕| 国产亚洲精品av在线| 国产精品人妻久久久久久| 国产综合懂色| 给我免费播放毛片高清在线观看| 中文字幕免费在线视频6| 国产高清三级在线| 免费看光身美女| 亚洲人成网站在线播放欧美日韩| 亚洲精品自拍成人| 日韩三级伦理在线观看| 99国产精品一区二区蜜桃av| a级毛色黄片| 高清在线视频一区二区三区 | 一进一出抽搐gif免费好疼| 亚洲精品日韩av片在线观看| 亚洲av中文av极速乱| 少妇人妻一区二区三区视频| 99久久中文字幕三级久久日本| 中国国产av一级| 国产成人影院久久av| 久久精品久久久久久久性| 中文在线观看免费www的网站| 在线国产一区二区在线| 成人毛片a级毛片在线播放| 国产高清激情床上av| 麻豆国产97在线/欧美| 欧美在线一区亚洲| 麻豆av噜噜一区二区三区| 国产综合懂色| 成人亚洲欧美一区二区av| av天堂在线播放| 美女国产视频在线观看| 午夜亚洲福利在线播放| 麻豆精品久久久久久蜜桃| 久久久a久久爽久久v久久| 日韩欧美在线乱码| 91久久精品电影网| 高清日韩中文字幕在线| 91久久精品国产一区二区成人| 成人性生交大片免费视频hd| 免费看日本二区| 大又大粗又爽又黄少妇毛片口| 亚洲最大成人中文| 神马国产精品三级电影在线观看| 中文精品一卡2卡3卡4更新| 欧美zozozo另类| a级毛色黄片| 黄色欧美视频在线观看| 尾随美女入室| 日本免费一区二区三区高清不卡| 亚洲综合色惰| 日本av手机在线免费观看| 久久午夜亚洲精品久久| 中文字幕制服av| 一级毛片aaaaaa免费看小| 大香蕉久久网| 日本成人三级电影网站| 国产精品无大码| 国产精品三级大全| 少妇猛男粗大的猛烈进出视频 | 自拍偷自拍亚洲精品老妇| 色哟哟哟哟哟哟| 国产三级在线视频| 欧美一区二区亚洲| 成人无遮挡网站| 尤物成人国产欧美一区二区三区| 看十八女毛片水多多多| 99久久中文字幕三级久久日本| 春色校园在线视频观看| 免费观看精品视频网站| 黄片无遮挡物在线观看| 日韩av在线大香蕉| 成人欧美大片| 一区福利在线观看| 可以在线观看的亚洲视频| 欧美日本视频| 国产精品,欧美在线| 日韩 亚洲 欧美在线| 国产精品99久久久久久久久| 国产一区二区三区av在线 | h日本视频在线播放| 美女国产视频在线观看| 国产美女午夜福利| 免费无遮挡裸体视频| 欧美bdsm另类| 乱系列少妇在线播放| 久久亚洲国产成人精品v| 亚洲成人精品中文字幕电影| 一本久久中文字幕| 搞女人的毛片| 久久人人爽人人爽人人片va| 中出人妻视频一区二区| 欧美精品一区二区大全| 卡戴珊不雅视频在线播放| 丝袜喷水一区| 亚洲成人精品中文字幕电影| 国产精品嫩草影院av在线观看| 18+在线观看网站| eeuss影院久久| 日韩高清综合在线| 欧美人与善性xxx| 国产精品久久久久久精品电影| 亚洲最大成人手机在线| 内射极品少妇av片p| 99热网站在线观看| 韩国av在线不卡| 老熟妇乱子伦视频在线观看| 哪里可以看免费的av片| 免费看日本二区| 国产精品久久久久久久电影| 国产黄色小视频在线观看| av黄色大香蕉| 国产真实伦视频高清在线观看| 午夜福利在线观看吧| 国产高清三级在线| 国内揄拍国产精品人妻在线| 非洲黑人性xxxx精品又粗又长| 亚洲国产精品合色在线| 综合色av麻豆| 午夜精品一区二区三区免费看| 一区二区三区高清视频在线| 亚洲婷婷狠狠爱综合网| 欧美+日韩+精品| 我的女老师完整版在线观看| 久久精品国产亚洲网站| 国内精品美女久久久久久| 偷拍熟女少妇极品色| 欧美高清性xxxxhd video| 在线a可以看的网站| 精品久久久久久久久av| 一卡2卡三卡四卡精品乱码亚洲| 欧美日韩一区二区视频在线观看视频在线 | 男女做爰动态图高潮gif福利片| 久久久色成人| 婷婷精品国产亚洲av| 丰满乱子伦码专区| 国产色爽女视频免费观看| 久久99精品国语久久久| 久久午夜亚洲精品久久| 国产精品国产高清国产av| 国产亚洲91精品色在线| 亚洲美女搞黄在线观看| 国产精品福利在线免费观看| 性插视频无遮挡在线免费观看| 可以在线观看毛片的网站| 欧美日韩精品成人综合77777| 99久国产av精品国产电影| 伊人久久精品亚洲午夜| 淫秽高清视频在线观看| 亚洲在线自拍视频| 久久精品国产亚洲av香蕉五月| 能在线免费观看的黄片| 久久精品91蜜桃| 超碰av人人做人人爽久久| 真实男女啪啪啪动态图| 欧美日本亚洲视频在线播放| 麻豆成人午夜福利视频| 女同久久另类99精品国产91| 成人午夜高清在线视频| 国产成人一区二区在线| 色综合站精品国产| 国产精品麻豆人妻色哟哟久久 | 久久精品人妻少妇| 日本一二三区视频观看| 波野结衣二区三区在线| 日本av手机在线免费观看| 国产成人精品一,二区 | 国产精品野战在线观看| 亚州av有码| 成人二区视频| 亚洲精品自拍成人| 毛片女人毛片| av天堂中文字幕网| 日本免费a在线| 亚洲欧美日韩高清在线视频| www.色视频.com| 久久精品国产亚洲网站| 在线观看av片永久免费下载| 毛片女人毛片| 熟女电影av网| 丝袜美腿在线中文| 一个人看视频在线观看www免费| 亚洲18禁久久av| 欧美另类亚洲清纯唯美| 亚洲激情五月婷婷啪啪| 免费观看的影片在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 一边亲一边摸免费视频| 国产精品久久电影中文字幕| 亚洲欧美日韩无卡精品| 少妇人妻精品综合一区二区 | 在线观看美女被高潮喷水网站| 国产日韩欧美在线精品| 精品一区二区免费观看| 亚洲精品乱码久久久v下载方式| 久久热精品热| 久久久国产成人免费| 网址你懂的国产日韩在线| 我要看日韩黄色一级片| 国产激情偷乱视频一区二区| av又黄又爽大尺度在线免费看 | 麻豆一二三区av精品| av免费在线看不卡| 国内久久婷婷六月综合欲色啪| 2021天堂中文幕一二区在线观| 久久午夜亚洲精品久久| 久久精品91蜜桃| 欧美成人a在线观看| 亚洲无线观看免费| 九九热线精品视视频播放| 久久久a久久爽久久v久久| 精华霜和精华液先用哪个| 亚洲人与动物交配视频| 日韩强制内射视频| 国产精品野战在线观看| 亚洲精品粉嫩美女一区| 国内精品久久久久精免费| 国产成人91sexporn| 久久99热6这里只有精品| 99久久久亚洲精品蜜臀av| 最后的刺客免费高清国语| 亚洲,欧美,日韩| 国产一区二区亚洲精品在线观看| 国产探花在线观看一区二区| 深夜精品福利| 国产色爽女视频免费观看| 国产高清激情床上av| 91精品国产九色| 日日啪夜夜撸| 久久精品夜色国产| 久久欧美精品欧美久久欧美| 乱系列少妇在线播放| 亚洲第一区二区三区不卡| 日韩av在线大香蕉| 免费无遮挡裸体视频| 亚洲av第一区精品v没综合| 亚洲人成网站在线观看播放| 精品久久国产蜜桃| 看片在线看免费视频| 国产成人精品婷婷| 国内精品久久久久精免费| 免费电影在线观看免费观看| 精品人妻一区二区三区麻豆| 久久久成人免费电影| 老司机福利观看| 国产精品福利在线免费观看| 免费看av在线观看网站| 黄片无遮挡物在线观看| 免费看日本二区| 亚洲乱码一区二区免费版| 欧美xxxx性猛交bbbb| 国产免费一级a男人的天堂| 男人的好看免费观看在线视频| 三级经典国产精品| 免费人成视频x8x8入口观看| 亚洲av二区三区四区| 日日摸夜夜添夜夜爱| 欧美+日韩+精品| 欧美性感艳星| 国产伦一二天堂av在线观看| 免费看美女性在线毛片视频| 国产精品99久久久久久久久| 国产高清有码在线观看视频| 九草在线视频观看| 欧美成人a在线观看| 少妇高潮的动态图| 91久久精品国产一区二区成人| 成人永久免费在线观看视频| 综合色丁香网| 免费看光身美女| 国产黄色小视频在线观看| 成人永久免费在线观看视频| 狠狠狠狠99中文字幕| 亚洲,欧美,日韩| 免费大片18禁| 内射极品少妇av片p| 午夜亚洲福利在线播放| 亚洲欧美日韩东京热| 免费看美女性在线毛片视频| 国产精品乱码一区二三区的特点| 国产高清有码在线观看视频| 网址你懂的国产日韩在线| 夜夜爽天天搞| av天堂中文字幕网| 国产亚洲精品久久久com| 高清午夜精品一区二区三区 | 大型黄色视频在线免费观看| 人妻系列 视频| 神马国产精品三级电影在线观看| 久久久久性生活片| 日本与韩国留学比较| 亚洲成a人片在线一区二区| 国产精品久久久久久亚洲av鲁大| 国产日本99.免费观看| 男女啪啪激烈高潮av片| 少妇熟女欧美另类| 亚洲欧美精品专区久久| 午夜爱爱视频在线播放| 欧美+亚洲+日韩+国产| 秋霞在线观看毛片| 久久99蜜桃精品久久| 99热这里只有是精品在线观看| 99久久久亚洲精品蜜臀av| 成人性生交大片免费视频hd| 久久精品人妻少妇| 少妇人妻精品综合一区二区 | 久久婷婷人人爽人人干人人爱| 午夜精品在线福利| 亚洲精华国产精华液的使用体验 | 韩国av在线不卡| 国产午夜精品论理片| 精品人妻视频免费看| 国产精华一区二区三区| 国产激情偷乱视频一区二区| 久久精品国产自在天天线| 亚洲四区av| av在线亚洲专区| 美女国产视频在线观看| 人妻系列 视频| 中文亚洲av片在线观看爽| 嫩草影院精品99| 观看免费一级毛片| 精品人妻视频免费看| 日本欧美国产在线视频| 三级经典国产精品| 国产成人影院久久av| 亚洲电影在线观看av| 国产淫片久久久久久久久| 久久精品人妻少妇| 久久精品国产99精品国产亚洲性色| 久久国产乱子免费精品| 日本三级黄在线观看| www.色视频.com| 国产又黄又爽又无遮挡在线| 久久人人爽人人爽人人片va| 亚洲欧美成人精品一区二区| 18禁在线播放成人免费| 欧美激情久久久久久爽电影| 欧美成人一区二区免费高清观看| 不卡一级毛片| 51国产日韩欧美| av在线观看视频网站免费| 蜜桃亚洲精品一区二区三区| 草草在线视频免费看| 国产在视频线在精品| 五月玫瑰六月丁香| 久久久久久久久中文| 观看美女的网站| 九草在线视频观看| 午夜福利视频1000在线观看| 欧美一区二区精品小视频在线| 久久精品夜色国产| 婷婷色综合大香蕉| 高清在线视频一区二区三区 | 一级二级三级毛片免费看| 亚洲在久久综合| 日韩av在线大香蕉| 国产午夜福利久久久久久| av女优亚洲男人天堂| 2022亚洲国产成人精品| 欧美成人免费av一区二区三区| www.色视频.com| 精品久久久噜噜| 综合色av麻豆| 99久久九九国产精品国产免费| 18禁裸乳无遮挡免费网站照片| 亚洲欧美精品自产自拍| 人妻系列 视频| 在线观看一区二区三区| 丝袜喷水一区| 亚洲精品国产成人久久av| 不卡视频在线观看欧美| 22中文网久久字幕| 男人和女人高潮做爰伦理| 22中文网久久字幕| 欧美成人精品欧美一级黄| 久久婷婷人人爽人人干人人爱| 国产精品麻豆人妻色哟哟久久 | av免费在线看不卡| 看片在线看免费视频| 日韩成人av中文字幕在线观看| 黄色一级大片看看| 亚洲欧美清纯卡通| 欧美又色又爽又黄视频| 国产精品蜜桃在线观看 | 精品久久久久久成人av| 久久精品夜色国产| 久久人人精品亚洲av| 成人亚洲精品av一区二区| 免费av毛片视频| 黄色欧美视频在线观看| 黄色视频,在线免费观看| 国产伦在线观看视频一区| 日韩av在线大香蕉| 12—13女人毛片做爰片一| 日本黄色片子视频| 高清日韩中文字幕在线| 国产激情偷乱视频一区二区| 嫩草影院新地址| 在线观看66精品国产| 国产老妇伦熟女老妇高清| 午夜福利在线在线| 日本三级黄在线观看| 国产黄a三级三级三级人| 亚洲性久久影院| 欧美成人a在线观看| 1000部很黄的大片| 国产v大片淫在线免费观看| 色吧在线观看| 少妇的逼水好多| 岛国在线免费视频观看| 久久精品91蜜桃| 日韩一区二区三区影片| 热99在线观看视频| 欧美最黄视频在线播放免费| 69人妻影院| 高清午夜精品一区二区三区 | 国产高清激情床上av| 国产在线精品亚洲第一网站| 中国国产av一级| 黄色一级大片看看| 国产真实乱freesex| 精品人妻视频免费看| 十八禁国产超污无遮挡网站| 欧美一区二区亚洲| 亚洲国产高清在线一区二区三| 韩国av在线不卡| 久久久成人免费电影| 天堂av国产一区二区熟女人妻| 69人妻影院| 亚洲一区二区三区色噜噜| 久久精品国产99精品国产亚洲性色| av免费在线看不卡| 国产午夜精品论理片| 99riav亚洲国产免费| 在线观看av片永久免费下载| 中文字幕免费在线视频6| 国产高清视频在线观看网站| 伦精品一区二区三区| 亚洲欧美清纯卡通| 欧美日韩在线观看h| 精品久久久久久久久久久久久| 亚洲无线在线观看| 变态另类成人亚洲欧美熟女| 美女 人体艺术 gogo| 哪个播放器可以免费观看大片| 久久久精品大字幕| 嫩草影院精品99| 女同久久另类99精品国产91| 国产亚洲欧美98| 丰满人妻一区二区三区视频av| 97超碰精品成人国产| 国模一区二区三区四区视频| 国产精品久久久久久亚洲av鲁大| 尤物成人国产欧美一区二区三区| 听说在线观看完整版免费高清| 国产精品1区2区在线观看.| 乱系列少妇在线播放| 日韩精品有码人妻一区| 国产av麻豆久久久久久久| 欧美激情在线99| 国产极品精品免费视频能看的| 蜜桃亚洲精品一区二区三区| 日韩三级伦理在线观看| 国产伦一二天堂av在线观看| 高清毛片免费观看视频网站| 国产真实伦视频高清在线观看| 国产成人a区在线观看| 免费观看人在逋| 不卡一级毛片| 精品99又大又爽又粗少妇毛片| 久久人人爽人人爽人人片va| 伦理电影大哥的女人| 人妻少妇偷人精品九色| eeuss影院久久| 国产精品野战在线观看| 国产精品久久视频播放| 国产三级中文精品| 亚洲不卡免费看| 色综合站精品国产| 国产亚洲精品久久久久久毛片| 在线免费观看不下载黄p国产| 亚洲av男天堂| 国产精品不卡视频一区二区| 久久久久九九精品影院| 亚洲欧美日韩无卡精品| 中文字幕久久专区| 99久久无色码亚洲精品果冻| 91午夜精品亚洲一区二区三区| 九九爱精品视频在线观看| 久久久精品欧美日韩精品| av.在线天堂| 91精品国产九色| 男女边吃奶边做爰视频| 国产精品.久久久| 国产麻豆成人av免费视频| 久久99热这里只有精品18| 欧美日本视频| 男人和女人高潮做爰伦理| 一区福利在线观看| 在线播放无遮挡| 日韩亚洲欧美综合| 少妇被粗大猛烈的视频| 国产伦精品一区二区三区四那| 国产成人a∨麻豆精品| 人人妻人人澡欧美一区二区| 国产色爽女视频免费观看| av视频在线观看入口| 精品人妻偷拍中文字幕| 老女人水多毛片| 综合色丁香网| 亚洲成人久久性| 亚洲久久久久久中文字幕| 成人漫画全彩无遮挡| av在线蜜桃| 亚洲久久久久久中文字幕| 深爱激情五月婷婷| 不卡一级毛片| 国产精品永久免费网站| or卡值多少钱| 天堂中文最新版在线下载 | 国产成人精品久久久久久| 性欧美人与动物交配| 国产国拍精品亚洲av在线观看| 99久久九九国产精品国产免费| a级毛色黄片| 精品午夜福利在线看|