龔文峰 曹凱華 周瑩

摘要：基于2006、2011和2016年的Landsat TM/OLI影像數(shù)據(jù)，反演農(nóng)牧過渡區(qū)植被覆蓋度。借助于RS和GIS技術(shù)，以黑龍江省甘南縣為研究對(duì)象，采用歸一化植被指數(shù)及像元二分模型，定量獲取研究區(qū)域植被覆蓋度，基于轉(zhuǎn)移矩陣模型定量分析10年間植被覆蓋度時(shí)空變化特征和演變進(jìn)程。結(jié)果表明，研究區(qū)植被覆蓋度整體狀況較好，2006、2011、2016年3期數(shù)據(jù)高覆蓋度和中高覆蓋度面積比例之和分別為87.29%、87.18%和87.10%。2006—2011年植被覆蓋度主要由高覆蓋度和中覆蓋度轉(zhuǎn)換為中高覆蓋度，轉(zhuǎn)換率分別為76.83%和55.49%，2011—2016年中高覆蓋度和高覆蓋度之間轉(zhuǎn)換強(qiáng)烈，相互間的轉(zhuǎn)換率分別為10.87%和83.00%。10年間，2006—2011、2006—2016年植被變化穩(wěn)定區(qū)的面積比例分別為64.32%和65.54%。2011—2016年較2006—2011年植被退化面積減少257 km2，植被改善面積增加199 km2，植被仍以退化為主。

關(guān)鍵詞：植被覆蓋度;農(nóng)牧交錯(cuò)帶;時(shí)空演變;甘南縣

中圖分類號(hào)：TP79? ? ? ? ?文獻(xiàn)標(biāo)識(shí)碼：A

文章編號(hào)：0439-8114（2019）12-0052-06

DOI：10.14088/j.cnki.issn0439-8114.2019.12.013? ? ? ? ? ?開放科學(xué)（資源服務(wù)）標(biāo)識(shí)碼（OSID）：

Abstract： In order to invert the vegetation coverage of the farming-pastoral transition zone by using the Landsat TM/OLI image data of 2006， 2011 and 2016 ，the Gannan county of Heilongjiang province was taken as the study region， based on the RS and GIS technology， the models of normalized difference vegetation index（NDVI） and dimidiate pixel were quantitatively obtained the vegetation coverage of the study area， and the transfer matrix model was used to analyze the temporal and spatial variation characteristics and evolution process of vegetation coverage within 10 years. The results showed that the overall coverage of vegetation cover in the study area was good. The sum of the data coverage and the medium-high coverage area of the three phases in 2006，2011，2016 was 87.29%，87.18% and 87.10%， respectively. In 2006—2011， the vegetation coverage was mainly converted from high coverage and middle coverage to medium and high coverage， and the conversion rates was 76.83% and 55.49%， respectively. In 2011—2016， the conversion between medium and high coverage and high coverage was strong， and the conversion rate was 10.87% and 83.00%， respectively. In the past 10 years， the area ratio of the vegetation change stable area was 64.32% and 65.54%， respectively. 2011—2016 compared with 2006—2011， the area of vegetation degradation decreased by 257 km2， while the area of vegetation improvement increased by 199 km2， and the vegetation was still dominated by degradation.

Key words： vegetation coverage; farming-pastoral transition zone; spatio-temporal evolution; Gannan county

收稿日期：2018-10-12

基金項(xiàng)目：國(guó)家科工局高分辨率對(duì)地重大觀測(cè)項(xiàng)目（民用部分：72-Y20A08-9002-16/17）;黑龍江省自然科學(xué)基金項(xiàng)目（QC2017036）

作者簡(jiǎn)介：龔文峰（1976-），男，河南南陽人，教授，博士，主要從事RS和GIS在資源與環(huán)境監(jiān)測(cè)上的應(yīng)用，（電話）13796674248（電子信箱）

gwf101@163.com。

