龔文峰,吳　娟,王笑峰,劉　濤,王鑫鑫,劉芳平

(黑龍江大學(xué) 水利電力學(xué)院，哈爾濱 150086)

?

基于RS與GIS密云區(qū)植被蓋度動(dòng)態(tài)演變及其預(yù)測(cè)研究

龔文峰,吳娟,王笑峰,劉濤,王鑫鑫,劉芳平

(黑龍江大學(xué) 水利電力學(xué)院，哈爾濱 150086)

基于RS和GIS技術(shù)，以密云區(qū)2001年和2010年兩期Landsat TM5遙感影像為主要數(shù)據(jù)源，在定量反演研究區(qū)歸一化植被指數(shù)(NDVI)基礎(chǔ)上，研究該區(qū)域植被覆蓋度并分析其時(shí)空演變，并基于馬爾科夫模型完成其植被蓋度動(dòng)態(tài)預(yù)測(cè)。結(jié)果表明：密云區(qū)整體植被狀況較好，中高覆蓋和高蓋度植被是研究區(qū)植被覆蓋的主體，占總面積89.39%，近10 a密云區(qū)植被蓋度總體呈上升趨勢(shì)，高蓋度面積增加448.87 km2，其它類(lèi)型的植被蓋度面積減少，中高蓋度植被面積減少最大，為401.01 km2，且植被蓋度重心總體向高覆蓋度移動(dòng)。基于馬爾科夫模型的預(yù)測(cè)發(fā)到2091年，各植被類(lèi)型的面積比例達(dá)到相對(duì)穩(wěn)定狀態(tài)，分別為0.48%、1.75%、2.68%、6.94%和88.15%，高蓋度植被占據(jù)絕對(duì)的優(yōu)勢(shì)地位。

RS；GIS；密云區(qū)；植被蓋度；馬爾科夫模型

植被是陸地生態(tài)系統(tǒng)的主要組成部分,也是其他生物生存的基礎(chǔ)，具有截流降雨、減緩徑流、防沙治沙、保持水土等功能[1]。植被蓋度(Vegetation coverage)通常是指植被冠層垂直投影面積占基準(zhǔn)地表單位面積的比例或百分?jǐn)?shù)[2]，是刻畫(huà)陸地表面植被數(shù)量的一個(gè)重要參數(shù)[3]。作為重要的生態(tài)氣候參數(shù)，植被覆蓋率是許多全球及區(qū)域氣候數(shù)值模型中所需的重要信息，也是描述生態(tài)系統(tǒng)的重要基礎(chǔ)數(shù)據(jù)[4]。因此，掌握植被覆蓋度的變化規(guī)律，對(duì)評(píng)價(jià)人類(lèi)生存環(huán)境質(zhì)量、調(diào)節(jié)生態(tài)過(guò)程具有重要的實(shí)際意義[5-7]。

然而傳統(tǒng)的估算植被蓋度方法主要是地面測(cè)量，包括目估法、采樣法、儀器法等[8-9]，這些方法不但費(fèi)時(shí)費(fèi)力，而且很難給出較大尺度上的宏觀植被分布信息[10]，但隨著空間信息技術(shù)飛速發(fā)展，尤其遙感技術(shù)的迅猛發(fā)展，監(jiān)測(cè)區(qū)域面積大、周期時(shí)間短和獲取植被信息快速準(zhǔn)確等這些遙感技術(shù)特征被廣泛用于城市植被監(jiān)測(cè)研究[11-12]。近年來(lái)，利用遙感的手段進(jìn)行植被蓋度動(dòng)態(tài)監(jiān)測(cè)日趨成熟，并開(kāi)展了基于遙感植被蓋度動(dòng)態(tài)監(jiān)測(cè)研究[13-16]，但基于Markov模型預(yù)測(cè)未來(lái)植被蓋度變化的研究相對(duì)較少。

