李 輝，周啟剛，焦 歡，王福海，國洪磊

（1.重慶工商大學(xué) 旅游與國土資源學(xué)院，重慶400067；2.重慶工商大學(xué) 融智學(xué)院，重慶400033）

植被覆蓋是全球及區(qū)域水文氣候模型模擬中所需要的重要信息，也是描述生態(tài)系統(tǒng)的重要基礎(chǔ)數(shù)據(jù)［1］。在城市化高速發(fā)展的現(xiàn)代社會(huì)，獲取地表植被覆蓋及其變化信息，對(duì)于揭示地表空間變化規(guī)律，探討變化的驅(qū)動(dòng)因子，研究區(qū)域生態(tài)環(huán)境以及優(yōu)化城市生態(tài)環(huán)境都具有重要意義［2］。

隨著“3S”技術(shù)的快速發(fā)展，通過遙感技術(shù)反演植被覆蓋逐漸成為城市生態(tài)環(huán)境監(jiān)測的重要手段［3－5］。目前國內(nèi)外有關(guān)城市植被覆蓋變化的遙感應(yīng)用研究方興未艾，國內(nèi)學(xué)者已利用不同傳感器、不同時(shí)相、不同分辨率的衛(wèi)星影像對(duì)一些城市的植被覆蓋情況進(jìn)行了專題研究?？梢钥闯?，此類研究大多集中在北京、上海、廣州等東部發(fā)達(dá)城市，對(duì)西部城市，特別是山地城市的相關(guān)研究相對(duì)欠缺［6－7］。

本文利用 ArcGIS 10.2和erdas imagine 2010軟件對(duì)重慶市都市區(qū)1988年、2000年、2002年以及2012年的4期LandsatTM.遙感數(shù)據(jù)進(jìn)行處理，提取植被覆蓋指數(shù)，監(jiān)測并分析了其植被空間覆蓋與動(dòng)態(tài)變化情況。

1 研究區(qū)概況

重慶市位于四川盆地的東南部，長江的上游，與湖北、湖南、貴州、四川、陜西等省接壤，地跨東經(jīng)105°11′—110°11′、北緯28°10′—32°13′的四川盆地與長江中下游平原的過渡地帶。轄區(qū)總面積為8.24萬平方km，2013年重慶市都市區(qū)人口已達(dá)808萬人，屬于中國面積最大的內(nèi)陸城市。有長江和嘉陵江兩江流經(jīng)區(qū)內(nèi)，氣候濕潤溫和。其都市區(qū)主要包括重慶都市九區(qū)，即：渝中區(qū)、大渡口區(qū)、江北區(qū)、南岸區(qū)、沙坪壩區(qū)、九龍坡區(qū)、北碚區(qū)、渝北區(qū)、巴南區(qū)；主要分布在長江沿線，四面環(huán)山，以丘陵，低山為主，地形起伏較大，平均高程為400m。研究區(qū)總面積為5 465 km2，占全市總面積的6.63%，截止2012年常駐人口為772.31萬人。研究區(qū)從南北面向長江嘉陵江河谷傾斜，起伏較大，呈現(xiàn)“一山一槽二嶺”的自然景觀，是典型的特大山地城市［8］。

重慶市作為新中國最年輕的直轄市，是西南地區(qū)經(jīng)濟(jì)發(fā)展的中心之一，對(duì)于帶動(dòng)西南片區(qū)經(jīng)濟(jì)飛速發(fā)展具有重要的意義。近年來，重慶市經(jīng)濟(jì)發(fā)展迅猛，城市建設(shè)用地?cái)U(kuò)展進(jìn)一步加快，城市生態(tài)環(huán)境問題日益突出。選取重慶山地都市區(qū)作為研究對(duì)象，對(duì)于城市建設(shè)用地?cái)U(kuò)展與地表特征參數(shù)變化研究具有理論與實(shí)踐意義。

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

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

