• <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雙精度浮點除法器
    日韩三级伦理在线观看| 日本爱情动作片www.在线观看| 精品少妇久久久久久888优播| 熟女人妻精品中文字幕| 国产淫语在线视频| 亚洲精品aⅴ在线观看| 草草在线视频免费看| 亚洲情色 制服丝袜| 亚洲欧美中文字幕日韩二区| 国产淫语在线视频| 日韩人妻精品一区2区三区| 国产极品粉嫩免费观看在线| 亚洲高清免费不卡视频| 香蕉精品网在线| 久久久欧美国产精品| 亚洲一码二码三码区别大吗| av免费在线看不卡| 咕卡用的链子| 国产精品麻豆人妻色哟哟久久| 国产男女超爽视频在线观看| 亚洲综合色惰| 女性被躁到高潮视频| 涩涩av久久男人的天堂| 中文字幕人妻丝袜制服| 女人精品久久久久毛片| 国语对白做爰xxxⅹ性视频网站| 国产爽快片一区二区三区| 亚洲精品美女久久久久99蜜臀 | 午夜老司机福利剧场| 亚洲精品国产av蜜桃| 精品人妻熟女毛片av久久网站| 欧美日韩成人在线一区二区| 观看美女的网站| 熟女人妻精品中文字幕| 妹子高潮喷水视频| 久久久久久久国产电影| 国产成人aa在线观看| 肉色欧美久久久久久久蜜桃| 在线观看免费视频网站a站| 国产精品99久久99久久久不卡 | 最近最新中文字幕免费大全7| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 18+在线观看网站| 午夜视频国产福利| 赤兔流量卡办理| 欧美激情国产日韩精品一区| 女性生殖器流出的白浆| 香蕉丝袜av| av网站免费在线观看视频| 亚洲av国产av综合av卡| 欧美精品一区二区免费开放| 精品国产一区二区久久| 久久久久视频综合| 免费观看性生交大片5| 大片免费播放器 马上看| 性色av一级| 国产伦理片在线播放av一区| 日本av免费视频播放| 天美传媒精品一区二区| 久久精品国产亚洲av天美| 伦理电影大哥的女人| 国产一区二区三区av在线| 欧美少妇被猛烈插入视频| 建设人人有责人人尽责人人享有的| 久久久久久久亚洲中文字幕| 男女啪啪激烈高潮av片| 天美传媒精品一区二区| 日日啪夜夜爽| 满18在线观看网站| 在线观看国产h片| 天堂俺去俺来也www色官网| 18禁国产床啪视频网站| 国产极品粉嫩免费观看在线| 色婷婷久久久亚洲欧美| 久久亚洲国产成人精品v| freevideosex欧美| 免费观看无遮挡的男女| 老女人水多毛片| 中文字幕另类日韩欧美亚洲嫩草| 一级a做视频免费观看| 国产精品一区www在线观看| 国产色婷婷99| 只有这里有精品99| 91精品伊人久久大香线蕉| 免费黄网站久久成人精品| 国产高清不卡午夜福利| 2018国产大陆天天弄谢| 欧美3d第一页| 久久精品国产亚洲av涩爱| 欧美日韩精品成人综合77777| 亚洲精品av麻豆狂野| 18禁观看日本| 免费观看在线日韩| 狂野欧美激情性xxxx在线观看| 欧美成人午夜精品| 亚洲国产日韩一区二区| 少妇的丰满在线观看| www日本在线高清视频| 亚洲三级黄色毛片| 黄片无遮挡物在线观看| 午夜91福利影院| 一个人免费看片子| 久久久亚洲精品成人影院| 国产日韩一区二区三区精品不卡| 内地一区二区视频在线| 亚洲国产精品一区二区三区在线| 黑人欧美特级aaaaaa片| av在线app专区| 天堂中文最新版在线下载| 久久久国产一区二区| 国产成人aa在线观看| 人妻一区二区av| 美女主播在线视频| 99久国产av精品国产电影| 免费人成在线观看视频色| 丝袜人妻中文字幕| 黄色毛片三级朝国网站| 欧美精品高潮呻吟av久久| 欧美精品一区二区大全| av在线播放精品| av黄色大香蕉| 久久99一区二区三区| 国国产精品蜜臀av免费| 欧美成人午夜免费资源| 男女边摸边吃奶| 黄片播放在线免费| 乱码一卡2卡4卡精品| 久久99一区二区三区| av天堂久久9| 欧美亚洲日本最大视频资源| 久久久国产精品麻豆| 十分钟在线观看高清视频www| 女性被躁到高潮视频| 97超碰精品成人国产| 国产精品一国产av| 国内精品宾馆在线| 中国美白少妇内射xxxbb| 欧美日韩精品成人综合77777| 一级毛片我不卡| 亚洲三级黄色毛片| 99热国产这里只有精品6| 精品卡一卡二卡四卡免费| 亚洲国产毛片av蜜桃av| 亚洲色图 男人天堂 中文字幕 | 亚洲性久久影院| 日韩精品有码人妻一区| 亚洲av免费高清在线观看| 考比视频在线观看| 如日韩欧美国产精品一区二区三区| 91精品伊人久久大香线蕉| 国产日韩欧美在线精品| 成人黄色视频免费在线看| av国产久精品久网站免费入址| 免费少妇av软件| 视频在线观看一区二区三区| 9191精品国产免费久久| 色5月婷婷丁香| 亚洲欧美成人综合另类久久久| 欧美人与性动交α欧美软件 | 国产成人免费无遮挡视频| 涩涩av久久男人的天堂| 如日韩欧美国产精品一区二区三区| 9191精品国产免费久久| 日韩成人伦理影院| 国产麻豆69| 亚洲国产精品成人久久小说| 伦理电影大哥的女人| 啦啦啦在线观看免费高清www| 99热这里只有是精品在线观看| 岛国毛片在线播放| 男女边吃奶边做爰视频| 久久久久久人人人人人| 色网站视频免费| 人人妻人人澡人人看| 国产av精品麻豆| 日日啪夜夜爽| 免费看不卡的av| av在线观看视频网站免费| 精品少妇久久久久久888优播| 久久精品国产亚洲av天美| 寂寞人妻少妇视频99o| 欧美国产精品一级二级三级| 激情视频va一区二区三区| 欧美激情极品国产一区二区三区 | 蜜桃国产av成人99| 夫妻性生交免费视频一级片| 国产一区亚洲一区在线观看| 人妻少妇偷人精品九色| 亚洲欧美成人精品一区二区| 婷婷色综合大香蕉| 午夜老司机福利剧场| 日韩一区二区视频免费看| 国产精品无大码| 99久久综合免费| 美女中出高潮动态图| 2022亚洲国产成人精品| 色哟哟·www| 免费日韩欧美在线观看| 国产精品久久久久久精品电影小说| 亚洲第一av免费看| 永久免费av网站大全| 国产黄频视频在线观看| 亚洲伊人久久精品综合| 91在线精品国自产拍蜜月| 新久久久久国产一级毛片| 少妇被粗大的猛进出69影院 | 国产男女内射视频| 婷婷色综合www| 97在线人人人人妻| 看十八女毛片水多多多| 制服诱惑二区| 国产成人av激情在线播放| 久久人人爽人人片av| 乱码一卡2卡4卡精品| 久久久久久久久久人人人人人人| 综合色丁香网| 两个人免费观看高清视频| 视频区图区小说| 欧美日韩av久久| 国产在线视频一区二区| 精品一区在线观看国产| 欧美97在线视频| 观看av在线不卡| 老司机亚洲免费影院| 日韩,欧美,国产一区二区三区| 少妇的逼好多水| 国产爽快片一区二区三区| 男女啪啪激烈高潮av片| 精品国产国语对白av| 蜜桃在线观看..| 久久久久久伊人网av| 久久精品久久久久久久性| 欧美日韩一区二区视频在线观看视频在线| 免费观看性生交大片5| 亚洲丝袜综合中文字幕| 成人国产av品久久久| 久久97久久精品| 国产高清三级在线| 精品久久久久久电影网| 丝袜在线中文字幕| 男女下面插进去视频免费观看 | 亚洲av日韩在线播放| 亚洲精品av麻豆狂野| av视频免费观看在线观看| 精品国产一区二区久久| 亚洲国产欧美日韩在线播放| 欧美日韩一区二区视频在线观看视频在线| 日韩精品有码人妻一区| 亚洲性久久影院| av福利片在线| 亚洲成人手机| 久久久国产欧美日韩av| 香蕉丝袜av| 一本—道久久a久久精品蜜桃钙片| 高清不卡的av网站| 天天躁夜夜躁狠狠躁躁| 熟女av电影| 国产精品无大码| 亚洲,一卡二卡三卡| 黄色 视频免费看| 亚洲成色77777| 菩萨蛮人人尽说江南好唐韦庄| 久久精品国产综合久久久 | 国产黄频视频在线观看| 高清不卡的av网站| av黄色大香蕉| 国产免费福利视频在线观看| 高清av免费在线| 日韩一区二区三区影片| 尾随美女入室| 婷婷色综合大香蕉| 成人无遮挡网站| 十八禁高潮呻吟视频| 欧美最新免费一区二区三区| 亚洲国产av新网站| 啦啦啦中文免费视频观看日本| 伊人久久国产一区二区| 在线看a的网站| av一本久久久久| 人妻少妇偷人精品九色| 少妇的逼水好多| 亚洲国产精品一区二区三区在线| av福利片在线| 夜夜骑夜夜射夜夜干| 中文字幕制服av| 看免费成人av毛片| 中文字幕人妻丝袜制服| 高清av免费在线| 三级国产精品片| 97在线视频观看| 99久久中文字幕三级久久日本| 久久精品久久精品一区二区三区| 中文字幕亚洲精品专区| 大话2 男鬼变身卡| 天美传媒精品一区二区| 免费黄频网站在线观看国产| 亚洲一区二区三区欧美精品| 看免费av毛片| 丁香六月天网| 午夜日本视频在线| 大片免费播放器 马上看| 一区二区日韩欧美中文字幕 | 国产极品天堂在线| 综合色丁香网| 国产精品蜜桃在线观看| 国产精品一区二区在线观看99| 亚洲一级一片aⅴ在线观看| av国产久精品久网站免费入址| 国产淫语在线视频| 亚洲色图 男人天堂 中文字幕 | 午夜影院在线不卡| 日韩一本色道免费dvd| 国产日韩一区二区三区精品不卡| 欧美日韩亚洲高清精品| 大陆偷拍与自拍| 国产精品女同一区二区软件| 新久久久久国产一级毛片| 两性夫妻黄色片 | 考比视频在线观看| 少妇人妻 视频| 国产 精品1| 成人亚洲精品一区在线观看| 亚洲欧洲日产国产| 99热全是精品| 久久久久久久久久人人人人人人| 有码 亚洲区| 欧美成人午夜精品| 亚洲久久久国产精品| 精品午夜福利在线看| 9热在线视频观看99| 日本黄色日本黄色录像| 18禁裸乳无遮挡动漫免费视频| 欧美成人午夜精品| 日韩熟女老妇一区二区性免费视频| av一本久久久久| 18禁动态无遮挡网站| 免费观看在线日韩| 两个人看的免费小视频| 韩国av在线不卡| 一级毛片我不卡| 91成人精品电影| 日韩av在线免费看完整版不卡| 日韩三级伦理在线观看| 成年人午夜在线观看视频| 久久精品aⅴ一区二区三区四区 | 国产成人精品久久久久久| 91久久精品国产一区二区三区| 欧美成人午夜精品| 狠狠婷婷综合久久久久久88av| 青青草视频在线视频观看| 97在线人人人人妻| 亚洲精品乱码久久久久久按摩| 日日爽夜夜爽网站| 成年女人在线观看亚洲视频| 精品国产一区二区久久| 老司机亚洲免费影院| 国产一区二区三区综合在线观看 | 九九爱精品视频在线观看| 欧美精品人与动牲交sv欧美| 免费大片18禁| 免费人成在线观看视频色| 中文字幕人妻熟女乱码| 国产黄色视频一区二区在线观看| 性高湖久久久久久久久免费观看| 爱豆传媒免费全集在线观看| 韩国精品一区二区三区 | 中国三级夫妇交换| 国产精品久久久久成人av| 黄色配什么色好看| 精品一区在线观看国产| 亚洲av在线观看美女高潮| 亚洲精品久久久久久婷婷小说| 国产免费视频播放在线视频| 宅男免费午夜| 高清av免费在线| 熟妇人妻不卡中文字幕| 亚洲av日韩在线播放| 视频中文字幕在线观看| 国产色婷婷99| 26uuu在线亚洲综合色| 美女内射精品一级片tv| 久热这里只有精品99| 最新中文字幕久久久久| 在线免费观看不下载黄p国产| 在线天堂最新版资源| 视频中文字幕在线观看| 精品国产一区二区三区四区第35| 26uuu在线亚洲综合色| 最近2019中文字幕mv第一页| 久久国产精品大桥未久av| 亚洲国产精品一区三区| 午夜影院在线不卡| 国产一区二区在线观看av| 精品一区二区三区视频在线| 性高湖久久久久久久久免费观看| 久久久精品94久久精品| 日本欧美视频一区| 视频在线观看一区二区三区| 少妇猛男粗大的猛烈进出视频| 熟女av电影| 久久久久视频综合| 三上悠亚av全集在线观看| 男女免费视频国产| 亚洲成人手机| 涩涩av久久男人的天堂| 午夜福利网站1000一区二区三区| 国产高清国产精品国产三级| 久久精品国产亚洲av天美| 国产有黄有色有爽视频| 男女边吃奶边做爰视频| 精品亚洲成国产av| 只有这里有精品99| videossex国产| 国产在线视频一区二区| 国产精品免费大片| 99久久人妻综合| 亚洲色图 男人天堂 中文字幕 | 丝袜喷水一区| 极品人妻少妇av视频| 中文字幕制服av| 国产精品国产三级专区第一集| 五月天丁香电影| 欧美激情国产日韩精品一区| 中文精品一卡2卡3卡4更新| av线在线观看网站| 久久久久精品久久久久真实原创| 天堂俺去俺来也www色官网| 亚洲精品视频女| 最近的中文字幕免费完整| 老司机影院毛片| 日本av手机在线免费观看| 男人添女人高潮全过程视频| 国产成人免费观看mmmm| 欧美日韩亚洲高清精品| 1024视频免费在线观看| 国产免费一区二区三区四区乱码| 亚洲欧美成人精品一区二区| 日日撸夜夜添| 国产精品蜜桃在线观看| 国产精品国产三级专区第一集| 日本欧美国产在线视频| 久久人妻熟女aⅴ| 亚洲婷婷狠狠爱综合网| 伦精品一区二区三区| 少妇的丰满在线观看| 国产一级毛片在线| 亚洲av中文av极速乱| 大香蕉97超碰在线| 亚洲国产日韩一区二区| 日韩成人av中文字幕在线观看| 国产精品国产av在线观看| 亚洲国产成人一精品久久久| 午夜日本视频在线| 久久 成人 亚洲| 国产一区有黄有色的免费视频| 国产一区二区在线观看av| 妹子高潮喷水视频| 视频在线观看一区二区三区| 日韩免费高清中文字幕av| 黄色配什么色好看| 免费av中文字幕在线| 国产xxxxx性猛交| 日本av免费视频播放| 99热6这里只有精品| 日韩人妻精品一区2区三区| 男人爽女人下面视频在线观看| 久久久久久久久久久免费av| 亚洲国产精品一区三区| 另类亚洲欧美激情| 亚洲国产精品专区欧美| 热re99久久精品国产66热6| 中文字幕人妻丝袜制服| 久久精品久久精品一区二区三区| 日韩视频在线欧美| 69精品国产乱码久久久| 性色av一级| 只有这里有精品99| 精品熟女少妇av免费看| 国产探花极品一区二区| 中文欧美无线码| 久久久精品区二区三区| 亚洲国产毛片av蜜桃av| 国产麻豆69| 国产一级毛片在线| 大香蕉久久成人网| 亚洲av免费高清在线观看| 大陆偷拍与自拍| av卡一久久| 亚洲精品美女久久久久99蜜臀 | 免费观看性生交大片5| 宅男免费午夜| 久久综合国产亚洲精品| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产精品无大码| 飞空精品影院首页| 草草在线视频免费看| 一个人免费看片子| 五月玫瑰六月丁香| 美女中出高潮动态图| 亚洲国产欧美在线一区| 欧美97在线视频| 色哟哟·www| 亚洲欧美色中文字幕在线| 少妇人妻 视频| 黑人巨大精品欧美一区二区蜜桃 | 午夜影院在线不卡| 少妇人妻 视频| 高清毛片免费看| 91aial.com中文字幕在线观看| 51国产日韩欧美| 欧美xxⅹ黑人| 99久久精品国产国产毛片| 欧美日韩一区二区视频在线观看视频在线| 9191精品国产免费久久| 国产精品一区二区在线观看99| 一级毛片我不卡| 亚洲国产精品一区二区三区在线| 亚洲精品日本国产第一区| 纵有疾风起免费观看全集完整版| 高清欧美精品videossex| 高清黄色对白视频在线免费看| 欧美 亚洲 国产 日韩一| 最近中文字幕2019免费版| 大片免费播放器 马上看| 熟女电影av网| 曰老女人黄片| 人成视频在线观看免费观看| 丝袜脚勾引网站| 插逼视频在线观看| 一本久久精品| 制服丝袜香蕉在线| 九色亚洲精品在线播放| 亚洲av中文av极速乱| 黄网站色视频无遮挡免费观看| 中文欧美无线码| 国产无遮挡羞羞视频在线观看| 成人手机av| 在线亚洲精品国产二区图片欧美| 婷婷色综合大香蕉| 成年av动漫网址| 美女国产视频在线观看| 国产乱人偷精品视频| 草草在线视频免费看| 伊人久久国产一区二区| 黑人高潮一二区| 久久国产精品男人的天堂亚洲 | 国产色爽女视频免费观看| 亚洲成国产人片在线观看| 乱码一卡2卡4卡精品| 国产精品国产三级国产av玫瑰| 99久久人妻综合| 亚洲综合色惰| 高清毛片免费看| 少妇人妻久久综合中文| 成人综合一区亚洲| 午夜久久久在线观看| 国产av一区二区精品久久| 亚洲内射少妇av| 久久久欧美国产精品| 人人澡人人妻人| 色94色欧美一区二区| av不卡在线播放| av网站免费在线观看视频| 黄网站色视频无遮挡免费观看| av免费观看日本| 婷婷色av中文字幕| 日韩不卡一区二区三区视频在线| 涩涩av久久男人的天堂| 大码成人一级视频| 考比视频在线观看| 人人澡人人妻人| 色94色欧美一区二区| 日韩免费高清中文字幕av| 中文天堂在线官网| 亚洲精品第二区| 久久99蜜桃精品久久| 18禁动态无遮挡网站| 国产一区二区在线观看av| 久久久国产欧美日韩av| 九色成人免费人妻av| 另类亚洲欧美激情| 久久久久人妻精品一区果冻| 免费看光身美女| 大片免费播放器 马上看| 国产女主播在线喷水免费视频网站| www日本在线高清视频| 国产极品天堂在线| 国产精品蜜桃在线观看| 18禁裸乳无遮挡动漫免费视频| 午夜免费男女啪啪视频观看| 中国国产av一级| 亚洲精品美女久久av网站| 久久午夜福利片| 国产成人精品久久久久久| 亚洲五月色婷婷综合| 五月伊人婷婷丁香| 国产日韩欧美亚洲二区| 国产高清国产精品国产三级| 波野结衣二区三区在线| 男人爽女人下面视频在线观看| av播播在线观看一区| 欧美成人午夜免费资源| 一本色道久久久久久精品综合| 另类亚洲欧美激情| 精品少妇内射三级| 亚洲成人av在线免费| 久久国内精品自在自线图片| 男女啪啪激烈高潮av片| 亚洲国产成人一精品久久久| 亚洲成国产人片在线观看| 亚洲精华国产精华液的使用体验| 黄色毛片三级朝国网站| 边亲边吃奶的免费视频| 一二三四在线观看免费中文在 | 久久国产亚洲av麻豆专区|