李鈺溦， 賈坤， 魏香琴， 姚云軍， 孫俊， 牟麗秋

(1．北京師范大學(xué)地理學(xué)與遙感科學(xué)學(xué)院遙感科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室，北京 100875; 2．中國(guó)科學(xué)院遙感與數(shù)字地球研究所，北京 100101; 3．西林吉林業(yè)局綜合設(shè)計(jì)隊(duì)，大興安嶺 165300;4．大興安嶺神州北極木業(yè)有限責(zé)任公司宜家分公司，大興安嶺 165300)

中國(guó)北方地區(qū)植被覆蓋度遙感估算及其變化分析

李鈺溦1， 賈坤1， 魏香琴2， 姚云軍1， 孫俊3， 牟麗秋4

(1．北京師范大學(xué)地理學(xué)與遙感科學(xué)學(xué)院遙感科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室，北京 100875; 2．中國(guó)科學(xué)院遙感與數(shù)字地球研究所，北京 100101; 3．西林吉林業(yè)局綜合設(shè)計(jì)隊(duì)，大興安嶺 165300;4．大興安嶺神州北極木業(yè)有限責(zé)任公司宜家分公司，大興安嶺 165300)

為了分析中國(guó)北方地區(qū)2000年之后植被覆蓋度的時(shí)空分布及其變化，利用MODIS光譜反射率數(shù)據(jù)計(jì)算歸一化植被指數(shù)，采用像元二分模型對(duì)中國(guó)北方地區(qū)2000—2012年植被覆蓋度進(jìn)行定量估算，分析研究區(qū)13 a間植被覆蓋度的時(shí)空變化特征。研究結(jié)果表明：植被覆蓋度年內(nèi)變化特征體現(xiàn)在最大植被覆蓋度一般出現(xiàn)在7和8月份，與中國(guó)北方地區(qū)植被的生長(zhǎng)季相一致；整個(gè)中國(guó)北方地區(qū)年最大植被覆蓋度呈現(xiàn)緩慢增長(zhǎng)的趨勢(shì)，其增長(zhǎng)速率為每年0.2%；年最大植被覆蓋度變化的空間分布具有較大差異，其中東北、華北和黃土高原等三北防護(hù)林工程建設(shè)區(qū)的年最大植被覆蓋度有較明顯的增長(zhǎng)。

像元二分模型；中國(guó)北方；植被覆蓋度；時(shí)空變化特征

0 引言

植被覆蓋度(fractional vegetation cover，F(xiàn)VC)通常定義為植被(包括葉、莖及枝)在地面的垂直投影面積占統(tǒng)計(jì)區(qū)總面積的百分比[1]，是刻畫地表植被覆蓋狀況的重要參數(shù)[2]。植被覆蓋度及其變化特征是區(qū)域生態(tài)系統(tǒng)環(huán)境變化的重要指示[3]，在水文、生態(tài)及全球變化等領(lǐng)域研究中有重要意義[4]。

植被覆蓋度信息的獲取主要有地面測(cè)量和遙感估算2種方法。根據(jù)測(cè)量原理，地面測(cè)量方法主要分為目估法、采樣法、儀器測(cè)量法及模型法等[5]。地面測(cè)量方法的結(jié)果相對(duì)準(zhǔn)確、客觀；但費(fèi)時(shí)、費(fèi)力，而且只適用于小尺度范圍，很難在較大空間尺度上連續(xù)進(jìn)行。多光譜、多時(shí)相及多尺度的遙感數(shù)據(jù)為大區(qū)域甚至全球尺度連續(xù)監(jiān)測(cè)植被覆蓋度提供了可能。利用遙感數(shù)據(jù)獲取不同時(shí)空尺度的植被覆蓋度信息，已經(jīng)成為估算區(qū)域植被覆蓋度的主要技術(shù)手段。