本研究應(yīng)用的主要數(shù)據(jù)為研究區(qū)的數(shù)字高程圖和遙感數(shù)據(jù)。其中遙感數(shù)據(jù)為LandsatTM/ETM影像，及獲取年份依次為1988年、2000年、2002年、2012年。所有影像數(shù)據(jù)采用Albers投影，中央經(jīng)線為東經(jīng)105°，雙標(biāo)準(zhǔn)緯線分別為北緯25°和47°。LandsatTM 影像數(shù)據(jù)1，2，3，4，5，7波段的空間分辨率均為30m，6波段空間分辨率為120m。

2.2 研究方法

2.2.1 NDVI的表達(dá)式：

式中：p（nir）——TM4（近紅外）波段的亮度值；p（red）——TM3（可見光紅光）波段的亮度值。

對(duì)于Landsat TM數(shù)據(jù)而言，歸一化植被指數(shù)通常可表示為：

式中：R3，R4——TM影像波段3和波段4的反射率。

利用式（3）可得到1988—2000年和2000—2012年研究區(qū)歸一化植被指數(shù)變化分布圖。

式中：NDVIC——前后兩個(gè)目標(biāo)年份之間的NDVI變化值；NDVIn——第n年份研究區(qū)的NDVI值；NDVIn－1——前一個(gè)年份研究區(qū)的NDVI值。

將上式的計(jì)算結(jié)果進(jìn)行標(biāo)準(zhǔn)化分級(jí)，分別化分為負(fù)向變化、正向變化和無顯著變化［9－10］：

式中：n——像元數(shù)量，根據(jù)公式（4），（5）和（6），可以對(duì)重慶都市區(qū)不同研究時(shí)段內(nèi)的植被覆蓋情況進(jìn)行總體分析。

2.2.2 NDVI空間分布重心遷移模型 分析NDVI的重心變化能具體反映出NDVI在研究年限內(nèi)的空間分布變化，其研究的計(jì)算模型如下［9－10］：

式中：Xt，Yt——第t年 NDVI的重心坐標(biāo)；cti——第i塊NDVI斑塊的面積；Xi，Yi——第i塊NDVI斑塊的重心坐標(biāo)。

3 結(jié)果與分析

3.1 重慶市都市區(qū)NDVI時(shí)間變化分析

3.1.1 重慶市都市區(qū)NDVI總體分析 重慶市都市區(qū)NDVI在時(shí)間序列上呈現(xiàn)出明顯的差異。從1988—2012年，NDVI指數(shù)較低的區(qū)域面積逐漸擴(kuò)大，特別是渝北區(qū)、南岸區(qū)、九龍坡區(qū)、北碚區(qū)和巴南區(qū)增速最快。這主要是由于重慶直轄以來，這幾個(gè)區(qū)的經(jīng)濟(jì)發(fā)展迅速，建設(shè)用地面積急劇增加，城市綠化的速度跟不上城市建設(shè)的速度，所以導(dǎo)致NVDI指數(shù)降低。此外，NDVI指數(shù)較高的區(qū)域基本保持穩(wěn)定，只有局部地區(qū)產(chǎn)生變化。

通過數(shù)據(jù)分析可以得到1988年、2000年、2002年以及2012年研究區(qū)平均NDVI分別為0.2，0.13，0.32，0.3；由此可知，1988—2012年重慶市都市區(qū)平均NDVI水平變化十分明顯，都市區(qū)平均NDVI分布在時(shí)間上呈現(xiàn)出明顯的階段性。從1988—2002年，平均NDVI值變化十分明顯，特別是渝北區(qū)、南岸區(qū)、九龍坡區(qū)、北碚區(qū)和巴南區(qū)NDVI變化最快。這主要是由于直轄以來，這幾個(gè)區(qū)的經(jīng)濟(jì)發(fā)展迅速，建設(shè)用地面積急劇增加，城市綠化的速度跟不上建設(shè)用地?cái)U(kuò)張的速度，導(dǎo)致低NVDI區(qū)域面積不斷增加；同時(shí)在城市發(fā)展的過程中，由于生態(tài)觀念的加強(qiáng)，都市區(qū)對(duì)山林的管制力度也大幅增加。從2002—2012年，可以看出都市區(qū)平均NDVI區(qū)域基本保持穩(wěn)定，只有局部區(qū)域產(chǎn)生變化。

