周志強(qiáng)，曾 源，張 磊，杜 鑫，吳炳方

南水北調(diào)中線水源區(qū)植被覆蓋度遙感監(jiān)測(cè)分析

周志強(qiáng)1，2，曾 源1，張 磊1，杜 鑫1，吳炳方1

(1.中國(guó)科學(xué)院遙感應(yīng)用研究所，北京 100101;2.西南林業(yè)大學(xué)，昆明 650224)

南水北調(diào)中線工程是我國(guó)大規(guī)?？缌饔蛘{(diào)水工程的一部分，開(kāi)展該區(qū)域植被覆蓋度變化的研究與分析，對(duì)于保護(hù)該區(qū)域的生態(tài)環(huán)境及水質(zhì)具有重要意義。該文以2000年和2009年兩期遙感圖像為本底數(shù)據(jù)，利用基于NDVI的像元二分模型對(duì)南水北調(diào)中線水源區(qū)的植被覆蓋度進(jìn)行了估算，并分析了該區(qū)植被覆蓋度的時(shí)空變化特征。結(jié)果表明:2000年該水源區(qū)植被覆蓋度的平均值為67.5%，2009年的平均值達(dá)到72%，植被覆蓋度總體呈增長(zhǎng)趨勢(shì);植被覆蓋度增幅的空間特征表現(xiàn)為水源區(qū)中部地區(qū)高，東西部地區(qū)相對(duì)較低;在不同植被類(lèi)型中，落葉針葉林的覆蓋度平均值增幅最大，草地覆蓋度增幅最小;位于水源區(qū)的大多數(shù)縣(市)的植被覆蓋度在近十年來(lái)都有不同程度的增加，其中柞水縣的植被覆蓋度平均值增長(zhǎng)幅度最大，這與國(guó)家實(shí)施退耕還林、封山育林、基本農(nóng)田建設(shè)等政策有關(guān)。

植被覆蓋度;南水北調(diào)中線水源區(qū);遙感監(jiān)測(cè)

0 引言

植被覆蓋度是指一個(gè)地區(qū)所有植物(包括葉、莖、枝)在地面的垂直投影面積占總面積的百分比［1－4］。在土地沙漠化評(píng)價(jià)、水土流失監(jiān)測(cè)和分布式水文模型建立中，植被覆蓋度是重要的輸入?yún)?shù)，也是衡量一個(gè)地區(qū)生態(tài)環(huán)境好壞的重要指標(biāo)之一。因此獲取地表植被蓋度及其變化信息，對(duì)于揭示地表生態(tài)環(huán)境中各要素的空間變化規(guī)律，探討生態(tài)環(huán)境變化的驅(qū)動(dòng)因子，分析評(píng)價(jià)區(qū)域生態(tài)環(huán)境質(zhì)量都具有重要的現(xiàn)實(shí)意義［5］。

近年來(lái)，國(guó)內(nèi)外學(xué)者在植被覆蓋度估測(cè)方面已開(kāi)展了大量研究，其方法可分為地表實(shí)測(cè)和遙感監(jiān)測(cè)兩類(lèi)。根據(jù)測(cè)量原理，地表實(shí)測(cè)分為目估法、采樣法、儀器法和模型法［6－7］，但這些方法費(fèi)時(shí)費(fèi)力，且對(duì)于大面積實(shí)測(cè)來(lái)說(shuō)，更加不現(xiàn)實(shí)，因而遙感監(jiān)測(cè)已逐漸成為估算植被覆蓋度的主要手段［8］。根據(jù)模型建立的原理，遙感監(jiān)測(cè)方法分為物理模型法和統(tǒng)計(jì)模型法，前者由于變量多且難以測(cè)量，而會(huì)影響植被覆蓋度的提取精度;后者應(yīng)用最廣的方法是回歸模型法、像元分解法和植被指數(shù)法，回歸模型法只適用于特定區(qū)域與特定的植被類(lèi)型，一般不易推廣［9］，而像元分解法和植被指數(shù)法則不需要建立回歸模型［10］，因而對(duì)地表實(shí)測(cè)數(shù)據(jù)的依賴較小，是一種從宏觀上監(jiān)測(cè)特定區(qū)域植被治理情況的方法之一［11］，其相對(duì)于回歸模型法更具有實(shí)用價(jià)值，經(jīng)驗(yàn)證后可推廣應(yīng)用到大范圍地區(qū)的植被監(jiān)測(cè)，形成通用的植被覆蓋度計(jì)算方法［12］。