密云區(qū)境內(nèi)的密云水庫(kù)是北京市最大的飲用水源基地，它保障著首都1/ 3 的供水需求。然而隨著城市化進(jìn)程的加速和經(jīng)濟(jì)的快速發(fā)展，近年來(lái)人們對(duì)其周邊土地不合理的開(kāi)發(fā)利用以及采礦、筑路、盲目擴(kuò)大旅游業(yè)等，導(dǎo)致了庫(kù)區(qū)周?chē)?，尤其是?kù)東北、北部地區(qū)的水土流失嚴(yán)重[17]。密云區(qū)作為首都重要飲用水源地和生態(tài)涵養(yǎng)發(fā)展區(qū)，研究該區(qū)域植被蓋度及其變化信息，對(duì)于評(píng)價(jià)區(qū)域生態(tài)環(huán)境等具有重要現(xiàn)實(shí)意義[18-19]。基于此，本文在RS和GIS技術(shù)的支持下，以密云區(qū)2001年和2010年Landsat TM遙感影像為重要數(shù)據(jù)源，定量反演歸一化差異植被指數(shù)(Normalized Difference Vegetation Index -NDVI)，并計(jì)算分析植被蓋度時(shí)空演變規(guī)律，并基于Markov模型預(yù)測(cè)探究未來(lái)植被蓋度演變趨勢(shì)，旨在為該研究區(qū)域生態(tài)環(huán)境保護(hù)和土地利用結(jié)構(gòu)優(yōu)化提供技術(shù)參考。

1　研究區(qū)域概況

密云區(qū)位于北京市東北部，地理位置為E:116°39′33″-117°30′25″，N:40°13′7″-40°47′57″，地處華北平原與蒙古高原的過(guò)渡地帶，屬燕山山地與華北平原交接地，密云區(qū)屬于山區(qū)，東、西、北三面環(huán)山，中部低緩，西南開(kāi)闊，地勢(shì)自北向西南傾斜，呈簸箕形。密云區(qū)屬暖溫帶半濕潤(rùn)大陸性季風(fēng)氣候區(qū)，年平均氣溫10.8 ℃，年太陽(yáng)總輻射量約為135.0 kcal/cm2，年平均降水量為661.3 mm，平均降水日數(shù)為75 d，無(wú)霜期平均為177 d。

2　數(shù)據(jù)處理與研究方法

2.1數(shù)據(jù)來(lái)源及處理

以研究區(qū)1∶50 000的地形圖作為基礎(chǔ)地理參照數(shù)據(jù)，基于遙感處理軟件ENVI 5.1支持下，對(duì)2001年和2010年9月的Landsat TM遙感影像進(jìn)行幾何精校正和自動(dòng)匹配處理，保證其幾何精度，且影像投影方式為橫軸墨卡托投影UTM(Zone48N)。在此基礎(chǔ)上，完成遙感影像的大氣校正、圖像掩膜裁剪和去云及陰影等預(yù)處理。

2.2植被覆蓋提取

歸一化植被指數(shù)(NDVI)的計(jì)算公式為：

(1)

式中CH3為T(mén)M影像第三通道的反射值；CH4為T(mén)M影像第四通道的反射值[22]。利用植被蓋度與歸一化植被指數(shù)之間關(guān)系估算植被蓋度，如下式：

(2)

式中NDVI為所求像元的歸一化植被指數(shù)；NDVImin、NDVImax分別為非植被覆蓋部分(裸地和未利用地)和植被覆蓋部分(林地、草地等)歸一化植被指數(shù)值的最小值和最大值。

2.3植被覆蓋的分級(jí)

