基于NDVI的云南省植被覆被變化趨勢分析

劉珊珊　王建雄　牛超杰

摘要：選取云南省為研究區域，利用2001-2015年MODIS中國植被指數合成產品Tiff遙感圖像，通過Arcgis軟件對遙感圖像進行分區統計，得到研究區域各年每月的植被歸一化指數（NDVI）。在此基礎上進行統計分析，結合當地的生長季與種植制度，分析研究區域植被指數的年際、季節及月變化的特點。結果表明，在2001-2015年，云南省植被指數有增加的趨勢，表明在2001-2015年云南省植被覆蓋度基本保持穩定或略有增加的趨勢；2001-2015年云南省春季和冬季的NDVI整體呈增加趨勢，冬季的增速大于春季，夏季NDVI整體呈減少趨勢，而秋季處于平穩狀態，說明云南省2001-2015年春、冬兩季植被覆蓋度不斷增加；云南省NDVI月變化差異較為明顯，與當地生長季和種植制度密切相關，說明植被指數的月變化特征主要受種植作物的影響。

關鍵詞：NDVI；變化特征；趨勢分析；云南省

中圖分類號：P207；Q948 文獻標識碼：A 文章編號：0439-8114（2017）11-2037-04

DOI：10.14088/j.cnki.issn0439-8114.2017.11.009

Abstract： Taking Yunnan as the study area，using the Tiff remote sensing image of MODIS vegetation index composite products from 2001 to 2015，the remote sensing images were divided into districts by the Arcgis software，and the monthly NDVI values of each year in the study area were obtained. Based on the analysis，the characteristics of annual，seasonal and monthly changes of regional vegetation index were analyzed combination with local growing season and cropping system. The results showed that，during 2001 to 2015，there was an increasing trend of the vegetation index in Yunnan，which indicated that the vegetation coverage in Yunnan remained stable or increased slightly from 2001 to 2015. During the period of 2001 to 2015，the NDVI values of spring and winter in Yunnan were increasing in total， and the growth rate of winter was higher than that of spring. The NDVI value was decreasing in summer and the autumn was stable， which meant that the vegetation coverage of spring and winter from 2001 to 2015 was increasing. The monthly variation of the vegetation index was relatively obvious， which was closely related to the local growing season and cropping system， which indicated that the monthly change of vegetation index was mainly affected by the crop.

Key words： NDVI； variation characteristics； trend analysis； Yunnan Province

植被是聯結土壤、大氣和水分等要素的自然“紐帶”，具有明顯的年際變化和季節變化的特點，在一定程度上能代表土地覆蓋的變化[1，2]。植被指數是將遙感地物光譜資料經數學方法處理，以反映植被狀況的特征量[3]。在眾多植被指數中最為常用的是歸一化植被指數（NDVI），它能夠精確地反映植被綠度、光合作用強度、植被代謝強度及其季節和年際變化，是表征植物生長、植被覆蓋、生長狀況、生物量等的重要指標[4]，是監測植被變化的有效參數[5，6]，因此，在各地大尺度的植被動態監測、作物長勢監測、自然災害監測和作物產量預測等方面得到廣泛應用[7-10]。近年來，中國學者已經對很多地區的植被指數時空變化特征進行了深入的研究，但是不同的地域、不一樣的研究尺度、不同時間尺度、不同的數據源，植被指數的時空變化特征都會有所差異[1]。研究表明，中國西北地區20世紀80年代初期至1994年前后植被呈增加趨勢，而1994-2000年植被有退化趨勢[11]；而長江三角洲、珠江三角洲地區的NDVI有減少的趨勢[12]；秦嶺地區植被覆蓋度在1982-2007年呈降低趨勢[13]；貴州省植被覆蓋在1982-2003年有增加趨勢[14]。

云南省山地多、壩區極少的現狀使當地耕地少，林地多，復雜多樣的氣候條件以及地形地貌，使得云南省植被類型多樣。利用MODIS/NDVI數據集，使用Arcgis軟件對影像數據進行分區統計處理，通過一元線性回歸趨勢定量分析2001-2015年云南省植被年際以及月變化情況，進一步探討研究區域植被指數與當地種植制度及植被生長季的內在關聯。為進一步深入研究云南省生態環境現狀提供基礎參考資料。

1 研究區域概況與數據處理

1.1 研究區域概況