南水北調(diào)中線工程是保持或恢復(fù)水源區(qū)生態(tài)環(huán)境以及促進(jìn)社會(huì)經(jīng)濟(jì)持續(xù)發(fā)展的動(dòng)力［13］。我國(guó)學(xué)者針對(duì)該區(qū)域的生態(tài)、水土流失等也已進(jìn)行了大量研究:李思悅通過(guò)分析該區(qū)面臨的水質(zhì)、植被破壞、水土流失等主要生態(tài)環(huán)境問(wèn)題，提出了以異齡復(fù)層的喬灌混交林和草相結(jié)合為主的植被恢復(fù)、庫(kù)岸帶生態(tài)屏障建設(shè)的方法來(lái)改善庫(kù)區(qū)脆弱的生態(tài)環(huán)境［14］;朱桂香針對(duì)該區(qū)水土流失嚴(yán)重問(wèn)題，提出建立政府行為與市場(chǎng)經(jīng)濟(jì)相結(jié)合的運(yùn)行機(jī)制、封禁區(qū)群眾生產(chǎn)生活補(bǔ)貼機(jī)制［15］;金蓉玲和沈大軍等對(duì)該區(qū)調(diào)水進(jìn)行了水資源評(píng)價(jià)，并分析了對(duì)漢江中下游的水文情勢(shì)及用水帶來(lái)的影響［16－17］。上述研究針對(duì)水源區(qū)植被的時(shí)空變化分析的較少，特別是在植被自然生長(zhǎng)特征和人為因素影響方面。由于植被覆蓋度可以反映水源區(qū)的植被生長(zhǎng)狀況，因此通過(guò)該區(qū)調(diào)水前后植被覆蓋度的比較，可全面分析南水北調(diào)工程對(duì)本地區(qū)植被生長(zhǎng)狀況的影響。沈澤昊也指出今后在該區(qū)生態(tài)保護(hù)方面要特別考慮植被(特別是河岸帶植被)恢復(fù)等問(wèn)題［18］。本文在前人工作的基礎(chǔ)上，以南水北調(diào)中線水源區(qū)為研究對(duì)象，以Landsat TM圖像和國(guó)產(chǎn)環(huán)境災(zāi)害衛(wèi)星HJ－1A/1B CCD圖像為遙感數(shù)據(jù)，利用像元二分法和植被指數(shù)法對(duì)2000年和2009年的植被覆蓋度進(jìn)行估算，通過(guò)分析該區(qū)域植被覆蓋度的時(shí)空變化特征，以掌握水源區(qū)植被保護(hù)措施的實(shí)施情況，為進(jìn)行生態(tài)環(huán)境治理政策的制定提供遙感信息支持。

1 研究區(qū)概況

南水北調(diào)中線水源區(qū)范圍為 N 31°20'～34°10'，E 106°～112°，位置如圖1所示。集水面積約為9.45萬(wàn)km2，涉及6省37個(gè)縣(市)，其中長(zhǎng)安縣、巫溪縣及巫山縣面積非常小，在統(tǒng)計(jì)中并入周邊各縣。

圖1 南水北調(diào)中線水源區(qū)位置Fig.1 Regional location of water source area of the middle route of projects to divert water from the south to the north

研究區(qū)內(nèi)地貌形態(tài)多樣，特點(diǎn)主要是以中、低山和丘陵為主，一般高程在500～2 000 m之間，坡度陡、切割深、高差大;水源區(qū)氣候?yàn)楸眮啛釒Ъ撅L(fēng)氣候，多年平均氣溫12.2℃，多年平均降水量880 mm，緯向地帶性降水特征明顯。丹江口水庫(kù)是南水北調(diào)中線工程起點(diǎn)，也是漢江上游與中下游的分界點(diǎn)，其庫(kù)容量大，水體稀釋自凈能力強(qiáng)，水質(zhì)較好。

2 研究方法

2.1 數(shù)據(jù)采集與預(yù)處理