3.1.2 重慶市都市區(qū)分區(qū)NDVI變化分析 利用ArcGIS 10.2以及Erdas 8.7軟件的相關(guān)功能，可以準(zhǔn)確的得到1988—2012年重慶都市各區(qū)的NDVI平均值，再利用相關(guān)統(tǒng)計(jì)軟件即可得到重慶市都市各區(qū)在研究年限內(nèi)的變化趨勢圖。

從圖1可以明顯地看出，重慶市都市各區(qū)在1988—2012年平均NDVI的變化趨勢，總體而言，重慶都市各區(qū)在研究年限內(nèi)平均NDVI變化幅度均比較明顯；相比之下，沙坪壩區(qū)和渝北區(qū)的變化比較小。各區(qū)從1988—2012年平均NDVI呈現(xiàn)上下浮動(dòng)趨勢，其中渝中區(qū)、南岸區(qū)、大渡口區(qū)以及江北區(qū)變化較大，而渝中區(qū)平均NDVI值最低，且下降速度較快。九龍坡區(qū)、渝北區(qū)以及巴南區(qū)的平均NDVI指數(shù)始終保持較高水平，說明這幾個(gè)區(qū)的生態(tài)環(huán)境較好；另一方面，渝中區(qū)以及南岸區(qū)在2012年平均NDVI指數(shù)有所上升，可以分析出幾個(gè)區(qū)的NDVI呈逐步恢復(fù)的狀態(tài)。

圖1 1988－2012年重慶市都市各區(qū)平均NDVI變化趨勢

此外，由于重慶市都市區(qū)各區(qū)所處的地理位置不同，比如有些區(qū)處于地勢高區(qū)域，有些處于低地勢范圍。其次，有些區(qū)的山脈較多，有些區(qū)地勢平坦。因此，研究海拔與坡度對(duì)各區(qū)NDVI變化的影響具有重要意義。

3.1.3 重慶市都市區(qū)NDVI分級(jí)分析 為了更加科學(xué)和直觀的對(duì)比重慶市都市區(qū)1988年、2000年、2002年以及2012年的NDVI數(shù)據(jù)以及變化趨勢；在ArcGIS 10.2軟件的支持下，文章以自然斷點(diǎn)分類法將重慶市都市NDVI指數(shù)分為5個(gè)等級(jí)，分別是第一級(jí)（低）、第二級(jí)（較低）、第三級(jí)（中）、第四級(jí)（較高）、第五級(jí)（高），以便更深層次的揭示出重慶市都市區(qū)NDVI變化規(guī)律，具體分級(jí)數(shù)據(jù)見表1。

表1 1988－2012年都市區(qū)不同NDVI級(jí)別所占面積與比例

從表1中可以看出，重慶市都市區(qū)1988—2012年都市區(qū)不同NDVI級(jí)別所占面積與比例呈現(xiàn)一定的規(guī)律性；都市區(qū)NDVI第一級(jí)（低）水平面積比例呈現(xiàn)先增加后減少的趨勢，第二級(jí)（較低）水平面積比例總體呈現(xiàn)上升趨勢。另一方面，都市區(qū)NDVI第三級(jí)（中）水平、第四級(jí)（較高）水平以及第五級(jí)（高）水平所占比例在1988—2012年均較大，但是，隨著時(shí)間的推移，這三個(gè)水平的NDVI所占面積和比例呈現(xiàn)出上下浮動(dòng)的變化趨勢，尤其是第四等級(jí)與第五等級(jí)在總體上呈現(xiàn)下降趨勢。通過查閱重慶市都市區(qū)城市自然概括以及城市發(fā)展數(shù)據(jù)不難發(fā)現(xiàn)，由于重慶市都市區(qū)處于山區(qū)，研究區(qū)內(nèi)包含四條天然山脈，且植被覆蓋度保持良好，而重慶市城市的快速發(fā)展和建設(shè)用地的不斷擴(kuò)展以及相關(guān)政策的實(shí)施也影響到不同NDVI水平在區(qū)域范圍的不斷變化。

