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

    基于Python的Landsat8OLI影像建設用地自動識別與提取

    2018-01-19 11:35:34吳亞玲黃義忠彭秋志
    軟件導刊 2018年10期
    關鍵詞:建設用地遙感影像自動識別

    吳亞玲 黃義忠 彭秋志

    摘要:為了提高從海量遙感影像中解譯建設用地的效率,設計一種基于Python從Landsat8 OLI遙感影像中自動識別與提取建設用地的軟件系統(tǒng)。利用時域中值濾波去除云噪聲,并通過z值標準化獲得可適用統(tǒng)一閾值劃分的穩(wěn)定像元;結合多種地物指數(shù)開展閾值劃分以初步識別可能的建設用地像元;結合像元鄰域組合關系進一步精化提取結果,并自動計算對應的Kappa系數(shù)。研究表明,與傳統(tǒng)基于大型專業(yè)遙感影像解譯平臺相比,該系統(tǒng)不僅能保持較高的解譯精度,而且具有體積小、速度快、自動化程度高、可遷移性和可擴展性強等特點。

    關鍵詞:Python;遙感影像;自動識別;建設用地

    DOIDOI:10.11907/rjdk.181322

    中圖分類號:TP319

    文獻標識碼:A 文章編號:1672-7800(2018)010-0161-04

    英文摘要Abstract:In order to improve the efficiency of interpreting construction land from massive remote sensing images, the software system based on Python which can automatically identify and extract construction land from Landsat 8 OLI remote sensing images was designed. Firstly, median filtration was used to remove cloud interference, the stable pixels of uniform threshold value partitioning that was obtained through normalization of z value. Then,the threshold partitioning combined with various features indices are employed to identify possible pixels of construction land. Finally, the meta-neighborhood combination relationship is further refined to extract the result, and the corresponding Kappa coefficient is automatically calculated. The results show that the proposed system can not only maintain high accuracy of interpretation, but also has the advantages of small size, high speed and high degree of automation compared with the traditional method of recognition and extraction of construction land based on large professional remote sensing images interpretation system.

    英文關鍵詞Key Words:Python;remote sensing image;automatic identification;construction land

    0 引言

    社會經(jīng)濟的持續(xù)發(fā)展以及城市化進程中伴隨的城市人口、社會發(fā)展要素重組等,對城市建設用地利用方式產(chǎn)生了巨大影響[1-2]。城市建設用地動態(tài)監(jiān)測空間基礎數(shù)據(jù)中遙感影像分辨率的提高、分類方法的改進,對實現(xiàn)城市建設用地的高效利用具有重要意義[3]。

    隨著遙感影像大范圍普及,高效、快捷、自動化地獲取建設用地信息已經(jīng)成為一種趨勢[3-4]。近年來,為了獲得更高的分類精度,眾多從遙感影像中提取建設用地信息的方法相繼被提出。這些方法就分類單元而言可分為基于像元的分類和基于對象的分類[5],就分類途徑而言可分為監(jiān)督分類和非監(jiān)督分類[6],就分類結果表現(xiàn)形式而言可分為確定性分類和模糊分類[7-8],就分類所用信息類別而言可分為僅基于光譜信息的分類和以光譜信息為主的多信息組合分類等[9-10]。然而,日益多樣的方法雖然帶來更多選擇可能,但同時也增加了學習成本、增大了選擇困難。在實際工程應用中,不僅需要考慮解譯精度,而且需要兼顧簡單易學、便捷高效、成本低廉等要求。因此,有必要針對特定的專題解譯需求,開發(fā)出一系列自動化軟件系統(tǒng),以減輕對大型專業(yè)遙感影像處理平臺的依賴,大幅降低成本投入。

    本文借助面向對象的解釋型計算機程序設計語言——Python,綜合運用多波段指數(shù)法[11]和歸一化閾值控制法[12 -14],提出一種節(jié)約時間、更新速度快的建設用地識別與提取技術,以期自動化和批量化從海量遙感影像中解譯建設用地。

    1 系統(tǒng)原理與流程設計

    建設用地自動識別與提取系統(tǒng)流程如圖1所示?;驹硎牵米詣踊募闅v和迭代方法串聯(lián)起每一個解譯步驟,實現(xiàn)“進去的是原始數(shù)據(jù),出來的是最終產(chǎn)品”之目的。該系統(tǒng)主要包括6個關鍵解譯步驟:①對下載的原始遙感影像進行數(shù)據(jù)預處理,以便實現(xiàn)數(shù)據(jù)輸入自動化,主要模塊包括解壓縮文件、重命名文件、提取關鍵波段文件和獲取影像最小重疊區(qū);②疊合同一時期多期影像,利用時域中值濾波法實現(xiàn)去云降噪;③通過z值標準化獲得可適用統(tǒng)一閾值劃分的穩(wěn)定像元;④結合多種已被證明效果良好且穩(wěn)定的地物指數(shù)開展閾值劃分,以初步識別可能的建設用地像元;⑤結合像元鄰域組合關系進一步精化提取結果;⑥自動計算Kappa系數(shù),并根據(jù)該系數(shù)對前兩步進行迭代式調整,保證解譯結果的質量。

    2 系統(tǒng)實現(xiàn)

    系統(tǒng)實現(xiàn)基于Python語言進行編碼,數(shù)據(jù)預處理和質量控制部分通過自行編碼實現(xiàn),主要空間分析函數(shù)調用自ArcGIS中的空間分析模塊。

    2.1 數(shù)據(jù)預處理

    (1)對下載的原始數(shù)據(jù)進行批量解壓。對海量原始影像壓縮文件夾進行遍歷,識別并逐次解壓。文件解壓關鍵代碼如下:

    targzfile.extract(file, unzipedfolder) #從壓縮文件targzfile中把解壓出的文件file放到文件夾unzipedfolder里

    (2)重命名。在系統(tǒng)中,為了后續(xù)流程的完整性,對已解壓文件名不統(tǒng)一的波段重新命名。文件重命名關鍵代碼如下:

    arcpy.Rename_management(name_fr, name_to) #從將文件name_fr重命名為name_to

    (3)提取關鍵波段。海量影像的所有波段中,依據(jù)建設用地識別方法,對波段進行篩選,單獨提取出關鍵的波段數(shù),并進行另存。

    (4)獲取最小重疊區(qū)。去除空值和影像邊界,其留下的最小邊界作為后面研究范圍。以關鍵波段中某一波段為代表,利用Python腳本對其影像進行批量化掩摸,計算最小重疊區(qū)。關鍵代碼如下:

    arcpy.gp.CellStatistics_sa(in, out, "MINIMUM"; "DATA") #獲取所有輸入圖層的最小值

    剔除外圍以0值表示的空值黑邊,只保留有值區(qū),裁剪出最小重疊區(qū)。關鍵代碼如下:

    arcpy.gp.Con_sa(in , 1, out, "", "Value > 0") #剔除外圍黑邊

    2.2 云噪聲去除

    為了消除時相干擾,去除云噪聲,需進行時域中值濾波。其原理是云的反射率在各波段都位于異常極大值區(qū),而云影大多位于異常極小值區(qū),取中值具有很好的去云和去云影效果。關鍵代碼如下:

    arcpy.gp.CellStatistics_sa(in, out, "MEDIAN", "DATA") #獲取所有輸入圖層的中值

    2.3 穩(wěn)定像元獲取

    不同影像中波段值存在時相差異,需對生成的波段中值進行z值標準化,并計算影像的z值標準差,再以0.5倍標準差為閾值,替換不穩(wěn)定像元,最后反z值標準化,即基于統(tǒng)一的平均值和標準差系數(shù)反算回去,獲取每個單波段的穩(wěn)定像元圖層。該步驟能保證在后續(xù)步驟中影像均采用相同閾值進行分割。關鍵代碼如下:

    MEAN = arcpy.GetRasterProperties_management(in, "MEAN", "Band_1") #計算多期影像平均值

    STD = arcpy.GetRasterProperties_management(in, "STD", "Band_1") #計算多期影像標準差

    arcpy.gp.Minus_sa(in, MEAN, numerator) #獲得每期影像與多期影像平均值之差,以作為標準化的分子numerator

    arcpy.gp.Divide_sa(numerator, STD, out) #獲得每期影像與numerator之比,得到該期影像的標準化值

    arcpy.gp.Con_sa(in ,1, out, "", "Value < 0.5") #獲得均值兩側0.5倍標準差范圍內的像元,即穩(wěn)定像元

    2.4 建設用地初步識別

    建設用地提取過程實質是遙感影像信息的處理過程[15] ,遙感技術手段具有客觀、宏觀、快速、準確、動態(tài)等特點[16]。不同類型的地物反射光譜強度在不同波段上具有不一致性,分析和找出建設用地與其它背景地物間的差異,充分利用遙感影像光譜信息、時相信息、地形信息和鄰域信息,逐步排除確定性無關地物。運用波段差異提取出其特征信息的NDVI(歸一化植被指數(shù))[17]、MNDWI(歸一化水體指數(shù))[18],其值介于(-1,1)之間,經(jīng)過其處理后比值運算用于消除地形差異的影響。

    通過地物指數(shù)計算閾值分割方法,在已選用的地物指數(shù)參數(shù)中,構建各參數(shù)的閾值范圍,如:NDVI、MNDWI。同時在提取城市建設用地等地物中,僅依據(jù)指數(shù)去除植被和水體是完全不夠的,在系統(tǒng)中根據(jù)每兩個波段間的歸一化差值取閾值去除高亮地物(旱裸土、部分灘涂、部分高亮人工地面、部分高亮屋頂)、水體及陰坡植被、裸地及植被、海島湖泊及灘涂等其它不屬于建設用地的地物類型。其中地物參數(shù)選取和閾值分割均在Python中批量化處理,用到的波段組合原理為:

    NDVI=NIR-RedNIR+Red(1)

    MNDWI=Green-SWIR1Green+SWIR1(2)

    D73=SWIR2×Red(3)

    ND75=SWIR2-SWIR1SWIR2-SWIR1(4)

    ND42=NIR-GreenNIR+Green(5)

    ND75ND42=ND75-ND42(6)

    式(1)-式(6)均表示在Python中計算用到的指數(shù)。針對實驗區(qū)現(xiàn)狀,反復試驗后,優(yōu)先構建出一套完整的建設用地提取流程,用來提取建設用地信息。其中,部分關鍵代碼如下:

    (1)根據(jù)波段計算出NDVI、MNDWI指數(shù)。

    (2)采用歸一化閾值法確定植被的閾值范圍為0.25,粗略去除植被。

    arcpy.gp.Con_sa(MNDWI ,1, Pout, 0, 'Value > 0.25') #獲取像元大于0.25的NDVI,用于識別植被

    (3)采用步驟(2)中的歸一化閾值法,確定水體的閾值范圍為0.1,去除水體。

    arcpy.gp.Con_sa(MNDWI ,1, Pout, 0, 'Value > 0.10') #獲取像元大于0.1的MNDWI,用于識別水體

    (4)在前3步的基礎上,通過光譜和紋理分析,去除試驗區(qū)內比較明顯的植被、水體以后,通過光譜差異和借用實驗區(qū)范圍內的Google Earth影像對比,進一步識別出需要識別去除的高亮地物(含裸土、部分灘涂、部分高亮人工地面、部分高亮屋頂)、水體及偏陰坡植被、裸地及植被。

    arcpy.gp.Con_sa(D73 , 1, Pout, 0, 'Value > 0.06') #獲取像元大于0.06的D73,用于識別高亮地物

    arcpy.gp.Con_sa(D73 ,1, Pout,0, 'Value < 0.01') #獲取像元小于0.01的D73,用于識別水體及偏陰坡植被

    arcpy.gp.Con_sa(ND75ND42 ,1, Pout, 0, 'Value < -0.40') #獲取像元大于-0.04的ND75ND42,識別裸地及植被

    2.5 建設用地精化提取

    初步識別可能的建設用地后,系統(tǒng)用于精化提取建設用地的方法主要是采用空間鄰域關系計算。分別對水體、植被的邊界像元通過鄰域均值計算后,采用條件函數(shù)逐一判別消除其它地物。其關鍵代碼如下:

    arcpy.gp.Con_sa(Vm , 1, Vm01 , 0, "VALUE > 0.25") #獲取鄰域內均值Vm超過0.25的植被

    arcpy.gp.Con_sa(Wm , 1, Wm01 , 0, "VALUE > 0.25") ##獲取鄰域內均值Wm超過0.25的水體

    arcpy.gp.CellStatistics_sa(Vm01&Wm01;, Vm01sumWm01, "SUM", "DATA") #獲取鄰域內水體和植被Vm01& Wm01均超過0.25的像元的和

    arcpy.gp.Con_sa(Vm01sumWm01 ,1 , Pout, 0, "VALUE = 2") #獲取上一步中VALUE = 2的值

    2.6 質量控制

    對遙感分類精度進行評價,需要參考數(shù)據(jù)作為評價基準。系統(tǒng)選擇Google Earth影像作為參考數(shù)據(jù),在針對完整流程提取建設用地的基礎上,于所選實驗區(qū)內隨機生成1 000個參考點,對建設用地提取結果進行精度評價。采用的評價指標為Kappa系數(shù),其計算在Python中進行。

    3 實例驗證

    3.1 實驗區(qū)選擇與數(shù)據(jù)下載

    實驗選取昆明中心城區(qū)(五華區(qū)、盤龍區(qū)、西山區(qū)、官渡區(qū)、呈貢新區(qū)和空港經(jīng)濟區(qū))作為實驗區(qū),如圖2所示。采用2015年的Landsat 8 OLI遙感影像數(shù)據(jù),數(shù)據(jù)下載自https://glovis.usgs.gov/。為保障影像質量,避免云量干擾,下載的影像數(shù)據(jù)云量均控制在10%以內。對地面同一地點,下載同時區(qū)影像數(shù)量確保大于或等于10,本文實驗區(qū)所在的影像行列號為129043。以上條件限制主要為了在后續(xù)自動提取建設用地過程中,能利用盡量多的時相信息消除云干擾影響,減少季相差異,從而提高分類精度。其Landsat 8 OLI遙感影像數(shù)據(jù)均經(jīng)過Level1標準地形校正,即經(jīng)過系列輻射校正、地面控制點幾何校正和DEM地形校正,并以單波段輸出格式為GeoTIFF的形式封裝在壓縮包文件內,取樣方式均為三次卷積,地圖投影均采用UTM-WGS84投影。本文中運用到的波段為:Band3(Green)、Band4(Red)、Band5(NIR)、Band6(SWIR1)、Band7(SWIR2)5個波段,分辨率均為30m。

    3.2 建設用地提取結果

    圖2顯示了研究區(qū)Landsat 8 OLI影像的RGB波段合成效果,圖3是利用本文所述方法自動提取的建設用地效果。直觀對比,本文所述方法能夠很好地將建設用地與其它地物區(qū)別開來,特別是能很好地剔除耕地和灘涂;就精度定量評估而言,建設用地提取結果的Kappa系數(shù)為0.812 5,與大多數(shù)現(xiàn)有方法提取建設用地結果的精度水平相當,可以滿足大部分相關研究和應用的精度要求。

    4 結語

    本文在回顧和總結遙感影像建設用地提取研究的基礎上,介紹了一種基于Python從Landsat8 OLI遙感影像中自動識別與提取建設用地的技術方法。與其它建設用地提取方法相比,其可以在大范圍遙感影像覆蓋區(qū)域實時、高效、準確地自動提取建設用地信息,為快速監(jiān)測和評價建設用地變化提供新的途徑。在后續(xù)研究中,可將本文研究方法延伸到其它地物的獲取、監(jiān)測等領域。

    參考文獻:

    [1] 袁麗麗.城市化進程中城市土地可持續(xù)利用研究[D].武漢:華中農(nóng)業(yè)大學, 2005.

    [2] 馬兵.揚州市城市化進程中城鄉(xiāng)建設用地動態(tài)變化研究[D].南京:南京農(nóng)業(yè)大學, 2009.

    [3] 王宏勝,李永樹,吳璽,等.結合空間分析的面向對象無人機影像土地利用分類[J]. 測繪工程, 2018(2):57-61.

    [4] 何少林,徐京華,張帥毅.面向對象的多尺度無人機影像土地利用信息提取[J].國土資源遙感,2013,25(2):107-112.

    [5] 張洋華,劉慧平,楊曉彤.基于像元與對象分類的景觀指數(shù)差異性分析[J].遙感技術與應用,2016,31(1):119-125.

    [6] 趙春霞,錢樂祥.遙感影像監(jiān)督分類與非監(jiān)督分類的比較[J].河南大學學報:自然科學版,2004,34(3):90-93.

    [7] 張卉,張超,趙冬玲.特征權重優(yōu)化高分辨率遙感影像模糊分類研究[J].遙感信息,2010(1):94-98.

    [8] 劉艷芳,蘭澤英,劉洋,等.基于混合熵模型的遙感分類不確定性的多尺度評價方法研究[J].測繪學報,2009,38(1):82-87.

    [9] 陳晨,張友靜.基于多尺度紋理和光譜信息的SVM分類研究[J].測繪科學,2009,34(1):29-31.

    [10] 厲小潤,朱潔爾,王晶,等.組合核支持向量機高光譜圖像分類[J].浙江大學學報:工學版,2013,47(8):1403-1410.

    [11] 樊彥國,李翔宇,張磊,等.基于多波段分析的無閾值自動光譜角制圖分類法[J].地理與地理信息科學,2010, 26(2):38-41.

    [12] VIOVY N, ARINO O, BELWARD A S. The best index slope ex-traction(BISE):a method for reducing noise in NDVI time series[J].International Journal of Remote Sensing,1992,13:1585-1590.

    [13] 徐影,丁一匯,趙宗慈.人類活動引起的我國西北地區(qū)21世紀溫度和降水變化情景分析[J].冰川凍土,2003,25(3):327-330.

    [14] 付忠良.一些新的圖像閾值選取方法[J].計算機應用,2000,20(10):15-17.

    [15] 李愛民.基于遙感影像的城市建成區(qū)擴張與用地規(guī)模研究[D].鄭州:解放軍信息工程大學,2009.

    [16] 邢雅娟,劉東升,王鵬新.遙感信息與作物生長模型的耦合應用研究進展[J].地球科學進展,2009,24(4):444-451.

    [17] DEFRIES R, TOWNSHEND J. NDVI-derived land cover classification at global scales[J]. International Journal of Remote Sensing, 1994, 15(17):3567-3586.

    [18] 徐涵秋.福清市城鎮(zhèn)空間擴展規(guī)律及其驅動機制分析[J].遙感技術與應用, 2002, 17(2):86-92.

    (責任編輯:何 麗)

    猜你喜歡
    建設用地遙感影像自動識別
    自動識別系統(tǒng)
    特別健康(2018年3期)2018-07-04 00:40:18
    遙感影像資料在海圖制圖中的應用研究
    航海(2016年6期)2017-01-09 11:28:24
    論建設用地適宜性與國土資源規(guī)劃
    耕地保護約束下建設用地合理發(fā)展措施
    安徽省經(jīng)濟增長與建設用地擴張的脫鉤分析
    金屬垃圾自動識別回收箱
    遙感數(shù)字圖像處理課程實驗綜述
    科技視界(2016年13期)2016-06-13 12:01:14
    臥龍湖濕地遙感調查和演變分析
    基于IEC61850的配網(wǎng)終端自動識別技術
    電測與儀表(2016年6期)2016-04-11 12:06:38
    城鄉(xiāng)存量建設用地盤活再開發(fā)利益相關者訴求均衡分析及應對
    商(2016年1期)2016-03-03 10:51:18
    国产亚洲欧美精品永久| 在线观看一区二区三区激情| 亚洲精品久久午夜乱码| 国产激情久久老熟女| 好男人视频免费观看在线| 亚洲精品中文字幕在线视频| 国产精品久久久久久av不卡| 久久久久人妻精品一区果冻| 女的被弄到高潮叫床怎么办| 一区二区三区乱码不卡18| 极品人妻少妇av视频| 国产极品天堂在线| 这个男人来自地球电影免费观看 | 久久 成人 亚洲| 国产69精品久久久久777片| 国产日韩欧美视频二区| 在线观看免费视频网站a站| 在线观看一区二区三区激情| 欧美精品一区二区大全| 欧美另类一区| 亚洲欧美日韩另类电影网站| 日韩在线高清观看一区二区三区| videosex国产| 91aial.com中文字幕在线观看| 最近最新中文字幕免费大全7| 国产精品免费大片| 亚洲五月色婷婷综合| 成人综合一区亚洲| 亚洲精品自拍成人| 青青草视频在线视频观看| 麻豆精品久久久久久蜜桃| 欧美变态另类bdsm刘玥| 婷婷色麻豆天堂久久| 有码 亚洲区| 老司机影院成人| 国产高清三级在线| 成人毛片a级毛片在线播放| 亚洲精品,欧美精品| 丝袜人妻中文字幕| 国产精品久久久久久久久免| 晚上一个人看的免费电影| 国产精品久久久久久久久免| 黄色 视频免费看| 亚洲国产成人一精品久久久| 激情视频va一区二区三区| 亚洲伊人色综图| 久久久精品免费免费高清| 午夜福利在线观看免费完整高清在| 亚洲av.av天堂| 好男人视频免费观看在线| 丝瓜视频免费看黄片| 国产精品.久久久| 国产一区二区在线观看日韩| 久久 成人 亚洲| 男人舔女人的私密视频| 一级毛片黄色毛片免费观看视频| 中文精品一卡2卡3卡4更新| 日本欧美国产在线视频| 另类亚洲欧美激情| 亚洲精品视频女| 伦理电影免费视频| 如何舔出高潮| 国产av码专区亚洲av| 亚洲伊人久久精品综合| 日本欧美国产在线视频| 80岁老熟妇乱子伦牲交| 国产永久视频网站| 国产综合精华液| 婷婷色av中文字幕| www.av在线官网国产| 如日韩欧美国产精品一区二区三区| 欧美性感艳星| 欧美成人午夜精品| 国产日韩欧美亚洲二区| 桃花免费在线播放| 欧美精品国产亚洲| 免费观看性生交大片5| 老司机影院成人| 999精品在线视频| 成人午夜精彩视频在线观看| 亚洲久久久国产精品| 欧美人与善性xxx| 久久99一区二区三区| 成人国产麻豆网| 国产精品一区二区在线不卡| 一本大道久久a久久精品| 国产成人欧美| 久热这里只有精品99| 欧美最新免费一区二区三区| 久久99精品国语久久久| 少妇的丰满在线观看| 青青草视频在线视频观看| 高清欧美精品videossex| 国产精品蜜桃在线观看| 久久久久久久亚洲中文字幕| a级片在线免费高清观看视频| 天堂中文最新版在线下载| 免费观看a级毛片全部| 亚洲欧洲日产国产| 久久国产亚洲av麻豆专区| 亚洲精品第二区| 成年动漫av网址| 岛国毛片在线播放| 日本免费在线观看一区| 一级,二级,三级黄色视频| av免费观看日本| 搡女人真爽免费视频火全软件| 黑人高潮一二区| 在线观看免费高清a一片| 伊人久久国产一区二区| 一级,二级,三级黄色视频| 两个人看的免费小视频| 高清在线视频一区二区三区| 亚洲精品aⅴ在线观看| 男女高潮啪啪啪动态图| 韩国精品一区二区三区 | 校园人妻丝袜中文字幕| 免费黄频网站在线观看国产| 2022亚洲国产成人精品| 美女内射精品一级片tv| 精品人妻在线不人妻| 男女边摸边吃奶| 中国国产av一级| 香蕉精品网在线| 亚洲中文av在线| 七月丁香在线播放| 99热6这里只有精品| 高清视频免费观看一区二区| 我的女老师完整版在线观看| 一边摸一边做爽爽视频免费| 少妇高潮的动态图| 亚洲国产欧美在线一区| 国产成人欧美| 国产成人精品无人区| 亚洲精品国产色婷婷电影| 蜜桃国产av成人99| 中文字幕人妻丝袜制服| 有码 亚洲区| 欧美3d第一页| 欧美xxⅹ黑人| 激情五月婷婷亚洲| 国产在线一区二区三区精| 亚洲第一区二区三区不卡| 欧美少妇被猛烈插入视频| 婷婷色av中文字幕| 国产成人精品无人区| 成人漫画全彩无遮挡| 中文字幕亚洲精品专区| 人成视频在线观看免费观看| av在线播放精品| videos熟女内射| 伊人亚洲综合成人网| 观看美女的网站| 老熟女久久久| 超色免费av| 久久国产亚洲av麻豆专区| 黄色 视频免费看| 精品亚洲成a人片在线观看| 国产毛片在线视频| 色网站视频免费| av免费在线看不卡| 亚洲色图综合在线观看| 国产乱来视频区| 国产亚洲av片在线观看秒播厂| 成人免费观看视频高清| 日产精品乱码卡一卡2卡三| 大码成人一级视频| 黄片播放在线免费| 国产精品熟女久久久久浪| 另类精品久久| 黑人巨大精品欧美一区二区蜜桃 | 在线观看三级黄色| 亚洲内射少妇av| 美女内射精品一级片tv| 国产精品熟女久久久久浪| 国产 精品1| 老熟女久久久| 在线天堂中文资源库| 成人国语在线视频| 国产69精品久久久久777片| 免费高清在线观看日韩| 国产片内射在线| 免费黄网站久久成人精品| 人成视频在线观看免费观看| 18禁在线无遮挡免费观看视频| 国产xxxxx性猛交| 香蕉国产在线看| 美国免费a级毛片| a级片在线免费高清观看视频| 日韩制服骚丝袜av| 激情视频va一区二区三区| 人体艺术视频欧美日本| 精品少妇久久久久久888优播| 日本91视频免费播放| 免费高清在线观看视频在线观看| 国产精品久久久久成人av| 亚洲,欧美精品.| 国产成人91sexporn| 一级黄片播放器| 久久精品人人爽人人爽视色| 人人澡人人妻人| 日本91视频免费播放| av在线观看视频网站免费| 亚洲婷婷狠狠爱综合网| 日韩不卡一区二区三区视频在线| 免费不卡的大黄色大毛片视频在线观看| 亚洲精品日本国产第一区| 国产精品不卡视频一区二区| 国产精品人妻久久久影院| 婷婷色av中文字幕| 制服丝袜香蕉在线| 日韩成人av中文字幕在线观看| 国产高清三级在线| 少妇人妻 视频| 日韩 亚洲 欧美在线| 亚洲中文av在线| 色哟哟·www| 在线亚洲精品国产二区图片欧美| 亚洲精品乱码久久久久久按摩| 亚洲欧洲精品一区二区精品久久久 | 欧美人与性动交α欧美精品济南到 | 99精国产麻豆久久婷婷| 一级毛片 在线播放| 99久国产av精品国产电影| h视频一区二区三区| 亚洲精品视频女| 在线观看美女被高潮喷水网站| 婷婷色av中文字幕| av有码第一页| 精品一区二区三区四区五区乱码 | 成人国语在线视频| 毛片一级片免费看久久久久| 亚洲人与动物交配视频| 日本欧美视频一区| 成人亚洲欧美一区二区av| 一级毛片 在线播放| 啦啦啦视频在线资源免费观看| 国产69精品久久久久777片| 18禁国产床啪视频网站| 国产极品天堂在线| 欧美少妇被猛烈插入视频| 国产精品一区www在线观看| 午夜影院在线不卡| 寂寞人妻少妇视频99o| 一级片'在线观看视频| 免费黄网站久久成人精品| 国产黄色免费在线视频| a级毛片黄视频| 久久久久久久久久久免费av| 亚洲精品久久午夜乱码| 国产精品人妻久久久久久| 亚洲成国产人片在线观看| 丝袜在线中文字幕| av在线观看视频网站免费| 欧美国产精品va在线观看不卡| 欧美日本中文国产一区发布| 国产精品一区二区在线不卡| 成人午夜精彩视频在线观看| 啦啦啦在线观看免费高清www| 欧美成人午夜免费资源| 少妇被粗大猛烈的视频| 三级国产精品片| 国产免费福利视频在线观看| av国产久精品久网站免费入址| 亚洲久久久国产精品| 亚洲欧洲日产国产| 男女无遮挡免费网站观看| 草草在线视频免费看| 久久久精品区二区三区| 亚洲精品久久久久久婷婷小说| 日本欧美国产在线视频| 国产精品久久久久久久电影| 国产精品久久久久久av不卡| 国产一区二区激情短视频 | 黄网站色视频无遮挡免费观看| 精品熟女少妇av免费看| 亚洲国产最新在线播放| 看十八女毛片水多多多| 啦啦啦中文免费视频观看日本| 精品人妻在线不人妻| 一本色道久久久久久精品综合| 女人被躁到高潮嗷嗷叫费观| 午夜av观看不卡| 97人妻天天添夜夜摸| 国产精品久久久久久久电影| 精品久久久精品久久久| 免费高清在线观看日韩| 精品少妇内射三级| 黄色一级大片看看| xxxhd国产人妻xxx| 午夜福利,免费看| 有码 亚洲区| 亚洲国产毛片av蜜桃av| 日日撸夜夜添| 国产日韩欧美亚洲二区| 色婷婷av一区二区三区视频| 日韩制服丝袜自拍偷拍| 日韩中字成人| 丝袜在线中文字幕| 国产熟女午夜一区二区三区| 在线观看www视频免费| 成人毛片a级毛片在线播放| 亚洲欧美清纯卡通| 国产熟女欧美一区二区| 美女内射精品一级片tv| 波野结衣二区三区在线| 久久精品久久精品一区二区三区| 免费看光身美女| 中文字幕最新亚洲高清| 免费看不卡的av| 午夜久久久在线观看| 久久久国产欧美日韩av| 18禁裸乳无遮挡动漫免费视频| 久久人人97超碰香蕉20202| 亚洲色图 男人天堂 中文字幕 | 最近最新中文字幕大全免费视频 | 亚洲av国产av综合av卡| 亚洲五月色婷婷综合| 亚洲精品aⅴ在线观看| 欧美人与性动交α欧美软件 | 国产精品一区www在线观看| 免费观看在线日韩| 国内精品宾馆在线| 寂寞人妻少妇视频99o| 一本久久精品| 欧美最新免费一区二区三区| 午夜91福利影院| 大香蕉久久网| 精品少妇内射三级| 成人亚洲欧美一区二区av| 午夜福利视频精品| 亚洲经典国产精华液单| 好男人视频免费观看在线| 女人久久www免费人成看片| 日本av手机在线免费观看| 久久亚洲国产成人精品v| 黑人猛操日本美女一级片| 免费人成在线观看视频色| av又黄又爽大尺度在线免费看| 极品少妇高潮喷水抽搐| 色吧在线观看| 午夜影院在线不卡| 国产成人精品婷婷| 日韩欧美一区视频在线观看| 欧美少妇被猛烈插入视频| 日韩一区二区视频免费看| 啦啦啦在线观看免费高清www| 国产精品久久久av美女十八| 成人国产av品久久久| 亚洲国产精品成人久久小说| 边亲边吃奶的免费视频| 久久精品夜色国产| 免费人妻精品一区二区三区视频| 丰满饥渴人妻一区二区三| 日韩av在线免费看完整版不卡| 亚洲国产精品成人久久小说| 亚洲激情五月婷婷啪啪| 侵犯人妻中文字幕一二三四区| 日韩成人伦理影院| 青青草视频在线视频观看| 国产毛片在线视频| 校园人妻丝袜中文字幕| 亚洲中文av在线| 国产1区2区3区精品| 黄色一级大片看看| 好男人视频免费观看在线| 欧美最新免费一区二区三区| 亚洲精品国产av成人精品| 九九在线视频观看精品| 我要看黄色一级片免费的| 中文字幕人妻丝袜制服| 中文乱码字字幕精品一区二区三区| 亚洲av日韩在线播放| 欧美国产精品va在线观看不卡| 寂寞人妻少妇视频99o| 老司机亚洲免费影院| 中文字幕人妻丝袜制服| 亚洲国产色片| 国产永久视频网站| 天堂俺去俺来也www色官网| 黄片无遮挡物在线观看| 国产欧美日韩一区二区三区在线| 国产国语露脸激情在线看| 18禁观看日本| 国产一区二区在线观看av| 久久精品国产a三级三级三级| 日韩制服骚丝袜av| 亚洲第一av免费看| 黑人猛操日本美女一级片| 亚洲人成网站在线观看播放| 人人澡人人妻人| 精品人妻熟女毛片av久久网站| 久久精品久久精品一区二区三区| 国产精品一区二区在线观看99| 18禁在线无遮挡免费观看视频| 国产男女超爽视频在线观看| 在现免费观看毛片| 国产免费一级a男人的天堂| 99久久精品国产国产毛片| 高清黄色对白视频在线免费看| 亚洲精品日本国产第一区| 成人毛片a级毛片在线播放| 日韩 亚洲 欧美在线| 汤姆久久久久久久影院中文字幕| 韩国高清视频一区二区三区| 26uuu在线亚洲综合色| 老熟女久久久| 国产精品人妻久久久影院| 成年动漫av网址| 精品一区在线观看国产| 宅男免费午夜| 一二三四中文在线观看免费高清| 插逼视频在线观看| 成年人午夜在线观看视频| 男女边摸边吃奶| 久久精品久久久久久久性| 69精品国产乱码久久久| 亚洲国产精品国产精品| 一级片'在线观看视频| 考比视频在线观看| 校园人妻丝袜中文字幕| 免费高清在线观看日韩| 久久久久网色| 亚洲精品,欧美精品| 如何舔出高潮| 成人国产av品久久久| 婷婷色av中文字幕| 精品一区二区三区四区五区乱码 | 精品久久久久久电影网| 成人综合一区亚洲| 成人漫画全彩无遮挡| 久久鲁丝午夜福利片| 在线观看免费高清a一片| 熟女电影av网| 丝袜脚勾引网站| 国产精品国产三级专区第一集| 国产成人欧美| 色吧在线观看| 欧美变态另类bdsm刘玥| 久久精品久久精品一区二区三区| 九九在线视频观看精品| 97人妻天天添夜夜摸| 只有这里有精品99| 韩国精品一区二区三区 | 久久av网站| 中文字幕av电影在线播放| 久久久久久久久久人人人人人人| 亚洲国产精品成人久久小说| 黄片无遮挡物在线观看| 精品亚洲乱码少妇综合久久| 亚洲综合色网址| 亚洲丝袜综合中文字幕| 久久人人爽人人爽人人片va| 国产亚洲精品第一综合不卡 | 婷婷色av中文字幕| 在线观看免费日韩欧美大片| 天天躁夜夜躁狠狠躁躁| 日韩中文字幕视频在线看片| 两个人免费观看高清视频| 天堂8中文在线网| 欧美 亚洲 国产 日韩一| 中文字幕另类日韩欧美亚洲嫩草| 久久精品夜色国产| 99热国产这里只有精品6| 国产亚洲av片在线观看秒播厂| 国产激情久久老熟女| 亚洲,欧美精品.| 久久韩国三级中文字幕| 寂寞人妻少妇视频99o| 欧美成人午夜精品| 亚洲欧美成人综合另类久久久| 国产精品久久久av美女十八| 国产日韩一区二区三区精品不卡| 日韩熟女老妇一区二区性免费视频| 大香蕉97超碰在线| 这个男人来自地球电影免费观看 | 黄网站色视频无遮挡免费观看| 最新中文字幕久久久久| 国产色婷婷99| 校园人妻丝袜中文字幕| 欧美日韩一区二区视频在线观看视频在线| 亚洲欧美日韩另类电影网站| 亚洲欧洲国产日韩| 亚洲三级黄色毛片| 精品少妇久久久久久888优播| 人体艺术视频欧美日本| 亚洲精品成人av观看孕妇| av福利片在线| 男人添女人高潮全过程视频| 人人澡人人妻人| 最近中文字幕2019免费版| 国产 精品1| 观看av在线不卡| 免费黄频网站在线观看国产| 成人18禁高潮啪啪吃奶动态图| 国产免费现黄频在线看| 97人妻天天添夜夜摸| 欧美日韩精品成人综合77777| 国产亚洲av片在线观看秒播厂| 91久久精品国产一区二区三区| 国产成人精品福利久久| 色婷婷久久久亚洲欧美| 国产精品无大码| 精品亚洲成a人片在线观看| 男女高潮啪啪啪动态图| 亚洲经典国产精华液单| 51国产日韩欧美| 26uuu在线亚洲综合色| 亚洲国产看品久久| 日本爱情动作片www.在线观看| 精品熟女少妇av免费看| 亚洲av综合色区一区| 国产精品国产三级国产av玫瑰| 18+在线观看网站| 精品少妇内射三级| 亚洲精品乱码久久久久久按摩| 久久久久久人妻| 国产精品免费大片| 五月开心婷婷网| 国产精品久久久久久精品电影小说| 亚洲天堂av无毛| 亚洲第一av免费看| 久久青草综合色| 寂寞人妻少妇视频99o| 婷婷色av中文字幕| 亚洲av福利一区| 免费看光身美女| 老司机影院成人| 日本色播在线视频| 国产xxxxx性猛交| 日本爱情动作片www.在线观看| 成年美女黄网站色视频大全免费| 成人二区视频| 久久人人爽人人爽人人片va| 岛国毛片在线播放| 熟女人妻精品中文字幕| 久久这里只有精品19| 亚洲人成网站在线观看播放| 99九九在线精品视频| 国产极品天堂在线| 免费黄频网站在线观看国产| 综合色丁香网| 国产成人精品无人区| 亚洲国产成人一精品久久久| 日本猛色少妇xxxxx猛交久久| 国产精品.久久久| 欧美xxxx性猛交bbbb| 国产极品粉嫩免费观看在线| 一区二区日韩欧美中文字幕 | 精品少妇黑人巨大在线播放| 国产老妇伦熟女老妇高清| 久久精品国产自在天天线| 秋霞伦理黄片| 好男人视频免费观看在线| 蜜臀久久99精品久久宅男| 免费黄频网站在线观看国产| 日本av手机在线免费观看| 国产女主播在线喷水免费视频网站| 纯流量卡能插随身wifi吗| 男女边吃奶边做爰视频| 亚洲,欧美,日韩| 女性被躁到高潮视频| 日韩中字成人| 成人午夜精彩视频在线观看| 丰满少妇做爰视频| 如日韩欧美国产精品一区二区三区| 国产又爽黄色视频| 日日摸夜夜添夜夜爱| 国产精品 国内视频| 国产乱人偷精品视频| 久久人人爽av亚洲精品天堂| 在线看a的网站| 内地一区二区视频在线| 国产片内射在线| 你懂的网址亚洲精品在线观看| 蜜桃国产av成人99| 一级a做视频免费观看| 国产男女超爽视频在线观看| 妹子高潮喷水视频| 一边亲一边摸免费视频| 日韩熟女老妇一区二区性免费视频| 免费黄色在线免费观看| 亚洲成国产人片在线观看| 久久鲁丝午夜福利片| 51国产日韩欧美| 成人国产av品久久久| 又黄又爽又刺激的免费视频.| 精品少妇久久久久久888优播| 欧美精品人与动牲交sv欧美| 亚洲精品一区蜜桃| 视频中文字幕在线观看| 午夜福利网站1000一区二区三区| av线在线观看网站| 在线观看美女被高潮喷水网站| 纵有疾风起免费观看全集完整版| av线在线观看网站| 综合色丁香网| 插逼视频在线观看| 欧美xxxx性猛交bbbb| 日韩中字成人| 亚洲第一av免费看| 日本av手机在线免费观看| 天天躁夜夜躁狠狠久久av| 日本91视频免费播放| 日日爽夜夜爽网站| 国产精品国产三级专区第一集| 久久国产精品男人的天堂亚洲 | 欧美精品一区二区大全| 十八禁高潮呻吟视频| 天美传媒精品一区二区| 男女啪啪激烈高潮av片| 国产精品麻豆人妻色哟哟久久| av国产精品久久久久影院| 视频在线观看一区二区三区| 亚洲欧美色中文字幕在线|