1)遙感數(shù)據(jù)源。HJ－1A/1B衛(wèi)星CCD遙感數(shù)據(jù)的時(shí)間分辨率高、覆蓋范圍大，識(shí)別季節(jié)性變化強(qiáng)、光譜易混淆的植被類(lèi)型的識(shí)別能力強(qiáng)，因此本研究采用2009年6月5日6景30 m空間分辨率的HJ－1A/1B衛(wèi)星CCD遙感數(shù)據(jù)(波譜范圍0.43～0.9 μm［19］)進(jìn)行植被覆蓋度監(jiān)測(cè)。為了比較不同年份植被覆蓋度的變化，本文還采用了與其相同空間分辨率的2000年6—9月10景Landsat TM遙感數(shù)據(jù)(波譜范圍0.45 ～2.35 μm)來(lái)進(jìn)行植被覆蓋度監(jiān)測(cè)。

2)遙感數(shù)據(jù)預(yù)處理。由于大氣的吸收和散射會(huì)引起遙感圖像的輻射失真，而衛(wèi)星傳感器系統(tǒng)和地面系統(tǒng)中的儀器在接收、轉(zhuǎn)換、傳送及處理圖像信息的過(guò)程中也會(huì)使圖像產(chǎn)生幾何位置的變化［20］，因此在預(yù)處理過(guò)程中本文采用ERDAS ATCOR2模型，對(duì)TM和HJ－1A/1B CCD兩期遙感圖像分別進(jìn)行了輻射校正，并采用二次多項(xiàng)式幾何糾正模型進(jìn)行幾何精糾正，誤差控制在0.5個(gè)像元內(nèi)。

3)野外樣本采集。野外調(diào)查中，采集具有代表性、均質(zhì)的植被樣本122個(gè)，利用GPS定位樣本坐標(biāo)，通過(guò)魚(yú)眼相機(jī)獲取植被信息，初步掌握了研究區(qū)內(nèi)植被的空間分布情況。

2.2 植被覆蓋度估算方法

2.2.1 像元二分模型

像元二分模型［21－23］是線性模型中常見(jiàn)模型之一。假定遙感圖像中一個(gè)像元對(duì)應(yīng)地面一定面積的光譜信息是由兩個(gè)組分因子(有植被覆蓋的地表因子和無(wú)植被覆蓋的地表因子)的光譜信息線性組合而來(lái)。每個(gè)組分對(duì)傳感器所獲得的光譜信息都有貢獻(xiàn)，它們各自在像元對(duì)應(yīng)的地面面積中所占的比例即為該組分因子所貢獻(xiàn)信息的權(quán)重，其中植被覆蓋地表部分所占的百分比即為該像元的植被覆蓋度［24］。

根據(jù)上述像元二分模型，可以得到公式

該式可理解為傳感器所獲得的信息S由完全植被覆蓋地表所貢獻(xiàn)的信息Sveg和無(wú)植被覆蓋地表所貢獻(xiàn)的信息Ssoil組成;fc為植被覆蓋度，fc和(1－fc)分別為Sveg和Ssoil在該像元中所貢獻(xiàn)信息的權(quán)重。

對(duì)式(1)進(jìn)行變換可得到

采用這個(gè)計(jì)算模型可以最大程度地降低大氣氣溶膠、土壤背景等因素對(duì)計(jì)算植被覆蓋度的影響，只要把Sveg和Ssoil這兩個(gè)參數(shù)確定了，就可以通過(guò)式(2)求得植被覆蓋度。

2.2.2 歸一化植被指數(shù)模型

利用歸一化植被指數(shù)(normalized difference vegetation index，NDVI)進(jìn)行植被覆蓋度的估算，不僅方法簡(jiǎn)單易行，而且有利于對(duì)區(qū)域內(nèi)植被的總體狀況做定量評(píng)估［25］，因此該指數(shù)常用來(lái)反映植被生長(zhǎng)狀況、植被覆蓋、生物量等信息，是反映生態(tài)環(huán)境狀況的重要指標(biāo)［26］。計(jì)算公式為

式中:NIR為近紅外波段像素值;RED為可見(jiàn)光紅波段像素值;NDVI為歸一化植被指數(shù)。

將式(3)代入式(2)，即得到以下基于NDVI的植被覆蓋度像元二分模型，即

式中:NDVIveg為完全植被覆蓋地表所貢獻(xiàn)的信息;NDVIsoil為無(wú)植被覆蓋地表所貢獻(xiàn)的信息。

