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

    安太堡露天礦區(qū)復墾地植被覆蓋度反演估算研究

    2016-12-19 23:46:16耀,周
    中南林業(yè)科技大學學報 2016年11期
    關(guān)鍵詞:排土場植被指數(shù)覆蓋度

    張 耀 ,周 偉

    (1.中國地質(zhì)大學(北京) 土地科學技術(shù)學院,北京 100083;2.國土資源部 土地整治重點實驗室,北京 100035)

    安太堡露天礦區(qū)復墾地植被覆蓋度反演估算研究

    張 耀1,2,周 偉1,2

    (1.中國地質(zhì)大學(北京) 土地科學技術(shù)學院,北京 100083;2.國土資源部 土地整治重點實驗室,北京 100035)

    長期露天開采造成了植被大面積破壞、土地壓占損毀等環(huán)境問題。利用遙感技術(shù)準確、快速監(jiān)測礦區(qū)植被覆蓋度,對植被恢復和土地復墾監(jiān)管具有重要意義。以山西平朔安太堡露天煤礦為研究對象,利用多時相TM數(shù)據(jù)結(jié)合實地調(diào)查數(shù)據(jù),采用主成分分析法,從19個植被指數(shù)因子中提取主成分,建立植被覆蓋度反演模型。結(jié)果表明:(1)前兩個主成分PC1和PC2的累積貢獻率為91.508%,建立復墾植被覆蓋度反演模型,經(jīng)檢驗R2=0.742,達到顯著相關(guān)水平。(2)1990~2010年,南、西排土場植被覆蓋度分別增加了20.88%和12.89%,西擴排土場和內(nèi)排土場植被覆蓋度分別降低了18.23%和7.59%,南、西排土場植被覆蓋度和植被覆蓋增加量都大于西擴排土場和內(nèi)排土場的。(3)南、西排土場植被復墾時間都超過了10年,而西擴排土場和內(nèi)排土場的復墾時間少于10年,說明復墾周期長有利于植被恢復,進行復墾工程對礦區(qū)生態(tài)環(huán)境的修復具有直接作用和重要意義。

    植被覆蓋度;植被指數(shù);復墾植被;TM數(shù)據(jù);露天礦

    露天礦開采損毀原始地貌的植被覆蓋、降低土壤養(yǎng)分,嚴重損壞礦區(qū)及周邊的生態(tài)系統(tǒng)[1]。礦區(qū)土地復墾是將采礦對環(huán)境損毀降至最小、恢復當?shù)厣鷳B(tài)平衡的重要手段,是礦山生產(chǎn)建設中的重要環(huán)節(jié)[2]。復墾植被作為復墾最重要步驟之一,其覆蓋度的動態(tài)變化是復墾植被健康狀況的直接表現(xiàn),是評價復墾效果比較直觀和有效的指標[3]。植被復墾是一個長期的工程,隨著復墾時間的推移,會出現(xiàn)植被物種的更替。因此,根據(jù)植被在生態(tài)環(huán)境中的演化特征,要求復墾植被的監(jiān)測具有動態(tài)性和長期性。持續(xù)的、長期的監(jiān)測復墾植被恢復信息成為土地復墾工程的關(guān)鍵[4]。

    傳統(tǒng)樣地采樣方式受時空尺度的限制,參考樣地的選擇有限,且地面監(jiān)測復墾植被的樣方難以根據(jù)復墾植被年限的增長而變化,很難完成長周期多尺度監(jiān)測的目標[5]。而遙感技術(shù)具有覆蓋廣、采樣周期短、可動態(tài)地掌握信息等特點,已成為動態(tài)監(jiān)測區(qū)域植被覆蓋度的重要手段[6]。其中,植被指數(shù)法是從遙感圖像上快速、動態(tài)獲取大范圍植被信息常用的經(jīng)濟有效的方法,通過獲取植被生長狀態(tài)、植被參數(shù)反演、植被時序變化監(jiān)測等信息,來監(jiān)測植被生長狀態(tài)、綠色植被活力等信息,是定量研究植被覆蓋特征的好方法[7]。植被覆蓋度與植被指數(shù)耦合關(guān)系研究有:Purevdorj等通過研究蒙古和日本草原,對NDVI、土壤調(diào)節(jié)植被指數(shù)SAVI等5個植被指數(shù)和植被覆蓋度的相關(guān)性和多項式回歸分析,得到NDVI、SAVI與植被覆蓋度相關(guān)性較高的結(jié)論[8]。黎良財?shù)扔嬎鉔DVI、差值植被指數(shù)DVI、SR等5種植被指數(shù)結(jié)合相機拍攝樣方提取的植被覆蓋度進行分析,建立回歸模型[9]。

    國內(nèi)外學者以長時間序列Landsat影像為數(shù)據(jù)源,利用植被指數(shù)、主成分分析等方法分別對不同時空尺度的礦區(qū)植被覆蓋變化進行了研究[10-12]。研究發(fā)現(xiàn):Erener等用Landsat數(shù)據(jù)NDVI、SR研究土耳其Seyit?mer煤礦,分析表明復墾地植被覆蓋總體呈上升趨勢[13]。Halounova用歸一化植被指數(shù)NDVI、簡單比值指數(shù)SR等8個植被指數(shù)對捷克波西米亞北露天礦的100個復墾森林類型進行對比研究,評價10年間復墾植被恢復情況[14]。Susmita等用NDVI和纓帽變換研究美國阿巴拉契亞南部高植被覆蓋的露天煤礦區(qū),分析得到森林經(jīng)過擾動與恢復轉(zhuǎn)變的閾值[15]。黎良財?shù)扔肗DVI像元二分模型對鉛鋅礦區(qū)植被覆蓋度進行了時空變化分析[16]。

    針對生態(tài)極其脆弱的黃土高原露天礦區(qū),人工擾動大,破碎化程度高,篩選適宜的植被指數(shù)來反演植被覆蓋度,是定量植被遙感的關(guān)鍵技術(shù)。然而,研究露天礦復墾區(qū)域植被覆蓋的研究不多,尤其是用遙感技術(shù)做長期定量監(jiān)測的研究較少。所以結(jié)合植被地面調(diào)查,包括復墾植被覆蓋度、物種以及群落結(jié)構(gòu),通過對比復墾植被和原地貌自然植被的這些指標,來評價復墾植被是否達到標準。因此,本文以平朔生態(tài)脆弱區(qū)為研究對象,利用1990~2010年多時相TM影像,通過反演礦區(qū)植被覆蓋度,對復墾植被和當?shù)刂脖坏纳L狀況進行長期、動態(tài)監(jiān)測,分析其時空變化情況,為礦區(qū)土地復墾監(jiān)管的信息化提供支持,為礦區(qū)環(huán)境治理提供決策支持。

    1 研究區(qū)概況

    平朔礦區(qū)地處黃土高原東部,與晉陜蒙“黑三角”接壤,位于山西省朔州市平魯區(qū),112.29o~112.46oE,39.40o~ 39.54oN(圖 1)。

    圖1 山西平朔位置與礦區(qū)示意Fig.1 Location of Pingshuo mine in Shanxi province

    本區(qū)地貌景觀以黃土低山丘陵臺地為主,海拔1 200~1 600 m,地勢東北高,西南低,地表水系屬于海河流域的永定河水系,地表水資源貧乏。本區(qū)屬于半干旱大陸性溫帶季風氣候區(qū),年平均氣溫6.0℃,年降水量450.0 mm,年蒸發(fā)量2 160.0 mm。本區(qū)屬于中國植被區(qū)劃中的VI溫帶草原區(qū)域VIA東部草原亞區(qū)域VIAii溫帶北部草原地帶VIAiia溫帶南部森林(草甸)草原亞地帶VIAiia-4晉北山地森林草原區(qū)。植被以旱生和旱中生草原植物區(qū)系成分為主,植被覆蓋率低,主要植被群落類型有長芒草、扁穗冰草、克氏針茅等耐旱植物。礦區(qū)農(nóng)田栽培作物屬一年一熟制,主要農(nóng)作物有玉米、谷子、莜麥等。本區(qū)主要地帶性土壤為栗鈣土、栗褐土、黃綿土。平朔礦區(qū)屬寧武煤田北部區(qū)域,煤炭儲量127.5億t,開采時間1986~2095年,持續(xù)時間約為 110 年。平朔煤礦有三個露天礦:安太堡露天礦,安家?guī)X露天礦,東露天礦。本文主要研究區(qū)域為安太堡露天礦,其開采與復墾基本情況見表1。

    2 數(shù)據(jù)獲取與研究方法

    2.1 復墾植被覆蓋度的野外調(diào)查

    復墾區(qū)域的植被覆蓋度獲取對于模型建立和反演非常重要。選擇安太堡南排土場的四塊樣地(SⅠ、SⅢ、SⅣ、SⅤ)為樹木植被覆蓋度調(diào)查樣地,樣地被國土資源部列為復墾示范基地,樣地基本情況見表2。野外調(diào)查時間為2010年8月,以10m×10m樣方為樹種編號單元,調(diào)查內(nèi)容包括樹木的樹種名稱、高度、冠幅、胸徑、物候期和地理位置等信息,并建立數(shù)據(jù)庫。由于需要跟TM影像匹配分析,需要對樣方進行重采樣,生成30m×30 m網(wǎng)格,每個網(wǎng)格在TM影像上為1個像元大小,每個網(wǎng)格內(nèi)有9個原樣方(3×3),重采樣后樣方數(shù)量為32個。野外調(diào)查樹木植被覆蓋度是樹木冠層垂直投影面積占樣地面積的百分比,其取值范圍為0~100%[17]。

    表1 安太堡排土場的開采與復墾基本情況Table 1 Dumps area and the time of mining and reclamation

    表2 野外調(diào)查樣地信息概況?Table 2 Description of the sample plots

    2.2 遙感數(shù)據(jù)源

    采用安太堡礦區(qū)處于植被生長季節(jié)的Landsat TM數(shù)據(jù)(美國地質(zhì)調(diào)查局提供)作為數(shù)據(jù)源,數(shù)據(jù)獲取日期為1990-08-22、2000-08-10、2000-09-11、2000-09-18、2010-07-05、2010-07-12、2010-09-14。遙感影像研究區(qū)域無云覆蓋,已經(jīng)過輻射定標、大氣校正和幾何校正,數(shù)據(jù)質(zhì)量較好。為了減少物候變化對植被生長狀況的影響,取數(shù)據(jù)每年植被生長季的不同月份的均值為該年數(shù)值。利用外業(yè)調(diào)查范圍作為邊界,裁剪出研究區(qū)域TM影像。

    2.3 植被指數(shù)的計算

    植被指數(shù)是一種無量綱指數(shù),是多個波段光譜反射率的線性或非線性組合,從而將多光譜中的重要光譜信息壓縮到一個植被指數(shù)通道,實現(xiàn)多光譜數(shù)據(jù)的降維。遙感圖像上的植被信息,是通過植物葉子和植被冠層的光譜特征反映的??梢姽饩G光波段(520~590 nm)對區(qū)分植物類別敏感,紅光波段(630~690 nm)對植被覆蓋度、植被生長狀況敏感。可見光譜段(380~760 nm)受葉子葉綠素含量控制,近紅外譜段(760~1 100 nm)受葉內(nèi)細胞結(jié)構(gòu)控制,短波紅外譜段(1 300~2 500 nm)受葉細胞內(nèi)水分含量控制[17]。

    目前,已研究發(fā)展了40多種不同的植被指數(shù),每一種植被指數(shù)創(chuàng)立都有明確的理論依據(jù)和應用目的。最常用的歸一化植被指數(shù)NDVI能定量表達植被覆蓋信息,適用于植被發(fā)育中期或中等覆蓋度的植被監(jiān)測。結(jié)合植被在短波紅外波段反射率受到水分強吸收的特點,歸一化水分指數(shù)NDWI來更有效的反演植被冠層水分含量[18]。由于受到土壤背景和大氣噪聲的影響,對不同覆蓋區(qū)的出現(xiàn)敏感度不同的情況,需要加入不同參數(shù)來消除其影響的植被指數(shù)。在降低土壤背景的影響方面,SAVI引入土壤調(diào)節(jié)參數(shù)L能夠隨植被密度變化自動調(diào)節(jié)來減小土壤的影響[19];增強型植被指數(shù)EVI通過加入藍波段增強植被信號,校正土壤背景和氣溶膠散射的影響,更為適用于植被茂密地區(qū)[20]。在減少大氣散射的影響方面,大氣阻抗植被指數(shù)ARVI通過歸一化藍、紅和近紅外波段的輻射來減少大氣效應的敏感程度[21]。本文選用礦區(qū)常用的植被指數(shù)見表3。

    表3 植被指數(shù)計算公式Table 3 Equation of vegetation indices for Landsat image

    2.4 纓帽變換

    纓帽變換也成Kauth-Thomas變換(KT),該變換是一種特殊的主成分分析,是基于多波段圖像的線型變換,是一種全球性的植被指數(shù)。KT變換可從TM圖像中分解出每個像元的土壤亮度指數(shù)(KTB,反映土壤信息)、綠度指數(shù)(KTG,反映生物量特征,與圖像上植被數(shù)量相關(guān))和濕度指數(shù)(KTW,與冠層和土壤濕度相關(guān))[22]。計算公式:

    式中,Rblue、Rgreen、Rred、Rnir分別為藍、綠、紅和近紅外波段的反射率,R1.6um為1.6 um的短波紅外反射率,R2.2um為2.2 um的短波紅外反射率。

    3 結(jié)果與分析

    3.1 復墾排土場植被覆蓋度估測模型構(gòu)建

    用SPSS軟件對樣地的19個植被指數(shù)因子進行相關(guān)性分析(表4),大多數(shù)植被指數(shù)與植被覆蓋度顯著相關(guān)(|Pearson| >0.7,P<0.05),BGI與植被覆蓋度不相關(guān)(P>0.05),結(jié)果說明因子間可能存在多重相關(guān)性。在復墾區(qū)域植被覆蓋度的定量評價中,需要采用多元線性回歸或逐步回歸方法建立植被覆蓋度與植被指數(shù)回歸模型。這種方法簡便易行,但回歸模型中的相關(guān)變量會影響模型精度,甚至出現(xiàn)病態(tài)模型,因此采用主成分分析法,篩選植被指數(shù)因子,生成主成分分量。由于主成分分量間具有相互正交的特性,多個主成分分量之間相互獨立,這樣既可保留不同植被指數(shù)信息的主要特征,又可避免變量之間多重相關(guān)。因此,本文采用主成分分析篩選植被指數(shù)的主成分,再建立主成分與植被蓋度的回歸模型,用于估測復墾地植被覆蓋度。

    對樣地的19個植被指數(shù)因子進行主成分分析(表4,圖2),由表4看出,前2個主成分的樣本方差已達到91.508%,反映了樣本的主要信息,之后的主成分分量反映的信息量較少。信息量主要由PC1和PC2集中反映,變量由多個因子降為2個綜合變量,通過主成分變換起到了保留樣本主要信息和降維作用。因子載荷分析表明,第一主成分PC1反映大部分信息,貢獻率為84.656%,與PC1顯著相關(guān)的因子有:AFRI、ARVI、DVI、EVI、KTG、KTW、NDBI、NDVI、NDWI、PVI、RDVI、RGI、SR、TVI、VARI;第二主成分PC2反映部分信息,貢獻率為6.852%,與PC2相關(guān)的因子有:與土壤相關(guān)因子SAVI和KTB,與大氣相關(guān)因子VARI。

    表4 相關(guān)性系數(shù)、公因子方差提取和主成分因子載荷、特征值及貢獻率Table 4 Coefficient, common factor variance extraction,factors loading value, eigenvalue and contributiveness

    圖2 旋轉(zhuǎn)空間的成分圖Fig 2 Composition diagram of rotation space

    以主成分為自變量,將野外調(diào)查得到的32個樣方的植被覆蓋度為因變量,進行多元線性回歸分析,建立植被覆蓋度估測模型。計算公式為:

    式中:Y為復墾地植被覆蓋度(%)。模型方差分析及線型回歸關(guān)系顯著性檢驗結(jié)果見圖3、表5?;谶b感植被指數(shù)的主成分和礦區(qū)復墾排土場植被覆蓋度具有顯著相關(guān)關(guān)系,相關(guān)系數(shù)R2=0.742,說明反演的效果較好。但低植被覆蓋度區(qū)域反演蓋度比實測覆蓋度高,這可能與兩者的尺度差異和實際測定的方法有關(guān)。高植被覆蓋區(qū)域點的離散度較高,說明高植被覆蓋條件下反演結(jié)果可能受到植被密度、植被光學厚度等參數(shù)的影響較大。

    圖3 反演植被覆蓋度與實測植被覆蓋度的關(guān)系Fig. 3 Relationship between fi eld measured vegetation fraction and modeled vegetation fraction

    表5 方差分析Table 5 Variance analysis

    3.2 植被覆蓋度反演結(jié)果分析

    2010年,基于主成分分析植被覆蓋度反演結(jié)果,見表6、圖4。2010年,南排土場反演植被覆蓋度為98.85%,這與礦區(qū)實地調(diào)查的大于90%結(jié)果一致;西排土場估測植被覆蓋度為99.81%,這與礦區(qū)實地調(diào)查該區(qū)域植被覆蓋度大于95%結(jié)果一致。這說明,基于主成分分析植被覆蓋度反演模型取得良好結(jié)果,可用于分析研究區(qū)植被覆蓋狀況。

    表6 1990~2010年各排土場植被覆蓋度估測值Table 6 Vegetation fraction of dumps from 1990 to 2010

    圖4 1990~2010年各排土場植被覆蓋度估測趨勢Fig.4 Vegetation fraction of dumps from 1990 to 2010

    由表6、圖4可知,2010年礦區(qū)西擴排土場、內(nèi)排土場反演植被覆蓋度分別為81.77%、79.87%,南排土場和西排土場植被覆蓋優(yōu)于西擴排土場和內(nèi)排土場,說明西擴、內(nèi)排土場植被復墾情況弱于南、西排土場復墾情況。這是由于南排土場從1992-1997年開始復墾,西排土場從1994~1997年開始復墾,截止2010年兩個排土場的復墾時間都超過了10年。南排土場植被覆蓋度從1990年的77.97%升至2010年的98.85%,20年增加了20.88%。西排土場植被覆蓋度從1990年的86.92%增至2010年的99.81%,20年增加了12.89%,說明兩地植被恢復情況較好,復墾工程基本完成。西擴排土場1993年開始開采,1990年植被覆蓋度與原地貌林地覆蓋度一致,為100%;2000年降至84.53%,到了2010年降至81.77%,20年植被覆蓋度降低了18.23%,西擴排土場2001~2008年開始復墾,由于復墾時間較短,植被覆蓋度的變化不顯著,需要繼續(xù)進行植被的恢復工程。內(nèi)排土場1997年開始開采,1990年植被覆蓋度為87.46%,低于原地貌林地覆蓋度,這是由于內(nèi)排土場于1987年安太堡煤礦開采之初就作為生產(chǎn)區(qū)域,以礦山建筑和運煤通道為主,植被覆蓋已遭到破壞;2000年降至76.20%,2010年增至79.87%,這說明復墾植被已使區(qū)域內(nèi)植被覆蓋增加,但由于復墾時間較短和開采還在進行的情況,該區(qū)域植被覆蓋度較低,需要持續(xù)復墾。

    4 結(jié)論與討論

    (1)以安太堡露天礦為研究對象,采用主成分分析法,篩選提取了19個植被指數(shù)因子的主成分,利用得到的主成分建立了礦區(qū)復墾植被覆蓋度反演模型,經(jīng)過方差分析及相關(guān)性檢驗,達到顯著相關(guān)水平,取得了良好精度。在與植被覆蓋度相關(guān)的多個植被指數(shù)因子中,利用主成分分析方法,篩選主成分,再利用主成分值與復墾林地植被覆蓋度建立估測模型,既可以保留多個植被指數(shù)的主要信息,又可以避免因子間共線性的問題,起到降維、簡化模型的作用,是提取與復墾地植被覆蓋度相關(guān)的植被指數(shù)因子、估測復墾地植被覆蓋度的一種有效方法。

    這種方法可操作性強、技術(shù)路線簡單,也可適用于不同分辨率的數(shù)據(jù),對黃土高原煤礦區(qū)植被的生長信息有良好的監(jiān)測能力,可用于實時、動態(tài)對礦區(qū)復墾進展做長期的監(jiān)測和評價,能滿足礦區(qū)植被覆蓋定量化的研究要求,可以為建立土地復墾監(jiān)管指標體系做基礎數(shù)據(jù)。

    然而,由于實測數(shù)據(jù)是在覆蓋度較高的南排土場,反演結(jié)果在低植被覆蓋區(qū)域精度較低,所以今后的研究中需要加入更多不同覆蓋度的實地驗證數(shù)據(jù)來提高模型精度。

    (2)截止2010年,南排土場從1992~1997年開始復墾,西排土場從1994~1997年開始復墾,兩個排土場的復墾時間都超過10年,1990~2010年南排土場、西排土場的植被覆蓋度分別增加了20.88%和12.89%;西擴排土場1993年開始開采,2001~2008年開始復墾,內(nèi)排土場1997至今進行復墾,兩個排土場的復墾時間較短,1990~2010年西擴排土場、內(nèi)排土場的植被覆蓋度分別降低了18.23%和7.59%,南排土場和西排土場植被覆蓋度和增加量都優(yōu)于西擴排土場和內(nèi)排土場,這說明復墾周期長有利于植被恢復,復墾工程的進行對礦區(qū)生態(tài)環(huán)境的修復具有直接作用和重要意義。對比原地貌植被,經(jīng)過人工復墾的植被長勢持續(xù)良好,整體優(yōu)于未擾動區(qū)的植被。

    (3)平朔煤礦從1990~2010年期間的20年,在不同時期做了大量的生態(tài)補償措施,在一定程度上減輕周邊環(huán)境的影響。平朔礦區(qū)在生態(tài)修復方面成效顯著,有很多成功經(jīng)驗值得其他礦區(qū)借鑒:礦區(qū)采用邊開采邊復墾的技術(shù),土地利用類型呈“原地貌-剝離區(qū)-采場-排土場-復墾區(qū)”周期性轉(zhuǎn)換模式,礦區(qū)持續(xù)穩(wěn)定的環(huán)境保護投入與卓有成效的土地復墾工作,保障了礦區(qū)的生態(tài)環(huán)境修復的持續(xù)與強化。復墾區(qū)人工植被面積增大,隨著時間推移形成了較穩(wěn)定的植被群落,對新增排土場需及時進行植被修復,優(yōu)化植物物種,有助于人工植被向自然植被的演替,有利于黃土高原生態(tài)環(huán)境修復。

    致謝:本文特別感謝印第安納州立大學的謝燕華博士的悉心指導與幫助。

    [1] 胡振琪,龍精華,王新靜.論煤礦區(qū)生態(tài)環(huán)境的自修復、自然修復和人工修復[J].煤炭學報,2014(8): 1751-1757.

    [2] Prakash A, Gupta R P. Land-use mapping and change detection in a coal mining area - a case study in the Jharia coal fi eld, India[J].International Journal of Remote Sensing, 1998,19(3): 391-410.

    [3] Zhang Yan, Yang Jianying, Wu Hai-long, et al. Dynamic changes in soil and vegetation during varying ecological-recovery conditions of abandoned mines in Beijing[J]. Ecological Engineering,2014,73: 676-683.

    [4] Alday Josu G., Santana Víctor M., Marrs Rob H., et al. Shrubinduced understory vegetation changes in reclaimed mine sites[J]. Ecological Engineering, 2014,73: 691-698.

    [5] 宋祥蘭,王蘭英,鄺先松,等.贛南廢棄稀土礦區(qū)植被恢復模式試驗[J].中南林業(yè)科技大學學報, 2015,35(6): 58-62.

    [6] 嚴恩萍,林 輝,黨永峰,等.2000-2012年京津風沙源治理區(qū)植被覆蓋時空演變特征[J].生態(tài)學報, 2014,34(17):5007-5020.

    [7] Ramoelo Abel, Dzikiti Sebinasi, van Deventer Heidi, et al.Potential to monitor plant stress using remote sensing tools[J].Journal of Arid Environments, 2015,113: 134-144.

    [8] Purevdorj T, Tateishi R, Ishiyama T, et al. Relationships between percent vegetation cover and vegetation indices[J]. International Journal Remote Sensing, 1998,19(18): 3519-3535.

    [9] 黎良財,張曉麗,鄧 利,等.廣西巖溶地區(qū)植被覆蓋度與植被指數(shù)耦合模型研究[J].中南林業(yè)科技大學學報, 2015,35(5): 28-34.

    [10] Bakr N, Weindorf D C, Bahnassy M H, et al. Monitoring land cover changes in a newly reclaimed area of Egypt using multi-temporal Landsat data[J]. Applied Geography, 2010,30(4):592-605.

    [11] Larsson H. Linear regression s for canopy cover estimation in Acacia woodlands using Landsat-TM, -MSS and SPOT HRV XS data[J] .International Journal Remote Sensing, 1993, 14(11):2129-2136.

    [12] Prakash A, Gupta R P. Land-use mapping and change detection in a coal mining area-a case study in the Jharia coal fi eld, India[J].International Journal of Remote Sensing, 1998,19(3): 391-410.

    [13] Erener Arzu. Remote sensing of vegetation health for reclaimed areas of Seyit?mer open cast coal mine[J]. International Journal of Coal Geology, 2011,86(1): 20-26.

    [14] Halounová L. Reclamation areas and their development studied by vegetation indices[J]. International Journal of Digital Earth,2008,1(1): 155-164.

    [15] Sen Susmita, Zipper Carl E, Wynne Randolph H, et al.Identifying revegetated mines as disturbance/ recovery trajectories using an interannual Landsat chronosequence[J].Photogrammetric Engineering & Remote Sensing, 2012,78(3):223-235.

    [16] 黎良財,鄧 利,曹 穎,等.基于NDVI像元二分模型的礦區(qū)植被覆蓋動態(tài)監(jiān)測[J].中南林業(yè)科技大學學報,2012,32(6): 18-23.

    [17] Gitelson A, Kaufman YJ, Stark R, et al. Novel algorithms for remote estimation of vegetation fraction[J].Remote Sensing of environment, 2002,80(01): 76-87.

    [18] Gu Yingxin, Hunt Eric, Wardlow Brian, et al. Evaluation of MODIS NDVI and NDWI for vegetation drought monitoring using Oklahoma Mesonet soil moisture data[J]. Geophysical Research Letters, 2008,35(22).

    [19] Huete Alfredo R. A soil-adjusted vegetation index (SAVI)[J].Remote sensing of environment, 1988,25(3): 295-309.

    [20] Li Zhongfei, Li Xiaobing, Wei Dandan, et al. An assessment of correlation on MODIS-NDVI and EVI with natural vegetation coverage in Northern Hebei Province, China[J]. Procedia Environmental Sciences, 2010,2: 964-969.

    [21] Kaufman J, Tanre D. Atmospherically resistant vegetation index(ARVI) for EOS-MODIS[J]. IEEE Transactions on Geoscience and Remote sensing, 1992,30(2): 261-270.

    [22] Crist E P, Kauth R J. The Tasseled Cap De-Mystified[J].Photogrammetric Engineering and Remote Sensing, 1986,52(1):81-86.

    Remote estimation of vegetation fraction for reclaimed areas of Antaibao opencast mine

    ZHANG Yao1,2, ZHOU Wei1,2
    (1. School of Land Science and Technology, China University of Geosciences(Beijing), Beijing 100083, China;2. Key Laboratory of Land Consolidation and Rehabilitation, Ministry of Land Resources, Beijing 100035, China)

    In an area like the Antaibao coal mine, where extensive and rapid opencast mining is going on continuously, long-term reclamation monitoring is important. Accurate information of vegetation fraction is important for assessing how mining activities affect the ecosystem in mining areas. The principal components (PCs) for vegetation indices (VIs) are obtained by principal components analysis (PCA). A remote sensing method based on vegetation cover monitoring and assessment by using Landsat data sets with the temporal coverage from 1990 to 2010 was presented and applied to the Antaibao opencast mine, Shanxi Province, China. To reduce the impact of phenology on remote sensing based vegetation monitoring, 7 Landsat data were selected during vegetation growing season(June~October). A total of 33 pixels were selected as the survey fi elds during vegetation growth season in August 2010. VF of arbores canopy in the selected plots was computed. We monitored reclaimed vegetation variation through analyzing time series data of 19 vegetation indices, which were usually applied to mine area monitoring, such as NDVI, EVI. PCA for the VIs, the accumulative ratio of contribution of the fi rst two PCs is 91.508%,the variance contribution rate of PC1 is 84.656%, and the variance contribution rate of PC2 is 6.852%.The VF model was set up by regression analysis of SPSS based on the PC1 and PC2. T test examination showed that the fi eld VF was correlated signi fi cantly to PC1 and PC2, the correlation coef fi cient R2 of the model was 0.742. The result con fi rmed that the PCA in the inversion of vegetation fraction worked effectively, which could be widely used. In addition, we concluded that the vegetation system of South Dump, West Dump and West Enlarged Dump had gradually transformed to the natural ecosystem from 1990~2015.The vegetation in Inner Dump and Inner Enlarged Dump should be well restored. Reclamation was a long-time and dynamic project,whose effect is hard to determine in short term, while remote sensing technology is an effective way to monitor the inter-annual variation of vegetation recovery.

    Vegetation fraction; Vegetation indices; Revegetation; TM data; Opencast mine

    S771.8

    A

    1673-923X(2016)11-0113-07

    10.14067/j.cnki.1673-923x.2016.11.020

    2016-04-14

    國家自然科學基金資助項目(41571508)

    張 耀,碩士研究生 通訊作者:周 偉,教授;E-mail:zhouw@cugb.edu.cn

    張 耀,周 偉.安太堡露天礦區(qū)復墾地植被覆蓋度反演估算研究[J].中南林業(yè)科技大學學報,2016, 36(11): 113-119.

    [本文編校:吳 毅]

    猜你喜歡
    排土場植被指數(shù)覆蓋度
    呼和浩特市和林格爾縣植被覆蓋度變化遙感監(jiān)測
    基于NDVI的晉州市植被覆蓋信息提取
    石寶鐵礦排土場穩(wěn)定性研究
    排土場的安全防護對策措施與事故分析
    低覆蓋度CO分子在Ni(110)面的吸附研究
    露天礦排土場邊坡穩(wěn)定性分析及評價機制
    AMSR_2微波植被指數(shù)在黃河流域的適用性對比與分析
    河南省冬小麥產(chǎn)量遙感監(jiān)測精度比較研究
    主要植被指數(shù)在生態(tài)環(huán)評中的作用
    西藏科技(2015年1期)2015-09-26 12:09:29
    基于MODIS數(shù)據(jù)的植被指數(shù)與植被覆蓋度關(guān)系研究
    亚洲aⅴ乱码一区二区在线播放| 精品久久久久久久久亚洲| 99精国产麻豆久久婷婷| 人妻制服诱惑在线中文字幕| 亚洲av中文字字幕乱码综合| 丝袜喷水一区| 老女人水多毛片| 尤物成人国产欧美一区二区三区| 男人狂女人下面高潮的视频| 嘟嘟电影网在线观看| 少妇人妻精品综合一区二区| 狂野欧美激情性bbbbbb| 亚洲人与动物交配视频| 日韩 亚洲 欧美在线| 国产美女午夜福利| videos熟女内射| 看免费成人av毛片| 亚洲av成人精品一二三区| 校园人妻丝袜中文字幕| 免费看日本二区| 在线播放无遮挡| 男的添女的下面高潮视频| 99热网站在线观看| 美女cb高潮喷水在线观看| 亚洲不卡免费看| 久久精品国产亚洲av涩爱| 在线观看av片永久免费下载| 九色成人免费人妻av| 欧美成人一区二区免费高清观看| 国产亚洲av片在线观看秒播厂| 国产精品.久久久| 亚洲av中文av极速乱| 日韩av免费高清视频| 三级国产精品片| 波多野结衣巨乳人妻| 纵有疾风起免费观看全集完整版| 久久人人爽人人爽人人片va| 男人和女人高潮做爰伦理| 精品人妻一区二区三区麻豆| 人妻夜夜爽99麻豆av| xxx大片免费视频| 一本久久精品| 老师上课跳d突然被开到最大视频| 亚洲色图av天堂| 看非洲黑人一级黄片| 久久综合国产亚洲精品| 欧美 日韩 精品 国产| 久久鲁丝午夜福利片| 欧美日韩视频精品一区| av在线播放精品| 免费看日本二区| 久久鲁丝午夜福利片| 欧美bdsm另类| 最后的刺客免费高清国语| 丝瓜视频免费看黄片| 伦理电影大哥的女人| 欧美日韩精品成人综合77777| 国产综合精华液| 一级黄片播放器| 国产 一区 欧美 日韩| 免费黄频网站在线观看国产| 搡女人真爽免费视频火全软件| 欧美最新免费一区二区三区| 天堂中文最新版在线下载 | 欧美日本视频| 欧美三级亚洲精品| 1000部很黄的大片| 国产午夜福利久久久久久| 欧美日韩视频精品一区| 国产一区亚洲一区在线观看| 久久午夜福利片| 狂野欧美激情性bbbbbb| 日韩一区二区三区影片| 久久久a久久爽久久v久久| 亚洲av中文字字幕乱码综合| 你懂的网址亚洲精品在线观看| 欧美激情久久久久久爽电影| 亚洲色图av天堂| 亚洲欧美一区二区三区黑人 | av在线app专区| 久久久a久久爽久久v久久| 精品久久久久久久人妻蜜臀av| 尾随美女入室| 久久久久久久久久久免费av| 22中文网久久字幕| 欧美日韩在线观看h| 大码成人一级视频| 国产精品精品国产色婷婷| 丝袜美腿在线中文| 97热精品久久久久久| 日韩一区二区三区影片| 免费黄网站久久成人精品| 国产成人a区在线观看| 在线精品无人区一区二区三 | 亚洲自拍偷在线| 综合色丁香网| 亚洲国产精品国产精品| 午夜福利网站1000一区二区三区| av在线天堂中文字幕| 如何舔出高潮| 精品国产露脸久久av麻豆| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 特级一级黄色大片| 看非洲黑人一级黄片| 欧美成人一区二区免费高清观看| 人妻 亚洲 视频| 久久久午夜欧美精品| 日本av手机在线免费观看| 国产一级毛片在线| 夜夜看夜夜爽夜夜摸| 亚洲内射少妇av| 国产精品国产三级国产专区5o| 免费av观看视频| 亚洲自拍偷在线| 卡戴珊不雅视频在线播放| 精品亚洲乱码少妇综合久久| 欧美精品一区二区大全| 校园人妻丝袜中文字幕| 久久精品国产a三级三级三级| 欧美日韩综合久久久久久| 国产精品国产三级国产av玫瑰| 久久久久精品性色| 岛国毛片在线播放| 777米奇影视久久| 亚洲精品成人av观看孕妇| 搞女人的毛片| 国产亚洲91精品色在线| 天堂俺去俺来也www色官网| 欧美激情久久久久久爽电影| 久久久国产一区二区| 免费黄色在线免费观看| 国产黄片美女视频| 激情五月婷婷亚洲| 亚洲国产精品国产精品| 22中文网久久字幕| 嫩草影院入口| 人妻 亚洲 视频| 91aial.com中文字幕在线观看| 久久97久久精品| 伊人久久国产一区二区| 中文天堂在线官网| 男男h啪啪无遮挡| 最近最新中文字幕大全电影3| 人妻夜夜爽99麻豆av| 舔av片在线| 国产爽快片一区二区三区| 精品国产一区二区三区久久久樱花 | 日韩,欧美,国产一区二区三区| 欧美bdsm另类| 最近的中文字幕免费完整| 久久亚洲国产成人精品v| 国产真实伦视频高清在线观看| 卡戴珊不雅视频在线播放| 亚洲人成网站高清观看| 亚洲欧美精品专区久久| 美女高潮的动态| 黄色一级大片看看| 少妇的逼好多水| 亚洲成人一二三区av| 伦理电影大哥的女人| 777米奇影视久久| 麻豆精品久久久久久蜜桃| 人妻制服诱惑在线中文字幕| 亚洲国产欧美人成| 男人狂女人下面高潮的视频| 看免费成人av毛片| av播播在线观看一区| 七月丁香在线播放| 欧美日韩视频高清一区二区三区二| 91精品伊人久久大香线蕉| 亚洲欧美中文字幕日韩二区| 精品一区二区三卡| 两个人的视频大全免费| 在线观看三级黄色| 在现免费观看毛片| 久热久热在线精品观看| 国产成人午夜福利电影在线观看| 永久网站在线| 精品人妻视频免费看| 成人特级av手机在线观看| 啦啦啦中文免费视频观看日本| 日韩亚洲欧美综合| 一级av片app| 日韩电影二区| 亚洲av二区三区四区| 日本与韩国留学比较| 久久久精品94久久精品| 18禁裸乳无遮挡免费网站照片| 亚洲国产成人一精品久久久| 中国美白少妇内射xxxbb| 永久网站在线| 中国三级夫妇交换| videossex国产| 日本三级黄在线观看| 久久久久久九九精品二区国产| 亚洲,一卡二卡三卡| 夜夜爽夜夜爽视频| 99热国产这里只有精品6| 在现免费观看毛片| av福利片在线观看| 一级毛片电影观看| 99热这里只有是精品50| 26uuu在线亚洲综合色| 亚洲精品日本国产第一区| 中文资源天堂在线| 各种免费的搞黄视频| 女人被狂操c到高潮| 国产午夜福利久久久久久| 综合色丁香网| 看黄色毛片网站| 亚洲欧美成人综合另类久久久| 成年人午夜在线观看视频| 免费大片黄手机在线观看| 亚洲自偷自拍三级| 亚洲成人一二三区av| 国产成人精品久久久久久| 国产精品女同一区二区软件| 大又大粗又爽又黄少妇毛片口| 欧美日韩视频精品一区| 亚洲精品第二区| 成人免费观看视频高清| 国产亚洲午夜精品一区二区久久 | 精品一区二区免费观看| 女人被狂操c到高潮| 亚洲欧美一区二区三区黑人 | 综合色av麻豆| 男人添女人高潮全过程视频| 高清日韩中文字幕在线| 香蕉精品网在线| 午夜福利在线观看免费完整高清在| 欧美一区二区亚洲| 国产午夜精品一二区理论片| 搡女人真爽免费视频火全软件| 国产精品一区www在线观看| 久久久久国产网址| 亚洲最大成人av| 99久久人妻综合| 午夜福利高清视频| 六月丁香七月| 欧美+日韩+精品| 中国三级夫妇交换| 内地一区二区视频在线| 久久99精品国语久久久| 亚洲欧美一区二区三区国产| 国产老妇女一区| 少妇裸体淫交视频免费看高清| 国产高清三级在线| 亚洲伊人久久精品综合| 又爽又黄无遮挡网站| av.在线天堂| 美女xxoo啪啪120秒动态图| 免费黄频网站在线观看国产| 最新中文字幕久久久久| 男人狂女人下面高潮的视频| 亚洲精品,欧美精品| 亚洲va在线va天堂va国产| 波多野结衣巨乳人妻| 丰满人妻一区二区三区视频av| 久久久久性生活片| 边亲边吃奶的免费视频| 春色校园在线视频观看| 国产成人精品福利久久| 18禁在线播放成人免费| 精品久久久久久久人妻蜜臀av| 菩萨蛮人人尽说江南好唐韦庄| 狂野欧美白嫩少妇大欣赏| 国产有黄有色有爽视频| 国产精品久久久久久av不卡| 欧美激情久久久久久爽电影| 色视频在线一区二区三区| 国产精品麻豆人妻色哟哟久久| 亚洲熟女精品中文字幕| 麻豆精品久久久久久蜜桃| 欧美三级亚洲精品| 18禁在线无遮挡免费观看视频| 国产精品99久久久久久久久| 亚洲精品国产色婷婷电影| 观看免费一级毛片| 99热网站在线观看| 国产乱人偷精品视频| 日本色播在线视频| 韩国av在线不卡| av.在线天堂| 久久韩国三级中文字幕| 亚洲欧美成人精品一区二区| 国产精品三级大全| 又黄又爽又刺激的免费视频.| 婷婷色综合大香蕉| 久久精品国产亚洲网站| 2022亚洲国产成人精品| 国产成人aa在线观看| 在线天堂最新版资源| 国产成人免费无遮挡视频| 欧美潮喷喷水| 国产成人免费观看mmmm| 别揉我奶头 嗯啊视频| 一区二区三区精品91| 国产老妇伦熟女老妇高清| 国产白丝娇喘喷水9色精品| 国产精品精品国产色婷婷| 秋霞伦理黄片| 男女边吃奶边做爰视频| 久久久久久久亚洲中文字幕| 菩萨蛮人人尽说江南好唐韦庄| 女人被狂操c到高潮| 神马国产精品三级电影在线观看| 国产中年淑女户外野战色| 日本黄大片高清| 大片免费播放器 马上看| 久久久久久久久久成人| 成人黄色视频免费在线看| 久久久久网色| av.在线天堂| 能在线免费看毛片的网站| 女人久久www免费人成看片| 青春草亚洲视频在线观看| 成人综合一区亚洲| 丰满人妻一区二区三区视频av| 久久精品国产亚洲av天美| 边亲边吃奶的免费视频| 亚洲av一区综合| 最近最新中文字幕大全电影3| 国产 一区精品| 国产精品99久久99久久久不卡 | 中国美白少妇内射xxxbb| 最近最新中文字幕免费大全7| 成年人午夜在线观看视频| 亚洲精华国产精华液的使用体验| 最近最新中文字幕免费大全7| 99热6这里只有精品| 亚洲国产高清在线一区二区三| 精品少妇久久久久久888优播| 日韩一区二区视频免费看| 男女边摸边吃奶| 青青草视频在线视频观看| av黄色大香蕉| 久久99蜜桃精品久久| 午夜福利视频精品| 久久这里有精品视频免费| 亚洲欧美成人综合另类久久久| 久久人人爽人人片av| 在线a可以看的网站| 成人漫画全彩无遮挡| 精品99又大又爽又粗少妇毛片| 成年女人看的毛片在线观看| 91在线精品国自产拍蜜月| 人妻 亚洲 视频| 九九久久精品国产亚洲av麻豆| 久久久久网色| 黄色日韩在线| 国产亚洲午夜精品一区二区久久 | av播播在线观看一区| 美女cb高潮喷水在线观看| 亚洲最大成人中文| 夫妻午夜视频| 亚洲精品影视一区二区三区av| 熟女av电影| 成人国产av品久久久| 日韩,欧美,国产一区二区三区| www.色视频.com| 狂野欧美激情性bbbbbb| www.av在线官网国产| 亚洲精品色激情综合| 狂野欧美激情性bbbbbb| h日本视频在线播放| 在线观看三级黄色| 超碰97精品在线观看| 一级毛片黄色毛片免费观看视频| 久久精品综合一区二区三区| 国产精品久久久久久久久免| 国产成人精品福利久久| 啦啦啦在线观看免费高清www| 99热6这里只有精品| 国产黄片美女视频| 亚洲最大成人av| 国产高清三级在线| 国产免费福利视频在线观看| 国产成人精品久久久久久| 嫩草影院新地址| 国产免费又黄又爽又色| 亚洲欧美成人综合另类久久久| 久久精品久久久久久久性| 免费av不卡在线播放| 一级片'在线观看视频| 麻豆成人午夜福利视频| 欧美老熟妇乱子伦牲交| 我的女老师完整版在线观看| 女的被弄到高潮叫床怎么办| 亚洲精品aⅴ在线观看| 女人十人毛片免费观看3o分钟| 久久99热6这里只有精品| 国产黄片美女视频| 国产亚洲91精品色在线| 国产v大片淫在线免费观看| 国产精品一区二区在线观看99| 成人特级av手机在线观看| 新久久久久国产一级毛片| 免费高清在线观看视频在线观看| 亚洲精品成人av观看孕妇| 自拍偷自拍亚洲精品老妇| 成人二区视频| 亚洲欧美成人综合另类久久久| 亚洲国产精品国产精品| 中文精品一卡2卡3卡4更新| 国产免费又黄又爽又色| 国产高清国产精品国产三级 | 亚洲成人av在线免费| 亚洲国产最新在线播放| 午夜福利网站1000一区二区三区| 在线免费观看不下载黄p国产| 国产黄色免费在线视频| 特大巨黑吊av在线直播| 亚洲怡红院男人天堂| 丝袜喷水一区| 成人特级av手机在线观看| 国产毛片在线视频| 欧美zozozo另类| 成人亚洲欧美一区二区av| 2021少妇久久久久久久久久久| 国内精品美女久久久久久| 中文乱码字字幕精品一区二区三区| 国产乱人视频| 国产精品伦人一区二区| 亚洲伊人久久精品综合| 国产淫语在线视频| 久久人人爽av亚洲精品天堂 | 日韩制服骚丝袜av| 国产成人a区在线观看| 久久人人爽人人爽人人片va| 欧美老熟妇乱子伦牲交| 久久韩国三级中文字幕| 亚洲欧美精品自产自拍| 久久久久性生活片| 免费播放大片免费观看视频在线观看| 真实男女啪啪啪动态图| 男女国产视频网站| 岛国毛片在线播放| 久久国内精品自在自线图片| 中文字幕av成人在线电影| 欧美丝袜亚洲另类| 又大又黄又爽视频免费| 一级爰片在线观看| av一本久久久久| 日韩国内少妇激情av| 久久韩国三级中文字幕| 三级国产精品欧美在线观看| 少妇熟女欧美另类| 日韩不卡一区二区三区视频在线| 国产成人精品一,二区| 2018国产大陆天天弄谢| 中国美白少妇内射xxxbb| 女人十人毛片免费观看3o分钟| 精品国产露脸久久av麻豆| 亚洲第一区二区三区不卡| 欧美 日韩 精品 国产| 国产av码专区亚洲av| 综合色丁香网| 十八禁网站网址无遮挡 | 毛片女人毛片| 成年女人看的毛片在线观看| 亚洲四区av| 欧美区成人在线视频| 日日摸夜夜添夜夜爱| 日本熟妇午夜| 777米奇影视久久| 久久99热6这里只有精品| 国产黄频视频在线观看| 久久久色成人| 中国美白少妇内射xxxbb| 一本色道久久久久久精品综合| 一级av片app| 国内揄拍国产精品人妻在线| 免费在线观看成人毛片| 日韩电影二区| 日韩av免费高清视频| 男女无遮挡免费网站观看| 精品久久国产蜜桃| 日韩欧美一区视频在线观看 | 国产毛片a区久久久久| 免费观看无遮挡的男女| h日本视频在线播放| 国语对白做爰xxxⅹ性视频网站| 国产有黄有色有爽视频| 国产精品熟女久久久久浪| 看十八女毛片水多多多| 亚洲国产最新在线播放| 你懂的网址亚洲精品在线观看| 精品久久久精品久久久| 日本免费在线观看一区| 热99国产精品久久久久久7| 三级国产精品欧美在线观看| 国产黄色视频一区二区在线观看| 久久久久久九九精品二区国产| 丝袜脚勾引网站| 三级国产精品欧美在线观看| 亚洲欧美中文字幕日韩二区| 亚洲av中文av极速乱| 性色av一级| 成人黄色视频免费在线看| 九九爱精品视频在线观看| 日韩制服骚丝袜av| 国产一区有黄有色的免费视频| 国产一区亚洲一区在线观看| 亚洲人成网站高清观看| www.色视频.com| 嘟嘟电影网在线观看| 国产精品av视频在线免费观看| 国产爽快片一区二区三区| 97人妻精品一区二区三区麻豆| 成人二区视频| 2021天堂中文幕一二区在线观| 秋霞在线观看毛片| 男女啪啪激烈高潮av片| 亚洲国产精品成人久久小说| 成年女人看的毛片在线观看| 国产精品熟女久久久久浪| 日本黄大片高清| 国产人妻一区二区三区在| 午夜福利网站1000一区二区三区| 亚洲欧美中文字幕日韩二区| 男女啪啪激烈高潮av片| 久久精品夜色国产| 2022亚洲国产成人精品| 日本黄色片子视频| 亚洲精品日韩av片在线观看| 狠狠精品人妻久久久久久综合| 精华霜和精华液先用哪个| 美女视频免费永久观看网站| 亚洲欧美日韩无卡精品| 女人被狂操c到高潮| eeuss影院久久| 男人添女人高潮全过程视频| 久久97久久精品| 国产高潮美女av| 狠狠精品人妻久久久久久综合| 日韩av免费高清视频| 黄色一级大片看看| 亚洲成人久久爱视频| 国产一区二区在线观看日韩| 秋霞伦理黄片| 亚洲精品影视一区二区三区av| 男女边摸边吃奶| 亚洲国产精品国产精品| 嫩草影院精品99| 国产高清不卡午夜福利| 丝袜脚勾引网站| 视频区图区小说| 欧美变态另类bdsm刘玥| av卡一久久| 神马国产精品三级电影在线观看| 国产日韩欧美亚洲二区| 欧美日韩视频高清一区二区三区二| 美女脱内裤让男人舔精品视频| 国产亚洲一区二区精品| 国产一区亚洲一区在线观看| 看免费成人av毛片| 乱系列少妇在线播放| 26uuu在线亚洲综合色| 久久久精品免费免费高清| 久久韩国三级中文字幕| 国产黄色免费在线视频| videos熟女内射| 视频中文字幕在线观看| 亚洲欧美精品自产自拍| 狠狠精品人妻久久久久久综合| 97在线视频观看| 欧美日韩综合久久久久久| 日本黄色片子视频| 国产成人精品婷婷| 国产黄片美女视频| a级毛色黄片| 99久久中文字幕三级久久日本| 国产精品99久久久久久久久| 国产淫片久久久久久久久| 深夜a级毛片| 一级毛片aaaaaa免费看小| 亚洲高清免费不卡视频| 国产黄频视频在线观看| 亚洲三级黄色毛片| 女人十人毛片免费观看3o分钟| 久久久亚洲精品成人影院| 亚洲图色成人| 嫩草影院精品99| 欧美xxⅹ黑人| 国产爽快片一区二区三区| 日韩强制内射视频| 国产av码专区亚洲av| 久热久热在线精品观看| 久久久久久久午夜电影| 全区人妻精品视频| 18禁动态无遮挡网站| 精品一区在线观看国产| 肉色欧美久久久久久久蜜桃 | 高清日韩中文字幕在线| 不卡视频在线观看欧美| 国产黄色视频一区二区在线观看| 国产成人免费观看mmmm| 欧美+日韩+精品| 人妻少妇偷人精品九色| 国产成人免费观看mmmm| 成人高潮视频无遮挡免费网站| 最近的中文字幕免费完整| 一本色道久久久久久精品综合| 岛国毛片在线播放| 免费在线观看成人毛片| 自拍欧美九色日韩亚洲蝌蚪91 | 中文乱码字字幕精品一区二区三区| 日韩av不卡免费在线播放| 老司机影院成人| 亚洲精品久久久久久婷婷小说| 国产精品一二三区在线看| 中文欧美无线码| 18禁动态无遮挡网站| 欧美精品国产亚洲| 国产成人午夜福利电影在线观看|