植被覆蓋度遙感估算方法主要有應(yīng)用較廣的回歸模型法、混合像元分解法及機(jī)器學(xué)習(xí)法[6-8]，以及有應(yīng)用較少的物理模型法[9]、光譜梯度差法[10]及森林冠層密度分級(jí)法(forest canopy density，F(xiàn)CD)[11]等。回歸模型法是對(duì)植被覆蓋度與遙感數(shù)據(jù)某光譜波段、光譜波段組合或利用光譜波段計(jì)算出的植被指數(shù)進(jìn)行回歸分析，并建立反演模型[12]。這種方法對(duì)于局部區(qū)域具有較好的精度，但是需要大量地面測(cè)量數(shù)據(jù)，受時(shí)空條件限制大，只能針對(duì)特定的區(qū)域和植被類型，大范圍地推廣使用具有不確定性?；旌舷裨纸饽Ｐ图僭O(shè)每個(gè)地表組分對(duì)傳感器所觀測(cè)到的光譜反射率信號(hào)都有貢獻(xiàn)， 通過(guò)求解植被組分占混合像元所有組分的比例估算植被覆蓋度[7,13-14]?；旌舷裨纸馑惴ǘ嗖捎镁€性光譜混合理論，其中像元二分模型是線性混合像元分解模型中最為廣泛應(yīng)用的植被覆蓋度遙感估算模型[14]。機(jī)器學(xué)習(xí)算法主要包括神經(jīng)網(wǎng)絡(luò)[15-16]、決策樹[17]及支持向量機(jī)[18-19]等。機(jī)器學(xué)習(xí)算法的步驟一般為確定訓(xùn)練樣本、訓(xùn)練模型和估算FVC。機(jī)器學(xué)習(xí)算法的關(guān)鍵在于要確保所選訓(xùn)練樣本的準(zhǔn)確性和代表性。地面實(shí)測(cè)數(shù)據(jù)雖然可以獲得較高精度的訓(xùn)練樣本，但是在代表性和全面性方面具有一定的限制。

本研究旨在應(yīng)用長(zhǎng)時(shí)間序列的遙感數(shù)據(jù)，估算中國(guó)北方區(qū)域植被覆蓋度，并分析其時(shí)空變化特征。研究采用像元二分模型，基于MODIS數(shù)據(jù)歸一化植被指數(shù)(normalized difference vegetation index，NDVI)估算了中國(guó)北方地區(qū)2000—2012年13 a間的植被覆蓋度，并分析其時(shí)空變化特征，為進(jìn)一步評(píng)價(jià)中國(guó)北方地區(qū)環(huán)境演化提供基礎(chǔ)數(shù)據(jù)。

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

1.1 研究區(qū)概況

我國(guó)北方地區(qū)曾經(jīng)是風(fēng)沙危害、水土流失比較嚴(yán)重的區(qū)域。1979年，國(guó)家決定在西北、華北平原北部及東北西部等地區(qū)建設(shè)大型防護(hù)林工程，并在近幾年實(shí)施了大規(guī)模的植被保護(hù)與恢復(fù)工程，如退耕還林還草工程、天然林保護(hù)工程及京津風(fēng)沙源治理工程等。中國(guó)北方地區(qū)植被覆蓋情況發(fā)生了巨大的變化。因此，分析中國(guó)北方區(qū)域長(zhǎng)時(shí)間序列的植被覆蓋度時(shí)空變化特征對(duì)于評(píng)價(jià)各種生態(tài)工程效益和分析環(huán)境演化特征具有重要意義。

本文研究廣義上的北方地區(qū)，即淮河—秦嶺及其西線的延伸及其以北地區(qū)(涉及東北地區(qū)、華北地區(qū)、以及青海省以外的西北地區(qū))。由于青藏高原地區(qū)氣候環(huán)境具有一定的獨(dú)特性，因此本研究不涉及青藏高原北部地區(qū)。故研究區(qū)域以省界為基礎(chǔ)，由東向西包括了黑龍江、吉林、遼寧、北京、天津、河北、山西、內(nèi)蒙古、陜西、寧夏、甘肅及新疆等省(市、自治區(qū))。用于驗(yàn)證MODIS植被覆蓋度反演精度的Landsat數(shù)據(jù)覆蓋了河北省北部，中心點(diǎn)坐標(biāo)為E115°32′52″，N41°45′58″。

1.2 數(shù)據(jù)源及其預(yù)處理