2.2.3 閾值的確定

NDVIsoil在理論上是不隨時(shí)間和空間改變的，然而由于受大氣影響和地表濕度條件的改變，以及粗糙度、土壤類(lèi)型、土壤顏色等條件的不同，NDVIsoil也會(huì)隨之發(fā)生變化，變化范圍一般在 －0.1～0.2之間［27］。應(yīng)用NDVI閾值監(jiān)測(cè)植被的限制在于區(qū)分不同的植被覆蓋類(lèi)型需要設(shè)置不同的閾值，這是因?yàn)榇蟛糠种脖桓采w類(lèi)型是不同植物類(lèi)型的混合體［28］，所以不能采用固定的 NDVIsoil和 NDVIveg值。本文根據(jù)NDVI數(shù)據(jù)頻率統(tǒng)計(jì)表，計(jì)算NDVI數(shù)據(jù)的頻率累積值，取累積頻率為2%的 NDVI值為NDVIsoil，取累積頻率為98%的 NDVI值為NDVIveg。由于TM和HJ－1A/1B CCD數(shù)據(jù)的空間分辨率均為30 m，并且NDVI值范圍相近，因此兩期數(shù)據(jù)采用的閾值相同。

2.2.4 精度驗(yàn)證與變化分析方法

1)精度驗(yàn)證方法。利用2009年122個(gè)植被覆蓋度野外調(diào)查樣點(diǎn)對(duì)2009年植被覆蓋度遙感監(jiān)測(cè)結(jié)果進(jìn)行精度驗(yàn)證，具體方法是將每個(gè)樣點(diǎn)根據(jù)葉面積指數(shù)儀LAI2000和魚(yú)眼鏡頭觀測(cè)獲得的100 m×100 m空間分辨率的平均植被覆蓋度，與相對(duì)應(yīng)的3像元×3像元遙感監(jiān)測(cè)植被覆蓋度進(jìn)行匹配，獲得植被覆蓋度的觀測(cè)值(即真實(shí)值)與遙感測(cè)量值之間的決策系數(shù)R2，以驗(yàn)證遙感估算精度。

2)變化分析方法。通過(guò)植被覆蓋度逐像元相減的方法，即2009年像元植被覆蓋度減去2000年相對(duì)應(yīng)像元植被覆蓋度，就可得到2000—2009年植被覆蓋度變化。同時(shí)分別對(duì)研究區(qū)不同植被類(lèi)型、不同縣(市)的植被覆蓋度變化結(jié)果進(jìn)行了分析。

3 結(jié)果與分析

3.1 植被覆蓋度估算結(jié)果

將研究區(qū)劃分為植被覆蓋區(qū)和非植被覆蓋區(qū)，利用式(3)分別獲得2000年與2009年南水北調(diào)中線水源區(qū)的NDVI數(shù)據(jù)。然后根據(jù)該數(shù)據(jù)以及相應(yīng)的NDVIsoil和NDVIveg閾值，通過(guò)式(4)得到2000年與2009年的植被覆蓋度遙感監(jiān)測(cè)結(jié)果(圖2)。

圖2 南水北調(diào)中線水源區(qū)植被覆蓋度遙感監(jiān)測(cè)結(jié)果Fig.2 Remote sensing monitoring results of fractional vegetation cover in water source area of the middle route of projects to divert water from the south to the north

2000年和2009年南水北調(diào)中線水源區(qū)植被覆蓋度的平均值分別為67.5%和72%。由2期監(jiān)測(cè)結(jié)果可見(jiàn):植被覆蓋度的平均值整體呈增長(zhǎng)趨勢(shì)，其增幅的空間特征表現(xiàn)為水源區(qū)中部地區(qū)高，東、西部相對(duì)較低，中、西部高山地帶增幅較顯著，水源區(qū)中低植被覆蓋度區(qū)域的面積明顯減少。

3.2 精度驗(yàn)證

野外調(diào)查結(jié)果表明，南水北調(diào)中線水源區(qū)內(nèi)，自然植被占85.9%，植被群落主要以落葉闊葉/針葉混交林為主。圖3是根據(jù)研究區(qū)122個(gè)植被覆蓋度的實(shí)測(cè)數(shù)據(jù)(即真實(shí)值)與遙感估算數(shù)據(jù)得到的散點(diǎn)圖。2種數(shù)據(jù)之間相關(guān)性的決策系數(shù)R2為0.816 1。