從表2中，我們可以明顯看出，重慶市都市區(qū)在研究年限內(nèi)，平均NDVI指數(shù)均發(fā)生了明顯變化，通過分析可以發(fā)現(xiàn)，這種變化呈現(xiàn)出一定的規(guī)律性。

此外，都市區(qū)NDVI第三級(jí)（中）水平、第四級(jí)（較高）水平以及第五級(jí)（高）水平均變化不大；這是由于都市所處地理位置的原因，都市區(qū)的四條山脈的植被覆蓋比較完好。此外，無論是1988—2000年，2000—2002年，還是2002—2012年；都市區(qū)NDVI第一級(jí)（低）水平和NDVI第二級(jí)（較低）水平在研究年限內(nèi)，均發(fā)生了十分顯著的變化。其中，都市區(qū)NDVI第一級(jí)（低）水平是變化最明顯的等級(jí)，在ArcGIS的技術(shù)支持下，可以分析出，NDVI第一級(jí)（低）水平所代表的區(qū)域類型主要是建設(shè)用地和河流；在重慶市城市化水平的逐年提高，特別是重慶直轄以來，城市迅速發(fā)展以及建設(shè)用地大規(guī)模擴(kuò)張的背景下，導(dǎo)致研究區(qū)NDVI第一級(jí)（低）水平變化十分明顯。

總體而言，隨著時(shí)間的推移，研究區(qū)不同級(jí)別NDVI指數(shù)所占的面積和比例不斷地發(fā)生著變化；因此，關(guān)注區(qū)域NDVI指數(shù)的變化有利于城市生態(tài)環(huán)境的保護(hù)和改善。

表2 1988－2012年都市區(qū)不同NDVI級(jí)別轉(zhuǎn)化面積及速率

表3 1988－2002年都市區(qū)不同NDVI級(jí)別面積轉(zhuǎn)移矩陣

表4 2002－2012年都市區(qū)不同NDVI級(jí)別面積轉(zhuǎn)移矩陣

通過表3，可以十分清晰的看出，1988—2002年重慶市都市區(qū)不同NDVI級(jí)別相互間轉(zhuǎn)移的情況。從總體來看，不同NDVI級(jí)別逐級(jí)轉(zhuǎn)移的情況比較明顯，而且轉(zhuǎn)移面積和轉(zhuǎn)移比例都比較大；跨級(jí)轉(zhuǎn)移的情況雖然存在，但是所轉(zhuǎn)移的面積和比例都相對(duì)較少。因此，1988—2002年重慶市都市區(qū)不同NDVI級(jí)別相互間轉(zhuǎn)移呈現(xiàn)出一定的規(guī)律性，研究區(qū)NDVI整體水平有所提升。

從上述都市區(qū)NDVI轉(zhuǎn)移的情況，可以大致分析產(chǎn)生這種變化的主要原因，在1988—2002年間，由于生態(tài)環(huán)境保護(hù)觀念的不斷深入，重慶市加強(qiáng)了對(duì)都市區(qū)山脈的開發(fā)管制與綠地保護(hù)，此外，不斷提升的都市區(qū)的綠化水平也是一個(gè)重要原因。

結(jié)合 表3和表4，我們可以明顯看出，一方面，從2002—2012年的總體情況而言，研究區(qū)NDVI指數(shù)變化不大，另一方面，在2002—2012年間，重慶市都市區(qū)各等級(jí)轉(zhuǎn)移情況對(duì)比1988—2002年較明顯。