根據(jù)《土壤侵蝕分類(lèi)分級(jí)標(biāo)準(zhǔn)》[20]，本文將植被蓋度劃分為5個(gè)等級(jí)，Ⅰ級(jí)(極低植被蓋度水域)：植被覆蓋度0～0.3，主要為水域；Ⅱ級(jí)(低蓋度植被區(qū)域)：植被覆蓋度0.3～0.45，主要是裸土、低產(chǎn)草地、灌木林地、撂荒地和居民用地等；Ⅲ級(jí)(中蓋度植被區(qū)域)：植被覆蓋度0.45～0.6，主要是低產(chǎn)草地、灌木林地和撂荒地；Ⅳ級(jí)(高蓋度植被區(qū)域)：植被覆蓋度0.6～0.75，主要是中高產(chǎn)草地、林地和部分農(nóng)田等，一般呈片狀或塊狀分布；Ⅴ級(jí)(高蓋度植被區(qū)域)：植被覆蓋度0.75～1，主要是密林地和優(yōu)良耕地等，呈現(xiàn)片狀或塊狀分布。在ArcGIS 9.3支持下，結(jié)合一類(lèi)森林清查資料及高分辨率Google Earth圖片，對(duì)分類(lèi)的結(jié)果就進(jìn)行了精度驗(yàn)證。

2.4馬爾科夫模型

馬爾科夫過(guò)程是一種遵循“無(wú)后效性原則”的特殊隨機(jī)運(yùn)動(dòng)過(guò)程，將時(shí)間序列看作一個(gè)隨機(jī)過(guò)程，能夠定量地揭示各植被蓋度類(lèi)型之間的相互轉(zhuǎn)化狀況，進(jìn)而預(yù)測(cè)植被蓋度類(lèi)型隨時(shí)空的變化情況。該模型關(guān)鍵是確定植被蓋度類(lèi)型之間轉(zhuǎn)化初始轉(zhuǎn)移概率矩陣，其數(shù)學(xué)表達(dá)式為[21-23]:

(3)

式中n為研究區(qū)植被蓋度類(lèi)型的數(shù)目;Pij為初期到末期時(shí)段內(nèi)由類(lèi)型i轉(zhuǎn)為類(lèi)型j的轉(zhuǎn)移概率，且滿足以下條件:

(4)

對(duì)植被蓋度類(lèi)型預(yù)測(cè)的數(shù)學(xué)模型為:

(5)

式中P(k)為預(yù)測(cè)對(duì)象在k時(shí)刻的狀態(tài)向量，P為轉(zhuǎn)移矩陣概率，P(k+1)表示預(yù)測(cè)對(duì)象在k+1時(shí)刻的狀態(tài)向量。

3　結(jié)果分析

3.1植被蓋度空間分布

基于RS和GIS獲取研究區(qū)域2001年和2010 年植被蓋度圖(圖1)。由圖1可見(jiàn)，研究區(qū)植被覆蓋的主體是中高覆蓋植被和高蓋度植被，符合研究區(qū)域的土地利用結(jié)構(gòu)特征，密云區(qū)的土地利用主要以林地和耕地為主，且森林覆蓋率約70%[24-25]。2001年，密云區(qū)低覆蓋度分布于密云西南部密云鎮(zhèn)、十里堡、河南寨等居民用地，中蓋度主要是耕地，分布在低覆蓋度居民用地周?chē)?，中高蓋度分布于密云區(qū)西北部古北口、新城子以及密云水庫(kù)周?chē)呱w度植被主要呈大小不一的片狀分布于西北部馮家峪、不老屯、石城，東部北莊和大城子，南部東邵渠。2010年，密云區(qū)極低蓋度(水域)主要是密云水庫(kù)以及河流。低蓋度主要是居民用地，主要分布于密云鎮(zhèn)、十里堡、河南寨，還有一小部分零星散落在各鄉(xiāng)鎮(zhèn)內(nèi)。中蓋度主要為耕地，主要呈包圍狀分布于低蓋度居民用地周?chē)?，中高蓋度面積較小，零星散落在高蓋度植被區(qū)域內(nèi)。高蓋度植被呈片狀分布于密云區(qū)北部山區(qū)、東部山區(qū)以及密云區(qū)南部東邵渠。密云區(qū)城區(qū)內(nèi)的植被蓋度由低蓋度范圍擴(kuò)大，說(shuō)明密云區(qū)城區(qū)植被蓋度受到城鎮(zhèn)化的推進(jìn)的影響。密云區(qū)西北部和東部山區(qū)植被覆蓋度均大幅度提升，植被蓋度重心向高蓋度植被轉(zhuǎn)移。