植被覆蓋度是指植被（包括葉、莖、枝）在地面的垂直投影面積占統(tǒng)計(jì)區(qū)總面積的百分比[1，2]，它反映了植被在水平方向上的密度情況，是描述生態(tài)系統(tǒng)的重要基礎(chǔ)數(shù)據(jù)[3]，也是控制土壤侵蝕的關(guān)鍵因素和指示生態(tài)系統(tǒng)質(zhì)量變化的一個(gè)重要指標(biāo)[4，5]。樣本估算法是提取植被信息常用的傳統(tǒng)方法[6]，該方法不僅成本高、誤差大，且不適用于大面積監(jiān)測(cè)植被生長(zhǎng)狀況。然而遙感技術(shù)的不斷發(fā)展為大區(qū)域植被覆蓋甚至全球的植被覆蓋監(jiān)測(cè)提供了可能?；谶b感數(shù)據(jù)計(jì)算出的植被指數(shù)可以直接反映地表植被狀況，且植被指數(shù)與植被覆蓋度具有密切的相關(guān)關(guān)系[7]，已被廣泛應(yīng)用于植被覆蓋度的定性、定量估算，逐漸成為相關(guān)領(lǐng)域研究的一種新趨勢(shì)[8]。目前，基于植被指數(shù)估算植被覆蓋度模型的研究主要集中在植被覆蓋度時(shí)空演變及其發(fā)展趨勢(shì)[9，10]、基于二分模型探究植被覆蓋度[11，12]、城市化進(jìn)程與植被覆蓋度的定量分析[13，14]、植被覆蓋度變化的影響因素分析[15]、植被覆蓋度與土地利用的分析[16]等，但這些研究多集中在流域、湖泊及城市區(qū)域，而在東北農(nóng)牧過渡地區(qū)的研究相對(duì)較少。

農(nóng)牧交錯(cuò)帶是中國(guó)北方農(nóng)業(yè)區(qū)與草原區(qū)之間存在的一種寬窄不一的過渡地帶，其生態(tài)系統(tǒng)結(jié)構(gòu)較為簡(jiǎn)單，生態(tài)環(huán)境十分脆弱且對(duì)全球變化高度敏感[17]，區(qū)域生態(tài)環(huán)境保護(hù)和可持續(xù)發(fā)展方式一直是該區(qū)域的研究主題?；诖?，以位于農(nóng)牧交錯(cuò)區(qū)域的甘南縣為研究對(duì)象，以2006—2016年Landsat TM和OLI為主要數(shù)據(jù)源，借助于RS和GIS技術(shù)，基于NDVI及像元二分模型定量反演研究區(qū)域的植被覆蓋度，借助轉(zhuǎn)移矩陣模型，研究分析甘南縣10年間植被覆蓋度時(shí)空變化特征，探求研究區(qū)植被覆蓋情況及其變化規(guī)律，以期為該區(qū)域草地保護(hù)、土地利用優(yōu)化和生態(tài)環(huán)境保護(hù)提供理論依據(jù)，對(duì)該區(qū)域的農(nóng)業(yè)發(fā)展、生態(tài)建設(shè)和可持續(xù)發(fā)展具有一定借鑒意義。

1? 研究區(qū)概況

甘南縣位于黑龍江省西部，嫩江中游右岸，地理坐標(biāo)為122°54′—124°28′E，47°35′—48°32′N，地處大興安嶺南麓與松嫩平原過渡地帶，西部、北部是丘陵漫崗區(qū)，南部和東南部是地勢(shì)平坦的開闊平原。地形西北高，東南低，海拔為160～380 m，屬于寒溫帶大陸性季風(fēng)氣候，四季冷暖干濕分明：春季風(fēng)大雨少，易于發(fā)生干旱;夏季高溫多雨;秋冬季急劇降溫，大地封凍。年均氣溫2.6 ℃，年均降水量為455.2 mm，無霜期在120～132 d。甘南縣屬半農(nóng)半牧區(qū)域，是國(guó)家扶貧開發(fā)工作重點(diǎn)縣之一。土壤類型以黑土和黑鈣土為主，適合多種作物生長(zhǎng)。甘南縣自然資源豐富，林木主要有松樹、楊樹、榆樹、果樹等。

2? 數(shù)據(jù)與方法

2.1? 數(shù)據(jù)來源