通過分析可以看出，2002—2012年重慶市都市區(qū)各等級(jí)NDVI較顯著地相互間轉(zhuǎn)移現(xiàn)象，其主要體現(xiàn)在高NDVI指數(shù)向較高NDVI指數(shù)轉(zhuǎn)移上，二者相互間轉(zhuǎn)移的比例甚至達(dá)到了90%；其次也體現(xiàn)在較高NDVI指數(shù)向中NDVI指數(shù)的轉(zhuǎn)移上。這說明，在研究年限內(nèi)，研究區(qū)高NDVI區(qū)域發(fā)生了較大變化，造成這種變化的原因是復(fù)雜多樣的；通過分析，我們可以發(fā)現(xiàn)，造成這種變化的主要原因還是由于重慶市自直轄以來，快速的城市化和建設(shè)用地的大規(guī)模擴(kuò)張。研究區(qū)在1988年，2000年，2011年城市建設(shè)用地面積分別達(dá)到17 145hm2，28 603hm2，80 396hm2。

3.2 重慶市都市區(qū)NDVI空間變化分析

利用 NDVI計(jì)算公式，在 ArcGIS 10.2以及Erdas 8.7軟件的支持下，對(duì)重慶市都市區(qū)1988年、2000年、2002年和2012年的NDVI進(jìn)行反演，獲得了都市區(qū)四年的NDVI分布圖（附圖9—10）。從分布圖中可以看出，都市區(qū)NDVI在不同年份的空間分布呈現(xiàn)出明顯的差異性，這種差異隨時(shí)間而改變而呈現(xiàn)出一定規(guī)律性。

通過都市區(qū)NDVI分布圖可知，研究區(qū)歸一化植被指數(shù)分布在空間上呈明顯的差異性，這種差異隨時(shí)間的變化呈現(xiàn)出一定的規(guī)律性，且這種規(guī)律性隨著時(shí)間的不斷推移，越發(fā)明顯，城市建設(shè)區(qū)的不斷擴(kuò)展使得大片區(qū)域的NDVI指數(shù)降低，另一方面，由于“四山管制”規(guī)劃的實(shí)施，在一定程度上緩解了都市區(qū)NDVI指數(shù)的快速下降。

3.2.1 NDVI水平空間變化分析 從NDVI指數(shù)的水平空間分布情況來看，研究區(qū)NDVI指數(shù)總體呈現(xiàn)降低。此外，低NDVI指數(shù)的區(qū)域除分布在城市建成區(qū)外，還集中在長江與嘉陵江流域；另一方面，山地丘陵區(qū)的植被覆蓋率與較平坦地區(qū)的植被覆蓋率呈明顯的差異，呈現(xiàn)出明顯的垂直地帶性。在山地丘陵地區(qū)，植被覆蓋率保持相對(duì)穩(wěn)定，且海拔越高，植被覆蓋率越高。這主要是由兩個(gè)方面的原因造成的，一方面是由于高海拔地區(qū)地勢陡峭，土地不利于開發(fā)和利用；另一方面，都市區(qū)各大山脈都劃入了限制建設(shè)區(qū)或禁止建設(shè)區(qū)，使得這些區(qū)域的植被保護(hù)較好。地勢平緩區(qū)域?yàn)橹貞c市農(nóng)田和園地的集中分布區(qū)，也是研究年限內(nèi)建設(shè)用地?cái)U(kuò)張的主要區(qū)域，也是植被覆蓋率變化最大的區(qū)域，這也說明建設(shè)用地對(duì)農(nóng)田和園地的大規(guī)模替代能顯著改變地表植被的覆蓋狀況。

在ArcGIS 10.2的支持下，利用處理好的研究區(qū)數(shù)據(jù)，得到了1988—2012年重慶市都市區(qū)各NDVI等級(jí)重心坐標(biāo)轉(zhuǎn)移趨勢圖以及重心轉(zhuǎn)移示意圖（圖2）。

從圖2，我們可以很直觀的看出1988—2012年都市區(qū)各等級(jí)NDVI重心轉(zhuǎn)移坐標(biāo)變化趨勢，這種變化趨勢呈現(xiàn)出一定的等級(jí)和年份規(guī)律性?？傮w上來說，各等級(jí)在研究年限內(nèi)重心變化相對(duì)較小。