圖1　密云區(qū)植被蓋度空間分布圖Fig.1　Spatial distribution of vegetation coverage in Miyun District

3.2植被蓋度動(dòng)態(tài)變化

植被覆蓋變化總的態(tài)勢(shì)見(jiàn)表1。由表1可見(jiàn)，研究區(qū)的中植被覆蓋的面積變化高覆蓋和高蓋度植被主要分布于西北部和東部，源于該區(qū)域?qū)儆谏絽^(qū)，受自然環(huán)境條件的影響，人為干擾相對(duì)較少，分布著大量的森林，故其植被蓋度相對(duì)較高，兩期的中高覆蓋和高蓋度植被面積之和分別為1 941.96 km2和1 989.83 km2，且所占的面積比例分別為87.23%和89.39%，在研究區(qū)域中占據(jù)處于絕對(duì)的優(yōu)勢(shì)地位，構(gòu)成該區(qū)域的主體植被蓋度，表明密云區(qū)的植被覆蓋較好且得到很大程度的保護(hù)，這與密云水庫(kù)作為北京市重要的飲用水源地有直接的關(guān)聯(lián)。高蓋度植被面積由2001年的1 284.92 km2增加到2010年的1 733.79 km2，增加448.87 km2，而其它類(lèi)型的植被蓋度面積都呈現(xiàn)一定程度的減少，其中中高蓋度減少幅度最大，由2001年657.04 km2減少到2010年256.04 km2，減少401 km2，所占比例由29.51%降至11.50%。其中，中蓋度和中高蓋度植減少的區(qū)域主要集中分布于密云水庫(kù)周?chē)兔茉茀^(qū)東北部山區(qū)以及潮河流域、白河流域，且中蓋度和中高蓋度植被主要轉(zhuǎn)換為高蓋度植被，主要源于“退耕還林”和“退耕還草”等生態(tài)工程的實(shí)施，致使該區(qū)域的植被蓋度提升，說(shuō)明區(qū)域生態(tài)環(huán)境質(zhì)量有所好轉(zhuǎn)，表明區(qū)域的天然林保護(hù)等生態(tài)工程已見(jiàn)成效。極低植被蓋度減少的面積為28.65 km2，減少區(qū)域主要為密云水庫(kù)面積，源于10 a間密云水庫(kù)入庫(kù)水量有所減少，水庫(kù)水位下降了17 m，水庫(kù)蓄水不及庫(kù)容總量的1/3?？傮w而言，植被覆蓋度年均變化量較大，高蓋度植被區(qū)域的增加量最大，中高蓋度植被區(qū)域的減少量最大，密云區(qū)植被蓋度總體呈上升趨勢(shì)，植被蓋度重心向高覆蓋度移動(dòng)。

3.3植被覆蓋模擬分析

3.3.1植被蓋度轉(zhuǎn)移概率

基于Arc GIS 9.3支持下，完成兩期植被蓋度數(shù)據(jù)的空間疊加分析，獲取其變化信息和屬性特征，構(gòu)建轉(zhuǎn)移概率矩陣(表2)，定量分析不同類(lèi)型植被蓋度間的轉(zhuǎn)移頻數(shù)，進(jìn)而探究不同植被蓋度之間的相互轉(zhuǎn)化特征。

表1　密云區(qū)植被覆蓋度面積統(tǒng)計(jì)

表2　密云區(qū)2001-2010年植被蓋度變化轉(zhuǎn)移矩陣