本文獲取的用于估算中國(guó)北方地區(qū)長(zhǎng)時(shí)間序列植被覆蓋度的MODIS光譜反射率數(shù)據(jù)已經(jīng)過(guò)大氣校正等預(yù)處理[20]。該數(shù)據(jù)處理綜合利用時(shí)間、空間和光譜信息，減少了云雪誤判和漏判，填補(bǔ)了數(shù)據(jù)的缺失，經(jīng)過(guò)時(shí)間序列的平滑，保證了時(shí)空覆蓋的完整性[20]。MODIS數(shù)據(jù)的空間分辨率為1 km，時(shí)間分辨率為8 d，時(shí)間范圍為2000—2012年，其中2000年缺少1月份數(shù)據(jù)，每月共有17個(gè)格網(wǎng)覆蓋研究區(qū)。本文MODIS數(shù)據(jù)處理主要包括：利用紅波段和近紅外光譜波段計(jì)算NDVI，按照月最大NDVI值對(duì)時(shí)間分辨率重采樣，獲得月分辨率的數(shù)據(jù)，并用MODIS重投影工具(MODIS reprojection tool，MRT)對(duì)數(shù)據(jù)進(jìn)行拼接，并轉(zhuǎn)換為經(jīng)緯度投影，最終裁剪出研究區(qū)域數(shù)據(jù)。為了檢驗(yàn)MODIS數(shù)據(jù)植被覆蓋度估算精度，本研究利用30 m空間分辨率的Landsat 5和Landsat 7衛(wèi)星多光譜數(shù)據(jù)計(jì)算的植被覆蓋度進(jìn)行精度驗(yàn)證。隨機(jī)選取軌道號(hào)為14031的Landsat數(shù)據(jù)作為驗(yàn)證樣區(qū)，選取該樣區(qū)4個(gè)季度的4景云覆蓋較少和質(zhì)量較好的數(shù)據(jù)作為MODIS數(shù)據(jù)植被覆蓋度驗(yàn)證的樣本數(shù)據(jù)，獲取日期分別是2000年4月5日、2001年7月4日、2006年9月20日及2006年10月6日。Landsat數(shù)據(jù)的預(yù)處理主要包括大氣校正和云雪掩模，并投影轉(zhuǎn)換為與MODIS一致的經(jīng)緯度投影。最后利用大氣校正后的紅波段和近紅外波段光譜反射率數(shù)據(jù)計(jì)算NDVI值，進(jìn)而采用像元二分模型計(jì)算植被覆蓋度。

2 研究方法

2.1 像元二分模型

像元二分模型假設(shè)像元由植被和非植被2種端元組成。李苗苗等[13]討論了不同植被指數(shù)的特點(diǎn)，指出NDVI在利用像元二分模型估算植被覆蓋度中具有更好的表現(xiàn)。因此，本文像元二分模型計(jì)算植被覆蓋度的公式為

(1)

式中：FVC為植被覆蓋度；NDVIsoil為完全是裸土或者無(wú)植被覆蓋像元(即陸地土壤像元)的NDVI；NDVIveg為完全由植被所覆蓋像元(即純植被像元)的NDVI。只要確定NDVIsoil和NDVIveg的值，就能計(jì)算每一個(gè)像元的植被覆蓋度。

對(duì)于純土壤像元的NDVI，理論上應(yīng)接近于0，但是由于受到環(huán)境因素等影響，NDVIsoil通常在-0.1～0.2之間變化。一般情況下，NDVIsoil與NDVIveg取給定置信度的置信區(qū)間內(nèi)的最大值與最小值，可在一定程度上消除遙感圖像噪聲所帶來(lái)的誤差。本研究中，NDVIsoil與NDVIveg的取值是統(tǒng)計(jì)所有155景MODIS NDVI數(shù)據(jù)，并結(jié)合研究區(qū)植被覆蓋的實(shí)際情況得出的。根據(jù)頻率統(tǒng)計(jì)表，NDVIsoil是選取區(qū)域內(nèi)累計(jì)頻率為5%的像元的NDVI值，NDVIveg為累計(jì)頻率為90%的NDVI像元值；在統(tǒng)計(jì)計(jì)算的過(guò)程中，當(dāng)像元值≤0時(shí)，認(rèn)為是水體，不參與統(tǒng)計(jì)計(jì)算。在像元二分模型中，NDVIsoil和NDVIveg取統(tǒng)一的值，使得長(zhǎng)時(shí)間序列的數(shù)據(jù)都在一個(gè)相同的基準(zhǔn)上估算植被覆蓋度，因而其變化趨勢(shì)分析更為可靠。當(dāng)累計(jì)頻率分別為5%和90%時(shí)，NDVIsoil和NDVIveg的取值分別為0.013 3和0.848。