云南省地處中國西南部，北依廣袤的亞洲大陸，南鄰廣闊的印度洋及太平洋，處于東南季風和西南季風的控制下，又受青藏高原區的影響，從而形成復雜多樣的自然地理環境。云南省氣候兼具山原氣候、低緯氣候以及季風氣候的特點。夏季平均溫度在19～22 ℃，冬季平均溫度在6～8 ℃。年溫差一般為10～15 ℃，一天的溫度變化是早涼，午熱，尤其是春、冬兩季，日溫差高達12～20 ℃。全省大部分地區年降水量在1 100 mm左右，南部部分地區可達1 600 mm以上。由于冬夏兩季受不同大氣環流的控制和影響，降水量在時域上和地域上的分配也是極不均勻的。6-8月降水量最多，約占全年降水量的60%，11月至次年4月的冬春季節為旱季，降水量僅占全年的10%～20%。云南省相對平緩的山區只占總面積的10%，大面積土地高低差參，縱橫起伏，一定范圍又有和緩的高原面。云南全省面積在1 km2以上的壩子共有1 445個，面積在100 km2以上的壩子有49個，云南最大的壩子在陸良縣，面積為771.99 km2[15]。

由于復雜多樣的氣候條件、地形地貌以及降雨量在時間上分布極不均勻等特點，云南省的植被類型、種植制度和種植作物極具典型特征。云南省分為大春和小春，主要種植作物有小麥、玉米、水稻、烤煙等，植物類型繁多，主要有銀杏、龍柏、鵝掌楸等。氣候的多變使得云南省的植被覆蓋度在月、季節上不斷發生變化，加上干旱災害時有發生，年際上的植被覆蓋度又有所不同。

1.2 數據與處理

1.2.1 數據來源 NDVI數據集是由地理空間數據云提供的MODIS中國植被指數合成產品，它是由MODND1D計算得到，計算方法為取月內每天最大值。坐標系為EPSG：4326（WGS84），空間分辨率為500 m，時間分辨率為月。這里使用2001-2015年的NDVI數據集。

NDVI=（NIR-R）/（NIR+R） （1）

式中，NIR為近紅外波段的反射值；R為紅光波段的反射值。

1.2.2 處理方法

1）基于GIS的Tiff圖像處理。由于最初獲取的是全國范圍的NDVI Tiff圖像，因此需要進一步利用Arcgis軟件添加云南省省界對全國的NDVI進行分區統計，從而獲取云南省每年每月的NDVI數值，從中整合研究區域研究時段內的NDVI數據，以便后期進行統計分析。

NDVI的數值區間為（-1，1），負值表示地面覆蓋為云、水、雪等，對可見光高反射；0表示有巖石或裸土等，NIR與R近似相等；正值表示有植被覆蓋，且NDVI隨覆蓋度增大而增大。

2）NDVI數據處理。

2.2 季節變化分析

根據公式（3）計算得到2001-2015年4個季節的NDVI，結果如表2所示。由表2可知，秋季和冬季的NDVI普遍高于春季和夏季，春夏季節的NDVI波動范圍為0.598 362～0.742 288，秋冬季節的NDVI波動范圍為0.644 697～0.817 897，說明在秋季和冬季的植被覆蓋度高于其他兩個季節。2001-2015年云南省春季的NDVI整體上處于不斷增加的趨勢，趨勢線斜率為0.005 4；夏季NDVI整體上呈減少趨勢，趨勢線斜率為-0.003 8；秋季的NDVI處于平穩狀態，冬季NDVI整體上呈增加趨勢，趨勢線的斜率為0.008 7，冬季的增速大于春季，說明云南省2001-2015年春、冬兩季植被覆蓋度不斷增加。

2.3 月變化分析

根據公式（3），計算得到研究時間段各月的NDVI均值，結果如表3所示。由表3可知，一年當中，云南省植被指數的最高值出現在10月，為0.782 317，最小值出現在4月，為0.635 850。從植被生長季的角度進行分析，4-8月NDVI明顯增加，這與云南省植被生長期和枯萎期相對應，4月進入大春，總體從5月上旬生長開始，8月達到峰值；9月植被指數略有降低，這是由于云南省9月大春作物成熟且被收割；但隨著10月云南省進入小春，9-10月的植被指數再次呈現增加趨勢；10-12月以及1-4月之間的時間段，植被指數呈現顯著下滑趨勢。整體來講，NDVI的變化符合云南省的種植制度。

3 結論與展望

3.1 結論

通過使用基于時間序列的歸一化植被指數影像數據，利用一元線性趨勢分析法對2001-2015年云南省植被指數的年際、季節及月變化特征進行定量分析，得到以下3點結論。

