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

    基于Landsat 8 OLI影像的山區(qū)園地提取方法

    2021-05-28 12:38:10王剛王加勝苗旺元陳波
    現(xiàn)代計(jì)算機(jī) 2021年10期
    關(guān)鍵詞:昭陽區(qū)園地紋理

    王剛,王加勝,苗旺元,陳波

    (1.云南師范大學(xué)信息學(xué)院,昆明650500;2.西部資源環(huán)境地理信息技術(shù)教育部工程研究中心,昆明650500)

    0 引言

    山區(qū)園地的提取在國內(nèi)國外的相關(guān)研究普遍比較少,利用遙感技術(shù)對(duì)山區(qū)園地的提取也比較難。近年來,果園的提取主要是用光譜特征和植被指數(shù)來操作的。例如羅等人[1]把環(huán)境衛(wèi)星遙感影像作為數(shù)據(jù)基礎(chǔ)影像,把NDVI、波段、高程和紋理信息等作為提取果園的特征,然后使用決策樹的分類模型對(duì)研究區(qū)的果園進(jìn)行遙感提取和分類。李等人[2]和Immitzer等人[3]綜合了光譜信息提取并且對(duì)其進(jìn)行分類。因?yàn)榈孛嫣卣飨鄬?duì)來說比較復(fù)雜,結(jié)合光譜和植被指數(shù)特征的果樹提取精度通常是比較低的,容易出現(xiàn)分類錯(cuò)誤和遺漏的現(xiàn)象[4-5]。為了提高分類精度,國內(nèi)外的學(xué)者運(yùn)用了光譜結(jié)合紋理特征的提取方法進(jìn)行果樹研究[6-7]。例如岳等人[8]把不同高精度分辨率的高分一號(hào)遙感數(shù)據(jù)作為基礎(chǔ)影像,把光譜和紋理特征作為提取特征,分別使用多種遙感分類模型對(duì)研究區(qū)主要栽培的果樹進(jìn)行遙感識(shí)別分類。PU等人[9]和Dian等人[10]把空間信息和光譜特征作為提取特征對(duì)果樹進(jìn)行遙感提取分類,證明了把空間信息和光譜特征結(jié)合起來可以很大地提高果樹分類的準(zhǔn)確性。宋等人[11]把QuickBird影像作為基礎(chǔ)遙感影像,先利用小波變換的方法對(duì)紋理信息進(jìn)行提取,然后使用隨機(jī)森林這個(gè)分類方法自動(dòng)提取出獼猴桃果園。郝等人[12]利用面向?qū)ο蟮母叻直媛视跋窠Y(jié)合紋理特征對(duì)耕地進(jìn)行提取。黃等人[13]從紋理、形狀和對(duì)象的角度來進(jìn)行提取和分析影像的特征,在這個(gè)基礎(chǔ)上,用三個(gè)分辨率比較高的對(duì)影像有實(shí)際應(yīng)用的案例來驗(yàn)證自己提出的方法。陳等人[14]利用高分辨率影像,把光譜和紋理特征作為提取特征數(shù)據(jù),通過比較選取最好的紋理特征組合,使用決策樹這個(gè)分類方法進(jìn)行分類。

    綜上所述,在現(xiàn)有的果園提取方法中,結(jié)合光譜與植被指數(shù)特征的果樹提取精度相對(duì)較低,然而后來提出的一種把光譜和紋理特征作為提取特征的方法,很好地提升了對(duì)遙感影像的分類精度。但是因?yàn)檫b感影像空間分辨率的限制和制約,對(duì)影像的分類精度來說仍然具有比較大的提高空間。為此,亟需開展基于高分辨影像、融合光譜、地形、紋理等多特征的園地提取研究。

    GEE(Google Earth Engine)是一個(gè)綜合科學(xué)分析和地理信息數(shù)據(jù)可視化的平臺(tái)。本身存儲(chǔ)了近40年來主要公開的遙感影像數(shù)據(jù)集,并每天不斷進(jìn)行更新,使得GEE用戶可以方便地提取、調(diào)用和分析海量的遙感大數(shù)據(jù)資源。與常規(guī)影像處理軟件相比,GEE可以實(shí)現(xiàn)更加快速、在線批量處理大數(shù)據(jù)等功能,對(duì)各種植被指數(shù)、水體指數(shù)具有快速且精確計(jì)算的特點(diǎn)。

    本文擬以云南省昭通市昭陽區(qū)為研究區(qū)域,以GEE為軟件平臺(tái),把Landsat 8 OLI遙感影像作為基礎(chǔ)的數(shù)據(jù)影像,地形圖等地理數(shù)據(jù)作為輔助,對(duì)園地的特征進(jìn)行提取,主要包括了光譜特征、紋理特征、地形特征和纓帽變換特征,然后采用隨機(jī)森林(Random For-est)這個(gè)比較常見的分類方法進(jìn)行分類并且采用決策樹(Decision Tree)的分類方法和支持向量機(jī)(Support Vector Machine)的分類方法進(jìn)行提取2019年以昭陽區(qū)為研究區(qū)的園地并加以對(duì)比。

    1 研究區(qū)與數(shù)據(jù)

    1.1 研究區(qū)概況

    昭陽區(qū)隸屬于昭通市,在昭通市的西南方向,金沙江的南邊,在27.1°N-27.6°N、103.1°E-103.9°E之間,位于云貴高原的西北部(如圖1所示)。東、南、西、北邊分別與貴州省的威寧縣、昭通市的魯?shù)榭h、四川省的金陽縣、昭通市的彝良縣,大關(guān)縣和永善縣三縣相鄰,是昭通市最富裕的縣區(qū),占地總面積大約2167平方公里,占昭通市總面積的9.4%。昭陽區(qū)的氣候類型屬于北緯高原大陸季風(fēng)氣候,春夏秋冬四季的溫差比較小,干季濕季差別比較明顯,垂直差異比較大。全年無霜期約為220天,相對(duì)北方來說比較長,全年平均氣溫約為11.6攝氏度,比較適合水果的生長,全年平均降水量約為735毫米,全年平均日照時(shí)間約為1902.02小時(shí),晝夜的溫差較大,日溫差也比較大。該地區(qū)以紅壤為主,多為山區(qū),山區(qū)占64.3%,境內(nèi)水系發(fā)達(dá)。獨(dú)特的水文氣候條件,造就了“昭通蘋果”等水果品牌。

    1.2 數(shù)據(jù)與預(yù)處理

    (1)Landsat 8 OLI數(shù)據(jù)。

    根據(jù)昭陽區(qū)果園生長期以及影像的質(zhì)量等因素,選取來自GEE(https://earthengine.google.com/)平臺(tái)的28景Landsat 8 OLI遙感影像數(shù)據(jù)。時(shí)間是從2019年2月1日到2019年10月31日,包括果樹的發(fā)芽期、開花期、幼果期、膨大期、著色期和成熟期,分辨率為30m。首先對(duì)在GEE上獲取的Landsat 8 OLI遙感衛(wèi)星數(shù)據(jù)進(jìn)行數(shù)據(jù)的預(yù)處理工作,由于GEE中的Landsat 8 OLI遙感數(shù)據(jù)已經(jīng)進(jìn)行了大氣校正的工作,所以后續(xù)的預(yù)處理主要包括影像去云,影像融合鑲嵌,影像裁剪等過程,然后進(jìn)行特征提取。

    圖1 研究區(qū)位置示意圖

    (2)SRTM數(shù)據(jù)。

    本文使用SRTMGL1_003數(shù)據(jù),空間分辨率30 m。運(yùn)用STRM數(shù)據(jù)進(jìn)行特征提取,包括DEM、坡度、坡向和陰影等信息。

    (3)樣本數(shù)據(jù)。

    選擇合適的樣本點(diǎn)是提高植被分類結(jié)果的可靠性以及精度的前提。根據(jù)昭陽區(qū)的地形圖,土地利用方式和分布,將研究區(qū)的主要的土地作物種植類型化分成6種,分別為水體、草地、耕地、林地、建筑用地和園地,這6種土地作物種植類型在Google Earth平臺(tái)本身提供的高分辨率影像上具有非常明顯的識(shí)別特征。通過Google Earth平臺(tái)提供的2019年左右的高分辨率影像,利用目視解譯的辦法得到所有土地類型的樣本數(shù)據(jù),收集得到的樣本點(diǎn)隨機(jī)且均勻分布于整個(gè)昭陽區(qū)。把目標(biāo)年份的樣本數(shù)據(jù)先進(jìn)行隨機(jī)的分配,比例按照7:3來分配,其中采集到的70%的樣本數(shù)據(jù)放入三種分類器進(jìn)行分類訓(xùn)練,而剩下的30%的樣本數(shù)據(jù)用來進(jìn)行混淆矩陣的精度評(píng)價(jià)。樣本數(shù)如表1所示。

    2 研究方法

    研究技術(shù)路線如圖2所示。首先要對(duì)Landsat 8 OLI遙感影像進(jìn)行預(yù)處理的工作,由于GEE上的Land-sat數(shù)據(jù)已經(jīng)進(jìn)行了大氣校正,所以在這里主要進(jìn)行其他的一些預(yù)處理,包括影像去云、影像融合鑲嵌和影像裁剪等一系列過程,然后結(jié)合地形圖作為輔助數(shù)據(jù)進(jìn)行特征信息提取,包括物候特征、地物光譜特征、紋理特征、地形特征和纓帽變換特征,最后采用RF、CART、CART這三種常見的分類模型進(jìn)行分類并且進(jìn)行精度評(píng)價(jià)。

    圖2 技術(shù)路線圖

    2.1 特征提取

    特征變量的選擇對(duì)遙感分類來說非常重要,把不同的特征變量有機(jī)的組合起來可以使遙感分類的準(zhǔn)確性有比較明顯的改善[15]。本文綜合利用光譜特征、紋理特征、纓帽變換特征和地形特征進(jìn)行果園的遙感提取。

    (1)光譜特征。

    本文的光譜特征主要使用了一些多光譜波段和一些光譜指數(shù)[16]來表示的,光譜指數(shù)發(fā)展到今天有很多,本文使用了一些方便區(qū)分園地的光譜指數(shù),分別是歸一化植被指數(shù)(NDVI)、歸一化濕度指數(shù)(NDMI)、調(diào)整土壤亮度的植被指數(shù)(SAVI)、比值植被指數(shù)(RVI)和差值植被指數(shù)(DVI)。計(jì)算公式如表2所示。

    表2 光譜指數(shù)公式

    其中L表示的是土壤調(diào)整系數(shù),當(dāng)L=0時(shí),SAVI=NDVI,在研究中,L=0.5,植覆蓋度處于中等水平。

    (2)紋理特征。

    本文使用的紋理特征主要是通過灰度共生矩陣(GLCM)[17]計(jì)算得到的,分別是二階矩、對(duì)比度、相關(guān)性、方差和熵,計(jì)算公式如表3所示。

    表3 紋理特征和公式

    (3)纓帽變換特征。

    纓帽變換是根據(jù)固定的變換矩陣將原始影像投影綜合變換到具有物理意義的亮度(Brightness)、綠度(Greenness)和濕度(Wetness)特征向量的三維特征空間,分別反映了裸土巖石、植被覆蓋度和水分信息。這個(gè)變換的過程達(dá)到了減少特征維數(shù)、增強(qiáng)影像信息的效果。其變換公式如下:

    式中,X、Y分別為變換前后的多光譜空間的像元矢量;C為變換矩陣;b表示的是避免出現(xiàn)負(fù)值所添加的常數(shù)。Landsat 8 OLI影像下纓帽變換矩陣系數(shù)如表4所示[18]。

    表4 Landsat 8 OLI影像下纓帽變換矩陣系數(shù)

    (4)地形特征。

    地形特征包括基于SRTM的數(shù)字化高程模型(DEM)及其衍生產(chǎn)品(坡度、坡向、陰影)[19],引入地形特征作為限定條件,增加部分限定條件以提升精度。

    2.2 分類方法

    隨機(jī)森林分類算法是一種包含多個(gè)決策樹(CART)的機(jī)器學(xué)習(xí)算法,已被大量應(yīng)用在許多領(lǐng)域,例如土地利用變化圖等,這種分類方法在分類特征比較多的情況下被認(rèn)為具有比較好的分類效果[20-21]。隨機(jī)森林主要包括2個(gè)階段,分別是訓(xùn)練和分類。在訓(xùn)練這個(gè)階段,首先對(duì)輸入的數(shù)據(jù)通過每棵樹自舉重的方法對(duì)采集樣本進(jìn)行采樣,隨機(jī)重復(fù)抽取N個(gè)新的樣本集合來替代原來的訓(xùn)練樣本集合,然后通過這個(gè)樣本集合生成幾個(gè)獨(dú)立的CART決策樹。在分類這個(gè)階段,通過算術(shù)平均計(jì)算所有樹產(chǎn)生各種類別的概率,由每棵樹的分類結(jié)果通過投票得到[22-23]。最后通過GEE這個(gè)二次開發(fā)平臺(tái),使用隨機(jī)森林分類算法對(duì)2019年目標(biāo)年份的所有輸入的特征集合來進(jìn)行分類得到分類結(jié)果。

    2.3 精度評(píng)價(jià)

    在分類精度檢驗(yàn)與評(píng)估中,遙感分類結(jié)果需要使用參數(shù)指標(biāo)來片判別分類結(jié)果與實(shí)際土地類型的差異程度。把混淆矩陣(也稱誤差矩陣)作為評(píng)價(jià)分類的指標(biāo)[24]來展示結(jié)果的好壞,包括總體分類精度(OA)、Kap-pa系數(shù)(KC)、錯(cuò)分誤差(CE)、漏分誤差(OE)、用戶精度(UA)和生產(chǎn)者精度(PA)等,遙感影像分類的準(zhǔn)確性就是用這些評(píng)價(jià)指標(biāo)從不同的方面來表現(xiàn)出來的[25-26]。本文運(yùn)用了混淆矩陣中總體分類精度與Kappa系數(shù)來衡量試驗(yàn)區(qū)分類結(jié)果質(zhì)量。用驗(yàn)證樣本來計(jì)算目標(biāo)年份分類結(jié)果的OA、KC、UA、PA,用這些結(jié)果的好壞來評(píng)估果園遙感提取的準(zhǔn)確性。評(píng)價(jià)指標(biāo)如表5所示。

    表5 評(píng)價(jià)指標(biāo)

    3 結(jié)果與分析

    3.1 不同分類方法結(jié)果的精度對(duì)比

    用相同的樣本集合進(jìn)行SVM和CART的分類來評(píng)估RF算法的分類性能的好壞。3種算法的分類結(jié)果如表6所示。

    目標(biāo)年份的果園提取結(jié)果利用分類精度評(píng)價(jià)參數(shù)來評(píng)估。由表6可知,RF分類的OA、Kappa系數(shù)、UA和PA高于相同特征下的SVM和CART分類結(jié)果,與CART分類相比OA、KAPPA系數(shù)、UA和PA分別提升2%、3%、7%和6%,與SVM分類相比OA、Kappa系數(shù)、UA和PA分別提升6%、8%、11%和15%,表明RF提取果園的效果和穩(wěn)定性比較好。

    表6 分類結(jié)果

    3.2 三種不同分類方法的果園空間分布對(duì)比

    圖3表示的是通過不同的分類算法獲得的昭陽區(qū)的土地利用分類圖。可以明顯的看出,SVM分類結(jié)果中果園和耕地草地之間的錯(cuò)分現(xiàn)象較為嚴(yán)重,CART分類結(jié)果中果園分布完整性較好,與其他地物類型的混分狀況有較大改善。由圖3可以看出,CART分類與RF分類結(jié)果相對(duì)來說比較相似,但是果園和耕地之間錯(cuò)分漏分的現(xiàn)象有所增加。綜合三種分類結(jié)果表明,采用RF分類算法提取的果園與其他兩種分類方案相比具有更好的連續(xù)性,影像的整體分類效果在果園的提取上有明顯的提高,這與表6中的定量分析相符。經(jīng)統(tǒng)計(jì),基于RF分類算法提取的果園面積為318.9hm2。2019年果園分布圖顯示,昭陽區(qū)果園主要集中分布與灑漁、永豐、舊圃、北閘、小龍洞等鎮(zhèn)。

    3.3 研究區(qū)果園空間分布

    經(jīng)統(tǒng)計(jì),基于隨機(jī)森林分類算法提取的果園面積為318.9hm2,占昭陽區(qū)總面積的14.7%。2019年果園分布圖(如圖4所示)表明,昭陽區(qū)果園主要集中分布于灑漁、永豐、舊圃、北閘、小龍洞等鎮(zhèn)。

    4 結(jié)語

    本文通過GEE平臺(tái),選用覆蓋昭陽區(qū)的2019年共28個(gè)時(shí)相的Landsat 8 OLI遙感衛(wèi)星影像并完成數(shù)據(jù)的預(yù)處理,構(gòu)建包括光譜、紋理、地形、纓帽特征在內(nèi)的分類特征集,使用隨機(jī)森林分類算法,對(duì)昭陽區(qū)進(jìn)行園地信息提取,最終得到了2019年昭通區(qū)果園空間分布情況,對(duì)研究區(qū)園地栽種面積進(jìn)行估算,提取精度較高。研究得到的主要結(jié)論如下:

    圖3 三種分類方法結(jié)果圖

    圖4 果園分布圖

    (1)通過比較RF、CART和SVM這三種分類器的分類性能,結(jié)果表明,在相同分類特征下,RF算法的分類結(jié)果最好,OA和Kappa系數(shù)比其他兩種分類器要高,與CART分類相比OA提高了2%,與SVM分類相比OA提高了6%。

    (2)根據(jù)RF分類方法提取的園地結(jié)果,制作了昭陽區(qū)的園地分布圖,結(jié)果顯示,昭陽區(qū)園地總面積為318.9 hm2,占昭陽區(qū)總面積的14.7%,RF分類的OA達(dá)到96%,Kappa系數(shù)達(dá)到0.95,分類效果相對(duì)其他兩種分類器來說比較好。

    (3)在昭陽區(qū)各鎮(zhèn)均分布有園地,其中,灑漁鎮(zhèn)園地分布最多,果業(yè)為該鎮(zhèn)農(nóng)業(yè)產(chǎn)業(yè)中的主導(dǎo)產(chǎn)業(yè)。

    但是,基于像元的分類將無法避免導(dǎo)致“椒鹽現(xiàn)象”的發(fā)生,如何確定最佳的影像分割閾值,如何更有效地選取訓(xùn)練區(qū)域樣本,以及如何更準(zhǔn)確地選擇分類指標(biāo)與規(guī)則,都需要在后面繼續(xù)積極探索。

    猜你喜歡
    昭陽區(qū)園地紋理
    互動(dòng)園地
    封面人物:昭通市昭陽區(qū)第五小學(xué)黨總支書記楊正聰
    昭陽區(qū)蘋果園雜草種類及管理方式調(diào)查
    基于BM3D的復(fù)雜紋理區(qū)域圖像去噪
    軟件(2020年3期)2020-04-20 01:45:18
    昭陽區(qū)總工會(huì):舉辦職工氣排球周末聯(lián)賽初賽
    使用紋理疊加添加藝術(shù)畫特效
    昭陽區(qū)創(chuàng)建人與自然和諧生態(tài)園林城市的途徑和方法
    TEXTURE ON TEXTURE質(zhì)地上的紋理
    Coco薇(2017年8期)2017-08-03 15:23:38
    藝術(shù)園地
    藝術(shù)園地
    91老司机精品| 免费人成视频x8x8入口观看| 在线观看www视频免费| 亚洲人成网站高清观看| 麻豆国产av国片精品| 欧美性猛交╳xxx乱大交人| 精品人妻1区二区| 亚洲最大成人中文| 成人亚洲精品av一区二区| 18禁黄网站禁片午夜丰满| 国产高清有码在线观看视频 | 美女扒开内裤让男人捅视频| 亚洲五月婷婷丁香| 俄罗斯特黄特色一大片| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲成人久久性| 亚洲av中文字字幕乱码综合 | 老司机靠b影院| 亚洲国产精品合色在线| 国产成人精品久久二区二区免费| 啦啦啦韩国在线观看视频| 可以免费在线观看a视频的电影网站| 一级作爱视频免费观看| 久久久久国内视频| 免费看美女性在线毛片视频| 亚洲男人天堂网一区| 日韩 欧美 亚洲 中文字幕| 人人妻,人人澡人人爽秒播| 精品少妇一区二区三区视频日本电影| 欧美最黄视频在线播放免费| 亚洲男人的天堂狠狠| 中出人妻视频一区二区| 中文字幕精品亚洲无线码一区 | e午夜精品久久久久久久| 久久久久国内视频| 亚洲aⅴ乱码一区二区在线播放 | 国产一卡二卡三卡精品| 亚洲国产精品sss在线观看| 此物有八面人人有两片| 国产一区二区三区视频了| 欧美成人午夜精品| 在线观看66精品国产| 色综合婷婷激情| 国产成人av激情在线播放| 国内久久婷婷六月综合欲色啪| 日韩一卡2卡3卡4卡2021年| 欧美成狂野欧美在线观看| 18美女黄网站色大片免费观看| 叶爱在线成人免费视频播放| 精品免费久久久久久久清纯| 90打野战视频偷拍视频| 欧美激情 高清一区二区三区| 欧美日本亚洲视频在线播放| 丝袜在线中文字幕| 国产精品国产高清国产av| 国产熟女xx| 欧美一级a爱片免费观看看 | 国产精品 国内视频| 国产精品自产拍在线观看55亚洲| 草草在线视频免费看| 欧美av亚洲av综合av国产av| 亚洲中文字幕日韩| 91在线观看av| 亚洲国产精品久久男人天堂| 欧美乱妇无乱码| 日本a在线网址| 90打野战视频偷拍视频| 欧美激情 高清一区二区三区| 丰满的人妻完整版| 久久国产精品男人的天堂亚洲| 亚洲欧美日韩无卡精品| av中文乱码字幕在线| 国产乱人伦免费视频| 可以在线观看毛片的网站| 国产人伦9x9x在线观看| 最新在线观看一区二区三区| av福利片在线| 国产成人影院久久av| 国产高清有码在线观看视频 | 国产精品永久免费网站| 国产区一区二久久| 国产麻豆成人av免费视频| 日本撒尿小便嘘嘘汇集6| 伦理电影免费视频| 黄色丝袜av网址大全| 十分钟在线观看高清视频www| av福利片在线| 亚洲色图 男人天堂 中文字幕| 大香蕉久久成人网| 丝袜人妻中文字幕| 大香蕉久久成人网| 精品欧美国产一区二区三| 丝袜人妻中文字幕| 嫁个100分男人电影在线观看| 久久性视频一级片| 一区二区三区精品91| 免费看十八禁软件| 日韩欧美免费精品| 国产伦人伦偷精品视频| 91成年电影在线观看| 欧美性猛交黑人性爽| 亚洲熟妇熟女久久| 免费电影在线观看免费观看| 国产蜜桃级精品一区二区三区| 丁香欧美五月| 中文字幕人成人乱码亚洲影| 最好的美女福利视频网| 国产成人欧美| 一区二区三区高清视频在线| 天天躁夜夜躁狠狠躁躁| 久久精品人妻少妇| 国产精品综合久久久久久久免费| 色哟哟哟哟哟哟| 中文亚洲av片在线观看爽| ponron亚洲| 欧美+亚洲+日韩+国产| 成人国产一区最新在线观看| 波多野结衣av一区二区av| 美女大奶头视频| 精华霜和精华液先用哪个| 精华霜和精华液先用哪个| 国产又色又爽无遮挡免费看| www日本黄色视频网| 少妇被粗大的猛进出69影院| 精品乱码久久久久久99久播| 亚洲av成人av| 成人免费观看视频高清| 日韩视频一区二区在线观看| 91成年电影在线观看| 国产黄色小视频在线观看| 欧美激情高清一区二区三区| 国产伦人伦偷精品视频| 午夜免费激情av| 久久香蕉激情| 1024手机看黄色片| 国产麻豆成人av免费视频| x7x7x7水蜜桃| 99久久国产精品久久久| 午夜老司机福利片| 女警被强在线播放| 桃色一区二区三区在线观看| 免费在线观看亚洲国产| 久久精品国产99精品国产亚洲性色| 精品久久蜜臀av无| √禁漫天堂资源中文www| 亚洲在线自拍视频| 日日爽夜夜爽网站| 岛国在线观看网站| 免费无遮挡裸体视频| 麻豆成人av在线观看| 精品国产亚洲在线| 波多野结衣高清无吗| 亚洲自拍偷在线| 欧美午夜高清在线| 天堂动漫精品| 午夜福利免费观看在线| 18禁美女被吸乳视频| 99国产精品一区二区蜜桃av| 国产精品二区激情视频| www国产在线视频色| 精品少妇一区二区三区视频日本电影| 久久久久久久久中文| 精品国产乱码久久久久久男人| 日韩三级视频一区二区三区| 国内毛片毛片毛片毛片毛片| 色播亚洲综合网| 欧美成狂野欧美在线观看| 亚洲人成电影免费在线| 婷婷丁香在线五月| 伦理电影免费视频| 91大片在线观看| 日日摸夜夜添夜夜添小说| 久久 成人 亚洲| 天天一区二区日本电影三级| 女同久久另类99精品国产91| 久久精品人妻少妇| 亚洲第一av免费看| 老汉色av国产亚洲站长工具| www.999成人在线观看| 国产精品自产拍在线观看55亚洲| 丝袜人妻中文字幕| 欧美丝袜亚洲另类 | 黄色片一级片一级黄色片| 俺也久久电影网| 男人操女人黄网站| 久久久国产成人精品二区| 久久午夜综合久久蜜桃| 亚洲国产欧美网| 在线观看免费视频日本深夜| 村上凉子中文字幕在线| 国产精品免费一区二区三区在线| a级毛片a级免费在线| 1024视频免费在线观看| 99久久精品国产亚洲精品| 亚洲国产中文字幕在线视频| 韩国av一区二区三区四区| 亚洲熟女毛片儿| 午夜福利在线在线| 亚洲av美国av| 久久天躁狠狠躁夜夜2o2o| 一进一出抽搐动态| 精品第一国产精品| 一本大道久久a久久精品| 精品少妇一区二区三区视频日本电影| 亚洲精品在线观看二区| 国产欧美日韩精品亚洲av| 男女之事视频高清在线观看| 在线观看免费日韩欧美大片| 国产精品亚洲av一区麻豆| 精品一区二区三区av网在线观看| 又黄又粗又硬又大视频| 国产高清视频在线播放一区| 18禁黄网站禁片免费观看直播| 免费电影在线观看免费观看| 黄色视频不卡| 国产高清有码在线观看视频 | 久久99热这里只有精品18| 欧美日韩乱码在线| 麻豆成人午夜福利视频| 国产单亲对白刺激| 人人澡人人妻人| 色婷婷久久久亚洲欧美| 此物有八面人人有两片| 久久婷婷人人爽人人干人人爱| 国产精品爽爽va在线观看网站 | 一级作爱视频免费观看| 国内久久婷婷六月综合欲色啪| 黑人操中国人逼视频| 国产视频一区二区在线看| 波多野结衣av一区二区av| 女生性感内裤真人,穿戴方法视频| 亚洲精华国产精华精| 精品日产1卡2卡| 一区二区三区国产精品乱码| 一级a爱视频在线免费观看| 搡老岳熟女国产| 午夜福利免费观看在线| 午夜精品在线福利| 黄色丝袜av网址大全| 又黄又粗又硬又大视频| 老汉色∧v一级毛片| 久久婷婷人人爽人人干人人爱| 亚洲性夜色夜夜综合| 人人妻人人澡人人看| 久久久久精品国产欧美久久久| 日韩国内少妇激情av| 国产高清视频在线播放一区| 午夜a级毛片| 18禁美女被吸乳视频| 久9热在线精品视频| 校园春色视频在线观看| 91在线观看av| 一区二区三区高清视频在线| 国产蜜桃级精品一区二区三区| av在线天堂中文字幕| 中文字幕精品亚洲无线码一区 | 亚洲色图av天堂| 自线自在国产av| 两个人免费观看高清视频| 精品久久久久久久毛片微露脸| 1024手机看黄色片| 久久久国产成人精品二区| 国产1区2区3区精品| 一边摸一边做爽爽视频免费| 久久婷婷人人爽人人干人人爱| 久久久久久大精品| 亚洲一卡2卡3卡4卡5卡精品中文| 成熟少妇高潮喷水视频| 男人的好看免费观看在线视频 | 美女 人体艺术 gogo| 午夜成年电影在线免费观看| 老司机深夜福利视频在线观看| 午夜影院日韩av| 久久香蕉国产精品| 亚洲精品中文字幕一二三四区| 99精品久久久久人妻精品| 在线播放国产精品三级| 欧美激情久久久久久爽电影| 午夜免费成人在线视频| АⅤ资源中文在线天堂| 99riav亚洲国产免费| 精品午夜福利视频在线观看一区| 人人妻人人看人人澡| 国产爱豆传媒在线观看 | 国产高清激情床上av| 国产黄色小视频在线观看| 99久久精品国产亚洲精品| 欧美日本视频| 老司机靠b影院| 日本精品一区二区三区蜜桃| 18禁裸乳无遮挡免费网站照片 | 婷婷精品国产亚洲av| 免费av毛片视频| 欧美zozozo另类| 老司机午夜福利在线观看视频| 亚洲av片天天在线观看| 欧美亚洲日本最大视频资源| 丰满的人妻完整版| 亚洲五月天丁香| 99热这里只有精品一区 | 久久久国产欧美日韩av| 亚洲av日韩精品久久久久久密| 欧美乱妇无乱码| 精品第一国产精品| 亚洲三区欧美一区| 亚洲av美国av| 嫁个100分男人电影在线观看| 黄色丝袜av网址大全| 欧美精品亚洲一区二区| 一a级毛片在线观看| 俄罗斯特黄特色一大片| 女人被狂操c到高潮| 久久精品亚洲精品国产色婷小说| 国产男靠女视频免费网站| 国产精品日韩av在线免费观看| 久9热在线精品视频| 日日夜夜操网爽| 大型黄色视频在线免费观看| 99在线人妻在线中文字幕| 好男人在线观看高清免费视频 | 国产精品久久久av美女十八| 丝袜美腿诱惑在线| 99久久精品国产亚洲精品| 看片在线看免费视频| 一进一出抽搐gif免费好疼| 国产av不卡久久| 欧美成狂野欧美在线观看| 50天的宝宝边吃奶边哭怎么回事| 亚洲中文字幕日韩| 中文字幕精品亚洲无线码一区 | 黑人操中国人逼视频| 午夜福利在线在线| 免费看美女性在线毛片视频| 十分钟在线观看高清视频www| 老司机深夜福利视频在线观看| 日韩精品免费视频一区二区三区| 免费搜索国产男女视频| 高清在线国产一区| 妹子高潮喷水视频| 19禁男女啪啪无遮挡网站| 久久这里只有精品19| 国产成人精品久久二区二区免费| 欧美日本亚洲视频在线播放| 黄频高清免费视频| 熟女少妇亚洲综合色aaa.| 日韩精品中文字幕看吧| 国产高清激情床上av| 久久精品成人免费网站| 午夜久久久久精精品| 丝袜在线中文字幕| 色综合站精品国产| 十分钟在线观看高清视频www| 亚洲最大成人中文| 国内揄拍国产精品人妻在线 | 亚洲国产中文字幕在线视频| 日本成人三级电影网站| 亚洲av美国av| 久久久精品国产亚洲av高清涩受| 国产1区2区3区精品| 久久久久亚洲av毛片大全| 精品国产乱码久久久久久男人| 国产日本99.免费观看| 757午夜福利合集在线观看| 欧美黑人精品巨大| 高清毛片免费观看视频网站| 国产aⅴ精品一区二区三区波| 亚洲 欧美 日韩 在线 免费| 久久精品影院6| 精品无人区乱码1区二区| 欧美性长视频在线观看| 草草在线视频免费看| 99精品久久久久人妻精品| 国产成人精品久久二区二区免费| 午夜福利18| 欧美乱妇无乱码| 真人做人爱边吃奶动态| 色在线成人网| www.www免费av| 国产精品野战在线观看| 国产成人系列免费观看| 中文字幕另类日韩欧美亚洲嫩草| 亚洲欧美精品综合久久99| 欧美精品亚洲一区二区| 天堂动漫精品| 午夜老司机福利片| 又黄又粗又硬又大视频| 黄色女人牲交| 男人操女人黄网站| 国产亚洲精品综合一区在线观看 | 老司机深夜福利视频在线观看| 午夜a级毛片| 又大又爽又粗| 亚洲av中文字字幕乱码综合 | 中文字幕精品亚洲无线码一区 | 一进一出抽搐gif免费好疼| 俄罗斯特黄特色一大片| 亚洲精品国产一区二区精华液| 999久久久精品免费观看国产| 女性生殖器流出的白浆| 婷婷丁香在线五月| 一级a爱视频在线免费观看| 久久精品影院6| 欧美色欧美亚洲另类二区| 午夜激情av网站| 在线免费观看的www视频| 久久九九热精品免费| 美国免费a级毛片| 男女做爰动态图高潮gif福利片| 桃红色精品国产亚洲av| 日日夜夜操网爽| 18禁黄网站禁片免费观看直播| 久久精品成人免费网站| 国产精品电影一区二区三区| 男女午夜视频在线观看| 18禁美女被吸乳视频| 两个人看的免费小视频| 妹子高潮喷水视频| 成人午夜高清在线视频 | 国产一区在线观看成人免费| 久久亚洲真实| 窝窝影院91人妻| www日本黄色视频网| 黄片大片在线免费观看| 在线国产一区二区在线| 美女午夜性视频免费| 国产精品久久视频播放| 脱女人内裤的视频| 中文字幕人成人乱码亚洲影| 亚洲色图 男人天堂 中文字幕| 免费在线观看视频国产中文字幕亚洲| 一进一出好大好爽视频| 久99久视频精品免费| 亚洲一区中文字幕在线| 日本精品一区二区三区蜜桃| 国产单亲对白刺激| 欧美亚洲日本最大视频资源| 黄频高清免费视频| 精品第一国产精品| 2021天堂中文幕一二区在线观 | 天天一区二区日本电影三级| 亚洲精品粉嫩美女一区| 桃色一区二区三区在线观看| 国产精品亚洲一级av第二区| 国产精品av久久久久免费| 精品久久久久久久久久久久久 | 变态另类丝袜制服| 校园春色视频在线观看| 亚洲 国产 在线| 亚洲精华国产精华精| 亚洲av片天天在线观看| 国产私拍福利视频在线观看| 欧美中文日本在线观看视频| 身体一侧抽搐| а√天堂www在线а√下载| 国产精品电影一区二区三区| 国产高清激情床上av| 女警被强在线播放| 亚洲三区欧美一区| 中文字幕另类日韩欧美亚洲嫩草| www.熟女人妻精品国产| 一本综合久久免费| 天堂√8在线中文| 午夜福利成人在线免费观看| 欧美+亚洲+日韩+国产| 久久久久九九精品影院| 不卡av一区二区三区| 久久久久久国产a免费观看| 国产高清videossex| 伊人久久大香线蕉亚洲五| 18禁裸乳无遮挡免费网站照片 | 丝袜在线中文字幕| 男女视频在线观看网站免费 | 精品高清国产在线一区| 一区二区日韩欧美中文字幕| bbb黄色大片| 亚洲人成网站在线播放欧美日韩| 91成年电影在线观看| 国产高清videossex| 成人三级黄色视频| 欧美日韩福利视频一区二区| 曰老女人黄片| 国产精品精品国产色婷婷| 怎么达到女性高潮| 国产蜜桃级精品一区二区三区| 两人在一起打扑克的视频| 精品一区二区三区视频在线观看免费| 男女之事视频高清在线观看| 免费在线观看日本一区| 久久精品国产清高在天天线| 国产亚洲欧美精品永久| 99国产极品粉嫩在线观看| 国产av又大| 一级黄色大片毛片| 中文字幕另类日韩欧美亚洲嫩草| 啦啦啦观看免费观看视频高清| 日韩欧美一区视频在线观看| 国产精品亚洲av一区麻豆| 欧美一区二区精品小视频在线| 一级片免费观看大全| 精品久久久久久久人妻蜜臀av| 在线永久观看黄色视频| 亚洲一区高清亚洲精品| 成人国产一区最新在线观看| 香蕉国产在线看| 日本成人三级电影网站| 狠狠狠狠99中文字幕| netflix在线观看网站| 99在线人妻在线中文字幕| 午夜久久久久精精品| 国产精品香港三级国产av潘金莲| aaaaa片日本免费| 成人特级黄色片久久久久久久| 两个人视频免费观看高清| 午夜日韩欧美国产| 宅男免费午夜| √禁漫天堂资源中文www| 18禁黄网站禁片免费观看直播| 亚洲欧美激情综合另类| 不卡一级毛片| 一区二区三区精品91| 成年版毛片免费区| 欧美zozozo另类| 欧美 亚洲 国产 日韩一| 精品久久久久久久毛片微露脸| 久久久久久久久免费视频了| 国产成人系列免费观看| 91国产中文字幕| av超薄肉色丝袜交足视频| 啦啦啦韩国在线观看视频| 两个人免费观看高清视频| 非洲黑人性xxxx精品又粗又长| 曰老女人黄片| a级毛片a级免费在线| 国产精品久久电影中文字幕| 免费在线观看影片大全网站| а√天堂www在线а√下载| 欧美久久黑人一区二区| 国产真人三级小视频在线观看| 亚洲人成伊人成综合网2020| 成人av一区二区三区在线看| 嫁个100分男人电影在线观看| 久久香蕉精品热| 亚洲黑人精品在线| 亚洲午夜精品一区,二区,三区| 两个人视频免费观看高清| 亚洲自偷自拍图片 自拍| 精品久久久久久久毛片微露脸| 一区二区三区国产精品乱码| 国产一区二区在线av高清观看| 精品一区二区三区av网在线观看| 久久99热这里只有精品18| www.熟女人妻精品国产| 他把我摸到了高潮在线观看| 午夜成年电影在线免费观看| 日本 av在线| 国产av一区在线观看免费| 一进一出好大好爽视频| 在线播放国产精品三级| 欧美黑人精品巨大| 久久久久国产一级毛片高清牌| 天天躁狠狠躁夜夜躁狠狠躁| 十分钟在线观看高清视频www| 国产三级在线视频| 国产视频内射| 成人三级黄色视频| 一卡2卡三卡四卡精品乱码亚洲| 亚洲真实伦在线观看| 国产亚洲欧美98| 99re在线观看精品视频| 成年女人毛片免费观看观看9| 久久久久亚洲av毛片大全| 很黄的视频免费| 欧美日本亚洲视频在线播放| 成年免费大片在线观看| 午夜福利18| 99久久综合精品五月天人人| 午夜福利欧美成人| 国产久久久一区二区三区| 久热这里只有精品99| 国内久久婷婷六月综合欲色啪| 少妇被粗大的猛进出69影院| 久久午夜综合久久蜜桃| 美女高潮喷水抽搐中文字幕| 免费观看精品视频网站| 琪琪午夜伦伦电影理论片6080| 波多野结衣巨乳人妻| 国产精品自产拍在线观看55亚洲| 一本大道久久a久久精品| а√天堂www在线а√下载| 757午夜福利合集在线观看| 亚洲国产精品成人综合色| 国产精品一区二区精品视频观看| 最近最新免费中文字幕在线| 侵犯人妻中文字幕一二三四区| 又大又爽又粗| 身体一侧抽搐| bbb黄色大片| 久久人人精品亚洲av| 成熟少妇高潮喷水视频| 一区二区日韩欧美中文字幕| 一本精品99久久精品77| 亚洲av成人av| 99国产精品99久久久久| 国内毛片毛片毛片毛片毛片| 激情在线观看视频在线高清| 久久天躁狠狠躁夜夜2o2o| 999久久久精品免费观看国产| 麻豆国产av国片精品| 精品电影一区二区在线| 精品一区二区三区四区五区乱码| 村上凉子中文字幕在线| 91麻豆av在线| 一区二区日韩欧美中文字幕| 欧美久久黑人一区二区| av有码第一页|