2.2 MODIS植被覆蓋度精度驗(yàn)證

利用像元二分模型對(duì)Landsat多光譜數(shù)據(jù)進(jìn)行植被覆蓋度計(jì)算，NDVIsoil與NDVIveg分別選取NDVI累積頻率5%和95%時(shí)的值[21]。對(duì)于每一景Landsat影像，找到其相同時(shí)相MODIS數(shù)據(jù)計(jì)算的植被覆蓋度(盡量選擇時(shí)相最接近的)，對(duì)MODIS數(shù)據(jù)計(jì)算的植被覆蓋度進(jìn)行精度檢驗(yàn)。按照Landsat數(shù)據(jù)中心點(diǎn)所在經(jīng)緯度坐標(biāo)，選擇MODIS數(shù)據(jù)的9像元×9像元，然后計(jì)算MODIS每一個(gè)像元對(duì)應(yīng)的所有Landsat數(shù)據(jù)像元的植被覆蓋度均值。

2.3 植被覆蓋度趨勢(shì)變化分析

采用線性回歸方法對(duì)年最大植被覆蓋度年際變化趨勢(shì)進(jìn)行分析[22]。逐像元計(jì)算年最大植被覆蓋度年際變化的最小二乘線性回歸方程的斜率，即

(2)

式中：n為研究階段的年數(shù)(本文n=13)；FVCi為第i年最大植被覆蓋度；NDVIi為第i年某像元的植被覆蓋度；θslope為研究區(qū)間內(nèi)植被覆蓋度變化斜率，反映了在研究期間植被覆蓋度的變化趨勢(shì)。當(dāng)θslope>0時(shí)，表示植被覆蓋度增加；當(dāng)θslope=0時(shí)表示植被覆蓋度沒有變化；當(dāng)θslope<0時(shí)，則表示植被覆蓋度降低。

3 結(jié)果與分析

3.1 中國(guó)北方地區(qū)植被覆蓋度估算

根據(jù)水利部2008年頒布的《土壤侵蝕分類分級(jí)標(biāo)準(zhǔn)》中植被覆蓋度分級(jí)標(biāo)準(zhǔn)，將植被覆蓋度劃分為5個(gè)等級(jí)：低植被覆蓋度[0,30%)、中低覆蓋度[30%,45%)、中等覆蓋度[45%,60%)、中高覆蓋度[60%,75%)和高覆蓋度[75%,100%]。利用MODIS NDVI數(shù)據(jù)和像元二分模型估算中國(guó)北方地區(qū)2000—2012年長(zhǎng)時(shí)間序列的植被覆蓋度，并進(jìn)一步利用年最大植被覆蓋度，合成得到中國(guó)北方地區(qū)多年最大植被覆蓋度圖，直觀反映北方地區(qū)植被覆蓋度空間分布特征，如圖1所示。

圖1 北方地區(qū)最大植被覆蓋度空間分布Fig.1 Spatial distribution of maximum FVC in the north region of China

從圖1可以看出，中國(guó)北方地區(qū)植被覆蓋度較大的地區(qū)主要在東北地區(qū)，由于該區(qū)域有大、小興安嶺等原始森林、遼闊的草原以及松-嫩-遼平原農(nóng)業(yè)區(qū)等，因此最大植被覆蓋度都處于高植被覆蓋度等級(jí)。華北平原和黃土高原南部等地區(qū)是我國(guó)重要的農(nóng)業(yè)產(chǎn)區(qū)，在農(nóng)作物生長(zhǎng)茂盛的夏季，植被覆蓋度也相對(duì)較高。

為定量評(píng)估MODIS數(shù)據(jù)估算植被覆蓋度的精度，采用Landsat數(shù)據(jù)計(jì)算的植被覆蓋度對(duì)MODIS數(shù)據(jù)估算結(jié)果進(jìn)行定量的精度評(píng)價(jià)，如圖2所示。

圖2 MODIS數(shù)據(jù)植被覆蓋度和其對(duì)應(yīng)的TM數(shù)據(jù)植被覆蓋度散點(diǎn)圖Fig.2 Scatter diagram of the MODIS FVC data with the corresponding TM FVC data