由表2可見(jiàn)，高蓋度和極低蓋度具有較高的保留率，分別為93.36%和68.63%，中高蓋度是高蓋度植被的最大轉(zhuǎn)入來(lái)源，轉(zhuǎn)移面積為489.79 km2，其次為中蓋度植被，其轉(zhuǎn)入率為19.52%，該結(jié)果表明隨著北京市“十二五”時(shí)期環(huán)境保護(hù)和建設(shè)規(guī)劃的實(shí)施，通過(guò)人工育林、封山育林、退耕還林還草和治理水土流失等生態(tài)措施，森林覆蓋率得到了提高，在一定程度上表明區(qū)域的生態(tài)環(huán)境質(zhì)量有所改善。極低蓋度主體是密云水庫(kù)，由于2000年以來(lái)水庫(kù)水位不斷下降，部分土地逐漸露出水面，庫(kù)區(qū)周邊的農(nóng)民開(kāi)始利用部分區(qū)域的土地進(jìn)行農(nóng)業(yè)生產(chǎn)，并形成了優(yōu)質(zhì)的耕地，致使16.59 km2的極低蓋度轉(zhuǎn)換為高蓋度植被，這種局部不合理的土地利用轉(zhuǎn)化對(duì)于密云水庫(kù)極為不利，今后要加大“退耕還林”和“退耕還草”的力度，促使區(qū)域土地利用的合理化。低蓋度和中蓋度的保留率分別為46.25%和27.82%，中蓋度植被是中高蓋度和高蓋度的主要轉(zhuǎn)入來(lái)源，轉(zhuǎn)移面積分別為43.75 km2和21.70 km2，轉(zhuǎn)入中高蓋度植被的面積最大，說(shuō)明植被蓋度由低等級(jí)向高等級(jí)轉(zhuǎn)變，在一定程度上表明該區(qū)域的植被蓋度得到一定程度的恢復(fù)，區(qū)域生態(tài)環(huán)境質(zhì)量得到了提升。此外，中蓋度也是低蓋度植被的主要轉(zhuǎn)入來(lái)源，轉(zhuǎn)移面積為14.02 km2，隨著密云區(qū)城鎮(zhèn)化的快速發(fā)展，人類(lèi)加大對(duì)中蓋度植被的干擾程度，導(dǎo)致其部分農(nóng)田被建成學(xué)校、道路、醫(yī)院等居民用地建筑[26],造成區(qū)域生態(tài)環(huán)境質(zhì)量有所惡化。低蓋度植被的轉(zhuǎn)入率和轉(zhuǎn)出率相當(dāng)，故其面積沒(méi)有發(fā)生顯著變化。中高蓋度保留率最低，僅為20.65%，74.56%。中高蓋度轉(zhuǎn)化為高蓋度，對(duì)高蓋度植被增加的貢獻(xiàn)最大，這種轉(zhuǎn)化主要發(fā)生在密云東北山區(qū)、潮河流域以及白河流域，源于該區(qū)域積極實(shí)施“退耕還林，退耕還草”等林業(yè)生態(tài)措施。

3.3.2植被蓋度動(dòng)態(tài)預(yù)測(cè)

以2010年各植被蓋度類(lèi)型的面積百分比為起始向量，基于2001-2010年時(shí)間段轉(zhuǎn)移概率構(gòu)建初始狀態(tài)矩陣(表2)，在MATLAB軟件的支持下，以9 a為時(shí)間間隔為步長(zhǎng)，運(yùn)用馬爾柯夫鏈模型，對(duì)密云區(qū)未來(lái)植被覆蓋狀態(tài)進(jìn)行動(dòng)態(tài)預(yù)測(cè)，預(yù)測(cè)2019年以后直到相對(duì)穩(wěn)定平衡狀態(tài)時(shí)研究區(qū)植被覆蓋類(lèi)型的變化趨勢(shì)。