選取2006、2011和2016年Landsat 5 TM和Landsat 8 OLI的遙感數(shù)據(jù)為主要數(shù)據(jù)源（http：//www.gscloud.cn/）。所選取3期影像的時(shí)間均為8月，正是植被長(zhǎng)勢(shì)旺盛的時(shí)期，該時(shí)期遙感影像既能較好地反映植被生長(zhǎng)狀況，同時(shí)云量又相對(duì)較少。此外，收集到研究區(qū)域1∶50 000的行政區(qū)劃圖、DEM、部分野外GPS調(diào)查數(shù)據(jù)和其他相關(guān)統(tǒng)計(jì)數(shù)據(jù)。

2.2? 數(shù)據(jù)處理

2.2.1? 數(shù)據(jù)預(yù)處理? 基于遙感軟件ENVI 5.3對(duì)3期數(shù)據(jù)進(jìn)行輻射定標(biāo)、幾何校正、快速大氣校正、鑲嵌、拼接等預(yù)處理，圖像的投影方式為橫軸墨卡托投影UTM（Zone51N），然后利用行政邊界完成遙感影像的裁剪。

2.2.2? 植被指數(shù)獲取? 歸一化植被指數(shù)（Normalized difference vegetation index，NDVI）是應(yīng)用最廣泛的一種植被指數(shù)[18]，NDVI的變化在一定程度上能代表地表覆蓋度變化[19]，計(jì)算公式如下：

式中，Bnir代表近紅外波段的反射值;Bred代表紅光波段的反射值;NDVI是植物生長(zhǎng)狀態(tài)及植被空間分布密度的最佳指示因子。

2.2.3? 植被覆蓋度獲取? 像元二分模型是目前植被覆蓋度反演的有效方法[20]，其原理為通過遙感傳感器所觀測(cè)到的信息S，可表達(dá)為由綠色植被成分所貢獻(xiàn)的信息Sveg和由土壤成分所貢獻(xiàn)的信息Ssoil這兩部分組成[21]，植被覆蓋度Fc就是像元中植被覆蓋所占的面積，計(jì)算公式如下：

將歸一化植被指數(shù)NDVI代入相元二分模型，得到植被覆蓋度轉(zhuǎn)化公式：

式中，NDVImin和NDVImax分別代表無植被覆蓋和全植被覆蓋的NDVI。結(jié)合本研究狀況，NDVImin和NDVImax分別取NDVI圖像直方圖5%處和95%處NDVI，并利用式（3）對(duì)植被覆蓋度進(jìn)行求解。

2.2.4? 植被覆蓋度分類? 按照水利部1996年頒布的《土壤侵蝕分類分級(jí)標(biāo)準(zhǔn)》及相關(guān)研究成果[22，23]，結(jié)合研究區(qū)的實(shí)際情況及一類森林清查資料、野外調(diào)查點(diǎn)數(shù)據(jù)及高分辨率Google Earth影像，對(duì)分類的結(jié)果進(jìn)行了精度驗(yàn)證[24]，在此基礎(chǔ)上，基于ArcGIS 10.3將植被覆蓋度進(jìn)行等級(jí)劃分，具體結(jié)果見表1。

2.2.5? 轉(zhuǎn)移矩陣模型? 轉(zhuǎn)移矩陣模型法來源于系統(tǒng)分析中對(duì)系統(tǒng)狀態(tài)與形態(tài)轉(zhuǎn)移的定量描述[25]。此模型能綜合考慮第i類植被覆蓋類型轉(zhuǎn)變?yōu)槠渌莍類植被覆蓋類型與其他非i類植被覆蓋類型在該研究時(shí)期內(nèi)由其他空間區(qū)位上同時(shí)轉(zhuǎn)變?yōu)榈趇類植被覆蓋類型的變化過程，可以更精細(xì)和準(zhǔn)確地測(cè)算植被覆蓋度變化的空間過程。其表達(dá)式如下：

式中，Sij表示植被覆蓋度由i類轉(zhuǎn)為j類的面積;n表示植被覆蓋度類型數(shù)目;i、j分別表示研究初期與末期的植被覆蓋度類型。

3? 結(jié)果與分析

3.1? 植被覆蓋度時(shí)間變化動(dòng)態(tài)分析