在精度評(píng)價(jià)過(guò)程中，對(duì)于Landsat數(shù)據(jù)上有云的像元進(jìn)行剔除。MODIS數(shù)據(jù)估算的植被覆蓋度與對(duì)應(yīng)的Landsat數(shù)據(jù)植被覆蓋度的線性回歸系數(shù)為0.967，常數(shù)項(xiàng)為0.034，R2為0.82，顯著性較為明顯；回歸方程計(jì)算殘差平方和為1.36，中誤差為0.06可以看出中誤差相對(duì)較小，表明MODIS數(shù)據(jù)估算植被覆蓋度和Landsat數(shù)據(jù)植被覆蓋度相差不大，具有較高的估算精度。

MODIS數(shù)據(jù)估算結(jié)果可以用來(lái)分析中國(guó)北方地區(qū)13 a來(lái)的植被覆蓋度的動(dòng)態(tài)變化特征和規(guī)律。但是，有些數(shù)據(jù)點(diǎn)偏差較大，可能與2種數(shù)據(jù)的成像時(shí)間有關(guān)，因?yàn)镸ODIS數(shù)據(jù)是按照1個(gè)月為時(shí)間分辨率計(jì)算的，目前還很難找到2種在時(shí)間上精確匹配的數(shù)據(jù)。

3.2 北方地區(qū)植被覆蓋度年內(nèi)變化規(guī)律

中國(guó)北方地區(qū)植被的季相變化明顯，因此植被覆蓋度在不同季節(jié)也有明顯變化。本文以我國(guó)東北地區(qū)多年植被覆蓋度月均值為例說(shuō)明1 a中植被覆蓋度的變化特點(diǎn)。圖3為我國(guó)東北地區(qū)1—12月份2000—2012年植被覆蓋度的對(duì)應(yīng)月份平均值的變化曲線。從中可以看出，植被覆蓋度最大月均值出現(xiàn)在8月份，7月次之，并且7—9月植被覆蓋度變化不大；最小月均值則出現(xiàn)在1—2月，4—7月為植被快速生長(zhǎng)期，進(jìn)入秋季之后植被開始落葉，在9—11月，植被覆蓋度迅速減小。統(tǒng)計(jì)每年的最大和最小月均值可以發(fā)現(xiàn)，植被覆蓋度的變化和植被季相變化基本一致。

圖3 東北地區(qū)植被覆蓋度月均值變化曲線Fig.3 Month mean FVC value of northeast region

對(duì)2012年不同季度的植被覆蓋度進(jìn)行分等級(jí)制圖，如圖4所示。

圖4 2012年不同季度植被覆蓋度等級(jí)分布圖

Fig.4 Grade distribution of FVC in different quarters in 2012

由圖4可以看出，第1季度植被覆蓋度普遍偏低，即使在東北的原始森林區(qū)，植被覆蓋度也處于中等水平，其他地區(qū)大多是低植被覆蓋度水平。第2季度是植被返青和快速生長(zhǎng)的季節(jié)，植被覆蓋度變化較為明顯。在第3季度中期(7—8月)植被覆蓋度達(dá)到最大，特別是東北地區(qū)、華北地區(qū)和黃土高原地區(qū)。內(nèi)蒙古西部、甘肅北部和新疆地區(qū)則變化不明顯，且大多處于低覆蓋度水平，這與土地覆蓋類型有關(guān)，比如內(nèi)蒙古西部多為荒漠地區(qū)。第3季度的植被覆蓋度明顯比其他3個(gè)季度大，第1季度最小，這與圖3中曲線變化特點(diǎn)相一致。年最大植被覆蓋度和第3季度的分布基本一致，大部分像元的植被覆蓋度在第3季度達(dá)到最大值，充分說(shuō)明植被覆蓋度的變化和植被的生長(zhǎng)季相一致，在生長(zhǎng)最為旺盛的第3季度植被覆蓋度最大。第4季度植被覆蓋度和第1季度相差不大。

3.3 北方地區(qū)植被覆蓋度年際變化趨勢(shì)

圖5是2000和2012年的年最大和年平均不同等級(jí)植被覆蓋度占區(qū)域面積的比例。對(duì)于年最大植被覆蓋度，2012年高等植被覆蓋度的比例相對(duì)2000年高出很多，低、中低、中和中高等級(jí)則相對(duì)2000年較低。在年平均植被覆蓋度中，2012年的中低、中和高等級(jí)植被覆蓋度都較2000年高。因此，從整體來(lái)看，北方地區(qū)2012年植被覆蓋度比2000年有較大幅度的升高。