2010年后的植被覆蓋變化趨勢(shì)是高蓋度以下的植被類(lèi)型面積均呈現(xiàn)減少的趨勢(shì)，而高蓋度面積呈現(xiàn)穩(wěn)步增加的變化趨勢(shì)，見(jiàn)表3。假若當(dāng)前的自然條件和國(guó)家生態(tài)建設(shè)政策不出現(xiàn)太大的變動(dòng)的情況下，研究區(qū)域各植被蓋度類(lèi)型在2091年達(dá)到相對(duì)平衡的穩(wěn)定狀態(tài)，即之后若干年該區(qū)域各植被蓋度比例保持相對(duì)穩(wěn)定，此時(shí)各植被覆蓋類(lèi)型比例分別為0.48%、1.75%、2.68%、6.94%和88.15%?？梢?jiàn)高蓋度所占比例最大，完全處于絕對(duì)的優(yōu)勢(shì)地位，說(shuō)明密云區(qū)在水源保護(hù)林建設(shè)工程、京津風(fēng)沙源治理工程和天然林保護(hù)等生態(tài)保護(hù)措施下，使得高蓋度面積得以穩(wěn)步增加，最終的預(yù)測(cè)結(jié)果不僅表明了該區(qū)域高蓋度植被的未來(lái)發(fā)展趨勢(shì)，更預(yù)測(cè)其對(duì)區(qū)域生態(tài)環(huán)境的動(dòng)態(tài)影響作用，通過(guò)發(fā)揮森林植被的水土保持和水源涵養(yǎng)功能來(lái)改善北京市供水的質(zhì)和量。極低蓋度下降速度最快，根據(jù)預(yù)測(cè)趨勢(shì)，密云水庫(kù)水域面積達(dá)到穩(wěn)定狀態(tài)時(shí)只有10.76 km2，從2010年到2091年面積減少44.68 km2，面積減少的最大。如果照此預(yù)測(cè)結(jié)果趨勢(shì)發(fā)展下去，將直接威脅到北京市的用水安全。然而近些年來(lái)，隨著北京經(jīng)濟(jì)的高速發(fā)展，人口數(shù)量的不斷膨脹，城市需水量的急劇增加，為解決日益突出的用水問(wèn)題，水庫(kù)流域內(nèi)對(duì)地下水進(jìn)行大量開(kāi)采，直接導(dǎo)致研究區(qū)域的地下水位不斷下降，密云水庫(kù)面積大幅縮小，庫(kù)容減少，直至2014年，南水北調(diào)中線工程正式通水后，每年補(bǔ)水密云水庫(kù)(6～10)×108m3，密云水庫(kù)的蓄水量逐漸遞增，水庫(kù)面積減小的情況得以改觀，北京市水資源短缺問(wèn)題將得以緩解。

表3　密云區(qū)各類(lèi)植被蓋度變化趨勢(shì)預(yù)測(cè)

4　結(jié)　論

基于RS和GIS支持下，以密云區(qū)作為研究對(duì)象，借助歸一化植被指數(shù)(NDVI)獲取植被覆蓋度狀況，基于ArcGIS9.3定量分析研究區(qū)域植被蓋度時(shí)空演變，并基于Markov模型預(yù)測(cè)探究未來(lái)植被蓋度演變趨勢(shì)，旨在為該研究區(qū)域?yàn)樗亮魇Х乐魏椭匾幢Ｗo(hù)提供理論依據(jù)。結(jié)果表明：中高覆蓋植被和高蓋度植被占據(jù)密云區(qū)的主體，面積比例分別為87.23%和89.39%。高蓋度面積增加，為448.87 km2，而其它類(lèi)型植被蓋度面積均有不同程度減少，中高蓋度植被面積減少最大，為401.01 km2。密云區(qū)植被蓋度總體呈上升趨勢(shì)，植被蓋度重心向高覆蓋度移動(dòng)。