圖3 2009年植被覆蓋度遙感監(jiān)測(cè)結(jié)果精度驗(yàn)證Fig.3 Verified the accuracy of remote sensing monitoring results of fractional vegetation cover in 2009

結(jié)果表明，應(yīng)用像元二分模型估算植被覆蓋度具有一定的準(zhǔn)確性，雖然不同的研究區(qū)需要重新判定模型的輸入閾值，但模型在區(qū)域尺度監(jiān)測(cè)上的確有很好的適用性。

3.3 植被覆蓋度變化分析

3.3.1 不同植被類(lèi)型覆蓋度比較

圖4為2000年與2009年水源區(qū)不同植被類(lèi)型的覆蓋度對(duì)比。

圖4 2000年與2009年不同植被類(lèi)型覆蓋度對(duì)比Fig.4 Comparing the fractional cover of different vegetation types in 2000 and 2009

圖中顯示，草地、灌木、混交林、常綠闊葉林、常綠針葉林、落葉闊葉林及落葉針葉林的植被覆蓋度均有不同程度的增加，這與該區(qū)域以原始自然植被為主、人類(lèi)活動(dòng)干擾得到較好控制的實(shí)際情況相符。水源區(qū)植被覆蓋度增加的原因，除了自然因素以外，國(guó)家采取的植被保護(hù)政策和相關(guān)措施在一定程度上改善了植被的生長(zhǎng)條件則是重要因素。落葉針葉林平均植被覆蓋度增長(zhǎng)最快，達(dá)15.64%，增幅為25.49%。這是由于亞熱帶、熱帶落葉針葉林，如巴山冷杉、太白冷杉等分布在秦嶺南坡中段與西部海拔1 700～1 900 m的高山地帶，人為干擾因素很少的緣故。水源區(qū)2009年草地平均植被覆蓋度較2000年增加了4.44%。這是由于草地多分布于河谷及山區(qū)，特別是坡度較大、土壤較薄的山坡上，因此自然生長(zhǎng)面積的增幅相對(duì)較小。

3.3.2 各縣(市)植被覆蓋度比較

為了多方位地掌握該區(qū)植被覆蓋度的動(dòng)態(tài)變化情況，本文對(duì)位于中心區(qū)域的37個(gè)縣(市)2000年和2009年的植被覆蓋度變化情況進(jìn)行了分析，制作了“2000年和2009年植被覆蓋度平均值對(duì)比圖”和“植被覆蓋度變化分析圖”(圖5，6)。

圖5 2000年與2009年各縣(市)植被覆蓋度平均值對(duì)比Fig.5 Comparing the average of fractional vegetation cover of various towns in 2000 and 2009

圖6 2000—2009年植被覆蓋度變化分析Fig.6 The analysis of fractional vegetation cover changes in 2000 and 2009

由圖5和圖6得出，水源區(qū)內(nèi)植被覆蓋度整體呈增長(zhǎng)趨勢(shì)，其中柞水縣平均植被覆蓋度增幅最大;作為丹江口水庫(kù)的水利樞紐和渠首丹江口市和淅川縣的植被覆蓋度也有不同程度的增加。以下詳細(xì)分析這4個(gè)縣(市)的具體情況。

在水源區(qū)37個(gè)縣市中，柞水縣的平均植被覆蓋度增幅最大(圖7(a))，由59.72%增長(zhǎng)到77.48%。雖然柞水縣地處秦嶺南麓，是一個(gè)“九山半水半分田”的土石山區(qū)，但自2000年實(shí)施退耕還林政策以來(lái)，該縣已經(jīng)退耕造林28萬(wàn)畝??1畝≈0.066 7 hm2，目前森林復(fù)蓋率達(dá)66%［29］，植被已得到較好的恢復(fù)。

圖7 植被覆蓋度變化遙感監(jiān)測(cè)結(jié)果Fig.7 Monitoring results of fractional vegetation cover changes by remote sensing