圖5 2000年和2012年的年最大(左)和年平均(右)植被覆蓋度不同等級(jí)比例Fig.5 Maximum(left) and mean annual(right) FVC in different grade proportion in 2000 and 2012

為分析2000—2012年中國(guó)北方地區(qū)的植被覆蓋度總體變化趨勢(shì)，分別統(tǒng)計(jì)每月區(qū)域均值，計(jì)算年最大植被覆蓋度的變化趨勢(shì)(圖6)。

圖6 2000—2012年植被覆蓋度變化曲線(月分辨率)Fig.6 Variation profile of FVC from 2000 to 2012 (month resolution)

由圖6可知，植被覆蓋度隨時(shí)間呈周期性變化，在不同的年份，其月平均值有一定的浮動(dòng)，但是整體呈增長(zhǎng)趨勢(shì)，增長(zhǎng)速率可達(dá)到每年0.2%。

提取2000—2012年的MODIS數(shù)據(jù)每年最大植被覆蓋度，然后計(jì)算植被覆蓋度的變化趨勢(shì)，并按照斜率大小，劃分不同的等級(jí)，得到2000—2012年的植被覆蓋度變化趨勢(shì)空間分布圖，如圖7所示。

圖7 2000—2012年植被覆蓋度變化趨勢(shì)空間分布Fig.7 Spatial distribution of FVC change trend

統(tǒng)計(jì)不同變化等級(jí)的面積比例，如表1所示。

表1 植被覆蓋度變化斜率不同等級(jí)范圍的面積比例Tab.1 Area percentage in different grade scale of FVC change slopes

從圖7和表1可以看出，在2000—2012年13 a間，61.2%的地區(qū)植被覆蓋度沒有發(fā)生明顯的變化，植被覆蓋度降低的區(qū)域主要分布在新疆北部和內(nèi)蒙古東部的一些地區(qū)。但是，在內(nèi)蒙古東部邊界(位于東北造林區(qū)內(nèi))、內(nèi)蒙古中南部和寧夏、陜西北部(黃土高原造林區(qū))、華北平原，北部地區(qū)植被覆蓋度增加趨勢(shì)非常明顯。

4 結(jié)論

1)北方地區(qū)植被的季相特征明顯，植被覆蓋度呈現(xiàn)明顯的年周期變化規(guī)律，每年最大植被覆蓋度一般出現(xiàn)在7—8月，最小值一般在1—2月。

2)在2000—2012年這13 a間，北方大部分地區(qū)植被覆蓋度變化趨勢(shì)相對(duì)穩(wěn)定，但具有明顯的空間分布差異；整個(gè)北方地區(qū)年最大植被覆蓋度呈緩慢增長(zhǎng)的趨勢(shì)，其增長(zhǎng)率為每年0.2%。

下一步需要結(jié)合氣候資料和地面調(diào)查，對(duì)造成植被覆蓋度變化的自然因素和人為因素進(jìn)行分析。

[1] 梁順林,李小文,王錦地,等.定量遙感:理念與算法[M].北京:科學(xué)出版社,2003:382. Liang S L,Li X W,Wang J D,et al.Advanced Remote Sensing:Concepts and Algorithm[M].Beijing:Science Press,2003:382.

[2] 程紅芳,章文波,陳鋒.植被覆蓋度的遙感估算方法研究進(jìn)展[J].國(guó)土資源遙感,2008,20(1):13-18.doi:10.6046/gtzyyg.2008.01.02. Cheng H F,Zhang W B,Chen F.Advances in researches on application of remote sensing method to estimating vegetation coverage[J].Remote Sensing for Land and Resources,2008,20(1):13-18.doi:10.6046/gtzyyg.2008.01.02.

[3] 邢著榮,馮幼貴,楊貴軍,等.基于遙感的植被覆蓋度估算方法述評(píng)[J].遙感技術(shù)與應(yīng)用,2009,24(6):849-854. Xing Z R,Feng Y G,Yang G J,et al.Method of estimating vegetation coverage based on remote sensing[J].Remote Sensing Technology and Application,2009,24(6):849-854.

[4] Zeng X B,Dickinson R E,Walker A,et al.Derivation and evaluation of global 1 km fractional vegetation cover data for land modeling[J].Journal of Applied Meteorology,2000,39(6):826-839.