基于RS和GIS獲取研究區(qū)域的植被覆蓋度和面積統(tǒng)計(jì)數(shù)據(jù)（圖1和表2）。由表2分析可知，2006、2011、2016年高覆蓋度和中高覆蓋度面積之和最大，分別為4 174.22、4 168.87和4 165.15 km2。其中，中高覆蓋度面積占比最大，分別為70.84%、77.19%、73.16%，呈現(xiàn)先增加后減少的趨勢(shì)，2011年面積達(dá)到最大，為3 691.31 km2。中高覆蓋度所對(duì)應(yīng)的土地利用類型為部分耕地、稀疏林地和高產(chǎn)草地，且該區(qū)域土地類型主要以耕地和草地為主，符合研究區(qū)的土地利用結(jié)構(gòu)特征。高覆蓋度植被的面積變化趨勢(shì)和中覆蓋度植被相反，呈現(xiàn)先減少后增加的趨勢(shì)，減少和增加的面積分別為309.08和189.06 km2。中覆蓋度和低覆蓋度面積變化較小，中覆蓋度植被面積呈波動(dòng)式增長(zhǎng)趨勢(shì)，面積先減少后增加，總體增加了58.56 km2，而低覆蓋度植被面積呈波動(dòng)式下降趨勢(shì)，面積先增加后減少，總體減少了11.84 km2。極低覆蓋度植被面積逐年減少，由2006年的183.56 km2減至2016年的145.90 km2。

總體來看，研究區(qū)植被覆蓋度整體狀況較好。10年間中高覆蓋和高覆蓋面積變化較明顯，中高覆蓋度面積增加了2.32%。高覆蓋度面積減少了2.51%。中覆蓋度、低覆蓋度和極低覆蓋度植被面積的變化量較小，基本處于小波動(dòng)變化狀態(tài)。其中，低覆蓋度植被面積變化最小，由2006年的3.00%變?yōu)?016年的2.75%，變化量?jī)H減少了0.25%。

3.2? 植被覆蓋度空間動(dòng)態(tài)變化分析

為了更真實(shí)地反映植被覆蓋度的空間變化特征，本研究選用轉(zhuǎn)移矩陣模型定量分析各植被類型空間動(dòng)態(tài)變化特征?；贏rcGIS的空間分析，將2006、2011、2016年3個(gè)數(shù)據(jù)疊加分析，并將10年間植被變化特征分為兩個(gè)階段，即2006—2011年為第一階段，2011—2016年為第二階段，在此基礎(chǔ)上獲取不同植被覆蓋度等級(jí)面積轉(zhuǎn)移數(shù)據(jù)（圖2）。