洋縣2009年平均植被覆蓋度與2000年相比整體變化不大(圖7(b))。洋縣地表多為低山丘陵，地表土質(zhì)松散裸露，耕地面積較大，是植被覆蓋度變化不明顯的因素。區(qū)水生植被較多，基本位于丹江口水庫(kù)庫(kù)尾地帶;自2002年退耕還林等生態(tài)工程建設(shè)以來(lái)，水源區(qū)的植被生長(zhǎng)得到了有效保護(hù)。從實(shí)地調(diào)查也可以看出，通過(guò)工程與生物措施相結(jié)合，沿丹江口水庫(kù)兩岸形成了較高的植被覆蓋度，不同的造林模式，也增加了庫(kù)區(qū)生態(tài)林涵養(yǎng)水源和保護(hù)水土的功能。

2009年淅川縣的植被覆蓋度雖比2000年有所提高，但平均植被覆蓋度較其他36個(gè)縣市來(lái)說(shuō)還相對(duì)較低。淅川縣位于南水北調(diào)中線水源區(qū)渠首，在水源區(qū)內(nèi)的面積為2 407.15 km2，其中永久水面和季節(jié)性水面 302.47 km2，耕地 710.89 km2，水體和耕地共占該縣水源區(qū)內(nèi)面積的42.1%，由于水面及耕地所占比例較大，加上其他建設(shè)用地的使用，均使得該縣植被覆蓋度處于較低水平。

綜上所述，雖然南水北調(diào)中線水源區(qū)植被覆蓋度總體上逐年呈增長(zhǎng)趨勢(shì)，但該區(qū)植被覆蓋的整體格局仍保持不變，林地植被覆蓋度雖保持穩(wěn)定增加，但隨著水源區(qū)人口的不斷增加以及城鎮(zhèn)的擴(kuò)張，水源區(qū)水域河岸帶的植被將受到人為因素的干擾，部分地區(qū)的植被覆蓋度可能會(huì)有降低的趨勢(shì)。

4 結(jié)論

本文基于像元分解法和植被指數(shù)法，利用2000年生長(zhǎng)季的TM圖像和2009年生長(zhǎng)季的HJ－1A/1B圖像估算了南水北調(diào)中線水源區(qū)的植被覆蓋度，并分析其時(shí)空變化特征，主要結(jié)果如下:

1)利用像元二分模型估算南水北調(diào)中線水源區(qū)的植被覆蓋度，總體精度達(dá)到81.61%，證明用該方法進(jìn)行水源區(qū)植被覆蓋度的估算是可行的。

2)2000年水源區(qū)的平均植被覆蓋度為67.5%，2009年達(dá)到72%，植被覆蓋度總體呈增高趨勢(shì)，其增幅的空間特征表現(xiàn)為該區(qū)中部地區(qū)高，東部和西部地區(qū)相對(duì)較低。

3)在不同的植被類(lèi)型中，落葉針葉林植被覆蓋度的增長(zhǎng)幅度最大，增幅為25.49%，草地覆蓋度的增長(zhǎng)幅度最小，平均只增長(zhǎng)了4.44%。位于水源區(qū)的大多數(shù)縣(市)的植被覆蓋度在近十年內(nèi)都有不同程度的提高，其中柞水縣的平均植被覆蓋度的增長(zhǎng)幅度最大。植被覆蓋度的普遍增加與國(guó)家實(shí)施的退耕還林還草、荒山造林、封山育林、基本農(nóng)田建設(shè)等政策有密切關(guān)系。

本文的不足之處在于沒(méi)有對(duì)非密度模型加以討論，也未確定植被的消光系數(shù)和葉面積指數(shù)等重要的植被特征參量。此外，對(duì)綜合條件相近的混合像元區(qū)域取相同的參數(shù)，也是需要改進(jìn)的方面。今后還可進(jìn)一步開(kāi)展2000—2009年長(zhǎng)時(shí)間序列植被覆蓋度的月動(dòng)態(tài)監(jiān)測(cè)，對(duì)植被覆蓋度和累積生物量等生態(tài)環(huán)境效應(yīng)參數(shù)進(jìn)行動(dòng)態(tài)變化分析，以及結(jié)合降雨量進(jìn)行相關(guān)關(guān)系分析等工作。

［1］ Purevdorj T S，Tateishi R，Ishiyama T，et al.Relationships Between Percent Vegetation Cover and Vegetation Indices［J］.International Journal of Remote Sensing，1998，19(18):3519 －3535.

