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

    基于高分影像的林地變化檢測(cè)和自動(dòng)分割算法探討

    2022-02-21 11:04:36付曉燕岳彩榮蔣之富
    山東林業(yè)科技 2022年6期
    關(guān)鍵詞:變化檢測(cè)工具欄格網(wǎng)

    付曉燕,岳彩榮,蔣之富

    (1.西南林業(yè)大學(xué),云南 昆明 650224;2.國家林業(yè)和草原局西南調(diào)查規(guī)劃院,云南 昆明 650031)

    森林資源動(dòng)態(tài)監(jiān)測(cè)已經(jīng)成為林業(yè)的重要工作,當(dāng)前森林資源變化檢測(cè)主要依靠目視判讀對(duì)比前后兩期遙感影像,根據(jù)像素值對(duì)比發(fā)現(xiàn)森林資源疑似變化區(qū)域。但遙感影像內(nèi)容豐富、覆蓋區(qū)域往往比較大,通過人眼識(shí)別、手工區(qū)劃變化區(qū)域,不僅工作量大、效率低,而且可能存在漏劃、錯(cuò)劃的情況[1]。目前,業(yè)內(nèi)還沒有非常有效的技術(shù)手段對(duì)遙感影像變化區(qū)域進(jìn)行自動(dòng)檢測(cè)、區(qū)劃。本文旨在通過研究多期遙感影像的變化特點(diǎn),提出自動(dòng)區(qū)劃遙感影像變化區(qū)域的算法、方法,并在林業(yè)最為常用的地理信息系統(tǒng)軟件ArcGIS上實(shí)現(xiàn)該算法,開發(fā)相應(yīng)的插件和軟件工具,以快速、較為準(zhǔn)確地通過兩期遙感影像對(duì)比自動(dòng)獲得監(jiān)測(cè)區(qū)的變化圖斑。

    1 研究方法

    利用ArcGIS插件開發(fā)技術(shù)和ArcObjects SDK,在研究現(xiàn)有的遙感影像變化檢測(cè)算法的基礎(chǔ)上,提出一種快速進(jìn)行遙感影像變化檢測(cè)和變化區(qū)域分割、提取算法,實(shí)現(xiàn)遙感影像變化檢測(cè)半自動(dòng)區(qū)劃工具、遙感影像變化檢測(cè)自動(dòng)區(qū)劃工具。

    2 研究內(nèi)容

    研究內(nèi)容主要涉及到ArcGIS插件開發(fā)技術(shù)、遙感影像數(shù)據(jù)管理模塊開發(fā)、兩期遙感影像自動(dòng)檢查和搜索算法、遙感影像變化檢測(cè)和變化區(qū)域分割提取算法、算法優(yōu)化等內(nèi)容。需要指出的是,本文所說的林地,是指由林業(yè)主管部門管理的、有植被覆蓋的土地。

    2.1 遙感影像數(shù)據(jù)管理

    在有限的時(shí)間和資源基礎(chǔ)上,完全從底層開發(fā)影像數(shù)據(jù)管理功能已不太現(xiàn)實(shí),采用現(xiàn)有的、成熟的軟件平臺(tái)可以將重點(diǎn)放在遙感影像變化檢測(cè)算法的設(shè)計(jì)和改進(jìn)上,忽略細(xì)枝末節(jié)。基于這個(gè)原因,作者采用目前業(yè)內(nèi)最為成熟的ArcGIS軟件,作為影像數(shù)據(jù)加載、顯示、管理的基本平臺(tái),在ArcGIS軟件上,采用插件開發(fā)技術(shù)、ArcObjects SDK進(jìn)行影像數(shù)據(jù)處理、變化檢測(cè)的定制開發(fā)。

    2.2 遙感影像自動(dòng)檢查和搜索

    目前現(xiàn)有的軟件平臺(tái)在進(jìn)行變化檢測(cè)時(shí),均需要人為確保前后兩期影像是否為同一個(gè)區(qū)域、坐標(biāo)系是否一致、遙感影像范圍是否一致,影像檢測(cè)前往往需要裁剪、坐標(biāo)系變換等工作,增加了人為工作量。為減輕人為工作量,筆者設(shè)計(jì)了影像區(qū)域檢測(cè)模塊,對(duì)影像的坐標(biāo)范圍、坐標(biāo)系統(tǒng)進(jìn)行變化檢測(cè),自動(dòng)提取前后兩期影像的公共范圍。

    2.3 遙感影像變化檢測(cè)和變化區(qū)域提取算法研究

    影像變化檢測(cè)的關(guān)鍵在于變化檢測(cè)和變化區(qū)域提取算法研究和設(shè)計(jì),現(xiàn)有的軟件平臺(tái)主要采用監(jiān)督分類、無監(jiān)督分類兩類算法,監(jiān)督分類需要訓(xùn)練樣本,人為干預(yù)地方很多,要實(shí)現(xiàn)自動(dòng)變化檢測(cè),這一方法基本可以排除。無監(jiān)督分類大多基于統(tǒng)計(jì)理論,如主成分分析、NDVI指數(shù)、相似像素聚類等?,F(xiàn)有這些方法都基于嚴(yán)格的數(shù)學(xué)模型,從理論上說,能較好地提取變化區(qū)域,但在實(shí)際使用過程中,提取效果總是差強(qiáng)人意。筆者認(rèn)為主要原因在于,這些算法往往是基于比較好的影像數(shù)據(jù),設(shè)計(jì)算法的人也是基于規(guī)整的影像數(shù)據(jù),沒有考慮實(shí)際使用過程中諸多因素,如前后兩期影像的一小部分存在陰影,兩期影像校正效果不理想,局部存在偏移,兩期影像范圍不一致等等,導(dǎo)致產(chǎn)生很多偏差。故有必要對(duì)現(xiàn)有算法進(jìn)行改進(jìn),在算法的檢測(cè)過程中,對(duì)出現(xiàn)的偏差進(jìn)行自動(dòng)排除。

    2.4 基于多線程下的優(yōu)化算法

    對(duì)于大范圍的高分辨率遙感影像,數(shù)據(jù)量往往很大,難以一次性對(duì)整幅遙感影像進(jìn)行分割,需要設(shè)計(jì)一種優(yōu)化算法。本文采用的是異步分割算法,該算法在邏輯上將整幅遙感影像劃分成多個(gè)格網(wǎng),分割過程是按格網(wǎng)進(jìn)行異步分割。格網(wǎng)只是邏輯劃分,無需進(jìn)行柵格影像切割,原因是ArcObjects開發(fā)的程序可同時(shí)對(duì)一幅遙感影像加載多次,這與多個(gè)ArcGIS進(jìn)程可以加載同一幅影像是相類似的。當(dāng)要分割一個(gè)格網(wǎng)區(qū)域時(shí),只需要指定格網(wǎng)的左下角坐標(biāo)、行數(shù)、列數(shù)和需要分割的影像文件位置,就可以指定一個(gè)分割范圍,然后工具啟動(dòng)多個(gè)線程對(duì)格網(wǎng)進(jìn)行異步分割。

    計(jì)算機(jī)處理數(shù)據(jù)時(shí)有單線程和多線程之分,單線程就是進(jìn)程內(nèi)用1個(gè)主要線程進(jìn)行數(shù)據(jù)處理,多線程是啟用多個(gè)線程同時(shí)進(jìn)行數(shù)據(jù)處理。為提升分割效率,遙感影像分割肯定需要啟用多線程技術(shù)?,F(xiàn)在的計(jì)算機(jī)處理器普遍在4核心以上,可以充分利用多核心的優(yōu)勢(shì),進(jìn)行異步分割。本文設(shè)計(jì)了一種算法,該算法在運(yùn)行過程中,將會(huì)自動(dòng)獲取計(jì)算機(jī)物理核心數(shù)和可用內(nèi)存,根據(jù)物理核心數(shù),分配多個(gè)線程同時(shí)進(jìn)行分割。對(duì)于4核心處理器計(jì)算機(jī),工具最多啟動(dòng)2個(gè)線程;對(duì)于8核心計(jì)算機(jī),工具最多啟動(dòng)6個(gè)分割線程;對(duì)于16核心以上處理器,工具最多啟動(dòng)14個(gè)分割線程。也就是說,對(duì)于一臺(tái)16核心處理器的計(jì)算機(jī),如果內(nèi)存充足,工具將會(huì)啟動(dòng)14個(gè)分割線程同時(shí)進(jìn)行分割。假設(shè)每個(gè)分割線程分割的格網(wǎng)大小是3500×3500個(gè)像素,高分遙感影像分辨率為2m,則每個(gè)分割線程分割的面積是7000×7000m,即49.0 km2,14個(gè)線程可以同時(shí)分割686.0 km2的區(qū)域,因?yàn)樾枰指畹倪b感影像往往不是方形的,因此實(shí)際分割的有效區(qū)域可能會(huì)更小一些。通過多線程進(jìn)行優(yōu)化,可以極大地提升運(yùn)行效率。

    2.5 遙感影像半自動(dòng)區(qū)劃模塊、自動(dòng)分割模塊實(shí)現(xiàn)

    在Visual Studio 2010開發(fā)工具上,應(yīng)用C#和.Net 4.0,結(jié)合ArcGIS插件開發(fā)技術(shù),筆者將以上算法編寫代碼,封裝成插件和對(duì)應(yīng)的軟件工具,以工具欄的形式鑲嵌在ArcMap主界面上,既能方便林業(yè)工作人員調(diào)用,也能充分利用現(xiàn)有ArcGIS工具,提升工作效率。

    3 關(guān)鍵算法

    3.1 前后兩期遙感影像自動(dòng)檢查和搜索技術(shù)研究

    在用分割算法分割前,需要檢查影像坐標(biāo)系統(tǒng)、空間位置等是否正確。工具設(shè)計(jì)了影像質(zhì)量檢測(cè)和搜索算法,能根據(jù)遙感影像拍攝時(shí)間、遙感影像像素值,自動(dòng)確定子區(qū)域劃分方法。

    3.2 變化區(qū)域搜索、像元聚類、自動(dòng)分割算法實(shí)現(xiàn)

    根據(jù)搜索到的前后兩期遙感影像,在邏輯上劃分成3500×3500個(gè)像素大小的格網(wǎng),在格網(wǎng)內(nèi),再劃分成2×2個(gè)像素的細(xì)小格網(wǎng),按4個(gè)像素大小的區(qū)域進(jìn)行像素比較,對(duì)細(xì)小格網(wǎng)內(nèi)像素值很接近(如像素R、G、B值在15%以內(nèi))的,由該區(qū)域開始往外聚類,聚類后面積大于一定值(如100m2)的區(qū)域,標(biāo)記本期疑似區(qū)域。對(duì)搜索到的前期遙感影像,按相同方法處理,由該區(qū)域進(jìn)行聚類,標(biāo)記為前期疑似區(qū)域。用本期疑似區(qū)域按空間位置選中前期疑似區(qū)域的結(jié)果,將選中區(qū)域調(diào)用ArcObjects的圖層相減后,進(jìn)行多部件打散,刪除面積小于200m2的面,得到疑似前期林地轉(zhuǎn)變成本期非林地圖斑結(jié)果;分割結(jié)果中的前期疑似非林地圖斑減本期疑似林地圖斑,得到前期非林地轉(zhuǎn)化成本期有林地圖斑結(jié)果。

    3.3 區(qū)劃圖斑邊界調(diào)整算法

    自動(dòng)分割完成后,分割邊界沿像素值邊界前進(jìn),可能存在分割邊界點(diǎn)過密、鋸齒狀明顯的問題,實(shí)際應(yīng)用仍然存在諸多不便,需要對(duì)圖斑邊界進(jìn)行抽稀、修正、去除面積過小的內(nèi)環(huán)(如內(nèi)環(huán)面積小于50m2),按林業(yè)經(jīng)營界線(如林班界)進(jìn)行進(jìn)一步分割,達(dá)到符合林業(yè)區(qū)劃標(biāo)準(zhǔn)的要求。

    3.4 變化區(qū)域自動(dòng)分割算法優(yōu)化

    大范圍的高分遙感影像,數(shù)據(jù)量極大,自動(dòng)勾繪運(yùn)行速度可能較慢,需要進(jìn)一步優(yōu)化,優(yōu)化主要從兩方面進(jìn)行。

    (1)盡量縮小自動(dòng)分割范圍。例如,在只需要監(jiān)測(cè)林地減少的情況,結(jié)合林業(yè)經(jīng)營資料,剔除非林地區(qū)域,僅僅在林地范圍內(nèi)進(jìn)行自動(dòng)勾繪;在只需要監(jiān)測(cè)林地增加的情況,結(jié)合林業(yè)經(jīng)營資料,剔除有林地區(qū)域;既要監(jiān)測(cè)林地增加、也要監(jiān)測(cè)林地減少的情況,盡量利用監(jiān)測(cè)區(qū)域的行政邊界,減少區(qū)劃范圍。

    (2)充分利用現(xiàn)有計(jì)算機(jī)的多核心、大內(nèi)存等優(yōu)勢(shì),提高運(yùn)算性能?,F(xiàn)有的計(jì)算機(jī)一般都有四核甚至六核心以上處理,在變化區(qū)域勾繪時(shí),充分利用多核心優(yōu)勢(shì),同時(shí)異步分割多個(gè)3500×3500個(gè)像素格網(wǎng),如在六核處理器的計(jì)算機(jī)上,通過多線程技術(shù),同時(shí)對(duì)4個(gè)3500×3500個(gè)像素格網(wǎng)進(jìn)行自動(dòng)分割,對(duì)于具有16核心處理器的工作站,同時(shí)對(duì)14個(gè)3500×3500個(gè)像素的格網(wǎng)進(jìn)行自動(dòng)分割。

    研究過程中發(fā)現(xiàn),格網(wǎng)的形狀最好為正方形,正方形方便程序處理和優(yōu)化。格網(wǎng)范圍也不能過大,格網(wǎng)大小不能超過3500×3500個(gè)像素但最好大于1000×1000個(gè)像,原因是ArcObjects開發(fā)的程序是基于32位進(jìn)程的,32位進(jìn)程占用的總內(nèi)存不能超過4GB,在分割過程中,需要存儲(chǔ)每個(gè)像素點(diǎn)的RGB分量、主要波段等信息,如果范圍過大,會(huì)導(dǎo)致內(nèi)存溢出。對(duì)于1m分辨率的高分遙感影像,3500×3500個(gè)像素相當(dāng)于3500×3500m的區(qū)域,即12.25 km2;對(duì)于2m分辨率的高分遙感影像,3500×3500個(gè)像素相當(dāng)于7000×7000米的區(qū)域,即49 km2,每個(gè)子區(qū)域的實(shí)際范圍并不小。需要說明的是,因?yàn)橐指畹倪b感影像往往不是方形的,還存在一些背景區(qū)域,在邊界上需要分割的有效影像會(huì)更小一些。格網(wǎng)如果劃分過小,會(huì)導(dǎo)致格網(wǎng)數(shù)量增加,降低數(shù)據(jù)處理效率,增加后期融合的工作量。對(duì)于2m的高分遙感影像,云南省的一個(gè)縣,格網(wǎng)劃分個(gè)數(shù)一般在40~80個(gè)之間。對(duì)于16核心的計(jì)算機(jī)來說,可最多啟用14個(gè)分割線程,每個(gè)線程分割時(shí)間大致在3~10分鐘,每個(gè)計(jì)算機(jī)物理核心平均運(yùn)行3到6輪可以分割一個(gè)縣,運(yùn)行時(shí)間視政區(qū)面積,大致在15~70分鐘可以完成一個(gè)縣的自動(dòng)分割。

    3.5 分割區(qū)域合并和消除

    因?yàn)樵诜指钸^程中,劃分了多個(gè)格網(wǎng),因此異步分割后還需要對(duì)分割結(jié)果進(jìn)行融合、合并操作,將按格網(wǎng)的分割結(jié)果融合成一個(gè)整體。融合操作主要采用ArcObjects的Merge算法進(jìn)行合并,合并后調(diào)用ArcObjects的GP工具的融合算法,按類別進(jìn)行融合和消除,并刪除面積小于200m2的細(xì)碎圖斑。在工具運(yùn)行過程中,可以調(diào)用ArcGIS的工具箱來進(jìn)行一些常見的空間分析運(yùn)算,這也是插件開發(fā)最大的優(yōu)點(diǎn)。

    4 運(yùn)行效果

    4.1 工具欄封裝效果

    為便于使用,作者將算法封裝成多個(gè)模塊,利用C#4.0程序設(shè)計(jì)語言和Visual Studio 2010開發(fā)軟件,編寫成2個(gè)ArcGIS插件“遙感影像變化監(jiān)測(cè)半自動(dòng)區(qū)劃工具”和“遙感影像變化監(jiān)測(cè)自動(dòng)區(qū)劃工具”,輔助判讀工具可以內(nèi)嵌在ArcMap 10.2以上的版本中,方便與ArcGIS的其他工具協(xié)同操作,封裝后的工具欄運(yùn)行效果見圖1,上面的工具為半自動(dòng)區(qū)劃工具,下面的工具為自動(dòng)分割工具,因工具欄過短,無法顯示所有的文字。

    圖1 工具欄運(yùn)行效果Figure 1 Segmentation toolbar running effects

    4.2 半自動(dòng)區(qū)劃工具運(yùn)行效果

    半自動(dòng)區(qū)劃原理類似于Photoshop中的魔法棒,仍然需要借助人力查找到某個(gè)變化區(qū)域,在變化區(qū)域點(diǎn)一下鼠標(biāo),或者繪制一條線段,工具利用分割算法自動(dòng)分割該變化區(qū)域,因?yàn)樾枰肆φ业阶兓瘏^(qū)域,因此稱為半自動(dòng)區(qū)劃工具。通過在屏幕上單擊某個(gè)區(qū)域,工具自動(dòng)將該區(qū)域相鄰的像素進(jìn)行聚類和勾繪,勾繪方法和效果如圖2所示。在圖2的左圖,找到變化區(qū)域后,手工點(diǎn)一個(gè)點(diǎn)或者繪制一條線段,工具將會(huì)自動(dòng)分割并將結(jié)果形成面。勾繪結(jié)果自動(dòng)平滑后,效果不亞于手工勾繪。繪制過程中,如果存在已有圖斑,工具還能自動(dòng)扣除已經(jīng)有的圖斑,如圖2右下角的一小塊圖斑,工具可以自動(dòng)扣除,省去了手工捕捉的繁瑣工作過程。

    圖2 半自動(dòng)區(qū)劃工具運(yùn)行效果Figure 2 Semi-automatic segmentation running effects

    4.3 自動(dòng)分割工具運(yùn)行效果

    因自動(dòng)分割工具并不需要手工勾繪圖斑邊界,因此自動(dòng)分割工具欄按鈕要比輔助判讀工具欄界面簡潔,自動(dòng)分割工具欄只有6個(gè)按鈕(如圖1下半部分工具欄)。該工具欄最重要的按鈕是最右邊的遙感影像自動(dòng)分割按鈕,點(diǎn)擊該按鈕后,自動(dòng)分割運(yùn)行界面如圖3所示。自動(dòng)分割時(shí)需要錄入本期、前期遙感影像,通過圖3界面左側(cè)的功能區(qū)進(jìn)行數(shù)據(jù)錄入,通過圖3右側(cè)的進(jìn)程監(jiān)控區(qū)域進(jìn)程分割進(jìn)程監(jiān)控。

    圖3 自動(dòng)分割工具運(yùn)行效果Figure 3 Automatic segmentation running window

    自動(dòng)分割運(yùn)行效果如圖4所示,圖4為云南省昆明市高分1號(hào)遙感影像,遙感影像分辨率為2 m,圖4中左圖為本期變化后的影像,圖4右圖為前期未發(fā)生變化的影像,黃色邊界為自動(dòng)分割獲取的變化區(qū)域。自動(dòng)區(qū)劃能獲取到變化區(qū)域,但是在邊界線方面仍然不是很精細(xì),存在人工進(jìn)一步調(diào)整的必要。

    圖4 自動(dòng)分割運(yùn)行效果Figure 4 Automatic segmentation running results

    5 結(jié)語

    本文設(shè)計(jì)的方法,優(yōu)點(diǎn)是采用多線程技術(shù),提取速度較快,特別是對(duì)變化較為明顯的區(qū)域,效果較好。自動(dòng)和半自動(dòng)勾繪的邊界非常接近專業(yè)人員的勾繪(如圖2所示)。分割算法不僅適用于高分辨率遙感影像,也適應(yīng)于無人機(jī)航拍的高清影像。算法為筆者自行編寫,可以隨時(shí)根據(jù)業(yè)務(wù)需要進(jìn)行定制和更改,可以根據(jù)衛(wèi)星影像情況進(jìn)行優(yōu)化。但是改進(jìn)現(xiàn)有算法對(duì)計(jì)算機(jī)技術(shù)和算法要求很高,需要做很多測(cè)試,筆者在實(shí)現(xiàn)過程中也發(fā)現(xiàn)一些不足的地方,亟待解決,主要體現(xiàn)在以下3個(gè)方面:

    (1)細(xì)長型的地物識(shí)別,可能會(huì)存在間斷的問題(如圖5所示)。如在林區(qū)開挖的小路,因路面寬度較小,受樹冠陰影影響較大,可能會(huì)存在間斷的問題,如何對(duì)其改進(jìn)也是技術(shù)上的一大難點(diǎn)。

    圖5 難以連續(xù)分割的狹長小路Figure 5 Narrow path that is difficult to be divided

    (2)遙感影像質(zhì)量較差時(shí)(如存在云、霧覆蓋)難以自動(dòng)區(qū)劃,需人工輔助區(qū)劃。云霧覆蓋導(dǎo)致的誤判,需要借助深度學(xué)習(xí)才能有所緩解,下一步研究往深度學(xué)習(xí)方向拓展。

    (3)兩期影像配準(zhǔn)時(shí)要比較準(zhǔn)確,不能出現(xiàn)同一位置前后兩期影像偏差較大的情況。雖然在研究過程中對(duì)編寫的算法進(jìn)行了一定的偏移容差限制,但同一地物前后兩期遙感影像配準(zhǔn)偏差太大時(shí),容易導(dǎo)致誤判。

    猜你喜歡
    變化檢測(cè)工具欄格網(wǎng)
    用于遙感圖像變化檢測(cè)的全尺度特征聚合網(wǎng)絡(luò)
    “玩轉(zhuǎn)”西沃白板
    基于多尺度紋理特征的SAR影像變化檢測(cè)
    實(shí)時(shí)電離層格網(wǎng)數(shù)據(jù)精度評(píng)估
    基于稀疏表示的視網(wǎng)膜圖像對(duì)變化檢測(cè)
    基于Landsat影像的黃豐橋林場(chǎng)森林變化檢測(cè)研究
    基于空間信息格網(wǎng)與BP神經(jīng)網(wǎng)絡(luò)的災(zāi)損快速評(píng)估系統(tǒng)
    設(shè)計(jì)一種帶工具欄和留言功能的記事本
    平均Helmert空間重力異常格網(wǎng)構(gòu)制方法
    基于位置服務(wù)的地理格網(wǎng)編碼設(shè)計(jì)
    草草在线视频免费看| 欧美丝袜亚洲另类 | 老司机深夜福利视频在线观看| 老熟妇乱子伦视频在线观看| 神马国产精品三级电影在线观看| 五月玫瑰六月丁香| 一夜夜www| 成人午夜高清在线视频| 成人无遮挡网站| 亚洲人成伊人成综合网2020| 国产综合懂色| h日本视频在线播放| 日韩欧美 国产精品| 国产一区二区三区av在线 | 老熟妇仑乱视频hdxx| 精品久久久久久,| 亚洲美女黄片视频| 免费av毛片视频| 久久久久久久午夜电影| 99在线视频只有这里精品首页| 蜜桃亚洲精品一区二区三区| 国产精品亚洲美女久久久| 国产女主播在线喷水免费视频网站 | 国产激情偷乱视频一区二区| 一级av片app| 中文字幕av在线有码专区| 成人美女网站在线观看视频| 久久久久久久亚洲中文字幕| av在线老鸭窝| 亚洲三级黄色毛片| 联通29元200g的流量卡| 赤兔流量卡办理| 国产探花极品一区二区| 国产精品久久久久久久电影| 免费人成在线观看视频色| 久久九九热精品免费| 免费av毛片视频| 久久国产乱子免费精品| 九九在线视频观看精品| 午夜精品一区二区三区免费看| 内地一区二区视频在线| 熟妇人妻久久中文字幕3abv| 成人毛片a级毛片在线播放| 黄片wwwwww| 亚洲av成人av| 亚洲精品久久国产高清桃花| 极品教师在线视频| 日本在线视频免费播放| 国产精品亚洲美女久久久| 12—13女人毛片做爰片一| 国产精品av视频在线免费观看| 我要看日韩黄色一级片| 三级毛片av免费| 国产精品永久免费网站| 亚洲va日本ⅴa欧美va伊人久久| 91久久精品电影网| 国产国拍精品亚洲av在线观看| 久久天躁狠狠躁夜夜2o2o| 最新在线观看一区二区三区| 蜜桃亚洲精品一区二区三区| 搞女人的毛片| 午夜福利18| 有码 亚洲区| 欧美性感艳星| 国产单亲对白刺激| 联通29元200g的流量卡| 又紧又爽又黄一区二区| 亚洲精品成人久久久久久| 琪琪午夜伦伦电影理论片6080| 欧美xxxx黑人xx丫x性爽| 国产视频内射| 18+在线观看网站| 中文亚洲av片在线观看爽| 国产一区二区在线观看日韩| 中文字幕熟女人妻在线| 亚洲av免费高清在线观看| 黄色丝袜av网址大全| 欧美绝顶高潮抽搐喷水| 人人妻,人人澡人人爽秒播| 国产一区二区在线av高清观看| 性色avwww在线观看| 一本久久中文字幕| 国产aⅴ精品一区二区三区波| 亚洲av二区三区四区| 18禁黄网站禁片免费观看直播| 乱人视频在线观看| bbb黄色大片| 国产在线精品亚洲第一网站| 国产视频一区二区在线看| 3wmmmm亚洲av在线观看| 亚洲欧美日韩卡通动漫| 非洲黑人性xxxx精品又粗又长| av黄色大香蕉| 亚洲人成网站在线播放欧美日韩| 日韩欧美一区二区三区在线观看| 日日干狠狠操夜夜爽| 中出人妻视频一区二区| 国产精品1区2区在线观看.| 精品国产三级普通话版| 午夜精品在线福利| 18+在线观看网站| 亚洲真实伦在线观看| 国内毛片毛片毛片毛片毛片| 麻豆精品久久久久久蜜桃| 中文字幕av在线有码专区| 嫩草影视91久久| 色吧在线观看| 亚洲av免费高清在线观看| 能在线免费观看的黄片| 欧美高清成人免费视频www| 国产免费一级a男人的天堂| 91麻豆精品激情在线观看国产| 九九热线精品视视频播放| 国产真实乱freesex| 精品一区二区三区av网在线观看| 熟女电影av网| 亚州av有码| 国模一区二区三区四区视频| av在线老鸭窝| 亚洲成人免费电影在线观看| 亚洲国产欧洲综合997久久,| 精品无人区乱码1区二区| 国产伦精品一区二区三区四那| 亚洲va日本ⅴa欧美va伊人久久| 99国产极品粉嫩在线观看| 天堂√8在线中文| 真人一进一出gif抽搐免费| 亚洲第一电影网av| 日本-黄色视频高清免费观看| 久久中文看片网| 欧美成人一区二区免费高清观看| 中出人妻视频一区二区| 搡老熟女国产l中国老女人| 男女啪啪激烈高潮av片| 日本免费a在线| 干丝袜人妻中文字幕| 内射极品少妇av片p| 深夜a级毛片| 亚洲美女视频黄频| 国产精品三级大全| 黄色视频,在线免费观看| 国产精品三级大全| 一个人免费在线观看电影| 99久久九九国产精品国产免费| 小说图片视频综合网站| 国产高潮美女av| 亚洲av中文av极速乱 | 搡老妇女老女人老熟妇| 亚洲欧美日韩卡通动漫| 亚洲人与动物交配视频| 欧美极品一区二区三区四区| 免费在线观看日本一区| 亚洲中文字幕一区二区三区有码在线看| 免费av毛片视频| 日本 欧美在线| 如何舔出高潮| 三级男女做爰猛烈吃奶摸视频| 日韩欧美一区二区三区在线观看| 国产精品爽爽va在线观看网站| 国产日本99.免费观看| 黄色丝袜av网址大全| 老司机深夜福利视频在线观看| 亚洲中文字幕一区二区三区有码在线看| 日韩欧美一区二区三区在线观看| 午夜爱爱视频在线播放| 美女被艹到高潮喷水动态| 色哟哟哟哟哟哟| АⅤ资源中文在线天堂| 在线看三级毛片| 午夜免费男女啪啪视频观看 | 蜜桃亚洲精品一区二区三区| 一个人看的www免费观看视频| 久久婷婷人人爽人人干人人爱| а√天堂www在线а√下载| 欧美日韩亚洲国产一区二区在线观看| 少妇的逼水好多| 97碰自拍视频| 日本爱情动作片www.在线观看 | 中文字幕人妻熟人妻熟丝袜美| 亚洲欧美日韩东京热| 国产三级在线视频| 一区二区三区四区激情视频 | 日韩欧美三级三区| 国产精品一区www在线观看 | 欧美性感艳星| 免费看日本二区| 国产精品亚洲美女久久久| av国产免费在线观看| 色播亚洲综合网| 日韩欧美在线二视频| 亚洲av五月六月丁香网| 99久久中文字幕三级久久日本| 亚洲av一区综合| 久久草成人影院| 少妇熟女aⅴ在线视频| 黄色女人牲交| 1024手机看黄色片| 午夜视频国产福利| 婷婷精品国产亚洲av在线| 国产精品爽爽va在线观看网站| 少妇人妻精品综合一区二区 | 简卡轻食公司| 国产毛片a区久久久久| 99久久精品热视频| 很黄的视频免费| 欧美日韩亚洲国产一区二区在线观看| 成人毛片a级毛片在线播放| 欧美成人性av电影在线观看| 窝窝影院91人妻| 男人舔女人下体高潮全视频| 亚洲国产色片| 国产亚洲精品综合一区在线观看| 日本 av在线| 久久久久国内视频| 免费av观看视频| 91av网一区二区| or卡值多少钱| 亚洲第一电影网av| 97碰自拍视频| 精品久久国产蜜桃| 婷婷丁香在线五月| 12—13女人毛片做爰片一| 中文字幕熟女人妻在线| 内地一区二区视频在线| 亚洲 国产 在线| 搡老妇女老女人老熟妇| 亚洲一区高清亚洲精品| 中文资源天堂在线| 美女黄网站色视频| 热99re8久久精品国产| 成年版毛片免费区| 九九热线精品视视频播放| 欧美国产日韩亚洲一区| 日本五十路高清| 嫩草影院入口| 久久久久久久久中文| 亚洲欧美日韩无卡精品| 国产淫片久久久久久久久| 在线播放无遮挡| 久久这里只有精品中国| 女人十人毛片免费观看3o分钟| 亚洲avbb在线观看| 成人二区视频| av在线亚洲专区| 国产精品免费一区二区三区在线| 国产单亲对白刺激| 成人av一区二区三区在线看| 啪啪无遮挡十八禁网站| 日韩一区二区视频免费看| 黄片wwwwww| 简卡轻食公司| 免费观看在线日韩| 亚洲黑人精品在线| 91在线观看av| 国产高清激情床上av| 三级毛片av免费| 日本成人三级电影网站| 国产精品女同一区二区软件 | 国产欧美日韩精品亚洲av| 免费看日本二区| 12—13女人毛片做爰片一| 69av精品久久久久久| 免费人成在线观看视频色| 婷婷色综合大香蕉| 精品一区二区三区视频在线| 国产av麻豆久久久久久久| 国产高清视频在线观看网站| 国产爱豆传媒在线观看| 搞女人的毛片| 日本色播在线视频| 搡老妇女老女人老熟妇| 国产久久久一区二区三区| 看十八女毛片水多多多| 日韩中字成人| 亚洲aⅴ乱码一区二区在线播放| 国产v大片淫在线免费观看| 特级一级黄色大片| 亚洲美女黄片视频| 少妇裸体淫交视频免费看高清| 又爽又黄无遮挡网站| 天天一区二区日本电影三级| 久久久久久久久大av| www.www免费av| 成人综合一区亚洲| 国产蜜桃级精品一区二区三区| 非洲黑人性xxxx精品又粗又长| 小说图片视频综合网站| av视频在线观看入口| 欧美潮喷喷水| 在线免费观看的www视频| 夜夜夜夜夜久久久久| 国产一级毛片七仙女欲春2| 免费观看人在逋| 精品一区二区三区av网在线观看| 欧美bdsm另类| 免费看av在线观看网站| 99久久久亚洲精品蜜臀av| 99久久无色码亚洲精品果冻| 欧美人与善性xxx| 哪里可以看免费的av片| 成年女人毛片免费观看观看9| 亚洲一区二区三区色噜噜| 亚洲黑人精品在线| 欧美日韩瑟瑟在线播放| АⅤ资源中文在线天堂| 久久国产精品人妻蜜桃| 少妇丰满av| 久久精品国产亚洲网站| 久久久久久九九精品二区国产| 亚洲精品影视一区二区三区av| 国产伦精品一区二区三区四那| 精品一区二区免费观看| 99久国产av精品| 中文字幕av成人在线电影| 欧美成人一区二区免费高清观看| 黄色配什么色好看| 国内毛片毛片毛片毛片毛片| netflix在线观看网站| 亚洲av美国av| 亚洲不卡免费看| 日本与韩国留学比较| 久久久国产成人免费| 少妇高潮的动态图| 在线观看美女被高潮喷水网站| 亚洲av中文字字幕乱码综合| 天堂√8在线中文| 国产精品永久免费网站| 91av网一区二区| 99热6这里只有精品| 国产精品久久久久久久电影| 在线观看午夜福利视频| 色综合色国产| 精品免费久久久久久久清纯| 国产亚洲精品av在线| 日本精品一区二区三区蜜桃| 九色国产91popny在线| 色噜噜av男人的天堂激情| aaaaa片日本免费| 永久网站在线| 国产精华一区二区三区| 九九在线视频观看精品| 一区二区三区免费毛片| 国产毛片a区久久久久| 亚洲性夜色夜夜综合| 免费人成在线观看视频色| 99热网站在线观看| 直男gayav资源| 精品人妻视频免费看| 乱系列少妇在线播放| 亚洲人成网站在线播放欧美日韩| 亚洲黑人精品在线| 网址你懂的国产日韩在线| 国内毛片毛片毛片毛片毛片| 91在线观看av| 欧美一区二区亚洲| a级一级毛片免费在线观看| 国产成年人精品一区二区| 色尼玛亚洲综合影院| 91av网一区二区| aaaaa片日本免费| 欧美一区二区精品小视频在线| 日本一二三区视频观看| 国产视频一区二区在线看| 91久久精品国产一区二区成人| 伊人久久精品亚洲午夜| 18禁在线播放成人免费| 欧美+日韩+精品| 琪琪午夜伦伦电影理论片6080| 精品一区二区三区av网在线观看| 免费av毛片视频| 亚洲久久久久久中文字幕| 亚洲成人精品中文字幕电影| 国产69精品久久久久777片| 大型黄色视频在线免费观看| 国产精品一区www在线观看 | 能在线免费观看的黄片| 国产爱豆传媒在线观看| 国产真实乱freesex| 国产精品亚洲一级av第二区| 毛片一级片免费看久久久久 | 少妇的逼好多水| 99久久成人亚洲精品观看| 少妇被粗大猛烈的视频| 一卡2卡三卡四卡精品乱码亚洲| 2021天堂中文幕一二区在线观| 天天躁日日操中文字幕| 天天一区二区日本电影三级| 深爱激情五月婷婷| 中文资源天堂在线| 嫩草影院新地址| 中国美白少妇内射xxxbb| 九九热线精品视视频播放| 免费av观看视频| 国内毛片毛片毛片毛片毛片| 亚洲va日本ⅴa欧美va伊人久久| 国产伦人伦偷精品视频| 伦理电影大哥的女人| 女人十人毛片免费观看3o分钟| 欧美绝顶高潮抽搐喷水| 嫩草影院入口| 在线a可以看的网站| 久久久久九九精品影院| 身体一侧抽搐| 国产人妻一区二区三区在| 国内毛片毛片毛片毛片毛片| 欧美激情久久久久久爽电影| 亚洲五月天丁香| 欧美日韩综合久久久久久 | 免费无遮挡裸体视频| 午夜福利在线观看吧| 亚洲av二区三区四区| 国产成人影院久久av| 欧美日韩黄片免| 一本一本综合久久| 此物有八面人人有两片| 男人狂女人下面高潮的视频| 九九久久精品国产亚洲av麻豆| 国内久久婷婷六月综合欲色啪| 国产伦精品一区二区三区视频9| 天堂影院成人在线观看| 国产男人的电影天堂91| 国产精品人妻久久久久久| 舔av片在线| 国产免费男女视频| 久久天躁狠狠躁夜夜2o2o| 一区二区三区激情视频| 欧美色视频一区免费| 大型黄色视频在线免费观看| 一级av片app| 日本 av在线| 国产大屁股一区二区在线视频| 日韩欧美三级三区| 小蜜桃在线观看免费完整版高清| 国产精品女同一区二区软件 | 好男人在线观看高清免费视频| 国产高清视频在线播放一区| 别揉我奶头~嗯~啊~动态视频| 在线a可以看的网站| 成人午夜高清在线视频| 亚洲av不卡在线观看| 成人国产综合亚洲| 亚洲国产欧美人成| 成熟少妇高潮喷水视频| 国产精品综合久久久久久久免费| 国产高清有码在线观看视频| 亚洲无线在线观看| 99热这里只有是精品50| 国产主播在线观看一区二区| 成人二区视频| 午夜福利在线在线| 麻豆成人午夜福利视频| 国产亚洲精品综合一区在线观看| 国产一区二区在线av高清观看| 免费看a级黄色片| 女生性感内裤真人,穿戴方法视频| 久久久久久久亚洲中文字幕| or卡值多少钱| 国产色婷婷99| 狂野欧美激情性xxxx在线观看| 亚洲精品乱码久久久v下载方式| 亚洲综合色惰| 亚洲美女搞黄在线观看 | 国产精品电影一区二区三区| 色在线成人网| 日韩欧美在线二视频| 少妇被粗大猛烈的视频| 国产 一区 欧美 日韩| 久久精品久久久久久噜噜老黄 | 欧美不卡视频在线免费观看| 亚洲最大成人av| 天堂网av新在线| 成人特级av手机在线观看| 精品国内亚洲2022精品成人| 超碰av人人做人人爽久久| 国产大屁股一区二区在线视频| 国产一区二区在线av高清观看| 男女那种视频在线观看| 亚洲综合色惰| 中文字幕高清在线视频| 波多野结衣高清作品| 国产aⅴ精品一区二区三区波| 神马国产精品三级电影在线观看| 黄色视频,在线免费观看| or卡值多少钱| 国产精品日韩av在线免费观看| 乱系列少妇在线播放| 午夜日韩欧美国产| 亚洲国产精品sss在线观看| 免费看美女性在线毛片视频| 中出人妻视频一区二区| 欧美日韩国产亚洲二区| 美女高潮喷水抽搐中文字幕| 别揉我奶头~嗯~啊~动态视频| 麻豆一二三区av精品| 深夜a级毛片| 久久精品国产自在天天线| 久久亚洲真实| 国产在线精品亚洲第一网站| 日本黄大片高清| 久久婷婷人人爽人人干人人爱| aaaaa片日本免费| 在线a可以看的网站| 琪琪午夜伦伦电影理论片6080| 在线播放国产精品三级| 永久网站在线| 国产一区二区亚洲精品在线观看| 极品教师在线视频| 特级一级黄色大片| 亚洲精品粉嫩美女一区| 亚洲成a人片在线一区二区| 国产伦精品一区二区三区四那| 亚洲成人中文字幕在线播放| 中国美女看黄片| 欧美精品啪啪一区二区三区| 国产亚洲91精品色在线| 桃色一区二区三区在线观看| 成年版毛片免费区| 他把我摸到了高潮在线观看| 一进一出抽搐动态| 欧美绝顶高潮抽搐喷水| 亚洲熟妇熟女久久| 最近最新中文字幕大全电影3| 亚洲精品一区av在线观看| 一区二区三区免费毛片| 日日夜夜操网爽| 国产成人福利小说| 校园春色视频在线观看| 1000部很黄的大片| 麻豆成人午夜福利视频| 国内精品一区二区在线观看| 村上凉子中文字幕在线| 日本 欧美在线| 中文字幕久久专区| 欧美日韩综合久久久久久 | 中文字幕av成人在线电影| 欧美成人性av电影在线观看| 国产午夜福利久久久久久| 亚洲最大成人手机在线| 色在线成人网| 性欧美人与动物交配| 国产精品福利在线免费观看| 久久久久久久精品吃奶| 国产成人一区二区在线| 欧美黑人欧美精品刺激| 成人美女网站在线观看视频| 亚洲国产欧洲综合997久久,| 亚洲国产精品合色在线| 亚洲五月天丁香| 久久久午夜欧美精品| 少妇人妻精品综合一区二区 | 国产精品1区2区在线观看.| 久久久国产成人免费| 亚洲成人久久性| 亚洲avbb在线观看| 他把我摸到了高潮在线观看| 尾随美女入室| 少妇裸体淫交视频免费看高清| 不卡一级毛片| 国产精品美女特级片免费视频播放器| 欧美bdsm另类| 精品久久久久久久久久久久久| 淫秽高清视频在线观看| 久久精品国产亚洲av涩爱 | 免费av毛片视频| 99精品在免费线老司机午夜| 色综合站精品国产| 看片在线看免费视频| 俄罗斯特黄特色一大片| 午夜日韩欧美国产| 亚洲七黄色美女视频| 色综合色国产| 在线免费观看不下载黄p国产 | 欧美在线一区亚洲| 日本成人三级电影网站| 美女 人体艺术 gogo| 精品久久久久久久久亚洲 | 欧美日韩精品成人综合77777| 久久久久久久亚洲中文字幕| 亚洲精品粉嫩美女一区| 在线观看66精品国产| 国产亚洲精品久久久com| www日本黄色视频网| 日本在线视频免费播放| 国产欧美日韩精品亚洲av| 日本三级黄在线观看| 亚洲美女黄片视频| 国产精品伦人一区二区| 99久久久亚洲精品蜜臀av| 日韩高清综合在线| 夜夜看夜夜爽夜夜摸| 乱人视频在线观看| av在线老鸭窝| 少妇的逼好多水| 国产亚洲精品久久久久久毛片| 露出奶头的视频| 亚洲成av人片在线播放无| 一a级毛片在线观看| 成人性生交大片免费视频hd| 日本成人三级电影网站| 国产精品久久视频播放| 国产午夜精品论理片| 国产精品爽爽va在线观看网站| 久久久久免费精品人妻一区二区| 国产欧美日韩精品亚洲av| 中文字幕免费在线视频6| 精品久久久久久久久亚洲 | 日韩欧美国产在线观看| 成人三级黄色视频| 亚洲avbb在线观看| 久久久午夜欧美精品| 国产淫片久久久久久久久| 日本黄色视频三级网站网址| 日本爱情动作片www.在线观看 | 亚洲午夜理论影院| 国产精品av视频在线免费观看|