由表3可知，2006—2011年，高覆蓋度植被面積轉(zhuǎn)換最大，為683.43 km2，且有97.53%轉(zhuǎn)化為中高覆蓋度植被，主要發(fā)生在中興鄉(xiāng)、阿倫河流域及東陽鎮(zhèn)東部，主要由密林地和優(yōu)質(zhì)耕地變?yōu)橐话愀?，該轉(zhuǎn)變?cè)从谠诮?jīng)濟(jì)利益的驅(qū)使下，當(dāng)?shù)鼐用駷榱俗非蟾蟮慕?jīng)濟(jì)效益，加大了對(duì)林地的破壞力度，區(qū)域生態(tài)環(huán)境趨于惡化，土地利用呈現(xiàn)出一定的不合理化;其次是中高覆蓋度植被，轉(zhuǎn)換面積為610.7 km2，其中366.93 km2轉(zhuǎn)換為高覆蓋度植被，主要發(fā)生在寶山鄉(xiāng)東部、查哈陽鄉(xiāng)中南部和甘南鎮(zhèn)南部，源于國(guó)家推行退耕還林（草）等生態(tài)措施，部分耕地變?yōu)槊芰值鼗蛘邇?yōu)質(zhì)耕地，區(qū)域生態(tài)環(huán)境趨于好轉(zhuǎn)。在一定程度上表明中高覆蓋度和高覆蓋度植被間的轉(zhuǎn)換是該時(shí)期的主要轉(zhuǎn)換類型。中覆蓋度植被轉(zhuǎn)為中高覆蓋度植被的面積為169.14 km2，發(fā)生地點(diǎn)在查哈陽鄉(xiāng)北部、興十四鎮(zhèn)東部和巨寶鎮(zhèn)西部區(qū)域，對(duì)應(yīng)的土地利用類型為草地轉(zhuǎn)為耕地，該轉(zhuǎn)變?cè)从谌祟悓?duì)草地的不合理開墾，且在開發(fā)過程中沒有顧及生態(tài)保護(hù)問題，加之受氣候條件的影響，局部出現(xiàn)了草地沙漠化的現(xiàn)象，區(qū)域生態(tài)環(huán)境不斷惡化。低覆蓋度植被轉(zhuǎn)變?yōu)橹懈采w度和中高覆蓋度的面積分別為34.33和53.52 km2，發(fā)生地點(diǎn)在諾敏河周邊、查哈陽鄉(xiāng)北部和巨寶鎮(zhèn)的西部，土地利用類型由荒地（灘涂地）轉(zhuǎn)為耕地，主要源于周邊農(nóng)民干擾，填湖開荒，使?jié)竦孛娣e減少，區(qū)域生態(tài)環(huán)境進(jìn)一步惡化。極低覆蓋度植被面積轉(zhuǎn)變最小，共有84.17 km2發(fā)生轉(zhuǎn)變，主要轉(zhuǎn)變?yōu)榈透采w度和中高覆蓋度，轉(zhuǎn)換面積分別為40.36和25.22 km2，發(fā)生地點(diǎn)在太平湖水庫(kù)、音河水庫(kù)、音河流域、諾敏河流域，水域轉(zhuǎn)為灘涂地（水田）是其主要轉(zhuǎn)換地類，除了受氣候影響河流水位下降外，在經(jīng)濟(jì)利益驅(qū)使下，周邊農(nóng)民在河流周圍圍湖造田，不合理的土地利用使局部生態(tài)環(huán)境不斷惡化，今后應(yīng)加大對(duì)水域濕地的保護(hù)，促使區(qū)域土地合理使用。該時(shí)期極低覆蓋度是低覆蓋度和中高覆蓋度的主要轉(zhuǎn)入來源，源于人類對(duì)耕地及優(yōu)質(zhì)耕地需求增多及受城鎮(zhèn)化進(jìn)程的加速影響。

由表4可知，2011—2016年，中高覆蓋度植被的轉(zhuǎn)出面積最大，為827.82 km2，其中553.28 km2轉(zhuǎn)化為高覆蓋度，主要零星分布在中興鄉(xiāng)、巨寶鎮(zhèn)和甘南鎮(zhèn)，對(duì)應(yīng)的地類變化是耕地轉(zhuǎn)為密林地，這主要得益于國(guó)家退耕還林、還草工程等生態(tài)措施推行，促使局部生態(tài)環(huán)境質(zhì)量好轉(zhuǎn)，在一定程度上表明積極生態(tài)政策已經(jīng)取得成效，區(qū)域生態(tài)環(huán)境惡化趨勢(shì)得到緩解。高覆蓋度植被和中覆蓋度植被轉(zhuǎn)化為中高覆蓋度植被面積分別為380.17和148.46 km2，是中高覆蓋度植被的主要轉(zhuǎn)入來源，發(fā)生地點(diǎn)主要在甘南鎮(zhèn)南部、雙河農(nóng)場(chǎng)、查哈陽農(nóng)場(chǎng)及查哈陽鄉(xiāng)西北部、甘南鎮(zhèn)局部，對(duì)應(yīng)的轉(zhuǎn)換地類是林地（草地）轉(zhuǎn)為耕地，源于當(dāng)?shù)剞r(nóng)民為了追求更多的經(jīng)濟(jì)利益，加大了對(duì)林地（草地）的干擾和破壞力度，出現(xiàn)了毀草開荒等現(xiàn)象，使得局部生態(tài)環(huán)境遭到一定程度破壞，人為因素在環(huán)境惡化過程中處于主導(dǎo)作用。對(duì)于低覆蓋度植被而言，其轉(zhuǎn)換面積相對(duì)較小，為135.9 km2，地點(diǎn)集中在太平鎮(zhèn)東部和阿倫河下游，土地轉(zhuǎn)換類型是灘涂地轉(zhuǎn)為耕地，在經(jīng)濟(jì)利益的驅(qū)使下，當(dāng)?shù)剞r(nóng)民加大了對(duì)灘涂地的改造利用，在造成部分濕地消失的同時(shí)，進(jìn)一步惡化了局部生態(tài)環(huán)境。極低覆蓋度植被轉(zhuǎn)換面積變化最小，僅為85.38 km2，其中44.58 km2轉(zhuǎn)化為低覆蓋度，主要發(fā)生在諾敏河流域、音河水庫(kù)西部及太平湖水庫(kù)北部，水域轉(zhuǎn)為耕地是該區(qū)域主要土地利用轉(zhuǎn)變方式，主要源于人為因素使水域面積減少并加大了對(duì)濕地的改造利用?？傮w來看，2011—2016年，中高覆蓋保留率最高，達(dá)到81.85%，低覆蓋度和極低覆蓋度植被轉(zhuǎn)入率略小于轉(zhuǎn)出率，變化不是很明顯。受國(guó)家林業(yè)生態(tài)等政策的影響，呈現(xiàn)出了中高覆蓋度向高覆蓋度轉(zhuǎn)變的局面。