[5] 江輝.基于遙感的植被覆蓋度估算及其動(dòng)態(tài)研究——以鄱陽(yáng)湖區(qū)為例[D].南昌:南昌大學(xué),2005. Jiang H.Study on Vegetation Coverage and its Dynamic Change by Remote Sesing:A Case Study of the Poyang Lake[D].Nanchang:Nanchang University,2005.

[6] Van de Voorde T,Vlaeminck J,Canters F.Comparing different approaches for mapping urban vegetation cover from Landsat ETM+data:A case study on brussels[J].Sensors,2008,8(6):3880-3902.

[7] Xiao J F,Moody A.A comparison of methods for estimating fractional green vegetation cover within a desert-to-upland transition zone in central New Mexico,USA[J].Moody A.Remote Sensing of Environment,2005,98(2/3):237-250.

[8] North P R J.Estimation of FAPAR,LAI,and vegetation fractional cover from ATSR-2 imagery[J].Remote Sensing of Environment,2002,80(1):114-121.

[9] Gutman G,Ignatov A.The derivation of the green vegetation fraction from NOAA/AVHRR data for use in numerical weather prediction models[J].International Journal of Remote Sensing,1998,19(8):1533-1543.

[10]唐世浩,朱啟疆,王錦地,等.三波段梯度差植被指數(shù)的理論基礎(chǔ)及其應(yīng)用[J].中國(guó)科學(xué):D輯,2003,33(11):1094-1102. Tang S H,Zhu Q J,Wang J D,et al.Principle and application of three-band gradient difference vegetation index[J].Science in China Series D:Earth Sciences,2005,48(2):241-249.

[11]賈坤,姚云軍,魏香琴,等.植被覆蓋度遙感估算研究進(jìn)展[J].地球科學(xué)進(jìn)展,2013,28(7):774-782. Jia K,Yao Y J,Wei X Q,et al.A review on fractional vegetation cover estimation using remote sensing[J].Advances in Earth Science,2013,28(7):774-782.

[12]宋莎.基于多源遙感數(shù)據(jù)的植被覆蓋度研究[D].成都:四川農(nóng)業(yè)大學(xué),2010. Song S.Retrieval of Vegetation Coverage Using Multi-sensor Remote Sensing Data[D].Chengdu:Sichuan Agricultural University,2010.

[13]李苗苗,吳炳芳,顏長(zhǎng)珍,等.密云水庫(kù)上游植被覆蓋度的遙感估算[J].資源科學(xué),2004,26(4):153-159. Li M M,Wu B F,Yan C Z,et al.Estimaion of vegetation fraction in the upper basin of Miyun Reservior by remote sensing[J].Resources Science,2004,26(4):153-159.

[14]陳晉,陳云浩,何春陽(yáng),等.基于土地覆蓋分類的植被覆蓋率估算亞像元模型與應(yīng)用[J].遙感學(xué)報(bào),2001,5(6):416-422,481. Chen J,Chen Y H,He C Y,et al.Sub-pixel model for vegetation fraction estimation based on land cover classification[J].Journal of Remote Sensing,2001,5(6):416-422,481.

[15]Camacho F,Cernicharo J,Lacaze R,et al.GEOV1:LAI,FAPAR essential climate variables and FCOVER global time series capitalizing over existing products.Part 2:Validation and intercomparison with reference products[J].Remote Sensing of Environment,2013,137:310-329.

[16]Baret F,Weiss M,Lacaze R,et al.GEOV1:LAI and FAPAR essential climate variables and FCOVER global time series capitalizing over existing products.Part1:Principles of development and production[J].Remote Sensing of Environment,2013,137:299-309.

[17]Hansen M C,DeFries R S,Townshend J R G,et al.Towards an operational MODIS continuous field of percent tree cover algorithm:examples using AVHRR and MODIS data[J].Remote Sensing of Environment,2002,83(1/2):303-319.

[18]Huang C Q,Song K,Kim S,et al.Use of a dark object concept and support vector machines to automate forest cover change analysis[J].Remote Sensing of Environment,2008,112(3):970-985.

[19]Su L H.Optimizing support vector machine learning for semi-arid vegetation mapping by using clustering analysis[J].ISPRS Journal of Photogrammetry and Remote Sensing,2009,64(4):407-413.

