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

    ENDSI增強型雪指數(shù)提取積雪研究

    2018-03-06 08:41:24龐海洋孔祥生汪麗麗錢永剛
    自然資源遙感 2018年1期
    關(guān)鍵詞:綠光積雪反射率

    龐海洋, 孔祥生, 汪麗麗, 錢永剛

    (1.魯東大學(xué)資源與環(huán)境工程學(xué)院,煙臺 264025; 2.中國科學(xué)院光電研究院定量遙感信息技術(shù)重點實驗室,北京 100094 )

    0 引言

    降雪是人類生活環(huán)境中常見的自然現(xiàn)象。冬季北半球積雪和冰蓋面積超過40%; 在南、北半球的高緯度和高海拔地區(qū),大范圍的地表常年被積雪覆蓋[1]。積雪是影響全球熱量平衡的關(guān)鍵氣候因子,其變化會對區(qū)域氣候產(chǎn)生重要影響[2],是氣候系統(tǒng)的一個重要組成部分。積雪也是一種重要的淡水資源[3],在世界上很多地區(qū)都在依靠冰、雪融水來維持基本生活和農(nóng)業(yè)生產(chǎn)[4]。在我國干旱半干旱地區(qū),積雪融水是重要的淡水補充形式,對當?shù)氐淖匀画h(huán)境和人類生活意義深遠。在全球變暖背景下,積雪、凍土對氣候變化具有高度敏感性和重要反饋作用,已成為全球變化研究的核心內(nèi)容和熱點之一[5]。積雪覆蓋的時空特征研究在氣象學(xué)和水文學(xué)中也發(fā)揮著重要作用,對積雪的研究有重要現(xiàn)實意義[6]。

    傳統(tǒng)利用氣象站或者野外調(diào)查來監(jiān)測積雪的方法難度大、時效性差,存在許多弊端,遙感技術(shù)以其大范圍、多時相等特點逐漸成為積雪觀測的主要手段[7-8]。利用遙感技術(shù)提取積雪信息的方法可歸納成4種: 閾值法像元統(tǒng)計、監(jiān)督和非監(jiān)督分類、歸一化差值雪指數(shù)(normalized difference snow index,NDSI)和混合像元分解法[9]。其中,NDSI法是基于雪對可見光與短波紅外波段的反射特性和反射差的相對大小的一種測量方法[10],是目前最常用、精度較高的積雪提取方法之一,最初是Hall等在1995年提出的 “SNOMAP”方法[11],隨后國內(nèi)外學(xué)者對NDSI進行了不同程度的研究。趙軍等[12]以瑪納斯上游地區(qū)為例對雪蓋指數(shù)提取范圍的不確定性進行了研究; 陳文倩等[13]基于NDSI-NDVI特征對積雪面積進行了反演; Satir[14]利用NDSI等雪指數(shù)對意大利Cilo山的冰雪進行了連續(xù)性監(jiān)測分析。然而,NDSI只利用可見光中的綠光波段(0.53~0.59 μm)和短波紅外波段(1.57~1.65 μm),在薄雪區(qū)敏感度較低,不能精確地區(qū)分雪與非雪[15]。2013年2月12日發(fā)射成功的Landsat8 衛(wèi)星搭載了陸地成像儀(operational land imager,OLI)增加了B1波段(0.433~0.453 μm),雪對該波段和藍光波段(0.45~0.5 μm)也相當敏感,加之其較高的空間分辨率及時間分辨率,為小尺度積雪監(jiān)測又增添了新數(shù)據(jù)[16]。故以Landsat8 OLI為數(shù)據(jù)源,基于積雪光譜特征,提出增強型雪指數(shù)(enhanced normalized difference snow index,ENDSI)方法,試圖為精確提取積雪與雪厚遙感定量反演提供技術(shù)支持。

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

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

    Landsat8攜帶OLI和熱紅外傳感器(thermal infrared sensor,TIRS)2個傳感器,成像幅寬為185 km。其中OLI 陸地成像儀數(shù)據(jù)有9個波段,多光譜波段空間分辨率為30 m,全色波段空間分辨率達到15 m。與傳統(tǒng)的TM和ETM+不同,Landsat8在原藍光波段之外新增了B1(0.433~0.453 μm)[17],主要用于對海岸線變化的監(jiān)測,而積雪在該波段敏感度也較高。故本文以無云的Landsat8 OLI影像為數(shù)據(jù)源,構(gòu)建ENDSI,提高積雪提取精度。

    1.2 研究區(qū)概況

    NDSI與雪厚呈正相關(guān)關(guān)系,在厚雪區(qū)NDSI數(shù)值較大,在薄雪區(qū)其敏感度降低,受混合像元中其他地物類型的影響,NDSI數(shù)值較小,利用閾值法區(qū)分雪與非雪有一定難度。為探究在不同地理環(huán)境下ENDSI的適用性及其對積雪厚度的敏感程度,提高遙感監(jiān)測積雪的精度,根據(jù)經(jīng)緯度分布和海拔高低差異狀況選擇我國北方5個區(qū)域作為研究區(qū)(表1、圖1)。

    表1 研究區(qū)遙感數(shù)據(jù)源Tab.1 Remote sensing data sources of research areas

    (a) 研究區(qū)1 (b) 研究區(qū)2 (c) 研究區(qū)3

    (d) 研究區(qū)4 (e) 研究區(qū)5

    (Landsat8 OLI B6(R),B3(G),B2(B)假彩色合成)

    圖1研究區(qū)位置

    Fig.1Locationsofresearchareas

    圖1中偏藍色區(qū)域為積雪,偏紅色區(qū)域為土壤,其中圖1(a)為山東半島地區(qū),位于我國東部,平均海拔在500 m以下,屬于溫帶季風(fēng)氣候,降雪量豐富; 圖1(b)和(c)分別為內(nèi)蒙古自治區(qū)赤峰地區(qū)和甘肅張掖地區(qū),位于我國第二階梯,平均海拔在1 000~2 000 m之間,屬溫帶大陸性氣候,緯度較高,冬季積雪較厚; 圖1(d)為西藏自治區(qū)昌都地區(qū),屬于第一階梯,平均海拔在4 000 m以上; 而圖1(e)作為驗證區(qū),地處山西省、內(nèi)蒙古自治區(qū)和河北省的交界地區(qū),從東西方向來看,位于我國中部第二、三階梯的交界地帶,平均海拔為1 343 m,相對西部地區(qū)而言海拔較低,較東部地區(qū)海拔較高,屬于半干旱地區(qū),大陸性氣候顯著,冬季受西伯利亞高壓影響,具有典型的積雪特征; 同時,該地區(qū)是我國煤炭的重要產(chǎn)地,污染嚴重,生態(tài)環(huán)境脆弱,環(huán)境變化一直備受關(guān)注。

    2 研究方法

    積雪在可見光波段(0.38~0.76 μm)有較高的反射率,而且與雪厚有著較好的相關(guān)性,當積雪厚度小于20 cm 時,積雪的反射率隨著雪厚的增加而增加,但在短紅外波段,不如可見光波段明顯[18]。為了便于描述積雪的特性,根據(jù)積雪在綠光波段的反射率大小,將積雪分為3種類型[19]: 反射率在[0.45,1]范圍為厚雪(混合像元中積雪占50%以上),有明顯的積雪光譜特征; 反射率在[0.22,0.45)范圍內(nèi)為薄雪(混合像元中積雪占25%~50%),主要受積雪、土壤綜合影響; 反射率在[0.1,0.22)為超薄雪區(qū),積雪的反射特征較弱,以土壤的光譜特征為主(混合像元中積雪所占比例小于25%),故將超薄雪與土壤合稱非雪區(qū)。利用4個具有不同地理特征的研究區(qū)1—4分析區(qū)分雪與非雪的合理閾值,并以研究區(qū)5為驗證區(qū),利用目視解譯和混淆矩陣方法進行精度驗證。

    2.1 積雪的光譜特征

    與其他地物相比,積雪的光譜特征明顯(圖2)。

    圖2 研究區(qū)主要地物光譜曲線Fig.2 Spectral curves of the main objects in the study area

    圖2中,新雪具有獨特的光譜特征: 在可見光范圍內(nèi)有0.8以上的高反射率,0.6 μm左右形成一個高反射峰; 在近紅外波段(0.845~0.885 μm)反射率開始減?。?在短波紅外波段1.6 μm和2.0 μm左右分別形成2個吸收谷,反射率在0.15以下,而在1.8 μm和2.3 μm處形成2個反射峰,反射率在0.4以下; 2.5 μm以后反射率逐漸接近于0。老雪和融化的積雪變化趨勢與新雪相似,但反射峰出現(xiàn)在0.3~1.3 μm,1.7~1.8 μm和2.2~2.3 μm處,強度比新雪有不同程度的下降,其中壓實凍結(jié)的冰雪反射率最低[20]。土壤、植被和水體的反射率在可見光波段反射率較低,大都在0.15以下,而在短波紅外波段植被和水體的反射率較低,土壤的反射率在0.3以上,大于新雪的反射率。

    2.2 增強型雪指數(shù)(ENDSI)

    NDSI利用了積雪在綠光波段反射率較高,在短波紅外波段反射率較低的獨特變化特征,當NDSI≥0.4,ρgreen≥0.1且ρnir≥0.11時,判定為積雪[11],計算公式為

    (1)

    式中ρgreen,ρswir和ρnir分別為積雪在綠光波段、短波紅外波段和近紅外波段的反射率。

    NDSI的值取決于積雪在綠光波段和短波紅外波段反射率差值的大小,當ρgreen-ρswir接近ρgreen時,精度較高; 當ρgreen-ρswir非常小或者NDSI<0.4時,則可能存在較大誤差。此外由于遙感衛(wèi)星空間分辨率和地表環(huán)境復(fù)雜度影響,在設(shè)定閾值提取積雪時,混合像元區(qū)易產(chǎn)生誤差。

    為提高混合像元積雪的識別精度,可通過某種變換手段或積雪指數(shù)方法增大雪區(qū)與非雪區(qū)的差異,進而提高積雪提取精度。AvaField-3光譜儀實測和OLI影像采集得到的積雪光譜都表明,積雪光譜滿足

    (2)

    式中ρblue violet和ρblue分別是OLI數(shù)據(jù)在藍紫光和藍光波段的反射率。對式(2)中的3個判別標準進行累加處理,得

    ρblue violet+ρblue+ρgreen-3ρswir>0。

    (3)

    NDSI能夠?qū)崿F(xiàn)積雪的提取,但在所提取的積雪區(qū)內(nèi)部對雪厚的劃分不夠敏感,只能將像元辨別為雪或非雪,滿足不了高精度的流域制圖及雪蓋提取要求[15]。而式(3)增大了雪與非雪的差異,一定程度上也增大了雪間差異,在式(3)基礎(chǔ)上,提出ENDSI為

    (4)

    式中a為調(diào)節(jié)系數(shù)。當ENDSI≥0.3,ρgreen≥0.11且ρnir≥0.11時判定為積雪。

    為使積雪指數(shù)對混合像元中薄雪敏感,適當調(diào)節(jié)式(4)中a,當03.7時,部分雪區(qū)ENDSI<0。在大量實驗和仿真模擬的基礎(chǔ)上得出,當a=3.7時最佳,得到優(yōu)化后的ENDSI,即

    (5)

    3 結(jié)果分析與驗證

    3.1 閾值分析

    對5景影像進行輻射定標、FLAASH大氣校正和裁剪等預(yù)處理后,分別計算NDSI和ENDSI,在[0,1]之間每隔0.1選取一個閾值,提取積雪的像元數(shù)量變化,并對4個研究區(qū)進行綜合分析(圖3)。

    (a) 研究區(qū)1積雪像元總量 (b) 研究區(qū)1積雪像元變化數(shù)量 (c) 研究區(qū)2積雪像元總量 (d) 研究區(qū)2積雪像元變化數(shù)量

    (e) 研究區(qū)3積雪像元總量 (f) 研究區(qū)3積雪像元變化數(shù)量 (g) 研究區(qū)4積雪像元總量 (h) 研究區(qū)4積雪像元變化數(shù)量

    圖3ENDSI和NDSI提取積雪像元總量及變化數(shù)量

    Fig.3NumbersandvariationofextractedsnowpixelsbyENDSIandNDSImethods

    從像元總量上來看,ENDSI和NDSI提取的積雪像元總量均隨著閾值的增大而減小,在閾值[0,0.45)區(qū)間內(nèi),ENDSI提取積雪像元數(shù)小于NDSI提取結(jié)果,總量差距隨閾值增大逐漸減小; 當閾值在0.45左右時,兩者提取積雪像元數(shù)相當; 在(0.45,1)區(qū)間內(nèi),ENDSI提取雪像元量逐漸大于NDSI提取結(jié)果; 當閾值為1時,兩者提取積雪像元數(shù)相等。從像元變化數(shù)量來看,隨著閾值的增大,ENDSI和NDSI提取積雪像元變化量減小幅度不同,閾值在[0,0.7]區(qū)間ENDSI提取積雪像元減小量小于NDSI提取結(jié)果; 在(0.7,1]區(qū)間兩者變化趨于相同,在[0.6,1]區(qū)間積雪反射率較高,但積雪純度及厚度不同反射率仍有差異,NDSI和ENDSI在閾值為0.8左右時,常出現(xiàn)積雪像元數(shù)量急劇減少的現(xiàn)象,說明積雪主要集中在[0.6,0.8],在[0.8,1]區(qū)間積雪像元量較少。NDSI在[0,0.6]區(qū)間像元變化量相對平穩(wěn),變化不大,不利于選取閾值區(qū)分雪與非雪。冬季北方降雪地區(qū)地物類型主要以雪和土壤為主,積雪像元數(shù)量隨閾值的變化,正是提取類型由土壤向雪過渡的體現(xiàn): 在[0,0.3)區(qū)間主要為非雪或超薄雪,光譜特征相似,像元變化量較平穩(wěn); 在[0.3,0.6)之間地物類型由土壤向積雪區(qū)過渡,地物光譜差異開始增大,逐漸突破非雪的界限,積雪像元減小量逐漸增加; 在[0.6,1]區(qū)間地物主要為積雪,光譜特征相似,積雪像元數(shù)量變化趨于平穩(wěn),但在[0.8,1]區(qū)間時常出現(xiàn)積雪像元減小量劇增或驟減現(xiàn)象,這是由積雪像元過渡集中而閾值間隔過大造成的。由非雪向雪區(qū)過渡的轉(zhuǎn)折點0.3成為ENDSI判定雪與非雪的分割點,地理特征不同的4個研究區(qū)變化特征相似,說明0.3作為閾值分割點有一定的穩(wěn)定性、普適性。

    3.2 ENDSI變化曲線

    在研究區(qū)5 的Landsat8 OLI遙感影像中選取裸土和積雪各5條光譜曲線(圖4),在藍紫光—綠光波段,裸土反射率在0.15以下,積雪反射率在0.2以上; 短波紅外波段,裸土反射率在0.15~0.3之間,積雪反射率在0.2以下。

    圖4 研究區(qū)5中積雪和裸土光譜曲線Fig.4 Spectral curves of the snow and soil in study area 5

    郝曉華等[19]測量已知積雪覆蓋比例下的積雪光譜特征,也發(fā)現(xiàn)積雪在0.35~0.8 μm范圍內(nèi)反射率隨積雪純度的增加而減小,1.4~1.8 μm范圍內(nèi)反射率與積雪純度成反比,雖然積雪下墊面性質(zhì)不同反射率會有所差異,但以上變化規(guī)律不變。通過調(diào)節(jié)式(4)中系數(shù)a,調(diào)控可見光與短波紅外波段差異的大小,可以調(diào)節(jié)積雪指數(shù)曲度,實現(xiàn)對薄雪敏感。

    根據(jù)式 (6)分別計算裸土和積雪在短波紅外波段反射率ρsoilswir和ρsnowswir[21],即

    (6)

    式中:Γ(λ)是傳感器的固有參數(shù),各個波段的響應(yīng)函數(shù)曲線;ρ(λ)是光譜儀測試的目標反射率數(shù)據(jù)。

    在該地區(qū)雪區(qū)選取5個樣本點(113°8′49.22″E,40°9′30″N; 113°3′7.69″E,39°6′23.42″N; 112°4′9.32″E,40°4′58.33″N; 112°4′1.66″E,40°7′32.22″N; 112°2′28.22″E,39°1′13.64″N)發(fā)現(xiàn)雪區(qū)OLI前3個波段反射率之和相當于3倍綠光波段反射率減去0.1的值; 在裸土區(qū)選取5個樣本點(112°1′54.18″E,39°9′54.87″N; 113°0′56.11″E,40°7′48.27″N; 113°3′29.98″E,40°9′29.04″N; 114°3′35.25″E,40°0′40.52″N; 113°0′23.11″E,39°9′14.67″N)發(fā)現(xiàn)裸土區(qū)OLI前3個波段反射率之和相當于3倍綠光波段反射率減去0.06的值。以綠光反射率為橫坐標,裸土和積雪的ENDSI和NDSI值為縱坐標,根據(jù)表2中公式(7)—(14)模擬出NDSI與ENDSI變化曲線(圖5)。

    表2 ENDSI線性模擬公式Tab.2 ENDSI linear simulation formula

    (a) 積雪和裸土ENDSI和NDSI下的變化曲線 (b) 不同抑制系數(shù)下的ENDSI曲線

    圖5NDSI和NDSI在綠光波段反射率的變化趨勢(由特定值畫出)

    Fig.5VariabletrendgraphofNDSIandENDSIingreenband(drawfromaspecificvalue)

    圖5(a)中A為超薄雪區(qū),B為薄雪區(qū),C為厚雪區(qū),積雪在綠光波段反射率越高代表雪層越厚,積雪光譜特征越明顯,反射率越低則雪層越稀薄,積雪光譜特征越不明顯,ENDSI和NDSI隨綠光波段反射率的增大而增大。與NDSI相比,ENDSI在厚雪區(qū)增大,在薄雪和超薄雪區(qū)減小,超薄雪區(qū)減小幅度大于薄雪區(qū),厚雪區(qū)隨反射率的增加增大越明顯,但總體增加幅度有限。薄雪區(qū)在NDSI中主要位于[0.2,0.5],在ENDSI中位于[0,0.5],相同的反射率范圍下閾值區(qū)間擴大,顯示ENDSI在薄雪區(qū)更加敏感; 在超薄雪區(qū)ENDSI全部抑制在0以下; 裸土在ENDSI中較NDSI大幅度減小,被抑制在-0.2以下。ENDSI有效增大雪與非雪的差異,對薄雪具敏感性,為進一步精細化提取積雪和雪厚,以及遙感定量反演打下了基礎(chǔ)。從圖5(b)可見,ENDSI對雪厚敏感度總體上強于NDSI,ENDSI隨抑制系數(shù)a的增大,在薄雪區(qū)敏感度逐漸增強。當a=2和a=2.5時,ENDSI在超薄雪敏感度高于薄雪區(qū)和厚雪區(qū); 當a=3.7和a=4時,ENDSI在超薄雪區(qū)為負值,對雪厚敏感度在薄雪區(qū)強于厚雪區(qū),但當a=4時,薄雪區(qū)部分ENDSI值也為負值。超薄雪區(qū)雪含量較低,雪指數(shù)在厚雪區(qū)敏感度較低,對薄雪區(qū)精細化提取時,a=3.7較為合理。

    3.3 積雪提取結(jié)果

    驗證區(qū)不同方法的積雪提取結(jié)果如圖6所示,積雪表現(xiàn)為偏藍色特征。

    (a) OLI B6(R),B3(G),B2(B)假彩色合成影像 (b) ENDSI提取 (c) NDSI提取

    (d) ENDSI單獨提取部分(e) ENDSI和NDSI提取的相同區(qū)

    圖6驗證區(qū)ENDSI和NDSI積雪提取結(jié)果

    Fig.6SnowcovermapsextractedbyENDSIandNDSIinstudyarea5

    總體上看,ENDSI和NDSI都基本能夠滿足積雪提取要求,積雪面積較小不足全影像的1/10,主要集中在驗證區(qū)的5個區(qū)域,分布較零散。其中西部有一湖泊,結(jié)冰后被積雪覆蓋,反映積雪的光譜特征; 此外,在河流、山體陰影和林地中也有積雪零星分布。在積雪較厚的區(qū)域顏色明顯偏藍,在積雪較薄的地區(qū)顏色偏暗,積雪特征并不明顯,提取積雪時很容易誤分為裸土。從細節(jié)上看,NDSI提取的積雪(圖6(c))與ENDSI(圖6(b))提取的積雪相比,顯得較松散稀疏,許多零散分布的積雪不能提取,中部雪區(qū)差別尤為明顯。ENDSI和NDSI提取的積雪重合度較高(圖6(e)),主要集中在厚雪區(qū),而在薄雪區(qū),差異較大。

    經(jīng)過統(tǒng)計,研究區(qū)總面積約為37 407.8 km2,NDSI閾值為0.4時提取的積雪面積為1 843.377 km2,占總面積的4.927 8%; 利用ENDSI提取的積雪面積為2 342.04 km2,占總面積的6.260 8%(表3)。圖6(d)中ENDSI判斷為積雪而NDSI判斷為非積雪的面積有498.795 3 km2,占總面積的1.333 4%。表3中所有數(shù)據(jù)都經(jīng)過掩模處理,ENDSI-NDSI表示ENDSI提取為積雪而NDSI提取為非積雪部分,NDSI-ENDSI表示NDSI提取為積雪而ENDSI提取為非積雪部分。

    表3 ENDSI和NDSI提取積雪像元統(tǒng)計Tab.3 Statistics of extracted snow pixels by NDSI and ENDSI

    3.4 精度分析

    利用混淆矩陣和目視解譯2種方法對ENDSI提取精度進行比較分析。對Landsat8 OLI影像監(jiān)督分類,分為雪和非雪2類,將分類結(jié)果作為真實積雪分布,分別與ENDSI和NDSI提取的雪區(qū)與非雪區(qū)做混淆矩陣進行對比分析,同時分析其總體精度(表4—5)。其中NDSI一般設(shè)定0.4為雪與非雪的有效閾值分割點[11],吳曉晨等[22]證明NDSI>0.4可以用來提取冰雪覆蓋。當NDSI降低閾值時可以增加雪提取面積,由于NDSI對雪厚敏感度較低,降低閾值容易造成閾值減小過量。

    表4 驗證區(qū)監(jiān)督分類與NDSI分類結(jié)果混淆矩陣Tab.4 Confusion matrix between validation area supervised classification and NDSI classification result

    表5 驗證區(qū)監(jiān)督分類與ENDSI分類結(jié)果混淆矩陣Tab.5 Confusion matrix between validation area supervised classification and ENDSI classification result

    利用NDSI和ENDSI提取積雪的分類總體精度分別為97.636 7%和97.937 1%,Kappa相關(guān)系數(shù)分別為0.778 5和0.825 7,與NDSI相比,ENDSI總體分類精度提高了0.3%,Kappa系數(shù)提高了0.05。

    為了進一步檢驗ENDSI對厚雪和薄雪的敏感性,在研究區(qū)5遙感影像上分別選取同時存在土壤、薄雪和厚雪的3個子區(qū)域a,b和c,以目視解譯的方法分析ENDSI和NDSI提取積雪的分布情況。圖7由B6(R),B3(G),B2(B)假彩色合成,積雪在可見光具有高反射率,在短波紅外強吸收的光譜特征[23],在遙感影像上表現(xiàn)為偏藍色,由于厚雪在B2和B3波段反射率遠大于薄雪,在影像中亮度較高,呈現(xiàn)亮藍色,薄雪區(qū)亮度較低,呈現(xiàn)暗藍色,裸土呈現(xiàn)偏紅色的特征。在a,b和c 3個子區(qū)中,裸土和積雪作為標志性地物類型,分布特征實現(xiàn)了由裸土到薄雪再到厚雪的過渡,其中厚雪一般為團聚狀,薄雪以厚雪為中心向四周擴展,擴展距離不等。

    (a) 子驗證區(qū)a (b) 子驗證區(qū)b (c) 子驗證區(qū)c

    圖7目視解譯驗證實例

    Fig.7Accuracyverificationexamplediagrams

    圖7中黑色方框區(qū)域是ENDSI判定為雪而NDSI判定為非雪部分,綠色方框區(qū)域為ENDSI和 NDSI同時提取的積雪區(qū)。圖中綠色方框部分主要集中在亮藍色的厚雪區(qū),可見ENDSI和NDSI都能夠?qū)崿F(xiàn)對厚雪的提取,兩者差異較??; 而ENDSI除與NDSI共同提取的厚雪(綠色方框部分)部分外,在薄雪區(qū)(黑色方框部分)同樣判定為雪,在圖7(a),(b)和(c)中薄雪區(qū)分別占總面積的19.73%,11.84%和15.48%,經(jīng)目視解譯這些區(qū)域與土壤有明顯差異,應(yīng)判定為雪,可以得出ENDSI對薄雪更加敏感,能夠有效提取薄雪,且提取精度有所提高。在圖7(a)和(b)中提取的厚雪和薄雪成團聚狀分布,而圖7(c)中厚雪和薄雪離散分布,總體特征來看薄雪仍以厚雪為中心向四周擴展,擴展距離不等,提取結(jié)果與目視解譯一致。

    4 結(jié)論

    本文在傳統(tǒng)NDSI利用可見光范圍內(nèi)單一綠光波段和短波紅外波段提取積雪基礎(chǔ)上,基于Landsat8 OLI數(shù)據(jù)源,對藍紫光、藍光、綠光和短波紅外波段綜合考慮,探究ENDSI提取積雪的可能性,結(jié)論如下:

    1)依據(jù)積雪在綠光波段的反射率將雪分為厚雪區(qū)、薄雪區(qū)和超薄雪區(qū),ENDSI對雪厚更加敏感,ENDSI值變化幅度強于NDSI。與NDSI相比,ENDSI值在厚雪區(qū)增大,超薄雪區(qū)和裸土區(qū)減小,有效增大雪與非雪的差異,有利于對積雪與雪厚的精確提取和遙感定量反演。

    2)通過對4個不同經(jīng)緯度和海拔高度特征的地區(qū)提取雪像元數(shù)量變化分析,發(fā)現(xiàn)0.3為ENDSI區(qū)分雪與非雪的有效分割點,穩(wěn)定性較高。

    3)以驗證區(qū)監(jiān)督分類得到積雪分布為真值分別與ENDSI和NDSI提取積雪結(jié)果做混淆矩陣,發(fā)現(xiàn)ENDSI比NDSI總體精度增加了0.3%,2組Kappa系數(shù)均在0.77以上,分類結(jié)果有較高的可信度。

    利用ENDSI提取積雪精度有所提高,特別是在薄雪區(qū)尤為敏感。今后可以進一步探究積雪厚度與ENDSI的相關(guān)性,研究利用遙感技術(shù)定量反演雪厚的方法; 海冰與雪光譜特征相似,可以繼續(xù)探究ENDSI在海冰提取方面的應(yīng)用。

    [1] Jones H G,Pomeroy J W,Walker D A,et al.Snow Ecology:An Interdisciplinary Examination of Snow-covered Ecosystems[M].Cambridge:Cambridge University Press,2001.

    [2] 白淑英,吳 奇,史建橋,等.青藏高原積雪深度時空分布與地形的關(guān)系[J].國土資源遙感,2015,27(4):171-178.doi:10.6046/gtzyyg.2015.04.26.

    Bai S Y,Wu Q,Shi J Q,et al.Relationship between the spatial and temporal distribution of snow depth and the terrain over the Tibetan Plateau[J].Remote Sensing for Land and Resources,2015,27(4):171-178.doi:10.6046/gtzyyg.2015.04.26.

    [3] 郝曉華,王 建,李弘毅.MODIS雪蓋制圖中NDSI閾值的檢驗——以祁連山中部山區(qū)為例[J].冰川凍土,2008,30(1):132-138.

    Hao X H,Wang J,Li H Y.Evaluation of the NDSI threshold value in mapping snow cover of MODIS:A case study of snow in the Middle Qilian Mountains[J].Journal of Glaciology and Geocryology,2008,30(1):132-138.

    [4] Xiao X M,Shen Z X,Qin X G.Assessing the potential of VEGETATION sensor data for mapping snow and ice cover:A normalized difference snow and ice index[J].International Journal of Remote Sensing,2001,22(13):2479-2487.

    [5] 王國亞,毛煒嶧,賀 斌,等.新疆阿勒泰地區(qū)積雪變化特征及其對凍土的影響[J].冰川凍土,2012,34(6):1293-1300.

    Wang G Y,Mao W Y,He B,et al.Changes in snow covers during 1961—2011 and its effects on frozen ground in Altay Region, Xinjiang[J].Journal of Glaciology and Geocryology,2012,34(6):1293-1300.

    [6] Kour R,Patel N,Krishna A P.Assessment of relationship between snow cover characteristics(SGI and SCI) and snow cover indices(NDSI and S3)[J].Earth Science Informatics,2015,8(2):317-326.

    [7] 于泓峰,張顯峰.光學(xué)與微波遙感的新疆積雪覆蓋變化分析[J].地球信息科學(xué)學(xué)報,2015,17(2):244-252.

    Yu H F,Zhang X F.Retrieval and analysis of snow-covered days in Xinjiang based on optical and microwave remote sensing data[J].Journal of Geo-Information Science,2015,17(2):244-252.

    [8] 紀 鵬,郭華東,張 露.近20年西昆侖地區(qū)冰川動態(tài)變化遙感研究[J].國土資源遙感,2013,25(1):93-98.doi:10.6046/gtzyyg.2013.01.17.

    Ji P,Guo H D,Zhang L.Remote sensing study of glacier dynamic change in West Kunlun Mountains in the past 20 years[J].Remote Sensing for Land and Resources,2013,25(1):93-98.doi:10.6046/gtzyyg.2013.01.17.

    [9] 彥立利,王 建.基于遙感的冰川信息提取方法研究進展[J].冰川凍土,2013,35(1):110-118.

    Yan L L,Wang J.Study of extracting glacier information from remote sensing[J].Journal of Glaciology and Geocryology,2013,35(1):110-118.

    [10] 惠鳳鳴,田慶久,李英成,等.基于MODIS數(shù)據(jù)的雪情分析研究[J].遙感信息,2004,19(4):35-38.

    Hui F M,Tian Q J,Li Y C,et al.Research on snow condition analysis based on MODIS data[J].Remote Sensing Information,2004,19(4):35-38.

    [11] Hall D K,Riggs G A,Salomonson V V.Development of methods for mapping global snow cover using moderate resolution imaging spectroradiometer data[J].Remote Sensing of Environment,1995,54(2):127-140.

    [12] 趙 軍,付杰文,付 鵬.雪蓋指數(shù)法提取積雪范圍信息的不確定性研究——以瑪納斯上游地區(qū)為例[J].遙感技術(shù)與應(yīng)用,2014,29(2):293-299.

    Zhao J,Fu J W,Fu P.Accuracy and uncertainty of snow information based on NDSI:A case study of upstream region of Manas River[J].Remote Sensing Technology and Application,2014,29(2):293-299.

    [13] 陳文倩,丁建麗,孫永猛,等.基于NDSI-NDVI特征空間的積雪面積反演研究[J].冰川凍土,2015,37(4):1059-1066.

    Chen W Q,Ding J L,Sun Y M,et al.Retrieval of snow cover area based on NDSI-NDVI feature space[J].Journal of Glaciology and Geocryology,2015,37(4):1059-1066.

    [14] Satir O.Comparing the satellite image transformation techniques for detecting and monitoring the continuous snow cover and glacier in Cilo Mountain chain Turkey[J].Ecological Indicators,2016,69:261-268.

    [15] 裴 歡,房世峰,覃志豪,等.基于遙感的新疆北疆積雪蓋度及雪深監(jiān)測[J].自然災(zāi)害學(xué)報,2008,17(5):52-57.

    Pei H,Fang S F,Qin Z H,et al.Remote sensing-based monitoring of coverage and depth of snow in northern Xinjiang[J].Journal of Natural Disasters,2008,17(5):52-57.

    [16] 趙 軍,陳愷悅,師銀芳.提高森林覆蓋區(qū)積雪提取精度的方法研究——以瑪納斯河流域為例[J].遙感技術(shù)與應(yīng)用,2015,30(6):1051-1058.

    Zhao J,Chen K Y,Shi Y F.Methods research to improve the extraction accuracy of snow under forest cover:A case study of Manas River basin[J].Remote Sensing Technology and Application,2015,30(6):1051-1058.

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

    Xu H Q,Tang F.Analysis of new characteristics of the first Landsat8 image and their eco-environmental significance[J].Acta Ecologica Sinica,2013,33(11):3249-3257.

    [18] 魏 玥.北疆區(qū)域積雪深度變化的遙感監(jiān)測研究[D].烏魯木齊:新疆師范大學(xué),2010.

    Wei Y.Remote Sensing Monitoring of Snow Depth Change in North Part of Xinjiang[D].Urumqi:Xinjiang Normal University,2010.

    [19] 郝曉華,王 杰,王 建,等.積雪混合像元光譜特征觀測及解混方法比較[J].光譜學(xué)與光譜分析,2012,32(10):2753-2758.

    Hao X H,Wang J,Wang J,et al.Observations of snow mixed pixel spectral characteristics using a ground-based spectral radiometer and comparing with unmixing algorithms[J].Spectroscopy and Spectral Analysis,2012,32(10):2753-2758.

    [20] 閃 旭,劉志輝,張 波.新疆軍塘湖流域融雪期不同積雪及雪被地物光譜反射率特征分析[J].安徽農(nóng)業(yè)科學(xué),2014,42(3):853-855,887.

    Shan X,Liu Z H,Zhang B.Study of spectrum reflectance characteristics of snow and snow-covered land surface objects in the melting-snow period[J].Journal of Anhui Agricultural Sciences,2014,42(3):853-855,887.

    [21] 程 熙,沈占鋒,駱劍承,等.利用地物波譜學(xué)習(xí)的遙感影像波段模擬方法[J].紅外與毫米波學(xué)報,2010,29(1):45-48,62.

    Cheng X,Shen Z F,Luo J C,et al.Method on simulating remote sensing image band by using ground-object spectral features study[J].Journal of Infrared and Millimeter Waves,2010,29(1):45-48,62.

    [22] 吳曉晨,孟令奎,張東映,等.冰雪遙感監(jiān)測方法綜述[J].水利信息化,2013(1):35-39.

    Wu X C,Meng L K,Zhang D Y,et al.Overview on methods of snow and ice remote sensing monitoring[J].Water Resources Informatization,2013(1):35-39.

    [23] 劉玉潔,鄭照軍,王麗波.我國西部地區(qū)冬季雪蓋遙感和變化分析[J].氣候與環(huán)境研究,2003,8(1):114-123.

    Liu Y J,Zheng Z J,Wang L B.Remote sensing on snow cover and variation analyzing in west of China[J].Climatic and Environmental Research,2003,8(1):114-123.

    猜你喜歡
    綠光積雪反射率
    影響Mini LED板油墨層反射率的因素
    近岸水體異源遙感反射率產(chǎn)品的融合方法研究
    海洋通報(2022年4期)2022-10-10 07:40:26
    具有顏色恒常性的光譜反射率重建
    我們
    大糧積雪 誰解老將廉頗心
    炎黃地理(2017年10期)2018-01-31 02:15:05
    化學(xué)腐蝕硅表面結(jié)構(gòu)反射率影響因素的研究*
    電子器件(2017年2期)2017-04-25 08:58:37
    積雪
    2000~2014年西藏高原積雪覆蓋時空變化
    機器 人
    渴望
    亚洲av不卡在线观看| av.在线天堂| 三级经典国产精品| 亚洲国产精品sss在线观看| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲av二区三区四区| 日韩一本色道免费dvd| 亚洲av.av天堂| 深爱激情五月婷婷| 亚洲av免费在线观看| 国产精品日韩av在线免费观看| 女人十人毛片免费观看3o分钟| 国产亚洲av嫩草精品影院| 女的被弄到高潮叫床怎么办| 国产毛片a区久久久久| 波多野结衣高清无吗| 久久婷婷人人爽人人干人人爱| 国产精品女同一区二区软件| 亚洲一区高清亚洲精品| 伊人久久精品亚洲午夜| 少妇人妻一区二区三区视频| 久久综合国产亚洲精品| 成人二区视频| 日韩国内少妇激情av| 国产三级在线视频| 日韩欧美免费精品| 又爽又黄无遮挡网站| 亚洲欧美成人精品一区二区| av中文乱码字幕在线| 日本精品一区二区三区蜜桃| 禁无遮挡网站| 久久久精品大字幕| 中文亚洲av片在线观看爽| 亚洲第一区二区三区不卡| 久久韩国三级中文字幕| www.色视频.com| 搡女人真爽免费视频火全软件 | 日韩人妻高清精品专区| 不卡视频在线观看欧美| 久久久久国产精品人妻aⅴ院| 日韩成人伦理影院| 久久亚洲国产成人精品v| 国产精品一区二区免费欧美| 精品久久久久久久久久免费视频| 亚洲在线自拍视频| 联通29元200g的流量卡| 国产精品电影一区二区三区| 天堂网av新在线| 国产伦在线观看视频一区| 精品不卡国产一区二区三区| 男人舔奶头视频| 91麻豆精品激情在线观看国产| 亚洲精品色激情综合| 色吧在线观看| 国产一区二区亚洲精品在线观看| 亚洲精品色激情综合| 国产亚洲91精品色在线| 变态另类丝袜制服| 综合色av麻豆| av视频在线观看入口| 俺也久久电影网| 国产精品伦人一区二区| 久久久精品大字幕| 我的老师免费观看完整版| 色吧在线观看| 日韩 亚洲 欧美在线| 99精品在免费线老司机午夜| 亚洲精品在线观看二区| 丝袜美腿在线中文| av在线老鸭窝| 22中文网久久字幕| 日韩强制内射视频| 国产单亲对白刺激| 99久久九九国产精品国产免费| 一级a爱片免费观看的视频| 免费大片18禁| 久久精品国产亚洲av天美| 三级经典国产精品| 亚洲av成人av| 99热全是精品| 成人高潮视频无遮挡免费网站| 少妇丰满av| 日日摸夜夜添夜夜爱| 成年av动漫网址| 国语自产精品视频在线第100页| 看非洲黑人一级黄片| 欧美在线一区亚洲| 精品一区二区三区av网在线观看| 日韩人妻高清精品专区| av在线播放精品| 成熟少妇高潮喷水视频| 婷婷亚洲欧美| 日本精品一区二区三区蜜桃| 亚洲成a人片在线一区二区| 在线播放无遮挡| 色在线成人网| 精品久久久久久久末码| 亚洲欧美日韩高清在线视频| 免费看a级黄色片| 偷拍熟女少妇极品色| 又爽又黄无遮挡网站| 国内精品一区二区在线观看| 日韩欧美三级三区| 麻豆精品久久久久久蜜桃| 午夜日韩欧美国产| 亚洲人成网站在线播放欧美日韩| 我的老师免费观看完整版| 91久久精品电影网| 我的女老师完整版在线观看| 偷拍熟女少妇极品色| 日本爱情动作片www.在线观看 | 天堂网av新在线| 国产色婷婷99| 97超碰精品成人国产| 久久精品国产亚洲av涩爱 | 国产精品人妻久久久久久| 91狼人影院| 人妻制服诱惑在线中文字幕| 日韩欧美 国产精品| 草草在线视频免费看| 两个人视频免费观看高清| 午夜激情福利司机影院| 内地一区二区视频在线| 国产激情偷乱视频一区二区| 国产精品不卡视频一区二区| 人妻夜夜爽99麻豆av| 久久精品国产自在天天线| av黄色大香蕉| 国产aⅴ精品一区二区三区波| 久久精品国产亚洲av涩爱 | 国产精品人妻久久久久久| 欧美+日韩+精品| 欧美最新免费一区二区三区| 国产免费一级a男人的天堂| 成人三级黄色视频| 波多野结衣高清作品| 免费看美女性在线毛片视频| 亚洲av免费高清在线观看| a级毛色黄片| 午夜久久久久精精品| 看免费成人av毛片| 精品午夜福利在线看| 久久精品夜色国产| 少妇高潮的动态图| 日本精品一区二区三区蜜桃| 亚洲一区高清亚洲精品| 免费观看在线日韩| 九九热线精品视视频播放| 最新在线观看一区二区三区| 全区人妻精品视频| 深爱激情五月婷婷| av黄色大香蕉| 亚洲精品在线观看二区| 最后的刺客免费高清国语| 99热精品在线国产| 日韩精品有码人妻一区| 成人av在线播放网站| 日日摸夜夜添夜夜添av毛片| 草草在线视频免费看| 99在线视频只有这里精品首页| 在线国产一区二区在线| 婷婷色综合大香蕉| av女优亚洲男人天堂| 亚洲无线观看免费| 国产成年人精品一区二区| 国产精品美女特级片免费视频播放器| 国产伦在线观看视频一区| 婷婷六月久久综合丁香| 国产成人aa在线观看| 亚洲欧美日韩东京热| 丰满的人妻完整版| 久久人人爽人人片av| 看免费成人av毛片| 亚洲av免费在线观看| 免费在线观看成人毛片| 插逼视频在线观看| 亚洲av一区综合| 久久久久精品国产欧美久久久| 白带黄色成豆腐渣| 久久久久久九九精品二区国产| 99久久精品热视频| 精品免费久久久久久久清纯| 国产精华一区二区三区| 人妻久久中文字幕网| 国产高清三级在线| 国产激情偷乱视频一区二区| 网址你懂的国产日韩在线| 日韩av不卡免费在线播放| 亚洲精品色激情综合| 一级毛片我不卡| 高清毛片免费看| 欧美激情久久久久久爽电影| eeuss影院久久| 久久国内精品自在自线图片| 久久亚洲国产成人精品v| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲18禁久久av| 桃色一区二区三区在线观看| 成年女人看的毛片在线观看| 成人午夜高清在线视频| 国产大屁股一区二区在线视频| 神马国产精品三级电影在线观看| 午夜免费激情av| 国产精品人妻久久久久久| 日韩制服骚丝袜av| 亚洲成a人片在线一区二区| 亚洲成人av在线免费| 有码 亚洲区| 国产v大片淫在线免费观看| 国产成人精品久久久久久| 亚洲av一区综合| 男女那种视频在线观看| 婷婷精品国产亚洲av在线| 久久综合国产亚洲精品| 一卡2卡三卡四卡精品乱码亚洲| av.在线天堂| 综合色丁香网| 人人妻人人看人人澡| 国产精品av视频在线免费观看| 欧美xxxx性猛交bbbb| 免费无遮挡裸体视频| 女生性感内裤真人,穿戴方法视频| 晚上一个人看的免费电影| 免费在线观看影片大全网站| 99久久无色码亚洲精品果冻| 国产免费男女视频| a级毛片a级免费在线| 人妻久久中文字幕网| 91久久精品电影网| 国产亚洲av嫩草精品影院| 欧美绝顶高潮抽搐喷水| a级毛片免费高清观看在线播放| 欧美国产日韩亚洲一区| 久久国产乱子免费精品| 秋霞在线观看毛片| 国产精品久久久久久久电影| 亚洲国产精品成人综合色| 亚洲丝袜综合中文字幕| 亚洲av美国av| 亚洲精品456在线播放app| 99久久久亚洲精品蜜臀av| 丝袜美腿在线中文| 久久精品国产亚洲av涩爱 | 欧美xxxx性猛交bbbb| av天堂中文字幕网| 亚洲av熟女| 欧美潮喷喷水| 五月玫瑰六月丁香| 久久精品国产鲁丝片午夜精品| 欧美日本亚洲视频在线播放| 春色校园在线视频观看| 成人美女网站在线观看视频| 亚洲色图av天堂| 精华霜和精华液先用哪个| 日日干狠狠操夜夜爽| 国产黄色视频一区二区在线观看 | 搡老岳熟女国产| 成人漫画全彩无遮挡| aaaaa片日本免费| 亚洲av成人精品一区久久| 亚洲色图av天堂| 97人妻精品一区二区三区麻豆| 亚洲最大成人av| 国产真实乱freesex| 日韩欧美国产在线观看| 亚洲av免费在线观看| 日本一二三区视频观看| 欧美潮喷喷水| 精品人妻视频免费看| 成人无遮挡网站| 国产欧美日韩精品一区二区| 国产精品无大码| 亚洲精品国产成人久久av| 亚洲欧美精品综合久久99| 免费av不卡在线播放| 亚洲18禁久久av| 一级毛片我不卡| 老师上课跳d突然被开到最大视频| 亚洲欧美精品自产自拍| 女生性感内裤真人,穿戴方法视频| 少妇熟女欧美另类| 国产在线男女| 久久精品人妻少妇| 国产精品一区二区三区四区久久| 日韩高清综合在线| 老女人水多毛片| 欧美区成人在线视频| 99热只有精品国产| 欧美色视频一区免费| 日韩在线高清观看一区二区三区| av在线观看视频网站免费| 一区二区三区四区激情视频 | 亚洲精品乱码久久久v下载方式| 欧美日韩一区二区视频在线观看视频在线 | 欧美日本视频| 综合色av麻豆| 中文字幕av在线有码专区| 久久精品国产清高在天天线| 亚洲不卡免费看| 无遮挡黄片免费观看| 久99久视频精品免费| 欧美三级亚洲精品| 午夜免费激情av| h日本视频在线播放| 一夜夜www| 人人妻人人看人人澡| 99在线人妻在线中文字幕| 老女人水多毛片| 亚洲美女视频黄频| 九九爱精品视频在线观看| 午夜影院日韩av| 亚洲欧美中文字幕日韩二区| 亚洲电影在线观看av| 99久久无色码亚洲精品果冻| h日本视频在线播放| 亚洲av第一区精品v没综合| 97超碰精品成人国产| 国产一区二区激情短视频| 亚洲av成人精品一区久久| 国产精品av视频在线免费观看| 91精品国产九色| 午夜福利在线在线| 国产黄色视频一区二区在线观看 | 国产精品1区2区在线观看.| 欧美性猛交╳xxx乱大交人| 欧美最黄视频在线播放免费| 国产精品,欧美在线| 欧洲精品卡2卡3卡4卡5卡区| 男人舔女人下体高潮全视频| 久久草成人影院| av天堂在线播放| 午夜福利在线观看吧| 91在线观看av| 国产av在哪里看| 免费人成视频x8x8入口观看| 亚洲第一电影网av| 国产精品久久久久久久电影| 黑人高潮一二区| 久久人人爽人人片av| 国产午夜精品论理片| 国产午夜精品久久久久久一区二区三区 | 少妇人妻一区二区三区视频| 欧美成人一区二区免费高清观看| 看片在线看免费视频| 久久这里只有精品中国| 在线观看午夜福利视频| 亚洲成av人片在线播放无| 精品一区二区免费观看| 桃色一区二区三区在线观看| 国产精品无大码| 干丝袜人妻中文字幕| 最好的美女福利视频网| 久久草成人影院| 美女大奶头视频| 啦啦啦啦在线视频资源| 亚洲精品日韩av片在线观看| 欧美另类亚洲清纯唯美| www.色视频.com| 国产乱人偷精品视频| 少妇丰满av| 精品久久久久久成人av| 99久久精品热视频| 久久人妻av系列| 丰满人妻一区二区三区视频av| 又爽又黄a免费视频| 免费av毛片视频| 直男gayav资源| 菩萨蛮人人尽说江南好唐韦庄 | 俺也久久电影网| 九九在线视频观看精品| 禁无遮挡网站| 成人毛片a级毛片在线播放| 国产人妻一区二区三区在| 国产精品一区二区性色av| 亚洲国产色片| 欧美+亚洲+日韩+国产| 午夜激情欧美在线| 国产午夜精品久久久久久一区二区三区 | 午夜免费男女啪啪视频观看 | 亚洲色图av天堂| 69av精品久久久久久| 淫妇啪啪啪对白视频| 国产一区二区在线观看日韩| 国产精品爽爽va在线观看网站| 搡女人真爽免费视频火全软件 | 人妻丰满熟妇av一区二区三区| 精华霜和精华液先用哪个| 成人av在线播放网站| 高清午夜精品一区二区三区 | 欧美成人一区二区免费高清观看| 可以在线观看毛片的网站| 色在线成人网| 日韩亚洲欧美综合| 国产黄片美女视频| 久久久久久大精品| 久久久久久久久久久丰满| 国产三级中文精品| 久久久久性生活片| 成人三级黄色视频| 精品日产1卡2卡| 看非洲黑人一级黄片| 欧美极品一区二区三区四区| 日韩精品有码人妻一区| 久久精品国产亚洲av天美| 久久婷婷人人爽人人干人人爱| 免费黄网站久久成人精品| 欧美日韩精品成人综合77777| 女人被狂操c到高潮| 国产av在哪里看| 男人的好看免费观看在线视频| 男女那种视频在线观看| 高清毛片免费看| 国产老妇女一区| 国国产精品蜜臀av免费| 精品久久久久久成人av| 免费观看在线日韩| 黄色视频,在线免费观看| 超碰av人人做人人爽久久| 国产单亲对白刺激| 在线a可以看的网站| 色综合站精品国产| 直男gayav资源| ponron亚洲| 免费搜索国产男女视频| 久久人妻av系列| 在线观看av片永久免费下载| 两个人视频免费观看高清| 国产伦精品一区二区三区四那| 深夜a级毛片| 国产淫片久久久久久久久| 久久亚洲国产成人精品v| 18禁黄网站禁片免费观看直播| 99久久精品国产国产毛片| 亚洲欧美日韩东京热| 亚洲最大成人av| 99热这里只有是精品50| 国内揄拍国产精品人妻在线| 日本黄色片子视频| 91久久精品电影网| 午夜免费激情av| 91狼人影院| 国产成年人精品一区二区| 欧洲精品卡2卡3卡4卡5卡区| 搡老岳熟女国产| 秋霞在线观看毛片| 1000部很黄的大片| 六月丁香七月| 99精品在免费线老司机午夜| 久久精品国产自在天天线| 人妻少妇偷人精品九色| 国产精品一区二区三区四区免费观看 | 露出奶头的视频| 内射极品少妇av片p| 欧美成人a在线观看| 国产极品精品免费视频能看的| 乱系列少妇在线播放| 亚洲美女黄片视频| 国内久久婷婷六月综合欲色啪| 国产一区亚洲一区在线观看| 亚洲人成网站在线观看播放| 免费观看人在逋| 亚洲精品一卡2卡三卡4卡5卡| 欧美中文日本在线观看视频| 搡老岳熟女国产| av黄色大香蕉| a级毛片免费高清观看在线播放| 国产免费男女视频| 亚洲va在线va天堂va国产| 我的老师免费观看完整版| aaaaa片日本免费| 午夜a级毛片| 成人美女网站在线观看视频| 99九九线精品视频在线观看视频| 免费av不卡在线播放| 精品一区二区三区视频在线| 国产 一区精品| 中文字幕人妻熟人妻熟丝袜美| 亚洲精品粉嫩美女一区| 亚洲欧美日韩东京热| 一本久久中文字幕| 日产精品乱码卡一卡2卡三| 亚洲人成网站在线播| 国产高清视频在线观看网站| av在线亚洲专区| 中国美女看黄片| 日本免费一区二区三区高清不卡| 女的被弄到高潮叫床怎么办| 亚洲熟妇中文字幕五十中出| 国产免费男女视频| 亚洲欧美清纯卡通| 2021天堂中文幕一二区在线观| 日本三级黄在线观看| 久久精品影院6| 国产精品亚洲美女久久久| 成人无遮挡网站| 嫩草影视91久久| 精品久久久久久成人av| 九色成人免费人妻av| 国产成人freesex在线 | 色在线成人网| 国产 一区 欧美 日韩| 久久精品综合一区二区三区| 国产精品国产高清国产av| 精品欧美国产一区二区三| 国产免费一级a男人的天堂| 亚洲国产精品成人久久小说 | 久久九九热精品免费| 亚洲经典国产精华液单| 久久韩国三级中文字幕| 波野结衣二区三区在线| 一区二区三区四区激情视频 | 91狼人影院| or卡值多少钱| 日韩在线高清观看一区二区三区| 淫秽高清视频在线观看| 日韩欧美在线乱码| 一级毛片电影观看 | 少妇人妻精品综合一区二区 | 97超级碰碰碰精品色视频在线观看| 免费观看在线日韩| 午夜福利在线观看吧| 99久久精品一区二区三区| 久久精品国产亚洲av香蕉五月| 国产精品久久电影中文字幕| 天天躁夜夜躁狠狠久久av| 亚洲精品影视一区二区三区av| 免费av不卡在线播放| 性色avwww在线观看| 亚洲中文字幕日韩| 国产精品美女特级片免费视频播放器| 最新中文字幕久久久久| 大又大粗又爽又黄少妇毛片口| 日韩成人av中文字幕在线观看 | 久久精品国产亚洲网站| 久久久久久伊人网av| 91在线观看av| 男人的好看免费观看在线视频| 国产三级在线视频| 又黄又爽又刺激的免费视频.| 免费大片18禁| 寂寞人妻少妇视频99o| 久久精品国产99精品国产亚洲性色| 我的女老师完整版在线观看| 亚洲性久久影院| 给我免费播放毛片高清在线观看| 午夜福利18| 淫妇啪啪啪对白视频| а√天堂www在线а√下载| 人妻夜夜爽99麻豆av| 成人亚洲精品av一区二区| 九色成人免费人妻av| 欧美不卡视频在线免费观看| 精品人妻熟女av久视频| 久久99热这里只有精品18| 国产精品人妻久久久久久| 午夜老司机福利剧场| 国产三级在线视频| 淫妇啪啪啪对白视频| 一区二区三区四区激情视频 | a级毛片免费高清观看在线播放| 国产精品1区2区在线观看.| 精品无人区乱码1区二区| 色吧在线观看| 欧美绝顶高潮抽搐喷水| 亚洲国产精品成人久久小说 | 18禁在线无遮挡免费观看视频 | 在线播放国产精品三级| 欧美性猛交╳xxx乱大交人| 免费观看在线日韩| 精品午夜福利视频在线观看一区| 女生性感内裤真人,穿戴方法视频| 人人妻人人澡人人爽人人夜夜 | 日本免费一区二区三区高清不卡| 亚洲图色成人| 欧美一区二区亚洲| 国产精品一区二区性色av| 99热这里只有精品一区| 成人特级av手机在线观看| 国产熟女欧美一区二区| 成人鲁丝片一二三区免费| 久久人人爽人人爽人人片va| 网址你懂的国产日韩在线| 国产日本99.免费观看| 亚洲一级一片aⅴ在线观看| 精品少妇黑人巨大在线播放 | 美女cb高潮喷水在线观看| 五月伊人婷婷丁香| 成人精品一区二区免费| 嫩草影院入口| 色在线成人网| 欧美日韩综合久久久久久| aaaaa片日本免费| 久久久久久久久久黄片| 日韩精品中文字幕看吧| 精品无人区乱码1区二区| 亚州av有码| 老熟妇乱子伦视频在线观看| 麻豆乱淫一区二区| 天堂动漫精品| 人人妻人人澡人人爽人人夜夜 | 日韩人妻高清精品专区| 一本精品99久久精品77| 亚洲七黄色美女视频| 此物有八面人人有两片| 又黄又爽又刺激的免费视频.| 99国产精品一区二区蜜桃av| 亚洲精品日韩在线中文字幕 | 国产免费一级a男人的天堂| 日本在线视频免费播放| 偷拍熟女少妇极品色| 国产精品综合久久久久久久免费| 又爽又黄a免费视频| 最近中文字幕高清免费大全6| 欧美+亚洲+日韩+国产| 99热这里只有是精品50|