3.3? 植被覆蓋度時(shí)空演變分析

基于ArcGIS 10.3的空間分析，根據(jù)式（5），完成3期數(shù)據(jù)的疊加分析，在獲取2006—2016年疊加數(shù)據(jù)的基礎(chǔ)上，并對(duì)植被覆蓋度等級(jí)變化進(jìn)行定量分析，以期有助于深究該區(qū)域植被覆蓋度的時(shí)空演變特征。

式中，T是疊加結(jié)果;T1、T2表示要參與疊加的數(shù)據(jù);n表示參與運(yùn)算的數(shù)據(jù)期數(shù)。對(duì)于運(yùn)算結(jié)果而言，如果前后植被覆蓋度等級(jí)未發(fā)生變化的則定義為穩(wěn)定區(qū)，植被覆蓋度等級(jí)由低等轉(zhuǎn)為高等則為改善區(qū)，覆蓋度等級(jí)由高等轉(zhuǎn)為低等則視為退化區(qū)，并依次計(jì)算各個(gè)植被變化類型面積，探究植被覆蓋度變化總體特征（圖3）。

由圖3可知，2006—2011年穩(wěn)定區(qū)占據(jù)一定的優(yōu)勢(shì)地位，面積占比為64.32%，主要分布在中興鄉(xiāng)西南部區(qū)域、巨寶鎮(zhèn)、興十四鎮(zhèn)大部分區(qū)域、平陽鎮(zhèn)東部區(qū)域、寶山鄉(xiāng)及查哈陽鄉(xiāng)北部區(qū)域。2011—2016年穩(wěn)定區(qū)占研究區(qū)面積的65.54%，分布區(qū)域沒有發(fā)生太明顯的改變，增加的穩(wěn)定區(qū)面積在長(zhǎng)山鄉(xiāng)、查哈陽鄉(xiāng)西北部，源于這些區(qū)域主要分布著大量的耕地，土地利用結(jié)構(gòu)沒有發(fā)生改變。2006—2011年改善區(qū)域主要零星分布在中興鄉(xiāng)、巨寶鎮(zhèn)西北部、東陽鎮(zhèn)南部、查哈陽鄉(xiāng)北部區(qū)域，其區(qū)域面積約為714 km2，2011—2016年改善區(qū)面積為913 km2，面積增加了199 km2，增加的區(qū)域主要集中在興十四鎮(zhèn)南部、平陽鎮(zhèn)北部、太平湖流域及阿倫河中下游區(qū)域，受地形和林業(yè)生態(tài)政策因素的影響，植被覆蓋度增加，區(qū)域生態(tài)環(huán)境有所好轉(zhuǎn)。2006—2011年退化區(qū)的面積為991 km2，而2011—2016年退化區(qū)的面積為734 km2，面積減少了257 km2。退化區(qū)主要分布在甘南鎮(zhèn)、雙河農(nóng)場(chǎng)和查哈陽農(nóng)場(chǎng)及部分河流附近，主要由于城鎮(zhèn)化進(jìn)程的加速，使得部分耕地（草地）轉(zhuǎn)化為建筑用地，加之受經(jīng)濟(jì)利益驅(qū)使，當(dāng)?shù)剞r(nóng)民加大了對(duì)草地和林地的干擾力度，使得草地面積退化，部分區(qū)域出現(xiàn)沙化現(xiàn)象，區(qū)域生態(tài)環(huán)境遭到一定程度的破壞。