研究區(qū)NDVI指數(shù)第一等級(jí)，這個(gè)等級(jí)在研究區(qū)5個(gè)NDVI等級(jí)中處于最低水平，主要代表的土地利用類型為建設(shè)用地、裸地以及河流水面；其集中在重慶市都市區(qū)兩江交匯處。因此，通過觀察該NDVI等級(jí)重心的變化，也可以在一定程度上間接反映出研究區(qū)建設(shè)用地的擴(kuò)張方向。

圖2 1988－2012年都市區(qū)各等級(jí)NDVI重心轉(zhuǎn)移坐標(biāo)變化趨勢

從圖2可以看出，在1988—2012年間，研究區(qū)NDVI指數(shù)第一等級(jí)的重心坐標(biāo)轉(zhuǎn)移并不明顯，其轉(zhuǎn)移范圍主要限定在狹小的區(qū)域。研究區(qū)NDVI指數(shù)第二等級(jí)僅此于NDVI指數(shù)第一等級(jí)，其重心變化也不明顯，但較NDVI指數(shù)第一等級(jí)重心變化來說相對(duì)較大；說明在1988—2012年間，NDVI指數(shù)第二等級(jí)所占面積和區(qū)域有一定程度的變化以及轉(zhuǎn)移。

研究區(qū)NDVI指數(shù)第三等級(jí)以及第四等級(jí)所處范圍的重心位置變化不明顯，但是也其重心位置也經(jīng)歷了“北遷”和“南移”的過程；但是從該兩個(gè)等級(jí)重心的總體轉(zhuǎn)移情況來看，其變化范圍不大，并且也限定在了一定范圍之內(nèi)。

研究區(qū)NDVI第五等級(jí)是NDVI指數(shù)最高的等級(jí)，其主要代表的土地利用類型是森林和耕地；從圖2我們可以明顯的看出，該等級(jí)重心位置變化相當(dāng)明顯，其重心位置由最初的“南移”變?yōu)樽詈蟮摹氨边w”，變化范圍比較大。通過分析可知，研究區(qū)NDVI第五等級(jí)重心變化如此明顯，可以說明該等級(jí)NDVI指數(shù)所處的范圍在研究年限內(nèi)發(fā)生了顯著地變化。這種變化的主要原因是由于都市區(qū)高植被覆蓋區(qū)不斷變化引起的。

3.2.2 NDVI垂直空間變化分析 在垂直空間分布上，研究區(qū)NDVI分布呈現(xiàn)出明顯的垂直地帶性（表5）。

在山地丘陵地區(qū)，NDVI指數(shù)較高，且保持相對(duì)穩(wěn)定狀態(tài)，隨著海拔的升高，NDVI也越高。通過分析可知，造成這種情況的原因主要是由兩個(gè)方面的原因造成的，一方面是由于高海拔地區(qū)地勢陡峭，一般不利于土地的開發(fā)和利用；另一方面，由于重慶市將都市區(qū)各大山脈都劃入了限制建設(shè)區(qū)或禁止建設(shè)區(qū)，使植被得以較好的保護(hù)。地勢平緩區(qū)域是人類活動(dòng)的主要區(qū)域，也是植被覆蓋率變化最大的區(qū)域，說明建設(shè)用地?cái)U(kuò)展對(duì)農(nóng)田和園地的大規(guī)模替代能顯著改變區(qū)域NDVI狀況。通過對(duì)研究區(qū)不同海拔地區(qū)與NDVI分布情況之間相互關(guān)系的分析可以看出，不同級(jí)別NDVI的分布情況與海拔差異存在一定的聯(lián)系。在低海拔地區(qū)，隨著時(shí)間的推移，較低NDVI在低海拔區(qū)域保持相對(duì)穩(wěn)定，較高NDVI面積比例在低海拔區(qū)域逐年遞減；在中海拔地區(qū)，由于該區(qū)是城市建設(shè)的主要分布區(qū)域，較高NDVI有所較少，較低NDVI有所增加；在高海拔地區(qū)，在環(huán)境保護(hù)意識(shí)的加強(qiáng)和新能源燃料的普及背景下，較高NDVI面積比例有所提升。