［2］ Gitelson A A，Kaufman Y J，Stark R，et al.Novel Algorithms for Remote Estimation of Vegetation Fraction［J］.Remote Sensing of Environment，2002，80(1):76 －87.

［3］ 周?chē)?guó)林，袁正科.常用林業(yè)技術(shù)術(shù)語(yǔ)［M］.長(zhǎng)沙:湖南科學(xué)技術(shù)出版社，1982.

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

［5］ 陳云浩，李曉兵，史培軍，等.北京海淀區(qū)植被覆蓋的遙感動(dòng)態(tài)研究［J］.植物生態(tài)學(xué)報(bào)，2001，25(5):588 －593.

［6］ 張?jiān)葡?，李曉兵，陳云?草地植被蓋度的多尺度遙感與實(shí)地測(cè)量方法綜述［J］.地球科學(xué)進(jìn)展，2003，18(1):85 －93.

［7］ 吳素業(yè).一年生植被覆蓋度的簡(jiǎn)化測(cè)試方法［J］.水土保持科技情報(bào)，1999(1):45 －47.

［8］ 唐世浩，朱啟疆，周宇宇，等.一種簡(jiǎn)單的估算植被覆蓋度和恢復(fù)背景信息的方法［J］.中國(guó)圖象圖形學(xué)報(bào)，2003，8(11):1304－1308.

［9］ 李苗苗，吳炳方，顏長(zhǎng)珍，等.密云水庫(kù)上游植被覆蓋度的遙感估算［J］.資源科學(xué)，2004，26(4):153 －159.

［10］程紅芳，章文波，陳 鋒.植被覆蓋度遙感估算方法研究進(jìn)展［J］.國(guó)土資源遙感，2008(1):13 －18.

［11］李 攀，胡德勇，趙文吉.北方農(nóng)牧交錯(cuò)帶植被覆蓋變化遙感監(jiān)測(cè)研究——以河北省沽源縣為例［J］.國(guó)土資源遙感，2010(2):113－117.

［12］ Qi J，Marsett R C，Moran M S，et al.Spatial and Temporal Dynamics of Vegetation in the San Pedro River Basin Area［J］.Agricultural and Forest Meteorology，2000，105(1 －3):55 －68.

［13］李善同，許新宜.南水北調(diào)與中國(guó)發(fā)展［M］.北京:經(jīng)濟(jì)科學(xué)出版社，2004.

［14］李思悅，張全發(fā).南水北調(diào)中線丹江口庫(kù)區(qū)主要生態(tài)環(huán)境問(wèn)題及植被恢復(fù)［J］.中國(guó)農(nóng)村水利水電，2008(3):1－4.

［15］朱桂香.南水北調(diào)中線河南水源區(qū)的水土流失與治理對(duì)策［J］.中國(guó)水土保持，2005(8):30 －31.

［16］金蓉玲，郭海晉.南水北調(diào)中線水源區(qū)水資源評(píng)價(jià)及丹江口入庫(kù)水量分析［J］.人民長(zhǎng)江，1993，24(11):7 －12.

［17］沈大軍，劉昌明.南水北調(diào)中線工程不同調(diào)水規(guī)模對(duì)漢江中下游影響分析［J］.地理學(xué)報(bào)，1998，53(4):341 －348.

［18］沈澤昊，張全發(fā)，岳 超，等.南水北調(diào)中線水源區(qū)土地利用/土地覆被的空間格局［J］.地理學(xué)報(bào)，2006，61(6):633 －644.

［19］李石華，角媛梅.環(huán)境與災(zāi)害監(jiān)測(cè)預(yù)報(bào)小衛(wèi)星A星CCD影像質(zhì)量評(píng)價(jià)［J］.紅外技術(shù)，2009，31(9):167 －172.

［20］林 輝，劉泰龍，李際平.遙感技術(shù)基礎(chǔ)教程［M］.長(zhǎng)沙:中南大學(xué)出版社，2002.

［21］ Leprieur C，Verstraete M M，Pinty B.Evaluation of the Performance of Various Vegetation Indices to Retrieve Vegetation Cover from AVHRR Data［J］.Remote Sensing Review，1994，10(4):265 －284.

［22］陳 晉，陳云浩，何春陽(yáng)，等.基于土地覆蓋分類(lèi)的植被覆蓋率估算亞像元模型與應(yīng)用［J］.遙感學(xué)報(bào)，2001，5(6):416 －422.