上述分析表明，2006—2016年，植被覆蓋度總體穩(wěn)定。10年間植被覆蓋度退化面積為1 725 km2，改善面積為1 627 km2，植被仍以退化為主，研究區(qū)生態(tài)環(huán)境已被破壞，今后需要加強(qiáng)退耕還林還草及封山育林等措施的實(shí)施，保護(hù)甘南縣生態(tài)環(huán)境健康發(fā)展。

4? 小結(jié)

基于RS和GIS支持，以2006、2011和2016年的Landsat TM和Landsat 8 OLI影像為主要數(shù)據(jù)源，以甘南縣為研究區(qū)，運(yùn)用歸一化植被指數(shù)及相元二分模型，獲取植被覆蓋度數(shù)據(jù)，并定量分析10年間植被覆蓋度時(shí)空變化特征，結(jié)論如下。

1）2006—2016年，極低覆蓋度面積逐年減少，減少了37.66 km2;低覆蓋度和中高覆蓋度面積都呈現(xiàn)先增加后減少的變化趨勢(shì)，其中，低覆蓋度減少了11.84 km2;中覆蓋度和中高覆蓋度分別增加了58.56和110.95 km2;高覆蓋度面積減少最多，為120.02 km2。

2）甘南縣植被覆蓋度整體較好，主要以中高覆蓋度和高覆蓋度為主，2006—2011年低覆蓋度、中高覆蓋度植被面積均減少，極低覆蓋度、中覆蓋度和高覆蓋度植被面積均增加，高覆蓋度面積變化最大，共增加309.08 km2。2011—2016年中覆蓋度和高覆蓋度植被面積減少，極低覆蓋度、低覆蓋度、中高覆蓋度植被面積均增加，中高覆蓋度面積變化最大，共增加192.77 km2。

3）10年間研究區(qū)植被覆蓋度總體穩(wěn)定，2006—2011、2006—2016年穩(wěn)定區(qū)域面積比例分別為64.32%和65.54%。2011—2016年較2006—2011年退化面積減少257 km2，改善面積增加199 km2。2006—2016年植被退化區(qū)域大于改善區(qū)域，研究區(qū)生態(tài)環(huán)境趨于惡化。

參考文獻(xiàn)：

[1] 章文波，符素華，劉寶元.目估法測(cè)量植被覆蓋度的精度分析[J].北京師范大學(xué)學(xué)報(bào)（自然科學(xué)版），2001，37（3）：402-408.

[2] 秦? 偉，朱清科，張學(xué)霞，等.植被覆蓋度及其測(cè)算方法研究進(jìn)展[J].西北農(nóng)林科技大學(xué)學(xué)報(bào)（自然科學(xué)版），2006，34（9）：163-170.

[3] 牛寶茹，劉俊蓉，王政偉.干旱半干旱地區(qū)植被覆蓋度遙感信息提取研究[J].武漢大學(xué)學(xué)報(bào)（信息科學(xué)版），2005，30（1）：27-30.

[4] 鄧? 飛，全占軍，于云江.20年來烏蘭木倫河流域植被蓋度變化及影響因素[J].水土保持研究，2011，18（3）：137-140，152.

[5] 許? 旭，李曉兵，梁涵瑋，等.內(nèi)蒙古溫帶草原區(qū)植被蓋度變化及其與氣象因子的關(guān)系[J].生態(tài)學(xué)報(bào)，2010，30（14）：3733-3743.

[6] 周兆葉，儲(chǔ)少林，王志偉，等.基于NDVI的植被覆蓋度的變化分析-以甘肅省張掖市甘州區(qū)為例[J].草業(yè)科學(xué)，2008，25（12）：23-29.

[7] 程紅芳，章文波，陳? 鋒.植被覆蓋度遙感估算方法研究進(jìn)展[J].國(guó)土資源遙感，2008（1）：13-18.