高蓋度和極低蓋度具有較高的保留率，分別為93.36%和68.63%；中高蓋度保留率最低，僅為20.65%，74.56%中高蓋度轉(zhuǎn)化為高蓋度；從2010年模擬到穩(wěn)定平衡狀態(tài)時(shí)，各植被覆蓋類(lèi)型比例分別為0.48%、1.75%、2.68%、6.94%和88.15%，高蓋度所占比例最高，極低蓋度下降速度最快。

致謝：中國(guó)西部環(huán)境與生態(tài)科學(xué)數(shù)據(jù)中心提供部分精度驗(yàn)證數(shù)據(jù)，謹(jǐn)表謝意。

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

[2]Purevdorj T, 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.

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

[4]Sellers P J, S O Los, C J Tucker, et al. A revised land surface paramctcrization(Sib2)for atmospheric gems. Part II: The generation of global fields of terrestrial biophysical parameters from satellite data[J].JoarnaL of Climate, 1996, 9: 706-737.

[5]龐吉林,張克斌,喬娜,等.基于RS的鹽池縣10年植被覆蓋度動(dòng)態(tài)變化研究[J].水土保持研究,2012(4):112-116.

[6]胡江波.不同植被恢復(fù)模式對(duì)土壤肥力的影響[J].河南農(nóng)業(yè)科學(xué),2007, 3:69-72.

[7]王國(guó)芳.基于遙感技術(shù)的植被覆蓋度的動(dòng)態(tài)監(jiān)測(cè)[J].山西農(nóng)業(yè)科學(xué),2015,43(5):592-595.

[8]馬娜,胡云鋒.基于遙感和像元二分模型的內(nèi)蒙古正藍(lán)旗植被覆蓋度格局和動(dòng)態(tài)變化[J].地理科學(xué),2012,32(2): 251- 256.

[9]趙英時(shí).遙感應(yīng)用分析原理與方法[M].北京:科學(xué)出版社,2013.

[10] 史世蓮,章文波,王國(guó)燕.北京地區(qū)植被蓋度提取及其分布變化研究[J].遙感技術(shù)與應(yīng)用,2014,29(5):866-872.

[11] Weng Q,Lu D.Schubring J Estimation of land surface temperature-vegetation abundance relationship for urban heat island studies[J].Remote Sensing of Environment,2004,89(4):467-483.

[12] Yuan F,Bauer M E.Comparison of impervious suface area and nomalized difference vegetation index as indicators of surface urban heat island effects in Landsat imagery[J].Remote Sensing of Environment,2007,106(3):375-386.

[13] 萬(wàn)紅梅,李霞,董道瑞.基于多源遙感數(shù)據(jù)的荒漠植被覆蓋度估測(cè)[J].應(yīng)用生態(tài)學(xué)報(bào),2012,23(12):3331-3337.

[14] 宋挺,黃君.基于LandsatETM+等影像的無(wú)錫市十年陸地生態(tài)環(huán)境遙感評(píng)價(jià)[J].環(huán)境監(jiān)測(cè)管理與技術(shù),2014, 26(5):28-32.

[15] 王曉峰,任志遠(yuǎn).基于RS和GIS的榆林市植被覆蓋度動(dòng)態(tài)變化研究[J].陜西師范大學(xué)學(xué)報(bào),2008(3):101-104.

[16] 劉亞玲,潘志華,范錦龍.陰山北麓地區(qū)植被覆蓋動(dòng)態(tài)時(shí)空分析[J].資源科學(xué),2005,27(4):168-174.

[17] 張有全,宮輝力,趙文吉,等.基于GIS和USLE的密云縣土壤侵蝕評(píng)價(jià)及空間特征研究[J].水土保持研究,2007(3):358-362.

[18] 馬超飛,馬建文,布和放斯?fàn)?USLE模型中植被覆蓋因子的遙感數(shù)據(jù)定量估算[J].水土保持學(xué)報(bào),2001,21(4):6-9.[19] 高清竹,段敏杰,李玉娥,等.1981—2006年藏西北地區(qū)草地植被蓋度動(dòng)態(tài)變化分析[J].中國(guó)農(nóng)業(yè)氣象,2010,31(4):582-585.