［23］ Zribi M，Le H S，Taconet O，et al.Derivation of Wild Vegetation Cover Density in Semi－ arid Regions:ERS2/SAR Evaluation［J］.International Journal of Remote Sensing，2003，24(6):1335 －1352.

［24］吳 云，曾 源，吳炳方，等.基于MODIS數(shù)據(jù)的三北防護(hù)林工程區(qū)植被覆蓋度提取與分析［J］.生態(tài)學(xué)雜志，2009，28(9):1712－1718.

［25］吳 云，曾 源，趙 炎，等.基于MODIS數(shù)據(jù)的海河流域植被覆蓋度估算及動(dòng)態(tài)變化分析［J］.資源科學(xué)，2010，32(7):1417－1424.

［26］郭芬芬，范建容，嚴(yán) 冬，等.基于像元二分模型的昌都縣植被蓋度遙感估算［J］.中國(guó)水土保持，2010(5):65－67.

［27］ Carlson T N，Ripley D A.On the Relation Between NDVI，F(xiàn)ractional Vegetation Cover，and Leaf Area Index［J］.Remote Sensing of Environment，1997，62(3):241 －252.

［28］王 宏，李曉兵，鶯 歌，等.基于NOAA NDVI的植被生長(zhǎng)季模擬方法研究［J］.地理科學(xué)進(jìn)展，2006，25(6):21 －32.

［29］陜西柞水縣鞏固退耕還林成果項(xiàng)目對(duì)我省的啟示［EB/OL］.2010 －07 －07.http://plan.hainan.gov.cn/pwnview.php?id=1250.

Remote Sensing Monitoring and Analysis of Fractional Vegetation Cover in the Water Source Area of the Middle Route of Projects to Divert Water from the South to the North

ZHOU Zhi－ qiang1，2，ZENG Yuan1，ZHANG Lei1，DU Xin1，WU Bing － fang1
(1.Institute of Remote Sensing Applications，Chinese Academy of Sciences，Beijing 100101，China;2.Southwest Forestry University，Kunming 650224，China)

The middle route of the projects to divert water from the south to the north is a part of the large－scale inter－basin water transfer projects in China.It is important to carry out the study and analysis of changes of regional fractional vegetation cover for the protection of ecological environment and water quality.In this paper，based on the data of the remote sensing images obtained in 2000 and 2009，the authors estimated fractional vegetation cover of the water source area by using the method of dimidiate pixel model from normalized difference vegetation index(NDVI)，and analyzed the temporal and spatial variation characteristics of the fractional vegetation cover.The main conclusions are as follows:The average fractional vegetation cover in the water source area was 67.5%in 2000，and reached 72%in 2009，with the average fractional vegetation cover being increased in the whole area.The spatial characteristics of increased fractional vegetation cover show that the increase in central region is relatively higher than that in the eastern and western regions;in different kinds of vegetation types，the deciduous conifer forest shows the largest average increase of the fractional vegetation cover，while the grassland shows the smallest increase;fractional vegetation cover has increased in different degrees in most towns of the water source area in the past 10 years，with the increase in Zhashui County being most apparent.This is attributed to the Chinese Government’s policies such as quitting cultivation and returning to forest，closing hillside to facilitate afforestation and farmland construction.

fractional vegetation cover;water source area of the middle route of projects to divert water from the south to the north;remote sensing monitoring

TP 79;X 835

A

1001－070X(2012)01－0070－07

10.6046/gtzyyg.2012.01.13

2011－05－12;

2011－06－08

中國(guó)科學(xué)院戰(zhàn)略性先導(dǎo)科技專(zhuān)項(xiàng)(編號(hào):XDA05050108)和國(guó)務(wù)院南水北調(diào)項(xiàng)目“利用遙感技術(shù)對(duì)南水北調(diào)中線水源區(qū)生態(tài)環(huán)境變化的監(jiān)測(cè)研究”共同資助。

周志強(qiáng)(1985－)，男，碩士研究生，主要從事生態(tài)遙感研究。E－mail:upgis@foxmail.com。

曾 源(1979－)，女，博士，副研究員，主要從事植被定量遙感研究。E－mail:yuanz@irsa.ac.cn。

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