總的來說，研究區(qū)NDVI指數(shù)不僅受高程因素的影響，也在人們環(huán)保意識(shí)與相關(guān)政策的變化下產(chǎn)生差異性。因此，在保護(hù)好現(xiàn)有高植被覆蓋度的同時(shí)，也要大力加強(qiáng)對(duì)低海拔區(qū)域植被覆蓋度的恢復(fù)與提升。

表5 不同年份不同海拔區(qū)域各NDVI等級(jí)分布面積比例 %

4 結(jié) 論

本文基于RS和GIS技術(shù)，對(duì)重慶市都市區(qū)1988—2012年間的NDVI在數(shù)量、結(jié)構(gòu)、空間上進(jìn)行分析，得到結(jié)論如下：

（1）在研究年限內(nèi)，無論是在時(shí)間還是空間上，研究區(qū)的NDVI指數(shù)都發(fā)生了十分顯著的變化，且這種變化呈現(xiàn)出一定的規(guī)律性。

（2）在時(shí)間變化上，研究區(qū)低NDVI指數(shù)和較低NDVI指數(shù)區(qū)域面積在不斷增加，而較高NDVI指數(shù)與高NDVI指數(shù)區(qū)域面積有所減少；尤其是在2002—2012年間，這種現(xiàn)象更加突出。

（3）在空間變化上，通過分析不同NDVI指數(shù)等級(jí)重心遷移軌跡，可以看出各等級(jí)的重心位置雖在不斷變化，但從整個(gè)研究時(shí)相內(nèi)，各等級(jí)的重心位置偏移位置相對(duì)較小，呈現(xiàn)出比較穩(wěn)定的狀態(tài)。

［1］ 陳穎彪，吳志峰，程炯，等.遙感與GIS支持下的城市綠地提取方法研究：以深圳市為例［J］.生態(tài)環(huán)境，2004，13（3）：362－364.

［2］ 丁艷梅，張繼賢，王堅(jiān)，等.基于TM數(shù)據(jù)的植被覆蓋度反演［J］.測繪科學(xué)，2006，31（1）：43－45.

［3］ 龔建周，夏北成.基于遙感影像的廣州市植被覆蓋度內(nèi)部結(jié)構(gòu)與時(shí)空變化［J］.植物資源與環(huán)境學(xué)報(bào)，2006，15（4）：25－29.

［4］ 韓貴鋒，徐建華，袁興中，等.1998—2001年重慶市都市區(qū)植被的時(shí)空變化［J］.生態(tài)學(xué)，2007，26（9）：1412－1417.

［5］ 田野，徐恒力，劉雙，等.基于歸一化植被指數(shù)變化分級(jí)的千山植被變化［J］.中國水土保持，2010（1）：25－28.

［6］ 樸世龍，方精云.最近18年來中國植被覆蓋的動(dòng)態(tài)變化［J］.第四紀(jì)研究，2001，21（4）：294－302.

［7］ 馬志勇，沈濤，張軍海，等.基于植被覆蓋度的植被變化分析［J］.測繪通報(bào)，2007（3）：45－48.

［8］ 李建輝，周啟剛，陳倩，等.基于RS和GIS的建設(shè)用地變化研究：以重慶市都市九區(qū)為例［J］.重慶工商大學(xué)學(xué)報(bào)，2012，29（1）：76－81.

［9］ 湯國安，楊昕.地理信息系統(tǒng)空間分析試驗(yàn)教程［M］.北京：科學(xué)出版社，2006.

［10］ 李苗苗.植被覆蓋度的遙感估算方法研究［D］.北京：中國科學(xué)研究院遙感應(yīng)用研究所，2003.

［11］ 王福海，馬澤忠，向敬華，等.基于TM的重慶山地都市區(qū)建設(shè)用地?cái)U(kuò)展與地表特征參數(shù)時(shí)空變化研究［J］.水土保持研究，2014，21（4）：107－113.