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

    基于決策樹方法的Landsat8 OLI影像紅樹林信息自動(dòng)提取

    2016-06-27 06:05:41張雪紅
    自然資源遙感 2016年2期
    關(guān)鍵詞:紅樹林決策樹

    張雪紅

    (南京信息工程大學(xué)地理與遙感學(xué)院,南京 210044)

    基于決策樹方法的Landsat8OLI影像紅樹林信息自動(dòng)提取

    張雪紅

    (南京信息工程大學(xué)地理與遙感學(xué)院,南京210044)

    摘要:基于廣西山口國家紅樹林生態(tài)自然保護(hù)區(qū)的Landsat 8 OLI 影像數(shù)據(jù),選用廣泛應(yīng)用于植被液態(tài)水含量反演的歸一化差值濕度指數(shù)(normalized difference moisture index,NDMI)和修正的歸一化差值池塘指數(shù)(modified normalized difference pond index, MNDPI)作為分類特征,運(yùn)用決策樹方法進(jìn)行紅樹林信息的自動(dòng)提取。研究結(jié)果表明: 紅樹林獨(dú)特的濱海濕地生境特點(diǎn),使其光譜同時(shí)包含植被和濕地信息; MNDPI和NDMI可分別反映可見光-近紅外波段反射率同短波紅外波段反射光譜的反差,可成功應(yīng)用于濕地植被信息的提取,能有效地將紅樹林同其他地物相區(qū)分; 采用Landsat8 OLI遙感數(shù)據(jù),并結(jié)合NDMI和MNDPI分類特征構(gòu)建的決策樹模型可有效地提取紅樹林信息,其錯(cuò)分率和漏分率都較低,分別為5.34%和1.69%。

    關(guān)鍵詞:紅樹林; NDMI; MNDPI; OLI; 決策樹

    0引言

    紅樹林是生長在熱帶、亞熱帶沿海潮間帶灘涂上特有的木本植物群落,屬常綠闊葉林,主要分布于淤泥深厚的海灣或河口鹽漬土壤上。紅樹林具有促淤固灘、防浪護(hù)堤、保持生物多樣性和凈化環(huán)境污染等作用[1-2]。然而由于紅樹林生存于海洋與陸地交錯(cuò)的生態(tài)脆弱帶,長期受到沿海不合理開發(fā)活動(dòng)和氣候變化的雙重威脅,過去的20a間全世界紅樹林的分布面積下降了36%[3]。為了加強(qiáng)紅樹林資源的管理和監(jiān)督,需要借助于遙感技術(shù)快速全面地獲得和掌握紅樹林時(shí)空分布信息。

    紅樹林的遙感識(shí)別與提取方法主要有監(jiān)督分類[4-5]、非監(jiān)督分類[6-7]、專家分類[8-9]、波段比值組合或植被指數(shù)法[4-5,7,10-11]、3S綜合法[12-13]等。在分類特征方面,通常采用基于原始反射特征以及歸一化植被指數(shù)(normalizeddifferencevegetationindex,NDVI)[4],(TM5-TM7)/(TM5+TM7)[10],TM3/TM5[11],TM5/TM4[11]和溫濕度指數(shù)[14]等指數(shù)或波段組合來區(qū)分紅樹林與非紅樹林。部分研究還結(jié)合了潮位、高程、海岸線等輔助地學(xué)數(shù)據(jù),以提高紅樹林的識(shí)別精度[15-17]。所采用的遙感影像數(shù)據(jù)根據(jù)研究目標(biāo)可分為兩大類: 宏觀動(dòng)態(tài)監(jiān)測類主要采用LandsatMSS[18],LandsatTM(ETM+)[10-12, 18-19],SPOT[12,20-22],ASTER[23]和HJ/CCD[24]等數(shù)據(jù); 另一類則從微觀尺度側(cè)重研究紅樹林的生態(tài)群落分布及演變規(guī)律,多采用航空影像(如CASI影像[11]、Hypmap影像[25])、IKONOS[26-28]及QuickBird[26,29]等高空間分辨率或高光譜影像數(shù)據(jù)。為了延續(xù)Landsat近40a的對(duì)地觀測,2013年2月12日,美國成功發(fā)射了Landsat8陸地衛(wèi)星,其搭載了OLI(operationallandimager)和TIRS(thermalinfraredsensor)2個(gè)有效載荷[30]。但是,到目前為止未見利用Landsat8OLI遙感數(shù)據(jù)來監(jiān)測紅樹林的相關(guān)文獻(xiàn)。

    歸一化差值濕度指數(shù)(normalizeddifferencemoistureindex,NDMI)[31]、歸一化差值水體指數(shù)(Normalizeddifferencewaterindex,NDWI)[32]能有效反映地物在近紅外和短波紅外波段反射光譜之間的反差,因此常用來反演植被中的液態(tài)水含量[31-35]。同樣,能夠有效描述可見光與短波紅外波段反射光譜之間反差的歸一化差值池塘指數(shù)(normalizeddifferencepondindex,NDPI)[36]常用來提取水體信息。由于紅樹林處于獨(dú)特的濱海濕地環(huán)境,使其光譜特征同時(shí)含有植被特征和濕地背景信息。為此,擬以廣西壯族自治區(qū)北海市山口國家紅樹林生態(tài)自然保護(hù)區(qū)為研究區(qū),基于Landsat8OLI影像數(shù)據(jù),嘗試結(jié)合NDMI和NDPI等光譜指數(shù),運(yùn)用決策樹方法來提取紅樹林,并以此來評(píng)價(jià)Landsat8OLI遙感數(shù)據(jù)對(duì)紅樹林的監(jiān)測能力和驗(yàn)證NDMI和NDPI等光譜指數(shù)在紅樹林遙感識(shí)別中的有效性。

    1研究區(qū)概況及數(shù)據(jù)源

    1.1研究區(qū)概況

    研究區(qū)位于廣西壯族自治區(qū)的北海市合浦縣東部的沙田半島。該半島東西兩側(cè)為山口國家紅樹林生態(tài)自然保護(hù)區(qū)。研究區(qū)所處的地理坐標(biāo)為E109°42′03″~109°45′36″,N21°29′00″~21°33′25″,東西長約6km,南北長約8km。研究區(qū)內(nèi)紅樹植物主要有木欖(Bruguieragymnorrhiza)、紅海欖(Rhizophorastylosa)、桐花樹(Aegicerascorniculatum)、白骨壤(Avicenniamarina)、海漆(Exoecariaagallocha)和秋茄(Kandeliacandel)等種類,其中連片的紅海欖純林在我國極為罕見。研究區(qū)屬于南亞熱帶季風(fēng)型海洋性氣候帶, 年平均氣溫22.9℃, 年平均降水量1 573.4mm[37]。區(qū)內(nèi)的地物類型主要包括桉樹人工林、陸生天然林、紅樹林、人工設(shè)施、光灘、灘涂以及水體等。

    1.2數(shù)據(jù)源

    研究所采用的數(shù)據(jù)包括Landsat8OLI影像、2013年10月22日和12月6日GoogleEarth高空間分辨率影像、地面控制點(diǎn)地理坐標(biāo)數(shù)據(jù)、紅樹林樣點(diǎn)的野外調(diào)查數(shù)據(jù)等。Landsat8衛(wèi)星包含OLI和TIRS這2個(gè)遙感器,OLI遙感器性能參數(shù)詳見表1。Landsat8OLI影像成像時(shí)間為2014年1月14日,軌道號(hào)為P124/R45,成像時(shí)的潮位為185cm(潮高基準(zhǔn)面在平均海平面下255cm)。2013年10月對(duì)該研究區(qū)進(jìn)行了野外調(diào)查,用GPS記錄了各樣點(diǎn)的地理位置,樣點(diǎn)數(shù)據(jù)總共25個(gè),其中地面控制點(diǎn)9個(gè)、紅樹林7個(gè)(注: 紅樹林樣點(diǎn)主要選取了保護(hù)區(qū)內(nèi)英羅保護(hù)站的紅樹林)、桉樹6個(gè)及木薯3個(gè),其空間分布詳見圖1。

    表1 Landsat8 OLI 波段設(shè)置及空間分辨率

    圖1野外調(diào)查樣點(diǎn)分布

    (全景圖為OLI標(biāo)準(zhǔn)假彩色圖像,局部放大圖為2009年

    10月16日GeoEye-1標(biāo)準(zhǔn)假彩色圖像)

    Fig.1Distributionoffieldsurveysites

    2處理方法

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

    數(shù)據(jù)預(yù)處理包括對(duì)Landsat8OLI影像的輻射校正和幾何糾正,通過ENVI4.8軟件實(shí)現(xiàn)。其中輻射校正包括輻射定標(biāo)和地表反射率的計(jì)算,首先利用影像元數(shù)據(jù)文件提供的定標(biāo)系數(shù)將原始DN值轉(zhuǎn)換成表觀反射率,然后采用改進(jìn)的暗目標(biāo)方法計(jì)算地表反射率[38-39]。幾何糾正保證誤差小于0.5個(gè)像元。影像圖采用UTM投影。

    2.2分類特征及決策樹方法

    短波紅外對(duì)液態(tài)水敏感,1.20μm為液態(tài)水的吸收帶之一[32]。有學(xué)者基于以上波段特征提出了NDMI[31]和NDWI[32]光譜指數(shù),分別定義為

    NDMI=(ρNIR-ρMIR)/(ρNIR+ρMIR) ,

    (1)

    NDWI=(ρ0.86-ρ1.24)/(ρ0.86+ρ1.24) ,

    (2)

    式中: ρNIR,ρMIR,ρ0.86和ρ1.24分別為近紅外、短波紅外、0.86μm和1.24μm波段反射率。

    Lacaux等[36]提出了NDPI來識(shí)別和評(píng)價(jià)池塘水體,其定義為

    NDPI=(ρMIR-ρG)/(ρMIR+ρG) ,

    (3)

    式中ρG和ρMIR分別為綠波段和短波紅外波段反射率。

    研究區(qū)大面積分布的紅海欖純林的光譜反射率明顯低于其他紅樹林品種,尤其是短波紅外波段反射率小于綠光波段。因此本文通過對(duì)NDPI做修正,用紅光波段反射率來替換NDPI中的綠光波段反射率,修正后的NDPI稱之謂MNDPI(ModifiedNDPI),在決策樹建模時(shí)用以替換NDPI,即

    MNDPI=(ρMIR-ρR)/(ρMIR+ρR) ,

    (4)

    式中ρR和ρMIR分別為紅波段和對(duì)應(yīng)短波紅外波段反射率。

    3結(jié)果與分析

    3.1典型地物光譜特征

    為了便于比較研究區(qū)內(nèi)典型地物的光譜特征,基于野外典型地物樣點(diǎn)的調(diào)查數(shù)據(jù),結(jié)合GoogleEarth高空間分辨率影像,通過目視解譯對(duì)區(qū)內(nèi)典型地物各選取了50~100個(gè)像元樣本(表2)。

    表2 典型地物光譜特征統(tǒng)計(jì)值

    圖2為研究區(qū)內(nèi)紅樹林、陸生植被、人工建筑、水體、灘涂、光灘、陸生植被與水體的混合像元、裸土

    圖2 研究區(qū)典型地物反射光譜曲線

    等典型地物平均反射光譜曲線,圖中的“混合像元”為陸生植被與水體的混合像元。陸生植被與紅樹林光譜曲線具有典型的植被反射率特征。紅樹林由于處于濱海濕地環(huán)境,其短波紅外反射率明顯低于陸生植被。陸生植被與水體的混合像元反射率光譜曲線形狀同植被極其相似,不同之處在于近紅外反射率顯著偏小。而灘涂由于生長著海草,且在影像數(shù)據(jù)成像時(shí)處于低潮,因此光譜曲線也反映出部分植被的特征,只是可見光區(qū)域反射率高于植被,而紅外區(qū)域反射率則明顯低于植被。同植被相比,人工建筑、光灘、裸土及水體等非植被地物類型的光譜沒有可見光波段植被的“兩谷一峰一陡坡”特征。

    基于典型地物選取的像元樣本數(shù)據(jù),分別統(tǒng)計(jì)了各地物類型的NDMI和MNDPI等光譜特征。綜合圖2和表2,可以得出:

    1)紅樹林、水體、灘涂、光灘的NDMI明顯大于其他地物類型,陸生植被以及陸生植被與水體的混合像元次之,人工建筑和裸土則最小。從而利用NDMI可以將紅樹林同人工建筑、裸土相區(qū)分,也能將絕大大部分紅樹林同陸生植被與水體的混合像元相區(qū)分,而紅樹林同陸生植被則有一定程度的重疊區(qū)域。

    2)陸生植被的MNDPI最大,紅樹林、人工建筑、裸土以及陸生植被與水體的混合像元次之,而水體、灘涂和光灘最小。因此,利用MNDPI可以將紅樹林同水體、灘涂和光灘,尤其是同陸生植被進(jìn)行區(qū)分。

    3.2決策樹模型的構(gòu)建

    采用決策樹分類方法來提取紅樹林空間分布信息時(shí),將研究區(qū)域分為紅樹林和非紅樹林兩類,并基于ENVI4.8軟件中決策樹分類模塊進(jìn)行紅樹林的識(shí)別。決策規(guī)則如圖3所示。

    圖3 紅樹林識(shí)別的決策樹模型

    規(guī)則①: 0≤MNDPI≤0.35,此規(guī)則將地物分成2類(一類為水體、陸生植被、灘涂、光灘,另一類為紅樹林、人工建筑、裸土、陸生植被與水體的混合像元); 規(guī)則②: NDMI≥0.42,此規(guī)則可將紅樹林和人工建筑、裸土、陸生植被與水體的混合像元有效區(qū)分。滿足此規(guī)則的為紅樹林,否則為人工建筑、裸土或陸生植被與水體的混合像元。文中MNDPI和NDMI的閾值是基于樣本數(shù)據(jù)經(jīng)過反復(fù)分析而確定的。

    3.3紅樹林識(shí)別結(jié)果及精度評(píng)價(jià)

    基于構(gòu)建的決策樹模式,對(duì)研究區(qū)內(nèi)紅樹林進(jìn)行識(shí)別試驗(yàn),其識(shí)別結(jié)果見圖4,其中紅色區(qū)域?yàn)榧t樹林分布區(qū)。底圖為進(jìn)行對(duì)比度拉伸后的OLI影像真彩色合成圖像。

    圖4 2014年山口紅樹林自然保護(hù)區(qū)紅樹林分布圖

    從圖像上隨機(jī)提取2 249個(gè)像元作為驗(yàn)證樣本用于進(jìn)行精度評(píng)價(jià),其中紅樹林415個(gè),陸生植被、水體、灘涂、光灘共763個(gè),人工建筑、陸生植被與水體的混合像元、裸土共1 071個(gè)。利用ENVI4.8軟件構(gòu)建混淆矩陣,并計(jì)算紅樹林識(shí)別精度。結(jié)果表明,錯(cuò)分率為5.34%,漏分率為1.69%。其中錯(cuò)分為紅樹林的主要為山區(qū)陰影植被,這是因?yàn)橹脖辉陉幱皡^(qū)各波段的反射率急劇減小,尤其是短波紅外波段減小幅度最大,從而導(dǎo)致MNDPI減小而NDMI增大; 另外還有少量沿海養(yǎng)殖水體同植被混合像元也極易錯(cuò)分為紅樹林。漏分的紅樹林主要為紅樹林斑塊邊緣且地勢(shì)略低的稀疏紅樹林,這些紅樹林多為先鋒樹種白骨壤,其樹葉顏色呈現(xiàn)為灰綠色,再加上潮汐的周期性浸淹,樹葉上常粘有灰色的潮土,導(dǎo)致其近紅外段反射率相對(duì)稠密紅樹林偏小,而短波紅外段反射率則略微偏大,于是NDMI總體偏小,故容易漏分。

    4結(jié)論

    基于Landsat8OLI影像數(shù)據(jù),以廣西山口國家紅樹林生態(tài)自然保護(hù)區(qū)為例,通過選擇NDMI和MNDPI等光譜指數(shù),采用決策樹方法進(jìn)行紅樹林信息的自動(dòng)提取研究,研究結(jié)果如下:

    1)紅樹林生長于特殊的濱海濕地環(huán)境,其光譜同時(shí)包含植被和濕地信息,而且其短波紅外波段反射率明顯低于陸地植被。

    2)NDMI和MNDPI指數(shù)能分別反映可見光-近紅外同短波紅外波段光譜的反差,應(yīng)用于紅樹林濕地植被信息提取可有效地將紅樹林同其他地物區(qū)分。其中NDMI可以將紅樹林同人工建筑、裸土分離,也可以將絕大部分紅樹林同陸生植被與水體的混合像元區(qū)分,而紅樹林同陸生植被會(huì)有一定程度的重疊區(qū)域; 陸生植被同紅樹林存在明顯的“同譜異物”現(xiàn)象,MNDPI不但可將紅樹林同水體、灘涂和光灘區(qū)分,而且還可以有效地將其同陸生植被分離。

    3)Landsat8OLI能有效監(jiān)測紅樹林的空間分布特征?;贚andsat8OLI影像數(shù)據(jù),結(jié)合NDMI和MNDPI分類特征,通過選取相應(yīng)的閾值并構(gòu)建決策樹模型,可以有效地提取紅樹林信息。

    為了進(jìn)一步提高紅樹林的識(shí)別精度,在后續(xù)的研究中,還需著重關(guān)注陰影植被、陸生植被與沿海養(yǎng)殖水體的混合像元同紅樹林誤分現(xiàn)象及稀疏紅樹林的漏分問題,進(jìn)一步完善紅樹林信息提取的方法流程。此外,本文中的分類特征閾值是否具有時(shí)空拓展性,還有待選用不同時(shí)間和不同區(qū)域的數(shù)據(jù)來進(jìn)一步驗(yàn)證或修正。

    志謝: 文中Landsat8OLI影像由中國科學(xué)院計(jì)算機(jī)網(wǎng)絡(luò)信息中心地理空間數(shù)據(jù)云 (http: //www.gscloud.cn)提供,衷心感謝各位評(píng)審專家所提的寶貴意見。

    參考文獻(xiàn)(References):

    [1]BlascoF,SaengerP,JanodetE.Mangrovesasindicatorsofcoastalchange[J].Catena,1996,27(3/4):167-178.

    [2]LiuK,LiX,ShiX,etal.MonitoringmangroveforestchangesusingremotesensingandGISdatawithdecision-treelearning[J].Wetlands,2008,28(2):336-346.

    [3]FAO.TheWorld’sMangroves1980—2005,FAOForestryPaperNo.153[R/OL].FAO:Rome,Italy,2007;Availableonline:ftp://ftp.fao.org/docrep/fao/010/a1427e/a1427e00.pdf(accessedon9October2014).

    [4]JensenJR,LinHY,YangXH,etal.Themeasurementofmangrovecharacteristicsinsouth-westFloridausingspotmultispectraldata[J].GeocartoInternational,1991,6(2):13-21.

    [5]RasolofoharinoroM,BlascoF,BellanMF,etal.AremotesensingbasedmethodologyformangrovestudiesinMadagascar[J].InternationalJournalofRemoteSensing,1998,19(10):1873-1886.

    [6]FerreiraMA,AndradeF,BandeiraSO,etal.Analysisofcoverchange(1995—2005)ofTanzania/Mozambiquetrans-boundarymangrovesusingLandsatimagery[J].AquaticConservation:MarineandFreshwaterEcosystems,2009,19(S1):S38-S45.

    [7]LeeTM,YehHC.Applyingremotesensingtechniquestomonitorshiftingwetlandvegetation:AcasestudyofDanshuiRiverestuarymangrovecommunities,Taiwan[J].EcologicalEngineering,2009,35(4):487-496.

    [8]黎夏,劉凱,王樹功.珠江口紅樹林濕地演變的遙感分析[J].地理學(xué)報(bào),2006,61(1):26-34.

    LiX,LiuK,WangSG.MangrovewetlandchangesinthePearlRiverestuaryusingremotesensing[J].ActaGeographicaSinica,2006,61(1):26-34.

    [9]張雪紅.基于知識(shí)與規(guī)則的紅樹林遙感信息提取[J].南京信息工程大學(xué)學(xué)報(bào),2011,3(4):341-345.

    ZhangXH.Remotesensinginformationextractionofmangrovebasedonknowledgeandrules[J].JournalofNanjingUniversityofInformationScience&Technology,2011,3(4):341-345.

    [10]ChaudhuryMU.DigitalanalysisofremotesensingdataformonitoringtheecologicalstatusofthemangroveforestsofSunderbansinBangladesh[C]//Proceedingsofthe23rdInternationalSymposiumonRemoteSensingoftheEnvironment.Bangkok,Thailand,1990,1:493-497.

    [11]GreenEP,ClarkCD,MumbyPJ,etal.Remotesensingtechniquesformangrovemapping[J].InternationalJournalofRemoteSensing,1998,19(5):935-956.

    [12]AschbacherJ,OfrenR,DelsolJP,etal.AnintegratedcomparativeapproachtomangrovevegetationmappingusingadvancedremotesensingandGIStechnologies:Preliminaryresults[J].Hydrobiologia,1995,295(1/3):285-294.

    [13]LongBG,SkewesTD.GISandremotesensingimprovesmangrovemapping[C]//FloreatWA.7thAustralasianRemoteSensingConference,Volume1.Melbourne:RemoteSensingandPhotogrammetryAssociationAustraliaLtd.,1994:545-551.

    [14]張雪紅,田慶久.利用溫濕度指數(shù)提高紅樹林遙感識(shí)別精度[J].國土資源遙感,2012,24(3):65-70.doi:10.6046/gtzyyg.2012.03.13.

    ZhangXH,TianQJ.Applicationofthetemperature-moistureindextotheimprovementofremotesensingidentificationaccuracyofmangrove[J].RemoteSensingforLandandResources,2012,24(3):65-70.doi:10.6046/gtzyyg.2012.03.13.

    [15]張雪紅,周杰,魏瑗瑗,等.結(jié)合潮位信息的紅樹林遙感識(shí)別[J].南京信息工程大學(xué)學(xué)報(bào),2013,5(6):501-507.

    ZhangXH,ZhouJ,WeiYY,etal.Remotesensingidentificationofmangroveforestcombinedtidallevelinformation[J].JournalofNanjingUniversityofInformationScienceandTechnology,2013,5(6):501-507.

    [16]劉凱,黎夏,王樹功,等.珠江口近20年紅樹林濕地的遙感動(dòng)態(tài)監(jiān)測[J].熱帶地理,2005,25(2):111-116.

    LiuK,LiX,WangSG,etal.MonitoringofthechangesofmangrovewetlandaroundtheZhujiangestuaryinthepasttwodecadesbyremotesensing[J].TropicalGeography,2005,25(2):111-116.

    [17]LiuK,LiX,ShiX,etal.MonitoringmangroveforestchangesusingremotesensingandGISdatawithdecision-treelearning[J].Wetlands,2008,28(2):336-346.

    [18]GiriC,PengraB,ZhuZL,etal.MonitoringmangroveforestdynamicsofthesundarbansinbangladeshandIndiausingmulti-temporalsatellitedatafrom1973to2000[J].Estuarine,CoastalandShelfScience,2007,73(1/2):91-100.

    [19]RakotomavoA,FromardF.DynamicsofmangroveforestsintheMangokyRiverDelta,Madagascar,undertheinfluenceofnaturalandhumanfactors[J].ForestEcologyandManagement,2010,259(6):1161-1169.

    [20]BlascoF,AizpuruM,GersC.DepletionofthemangrovesofcontinentalAsia[J].WetlandsEcologyandManagement,2001,9(3):255-266.

    [21]ThuPM,PopulusJ.StatusandchangesofmangroveforestinMekongdelta:casestudyinTraVinh,Vietnam[J].Estuarine,CoastalandShelfScience,2007,71(1/2):98-109.

    [22]ConcheddaG,DurieuxL,MayauxP.Anobject-basedmethodformappingandchangeanalysisinmangroveecosystems[J].ISPRSJournalofPhotogrammetryandRemoteSensing,2008,63(5):578-589.

    [23]VaiphasaC,SkidmoreAK,DeBoerWF.Apost-classifierformangrovemappingusingecologicaldata[J].ISPRSJournalofPhotogrammetryandRemoteSensing,2006,61(1):1-10.

    [24]李?yuàn)檴?田慶久.北部灣紅樹林的HJ星多光譜遙感提取模式[J].中國科學(xué):信息科學(xué),2011,41(S1):153-165.

    LiSS,TianQJ.MangroveextractionmodebasedonHJmultispectralremotesensinginBeibuGulf[J].ScientiaSinica:Informationis,2011,41(S1):153-165.

    [25]肖海燕,曾輝,昝啟杰,等.基于高光譜數(shù)據(jù)和專家決策法提取紅樹林群落類型信息[J].遙感學(xué)報(bào),2007,11(4):531-537.

    XiaoHY,ZengH,ZanQJ,etal.Decisiontreemodelinextractionofmangrovecommunityinformationusinghyperspectralimagedata[J].JournalofRemoteSensing,2007,11(4):531-537.

    [26]WangL,SousaWP,GongP,etal.ComparisonofIKONOSandQuickBirdimagesformappingmangrovespeciesontheCaribbeancoastofPanama[J].RemoteSensingofEnvironment,2004,91(3/4):432-440.

    [27]KovacsJM,WangJF,Flores-VerdugoF.MappingmangroveleafareaindexatthespecieslevelusingIKONOSandLAI-2000sensorsfortheAguaBravaLagoon,MexicanPacific[J].Estuarine,CoastalandShelfScience,2005,62(1/2):377-384.

    [28]ProisyC,CouteronP,FromardF.PredictingandmappingmangrovebiomassfromcanopygrainanalysisusingFourier-basedtexturalordinationofIKONOSimages[J].RemoteSensingofEnvironment,2007,109(3):379-392.

    [29]SalehMA.AssessmentofmangrovevegetationonAbuMinqarislandoftheRedSea[J].JournalofAridEnvironments,2007,68(2):331-336.

    [30]徐涵秋,唐菲.新一代Landsat系列衛(wèi)星:Landsat8遙感影像新增特征及其生態(tài)環(huán)境意義[J].生態(tài)學(xué)報(bào),2013,33(11):3249-3257.XuHQ,TangF.AnalysisofnewcharacteristicsofthefirstLandsat8imageandtheireco-environmentalsignificance[J].ActaEcologicaSinica,2013,33(11):3249-3257.

    [31]WilsonEH,SaderSA.DetectionofforestharvesttypeusingmultipledatesofLandsatTMimagery[J].RemoteSensingofEnvironment,2002,80(3):385-396.

    [32]GaoBC.NDWI-anormalizeddifferencewaterindexforremotesensingofvegetationliquidwaterfromspace[J].RemoteSensingofEnvironment,1996,58(3):257-266.

    [33]JacksonTJ,ChenDY,CoshM,etal.VegetationwatercontentmappingusingLandsatdataderivednormalizeddifferencewaterindexforcornandsoybeans[J].RemoteSensingofEnvironment,2004,92(4):475-482.

    [34]ChenDY,HuangJF,JacksonTJ.VegetationwatercontentestimationforcornandsoybeansusingspectralindicesderivedfromMODISnear-andshort-waveinfraredbands[J].RemoteSensingofEnvironment,2005,98(2/3):225-236.

    [35]HuangJF,ChenDY,CoshMH.Sub-pixelreflectanceunmixinginestimatingvegetationwatercontentanddrybiomassofcornandsoybeanscroplandusingnormalizeddifferencewaterindex(NDWI)fromsatellites[J].InternationalJournalofRemoteSensing,2009,30(8):2075-2104.

    [36]LacauxJP,TourreYM,VignollesC,etal.Classificationofpondsfromhigh-spatialresolutionremotesensing:Applicationtoriftvalleyfeverepidemicsinsenegal[J].RemoteSensingofEnvironment,2007,106(1):66-74.

    [37]梁士楚.廣西英羅灣紅樹植物群落的研究[J].植物生態(tài)學(xué)報(bào),1996,20(4):310-321.

    LiangSC.StudiesonthemangrovecommunitiesinYingluoBayofGuangxi[J].ActaPhytoecologicaSinica,1996,20(4):310-321.

    [38]ChavezJrPS.Animproveddark-objectsubtractiontechniqueforatmosphericscatteringcorrectionofmultispectraldata[J].RemoteSensingofEnvironment,1988,24(3):459-479.

    [39]ChavezPS.Image-basedatmosphericcorrections:Revisitedandimproved[J].PhotogrammetricEngineeringandRemoteSensing,1996,62(9):1025-1036.

    (責(zé)任編輯: 李瑜)

    DecisiontreealgorithmofautomaticallyextractingmangroveforestsinformationfromLandsat8OLIimagery

    ZHANGXuehong

    (School of Geography and Remote Sensing, Nanjing University of Information Science & Technology, Nanjing 210044, China)

    Abstract:NDMI (normalized difference moisture index) is widely used to assess and retrieve vegetation liquid water content. In this study, decision tree method was employed to automatically extract mangrove forests information combining the NDMI and MNDPI (modified normalized difference pond index), modified according to the mangrove characteristics, with Landsat8 OLI imagery acquired at Shankou mangrove national ecosystem nature reserve in Guangxi. The research results show that mangrove forests spectra consist of vegetation and wetland characteristics due to the unique near-shore coastal habitat of mangrove forests. MNDPI and NDMI can represent the spectral contrast between shortwave infrared region and visible region, near infrared region respectively. Therefore, the two spectral indices can successfully be employed to extract wetland vegetation and effectively discriminate mangrove forests from other land cover types. The decision tree method effectively extracted mangrove forests information by combining the classification features of MNDPI and NDMI and using Landsat8 OLI remotely sensed data. The commission error and omission error of mangrove forests were 5.34% and 1.69% respectively.

    Keywords:mangrove forests; NDMI; MNDPI; OLI; decision tree

    doi:10.6046/gtzyyg.2016.02.28

    收稿日期:2014-11-12;

    修訂日期:2015-03-02

    基金項(xiàng)目:國家自然科學(xué)基金項(xiàng)目“紅樹林冠層高光譜探測及群落類型識(shí)別研究”(編號(hào): 41201461)和江蘇政府留學(xué)獎(jiǎng)學(xué)金共同資助。

    中圖法分類號(hào):TP 751.1; TP 753

    文獻(xiàn)標(biāo)志碼:A

    文章編號(hào):1001-070X(2016)02-0182-06

    作者簡介:張雪紅(1980- ), 男 ,博士,副教授,主要從事植被生態(tài)遙感方面的研究。 Email: zxhbnu@126.com。

    引用格式: 張雪紅.基于決策樹方法的Landsat8OLI影像紅樹林信息自動(dòng)提取[J].國土資源遙感,2016,28(2):182-187.(ZhangXH.DecisiontreealgorithmofautomaticallyextractingmangroveforestsinformationfromLandsat8OLIimagery[J].RemoteSensingforLandandResources,2016,28(2):182-187.)

    猜你喜歡
    紅樹林決策樹
    藏著寶藏的紅樹林
    海岸衛(wèi)士——紅樹林
    幼兒園(2021年4期)2021-07-28 07:38:04
    紅樹林為何能在海水中生長
    一種針對(duì)不均衡數(shù)據(jù)集的SVM決策樹算法
    神奇的紅樹林
    決策樹和隨機(jī)森林方法在管理決策中的應(yīng)用
    電子制作(2018年16期)2018-09-26 03:27:06
    走過紅樹林
    歌海(2018年4期)2018-05-14 12:46:15
    基于改進(jìn)決策樹的故障診斷方法研究
    基于決策樹的出租車乘客出行目的識(shí)別
    基于決策樹的復(fù)雜電網(wǎng)多諧波源監(jiān)管
    電測與儀表(2016年2期)2016-04-12 00:24:40
    日韩大尺度精品在线看网址| 精品欧美国产一区二区三| 又粗又爽又猛毛片免费看| www.999成人在线观看| 制服诱惑二区| 熟妇人妻久久中文字幕3abv| 国产精品日韩av在线免费观看| 国产久久久一区二区三区| 欧美高清成人免费视频www| 成人午夜高清在线视频| 麻豆一二三区av精品| 亚洲av美国av| 大型黄色视频在线免费观看| 亚洲自拍偷在线| 久久亚洲精品不卡| 神马国产精品三级电影在线观看 | 亚洲中文字幕日韩| 女警被强在线播放| 国产不卡一卡二| 久久欧美精品欧美久久欧美| 男人舔女人下体高潮全视频| 久久久久久久久久黄片| 久久中文字幕一级| 精品久久久久久久毛片微露脸| 久久久久久久午夜电影| 日韩av在线大香蕉| 色播亚洲综合网| 国内毛片毛片毛片毛片毛片| 我要搜黄色片| 777久久人妻少妇嫩草av网站| 国产精品 国内视频| 长腿黑丝高跟| 91九色精品人成在线观看| 亚洲男人天堂网一区| aaaaa片日本免费| 午夜a级毛片| 日韩国内少妇激情av| 中文字幕精品亚洲无线码一区| 亚洲狠狠婷婷综合久久图片| 18美女黄网站色大片免费观看| 午夜精品在线福利| 成人三级做爰电影| 看免费av毛片| www日本在线高清视频| 国产精品99久久99久久久不卡| 精品一区二区三区视频在线观看免费| 日韩成人在线观看一区二区三区| 波多野结衣高清作品| 极品教师在线免费播放| 精品欧美国产一区二区三| 国产97色在线日韩免费| bbb黄色大片| 久久精品国产99精品国产亚洲性色| 伦理电影免费视频| 真人做人爱边吃奶动态| 日韩中文字幕欧美一区二区| 精华霜和精华液先用哪个| 美女午夜性视频免费| 黄片小视频在线播放| 看黄色毛片网站| 国内精品久久久久精免费| 国产高清视频在线播放一区| 久久久精品国产亚洲av高清涩受| 老司机深夜福利视频在线观看| 老熟妇仑乱视频hdxx| 亚洲欧美激情综合另类| √禁漫天堂资源中文www| 免费电影在线观看免费观看| 可以在线观看毛片的网站| 国产成人av教育| 麻豆久久精品国产亚洲av| 亚洲精品国产一区二区精华液| www日本黄色视频网| 欧美精品啪啪一区二区三区| 欧美国产日韩亚洲一区| 中亚洲国语对白在线视频| 99在线人妻在线中文字幕| 一级黄色大片毛片| 女人被狂操c到高潮| 免费看a级黄色片| 桃色一区二区三区在线观看| 蜜桃久久精品国产亚洲av| 每晚都被弄得嗷嗷叫到高潮| 最新在线观看一区二区三区| 天天一区二区日本电影三级| 中文字幕熟女人妻在线| 午夜激情av网站| 亚洲欧美日韩高清在线视频| 69av精品久久久久久| ponron亚洲| 日韩欧美在线二视频| 亚洲av日韩精品久久久久久密| 久久伊人香网站| 一夜夜www| 国产一区二区三区视频了| 午夜福利视频1000在线观看| 不卡av一区二区三区| 精品免费久久久久久久清纯| 亚洲成人久久爱视频| 久久性视频一级片| x7x7x7水蜜桃| 久久热在线av| 日韩中文字幕欧美一区二区| 国产亚洲精品第一综合不卡| 日韩三级视频一区二区三区| 这个男人来自地球电影免费观看| 色尼玛亚洲综合影院| 搡老熟女国产l中国老女人| 中文字幕高清在线视频| 黄色 视频免费看| 精品福利观看| 麻豆久久精品国产亚洲av| 日本黄大片高清| 国产熟女午夜一区二区三区| 天天添夜夜摸| 日本三级黄在线观看| 色哟哟哟哟哟哟| 亚洲男人的天堂狠狠| 最好的美女福利视频网| 老司机午夜福利在线观看视频| 女生性感内裤真人,穿戴方法视频| 国产成年人精品一区二区| 嫁个100分男人电影在线观看| 操出白浆在线播放| 99精品在免费线老司机午夜| 精品国产超薄肉色丝袜足j| 五月伊人婷婷丁香| 岛国在线免费视频观看| 成人亚洲精品av一区二区| 亚洲av成人一区二区三| 人人妻,人人澡人人爽秒播| 岛国视频午夜一区免费看| 欧美zozozo另类| 日本三级黄在线观看| 亚洲成a人片在线一区二区| 中出人妻视频一区二区| 国产av麻豆久久久久久久| 欧美日韩黄片免| 他把我摸到了高潮在线观看| 18禁裸乳无遮挡免费网站照片| 欧美日韩亚洲国产一区二区在线观看| 岛国在线观看网站| 在线观看一区二区三区| а√天堂www在线а√下载| 精品无人区乱码1区二区| 99精品在免费线老司机午夜| 三级国产精品欧美在线观看 | 99国产极品粉嫩在线观看| 国产一区二区在线观看日韩 | 欧美zozozo另类| 搡老妇女老女人老熟妇| x7x7x7水蜜桃| 50天的宝宝边吃奶边哭怎么回事| 熟妇人妻久久中文字幕3abv| 中文字幕人成人乱码亚洲影| 亚洲av熟女| 国产视频内射| 日韩国内少妇激情av| 日韩欧美精品v在线| 中文字幕熟女人妻在线| 三级毛片av免费| 精品高清国产在线一区| 波多野结衣高清无吗| 99在线人妻在线中文字幕| 一夜夜www| 亚洲成a人片在线一区二区| 日日爽夜夜爽网站| 每晚都被弄得嗷嗷叫到高潮| 亚洲熟妇中文字幕五十中出| 嫁个100分男人电影在线观看| 中文在线观看免费www的网站 | 亚洲成人国产一区在线观看| 桃红色精品国产亚洲av| 日本一本二区三区精品| 久久久久久免费高清国产稀缺| 国产av一区二区精品久久| 男女视频在线观看网站免费 | 在线观看免费视频日本深夜| 国产高清激情床上av| 69av精品久久久久久| 免费看日本二区| 高清在线国产一区| 脱女人内裤的视频| 日韩成人在线观看一区二区三区| 国产精品久久久久久人妻精品电影| 可以免费在线观看a视频的电影网站| 搡老熟女国产l中国老女人| 午夜福利在线在线| 国产黄a三级三级三级人| 一本一本综合久久| 九色国产91popny在线| 欧美乱码精品一区二区三区| 亚洲自拍偷在线| 亚洲激情在线av| 色播亚洲综合网| 人妻久久中文字幕网| 久久精品亚洲精品国产色婷小说| 日韩精品青青久久久久久| 午夜福利18| 成在线人永久免费视频| 国产69精品久久久久777片 | 国产日本99.免费观看| 亚洲 欧美一区二区三区| 日本免费一区二区三区高清不卡| 久久中文字幕人妻熟女| 狂野欧美白嫩少妇大欣赏| 99国产精品99久久久久| 国产亚洲精品久久久久久毛片| 亚洲电影在线观看av| 成人av在线播放网站| 久久天躁狠狠躁夜夜2o2o| 国产午夜精品久久久久久| 亚洲精品久久成人aⅴ小说| av天堂在线播放| 亚洲精品久久国产高清桃花| 啪啪无遮挡十八禁网站| 99国产综合亚洲精品| 一个人免费在线观看电影 | 国产精品久久视频播放| 可以在线观看的亚洲视频| 免费看日本二区| 国产精品永久免费网站| 99国产精品99久久久久| 成人高潮视频无遮挡免费网站| 不卡一级毛片| 村上凉子中文字幕在线| 免费电影在线观看免费观看| 一区福利在线观看| 天堂动漫精品| 国产亚洲欧美在线一区二区| 九色成人免费人妻av| 大型av网站在线播放| 久99久视频精品免费| 高清在线国产一区| 无限看片的www在线观看| 亚洲中文字幕一区二区三区有码在线看 | 999精品在线视频| 露出奶头的视频| 久久久久精品国产欧美久久久| 久热爱精品视频在线9| 日韩精品免费视频一区二区三区| 免费在线观看成人毛片| 亚洲电影在线观看av| 51午夜福利影视在线观看| 黑人欧美特级aaaaaa片| 色噜噜av男人的天堂激情| 99久久国产精品久久久| 精品一区二区三区视频在线观看免费| 最近最新免费中文字幕在线| 老司机深夜福利视频在线观看| 男女床上黄色一级片免费看| svipshipincom国产片| 亚洲av成人不卡在线观看播放网| 午夜精品在线福利| 国产精品亚洲一级av第二区| 99国产精品一区二区蜜桃av| 搡老熟女国产l中国老女人| 国产高清视频在线观看网站| 国产久久久一区二区三区| 99国产极品粉嫩在线观看| 嫁个100分男人电影在线观看| 特级一级黄色大片| 一个人免费在线观看的高清视频| 亚洲欧美日韩高清在线视频| 精品国产美女av久久久久小说| 桃色一区二区三区在线观看| 波多野结衣高清无吗| 成人手机av| 99国产精品一区二区三区| www.999成人在线观看| 精品久久久久久成人av| 国产野战对白在线观看| 十八禁人妻一区二区| 久久久久免费精品人妻一区二区| 亚洲精品国产精品久久久不卡| 日韩欧美在线二视频| 日本五十路高清| 久久午夜亚洲精品久久| 在线观看免费视频日本深夜| 色精品久久人妻99蜜桃| www日本在线高清视频| 欧美性猛交黑人性爽| 国产在线精品亚洲第一网站| 1024香蕉在线观看| 法律面前人人平等表现在哪些方面| 此物有八面人人有两片| 一进一出抽搐动态| 少妇粗大呻吟视频| 最近在线观看免费完整版| 小说图片视频综合网站| 69av精品久久久久久| 欧美乱码精品一区二区三区| 中文字幕熟女人妻在线| 人人妻人人看人人澡| 男女床上黄色一级片免费看| 欧美久久黑人一区二区| 丰满的人妻完整版| 午夜成年电影在线免费观看| 一区福利在线观看| 99热这里只有精品一区 | 亚洲国产欧美人成| 色尼玛亚洲综合影院| 一级a爱片免费观看的视频| 女警被强在线播放| 一进一出抽搐gif免费好疼| 欧美性猛交╳xxx乱大交人| 法律面前人人平等表现在哪些方面| 色哟哟哟哟哟哟| 黄频高清免费视频| 精品一区二区三区av网在线观看| 每晚都被弄得嗷嗷叫到高潮| 性色av乱码一区二区三区2| 日韩欧美国产在线观看| 亚洲乱码一区二区免费版| 亚洲国产欧美一区二区综合| 免费在线观看日本一区| 久久国产乱子伦精品免费另类| 久久性视频一级片| 99久久99久久久精品蜜桃| 啦啦啦免费观看视频1| 最近在线观看免费完整版| 国内精品一区二区在线观看| 欧美黑人欧美精品刺激| 夜夜夜夜夜久久久久| 欧美精品亚洲一区二区| 色综合站精品国产| 狂野欧美激情性xxxx| 亚洲国产高清在线一区二区三| 国产精品自产拍在线观看55亚洲| 国产久久久一区二区三区| 国产精品综合久久久久久久免费| 欧美日韩乱码在线| 国产激情久久老熟女| 日韩中文字幕欧美一区二区| 亚洲18禁久久av| 日本一区二区免费在线视频| 亚洲av成人一区二区三| avwww免费| 一本久久中文字幕| 亚洲成av人片在线播放无| 久久国产乱子伦精品免费另类| 日日干狠狠操夜夜爽| 色综合欧美亚洲国产小说| 精品国产亚洲在线| 极品教师在线免费播放| 亚洲国产精品久久男人天堂| 亚洲精品中文字幕在线视频| 久久久精品欧美日韩精品| 欧美日韩亚洲国产一区二区在线观看| 午夜两性在线视频| or卡值多少钱| 亚洲av电影不卡..在线观看| 国产午夜福利久久久久久| 国产一区二区三区在线臀色熟女| 亚洲精品国产一区二区精华液| 亚洲专区中文字幕在线| 久久天躁狠狠躁夜夜2o2o| 真人做人爱边吃奶动态| 久久亚洲精品不卡| 亚洲精品久久国产高清桃花| 岛国在线免费视频观看| 亚洲人与动物交配视频| 欧美日韩瑟瑟在线播放| 国产精品,欧美在线| 午夜精品久久久久久毛片777| 国产亚洲精品av在线| 国产不卡一卡二| 亚洲成av人片免费观看| 9191精品国产免费久久| 最好的美女福利视频网| 69av精品久久久久久| 亚洲天堂国产精品一区在线| 神马国产精品三级电影在线观看 | 一级毛片女人18水好多| 欧美日本亚洲视频在线播放| 久久天躁狠狠躁夜夜2o2o| 久久国产精品人妻蜜桃| 久久久久亚洲av毛片大全| 大型av网站在线播放| 国内精品久久久久精免费| 成人三级黄色视频| 男人的好看免费观看在线视频 | 长腿黑丝高跟| 免费看日本二区| 亚洲中文字幕一区二区三区有码在线看 | 欧美在线黄色| 国产成人av教育| 亚洲国产精品合色在线| 亚洲中文av在线| 级片在线观看| 99re在线观看精品视频| 国产精品 欧美亚洲| 少妇被粗大的猛进出69影院| 国产在线精品亚洲第一网站| 国产一区二区在线av高清观看| 看片在线看免费视频| 国产三级黄色录像| 亚洲一区二区三区不卡视频| cao死你这个sao货| 久久国产精品影院| 毛片女人毛片| 久久精品人妻少妇| 两人在一起打扑克的视频| 在线视频色国产色| 成人三级做爰电影| 在线永久观看黄色视频| 国产精品一区二区免费欧美| 国产精品日韩av在线免费观看| 精品高清国产在线一区| 性色av乱码一区二区三区2| 色老头精品视频在线观看| 丝袜美腿诱惑在线| 亚洲天堂国产精品一区在线| xxx96com| 亚洲专区字幕在线| 免费电影在线观看免费观看| 好男人在线观看高清免费视频| 一区二区三区国产精品乱码| 欧美av亚洲av综合av国产av| 国产精品久久视频播放| 黄色毛片三级朝国网站| 国产成人系列免费观看| 国产真实乱freesex| 在线观看舔阴道视频| 中文字幕av在线有码专区| 99热只有精品国产| 成年女人毛片免费观看观看9| 日本五十路高清| 人妻少妇偷人精品九色| 国语自产精品视频在线第100页| 乱人视频在线观看| 人人妻人人看人人澡| 亚洲激情五月婷婷啪啪| 日韩欧美三级三区| 晚上一个人看的免费电影| 人体艺术视频欧美日本| 69人妻影院| 国产伦在线观看视频一区| 国产一区二区亚洲精品在线观看| 十八禁国产超污无遮挡网站| 在线观看一区二区三区| a级毛色黄片| 麻豆久久精品国产亚洲av| 色播亚洲综合网| av在线老鸭窝| 欧美日韩精品成人综合77777| 51国产日韩欧美| 美女脱内裤让男人舔精品视频 | 国产精华一区二区三区| 欧美激情久久久久久爽电影| 亚洲最大成人av| 可以在线观看毛片的网站| av免费在线看不卡| 啦啦啦啦在线视频资源| 直男gayav资源| 1000部很黄的大片| 国内精品宾馆在线| 婷婷精品国产亚洲av| 最近中文字幕高清免费大全6| 国产午夜精品一二区理论片| 国产精品一区二区三区四区久久| 欧美3d第一页| 精品熟女少妇av免费看| 99国产精品一区二区蜜桃av| 亚洲国产日韩欧美精品在线观看| 亚洲美女视频黄频| 久久精品国产亚洲网站| 日本黄色视频三级网站网址| 天天躁日日操中文字幕| 国产精品99久久久久久久久| 国产一区二区亚洲精品在线观看| 国产精品麻豆人妻色哟哟久久 | 尤物成人国产欧美一区二区三区| 亚洲精品成人久久久久久| 韩国av在线不卡| 成人毛片a级毛片在线播放| 亚洲精品国产av成人精品| 亚洲经典国产精华液单| 岛国在线免费视频观看| 亚洲最大成人中文| 听说在线观看完整版免费高清| 精品国内亚洲2022精品成人| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 特级一级黄色大片| 99riav亚洲国产免费| 大型黄色视频在线免费观看| 欧美极品一区二区三区四区| 亚洲av免费在线观看| 亚洲欧美成人精品一区二区| 一区二区三区免费毛片| 亚洲激情五月婷婷啪啪| 欧美高清性xxxxhd video| 欧美精品国产亚洲| 国产成人一区二区在线| 日韩欧美精品免费久久| 神马国产精品三级电影在线观看| 亚洲精品乱码久久久v下载方式| 国产国拍精品亚洲av在线观看| 熟妇人妻久久中文字幕3abv| 国产精华一区二区三区| 村上凉子中文字幕在线| 国产又黄又爽又无遮挡在线| 听说在线观看完整版免费高清| 国产精品一区二区在线观看99 | 深夜a级毛片| 女的被弄到高潮叫床怎么办| 黄片wwwwww| 国产精品一区二区三区四区久久| 在线免费观看不下载黄p国产| 1000部很黄的大片| 亚洲欧美日韩高清专用| 九九在线视频观看精品| 能在线免费看毛片的网站| 色综合色国产| 嫩草影院入口| 夫妻性生交免费视频一级片| 精品久久久久久久末码| 小说图片视频综合网站| 午夜免费激情av| 高清午夜精品一区二区三区 | 麻豆成人av视频| 国产黄片视频在线免费观看| 久久精品国产亚洲av涩爱 | 搞女人的毛片| 最近中文字幕高清免费大全6| www.av在线官网国产| av卡一久久| 在线天堂最新版资源| 久久久色成人| 日韩精品青青久久久久久| 免费在线观看成人毛片| 欧美日韩在线观看h| 精品99又大又爽又粗少妇毛片| 国产伦精品一区二区三区视频9| 亚洲人成网站在线播| 久久精品国产清高在天天线| 成人美女网站在线观看视频| 在线免费观看的www视频| 亚洲三级黄色毛片| 欧美性猛交黑人性爽| 久久久久国产网址| 可以在线观看的亚洲视频| 精品国产三级普通话版| 人人妻人人看人人澡| 在线观看午夜福利视频| 欧美日本亚洲视频在线播放| 国产在线精品亚洲第一网站| 日本爱情动作片www.在线观看| 人妻夜夜爽99麻豆av| 人妻久久中文字幕网| 亚洲精品粉嫩美女一区| 看十八女毛片水多多多| 国产中年淑女户外野战色| 99热网站在线观看| 亚洲人成网站在线观看播放| 18禁裸乳无遮挡免费网站照片| 2021天堂中文幕一二区在线观| 亚洲国产精品国产精品| 欧美性感艳星| 三级国产精品欧美在线观看| 一卡2卡三卡四卡精品乱码亚洲| 国产一区亚洲一区在线观看| 午夜视频国产福利| 校园人妻丝袜中文字幕| 你懂的网址亚洲精品在线观看 | 麻豆成人av视频| 在现免费观看毛片| 亚洲精品成人久久久久久| 真实男女啪啪啪动态图| 天堂√8在线中文| 日韩欧美在线乱码| 国产成人影院久久av| 国产精品福利在线免费观看| 成年av动漫网址| 深夜精品福利| 国产精品人妻久久久久久| videossex国产| 亚洲成人精品中文字幕电影| 青春草视频在线免费观看| av卡一久久| 国产极品天堂在线| 观看免费一级毛片| av在线播放精品| 亚洲精品乱码久久久久久按摩| 国内久久婷婷六月综合欲色啪| 国产一区二区三区在线臀色熟女| 免费av不卡在线播放| 日本黄色片子视频| 哪里可以看免费的av片| av福利片在线观看| 性欧美人与动物交配| 床上黄色一级片| 在线观看66精品国产| 干丝袜人妻中文字幕| 亚洲精品粉嫩美女一区| 免费av不卡在线播放| 久久99热6这里只有精品| 长腿黑丝高跟| 变态另类成人亚洲欧美熟女| 又爽又黄无遮挡网站| 亚洲国产欧美人成| 全区人妻精品视频| 成人亚洲精品av一区二区| or卡值多少钱| 国产精品一区二区性色av| 亚洲精品乱码久久久v下载方式| 国国产精品蜜臀av免费| 亚洲人成网站在线播放欧美日韩| 干丝袜人妻中文字幕| 亚洲欧洲日产国产| 高清午夜精品一区二区三区 | 国产男人的电影天堂91| 国产精品1区2区在线观看.| 婷婷色av中文字幕| 成人特级av手机在线观看|