[20]梁順林,袁文平,肖青,等.全球陸表特征參量產(chǎn)品生成與應(yīng)用研究[J].中國(guó)科學(xué)院院刊,2013,28(s1):122-131. Liang S L,Yuan W P,Xiao Q,et al.Generation and applications of global land surface satellite(GLASS)products[J].Bulletin of Chinese Academy of Sciences,2013,28(s1):122-131.

[21]賈維花,廉麗姝,呂宜平,等.基于TM數(shù)據(jù)的黃河三角洲地區(qū)植被覆蓋度提取[J].地理信息世界,2012,10(6):62-66,74. Jia W H,Lian L S,Lv Y P,et al.The derivation of vegetation fraction based on TM data in Yellow River Delta[J].Geomatics World,2012,10(6):62-66,74.

[22]Duan H C,Yan C Z,Tsunekawa A,et al.Assessing vegetation dynamics in the Three-North Shelter Forest region of China using AVHRR NDVI data[J].Environmental Earth Sciences,2011,64(4):1011-1020.

(責(zé)任編輯：邢宇)

Fractional vegetation cover estimation in northern China and its change analysis

LI Yuwei1,JIA Kun1,WEI Xiangqin2,YAO Yunjun1,SUN Jun3,MOU Liqiu4

(1.StateKeyLaboratoryofRemoteSensingScience,SchoolofGeography,BeijingNormalUniversity,Beijing100875,China;2.InstituteofRemoteSensingandDigitalEarth,ChineseAcademyofSciences,Beijing100101,China; 3.ComprehensiveDesignTeamofXilinjiForestryBureau,DaHingganLing165300,China; 4.GreaterKhinganRangeShenzhouPolarWoodIndustryCo.,Ltd.Yijia,DaHingganLing165300,China)

The aim of this study is to analyze the spatial-temporal pattern and change of fractional vegetation cover (FVC) in northern China since 2000. FVC of northern China from 2000 to 2012 was estimated using dimidiate pixel model based on the normalized difference vegetation index (NDVI) calculated by MODIS spectral reflectance data. The FVC change trends and characteristics of the study area during the 13 years were analyzed. The inner annual FVC change trend in northern China indicated that the maximum FVC generally appeared in July or August,which was consistent with the vegetation growth season. The maximum annual FVC showed a slightly increase trend in the whole study area,and the annual increase rate was 0.2%. However,the spatial distribution of maximum annual FVC change trend had great differences. The typical regions of the Three-North Shelter Forest Region such as Northeast China,North China and Loess Plateau region had an obviously increase in maximum annual FVC.

dimidiate pixel model; northern China; fractional vegetation cover; spatial-temporal change characteristics

2014-01-21；

2014-04-28

國(guó)家自然科學(xué)基金項(xiàng)目“基于機(jī)器學(xué)習(xí)和融合算法的全球陸表植被覆蓋度遙感估算方法研究”(編號(hào)：41301353)、遙感科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室自由探索項(xiàng)目“多尺度植被覆蓋度觀測(cè)實(shí)驗(yàn)與遙感產(chǎn)品驗(yàn)證數(shù)據(jù)集生產(chǎn)”(編號(hào)：14ZY-06)、國(guó)家“863”計(jì)劃項(xiàng)目“全球生態(tài)系統(tǒng)與表面能量平衡特征參量生成與應(yīng)用”(編號(hào)：2013AA122801)和國(guó)家測(cè)繪地理信息局科技領(lǐng)軍人才科技資助項(xiàng)目共同資助。

10.6046/gtzyyg.2015.02.18

李鈺溦，賈坤，魏香琴，等.中國(guó)北方地區(qū)植被覆蓋度遙感估算及其變化分析[J].國(guó)土資源遙感,2015,27(2):112-117.(Li Y W,Jia K,Wei X Q,et al.Fractional vegetation cover estimation in northern China and its change analysis[J].Remote Sensing for Land and Resources,2015,27(2):112-117.)

TP 79; S 127

A

1001-070X(2015)02-0112-06

李鈺溦(1990-)，女，碩士研究生，主要從事定量遙感方面的研究。Email: yuwei_breeze@163.com。

賈坤(1983-)，男，博士，碩士生導(dǎo)師，主要從事定量遙感方面的研究。 Email: jiakun@bnu.edu.cn。