[8] 蓋永芹，李曉兵，張? 立，等.土地利用/覆被變化與植被蓋度的遙感監(jiān)測(cè)-以北京市密云縣為例[J].資源科學(xué)，2009，31（3）：523-529.

[9] 孫久虎，劉曉萌，李佑鋼，等.北運(yùn)河地區(qū)植被覆蓋的遙感估算及變化分析[J].水土保持研究，2006，13（6）：97-99.

[10] 蘇? 嫄，王志杰.基于遙感和GIS的陜南地區(qū)近20年植被覆蓋時(shí)空變化特征[J].水土保持研究，2018，25（1）：250-256.

[11] ROUSE J W，HAAS R H，SCHELL J A，et al. Monitoring vegetation systems in the great plains with ERTS[C].Washington DC：Proceedings of the Third ERTS symposium，1973.309-317.

[12] 郭芬芬，范建容，嚴(yán)? 冬，等.基于像元二分模型的昌都縣植被蓋度遙感估算[J].中國(guó)水土保持，2010（5）：65-67.

[13] 雷聲劍，張福平，燕玉超，等.黑河上游植被覆蓋時(shí)空變化特征及其未來趨勢(shì)[J].水土保持通報(bào)，2016，36（3）：159-164.

[14] 肖? 驍，李京忠，韓? 彬，等.東北老工業(yè)區(qū)植被覆蓋度時(shí)空特征及城市化關(guān)聯(lián)分析[J].生態(tài)科學(xué)，2017，36（6）：71-77.

[15] 王? 平，沈潤(rùn)平.基于CLM模型的植被覆蓋變化對(duì)黃土高原氣溫和降水的影響研究[J].科學(xué)技術(shù)與工程，2013，13（20）：5754-5760.

[16] 王? 童，何? 海，吳志勇.近30年來窟野河流域土地利用與植被覆蓋度變化分析[J].水電能源科學(xué)，2017，35（11）：128-130，83.

[17] 許幼霞，周? 旭，趙? 娟，等.半濕潤(rùn)半干旱過渡區(qū)洮河流域植被蓋度變化特征[J].水土保持通報(bào)，2016，36（6）：308-314.

[18] 張學(xué)霞，葛全勝，鄭景云.近50年北京植被對(duì)全球變暖的響應(yīng)及其時(shí)效-基于遙感數(shù)據(jù)和物候資料的分析[J].生態(tài)學(xué)雜志，2005，24（2）：123-130.

[19] 夏照華，張克斌，李? 瑞，等.基于NDVI的農(nóng)牧交錯(cuò)區(qū)植被覆蓋度變化研究-以寧夏鹽池縣為例[J].水土保持研究，2006，13（6）：178-181.

[20] PUREVDORJ T S，TATEISHI R，ISHIYAMA T，et al. Relationships between percent vegetation cover and vegetation indices[J].International journal of remote sensing，2010，19（18）：3519-3535.

[21] 趙泰安，彭道黎.利用多時(shí)相遙感影像分析密云縣植被動(dòng)態(tài)變化[J].東北林業(yè)大學(xué)學(xué)報(bào)，2013，41（12）：30-34.

[22] 陳? 晉，陳云浩，何春陽，等.基于土地覆蓋分類的植被覆蓋率估算亞像元模型與應(yīng)用[J].遙感學(xué)報(bào)，2001，5（6）：416-422.

[23] 龔文峰，吳? 娟，王笑峰，等.基于RS與GIS密云區(qū)植被蓋度動(dòng)態(tài)演變及其預(yù)測(cè)研究[J].黑龍江大學(xué)工程學(xué)報(bào)，2016，7（3）：25-31.

[24] 郝敬鋒，劉紅玉，李玉鳳，等.基于轉(zhuǎn)移矩陣模型的江蘇海濱濕地資源時(shí)空演變特征及驅(qū)動(dòng)機(jī)制分析[J].自然資源學(xué)報(bào)，2010， 25（11）：1918-1929.

[25] 劉玉林，韋瑞生，陳偉杰.甘南縣水土流失現(xiàn)狀及防治措施[J].黑龍江水利科技，2012，40（1）：274-275.