[20] GB/SL190-1996,土壤侵蝕分類(lèi)分級(jí)標(biāo)準(zhǔn)[S].

[21] 鄔建國(guó).景觀生態(tài)學(xué)-格局、過(guò)程、尺度與等級(jí)[M].北京:高等教育出版社，2000.

[22] 史培軍,宮鵬,李曉兵,等.土地利用/覆蓋變化研究的方法與實(shí)踐[M].北京:科學(xué)出版社,2000.

[23] 張葉,江曉波,邱楓.LUCC模型研究綜述[J].資源開(kāi)發(fā)與市場(chǎng),2006,22(4):311-314.

[24] 盧玉東,尹黎明,何丙輝,等.利用TM影像在土地利用/覆蓋遙感解澤波段選取研究[J].西南農(nóng)業(yè)大學(xué)學(xué)報(bào)：自然科學(xué)版,2005,27(4):479-486.

[25] 何春陽(yáng),陳晉,陳石浩,等.土地利用/覆蓋變化混合動(dòng)態(tài)監(jiān)測(cè)方法研究[J].自然資源學(xué)報(bào),2001,16(3):255-262.

[26] 陳濤,牛瑞卿,李平湘,等.密云水庫(kù)流域植被覆蓋度變化對(duì)輸沙量的影響[J].生態(tài)環(huán)境學(xué)報(bào),2010,19(1): 152-159.

Vegetation coverage dynamic evolution and its forecastbased on RS and GIS in Miyun District

GONG Wen-Feng,WU Juan,WANG Xiao-Feng,LIU Tao,WANG Xin-Xin，LIU Fang-Ping

(SchoolofHydraulicandElectric-power,HeilongjiangUniversity,Harbin150080,China)

Remote sensing images of Landsat TM5 in 2001 and 2010 were used as the main data source of Miyun District, Beijing city of China. Based on RS and GIS, normalized differential vegetation index (NDVI) was quantitatively inversed from the remote sensing data, which was used to calculate the vegetation coverage, and then the dynamic evolution and spatial and temporal variation characteristics of vegetation coverage were analyzed in this paper. On this basis, Markov model was used for the quantitative forecast of the future vegetation coverage of the study region. The results were shown that: the vegetation in Miyun District is in good condition. The main types of vegetation accounting for 89.39% of the total area are mid-high, high vegetation coverage. The increasing area of high vegetation coverage was 448.87 km2, while areas for other types of vegetation coverage were decreased, among which the mid-high vegetation coverage has the great decreasing area of 401.01 km2.The whole vegetation in study area in 10 years are increased and the main focus moved to high vegetation coverage. The vegetation coverage tendency after 2010 is that expect for high vegetation , other kinds of vegetation has decreased year by year. This trend will last for a longtime until it reaches a relatively stable state, with the proportion 0.48%,1.75%,2.68%,6.94% and 88.15% for all types vegetation coverage, respectively, meanwhile, the high vegetation coverage occupied the absolute dominate position.

RS;GIS; Miyun District; vegetation coverage; Markov model

10.13524/j.2095-008x.2016.03.037

2016-04-21

國(guó)家高技術(shù)研究發(fā)展計(jì)劃項(xiàng)目(2012AA102001);黑龍江省自然基金資助項(xiàng)目(D201410); 黑龍江省教育廳科學(xué)技術(shù)研究項(xiàng)目(1253153)

龔文峰(1976-),男，河南南陽(yáng)人,教授,博士,研究方向：RS和GIS在資源與環(huán)境監(jiān)測(cè)上的應(yīng)用，E-mail:gwf101@163.com。

X87

A

2095-008X(2016)03-0025-07

網(wǎng)絡(luò)出版地址：http:www.cnki.net/kcms/detail/23.1566.T.20160719.1132.006.html