1）在2001-2015年，云南省植被指數有增加的趨勢，表明在2001-2015年云南省植被覆蓋度基本保持穩定或略有增加的趨勢。

2）2001-2015年云南省春季和冬季的NDVI整體呈增加趨勢，冬季的增速大于春季，夏季NDVI呈減少趨勢，而秋季的NDVI處于平穩狀態，說明云南省2001-2015年春、冬兩季植被覆蓋度不斷增加。

3）云南省植被指數月變化差異較為明顯，與當地生長季和種植制度密切相關。NDVI處于增長的時間段與云南省作物種植期大春、小春相對應，說明NDVI的月變化特征主要受種植作物的影響。

3.2 展望

隨著傳感器的不斷進步，每天能夠獲得的遙感影像也不斷增多[16]，大家的焦點問題聚集在如何能夠及時地將這些海量的影像數據轉化為有用信息為我們所用。對大量的影像數據進行時間序列建模分析是有效的解決途徑之一。

運用時間序列對植被覆被變化進行分析能夠直觀準確地反映研究區域的生態環境，但能夠反映生態環境的要素復雜多樣，植被僅僅是其中主要的一個要素。在把生態環境放在首位的大背景下，為了保護我們的生存環境，不僅需要對植被覆被變化進行分析，還要從影響其變化的內在因素進行研究，包括自然因素（土壤、氣候等）和人為因素之間的相互影響。從更深層次來看，協調人地關系是對保護生態環境更高層次的要求，同植被要素一樣，其他自然要素和社會要素時序數據普遍存在，對這些多類型時間序列數據綜合分析，發現長時期變化的相互影響機理，對于增加植被覆被面積、改善生態環境無疑是一件具有挑戰性和研究價值的工作。

參考文獻

[1] 何彬方，馮 妍，吳文玉，等.安徽省近十年植被指數時空變化特征[J].生態學雜志，2010，29（10）：1912-1918.

[2] 吳月圓，徐天蜀，岳彩榮.基于MODIS/NDVI的云南省近十年植被動態監測分析[J].綠色科技，2013，4（10）：134-135.

[3] 解文歡，張有智，吳 黎.基于植被指數對望奎縣糧食作物產量預測方法的研究[J].黑龍江農業科學，2011，34（4）：113-115.

[4] TUCKER C J，FUNG I Y，KEELING C D，et al. Relationship between atmospheric CO2 variations and a satellite-derived vegetation index[J]. Nature，1986，319（6050）：195-199.

[5] 畢曉麗，王 輝，葛劍平.植被歸一化指數（NDVI）及氣候因子相關起伏型時間序列變化分析[J].應用生態學報，2005，16（2）：284-288.

[6] 張戈麗，歐陽華，張憲洲，等.基于生態地理分區的青藏高原植被覆被變化及其對氣候變化的響應[J].地理研究，2010，29（11）：2004-2016.

[7] 李月臣，宮 鵬，劉春霞，等.北方13省1982年～1999年植被變化及其與氣候因子的關系[J].資源科學，2006，28（2）：109-117.

[8] 王金亮，高 雁.云南省近20年植被動態變化遙感時序分析[J].云南地理環境研究，2010，22（6）：1-7.

[9] 鄭玉坤.多時相AVHRR-NDVI數據的時間序列分析及其在土地覆蓋分類中的應用[D].北京：中國科學院，2002.

[10] 張 翔，王 勇.NDVI對氣候因子的響應研究[J].地理空間信息，2015，12（6）：39-41，50.

[11] 徐興奎，林朝暉，薛 峰，等.氣象因子與地表植被生長相關性分析[J].生態學報，2003，23（2）：221-230.

[12] 方精云，樸世龍，賀金生，等.近20年來中國植被活動在增強[J].中國科學（C輯：生命科學），2003，33（6）：554-565.

[13] 孫 華，白紅英，張清雨，等.秦嶺南北地區植被覆蓋對區域環境變化的響應[J].環境科學學報，2009，29（12）：2635-2641.

[14] 鄭有飛，牛魯燕，吳榮軍，等.1982-2003年貴州省植被覆蓋變化及其對氣候變化的響應[J].生態學雜志，2009，28（9）：1773-1778.

[15] 廖麗君，魏 巍，鄭宏剛，等.縣級農用地質量分等成果年度更新方法研究[J].中國農業資源與區劃，2015，36（1）：51-62.

[16] 于 嶸.基于遙感時序數據的中國陸地植被覆蓋變化分析研究[D].北京：中國科學院，2006.