空間羅倫茲曲線在土地利用及其變化*研究中的應(yīng)用——以西藏拉薩市城關(guān)區(qū)為例

劉敏紅,何政偉,3,姜蘭,仇文俠,張船紅

(1.地質(zhì)災(zāi)害防治與地質(zhì)環(huán)境保護(hù)國(guó)家重點(diǎn)實(shí)驗(yàn)室,成都 610059;2.成都理工大學(xué)地球科學(xué)學(xué)院,成都610059;3.首都師范大學(xué)資源環(huán)境與地理信息系統(tǒng)北京市重點(diǎn)實(shí)驗(yàn)室,北京100037)

土地利用的分布和變化,不僅反映了人與自然的關(guān)系變化,而且反映了人與人之間社會(huì)關(guān)系的變化[1],是研究全球變化的重要組成因素[2-3]。20世紀(jì)70年代以來(lái),隨著遙感與GIS技術(shù)的發(fā)展,土地利用變化的研究方法由定性向定量發(fā)展[4],眾多學(xué)者使用各類(lèi)指數(shù)模型如變化率、轉(zhuǎn)移矩陣、動(dòng)態(tài)度等來(lái)研究土地利用的變化量及變化規(guī)律,而對(duì)土地利用的集中與分散程度研究甚少。本文以西藏拉薩市城關(guān)區(qū)為例,采用空間羅倫茲曲線來(lái)研究土地利用的集散程度及其變化。

地貌因素對(duì)土地利用的影響巨大。研究區(qū)作為典型的山地地區(qū),人類(lèi)對(duì)土地的利用受到諸如地勢(shì)、地貌、氣候等因素的影響,表現(xiàn)出中緯度地帶高山寬谷地區(qū)鮮明的土地利用特點(diǎn),與內(nèi)地平原山地丘陵區(qū)明顯不同的是土地利用表現(xiàn)出隨海拔高度而變化的立體格局[5],因此本文著眼于研究土地利用類(lèi)型在不同海拔分布的集中分散程度,對(duì)研究區(qū)今后制定合理的土地利用規(guī)劃和土地資源的有效利用具有重要參考價(jià)值。

1 研究區(qū)概況

以西藏拉薩市城關(guān)區(qū)為研究區(qū),該區(qū)地處雅魯藏布江支流拉薩河中游,位于東經(jīng) 91°00′-91°18′,北緯 29°31′-29°49′,南北長(zhǎng) 31 km,東西寬 28 km,總面積525 km2,是拉薩市唯一的市轄區(qū),是西藏政治、經(jīng)濟(jì)、文化的中心和交通樞紐,也是西藏人流、物流、信息流的中心。該區(qū)山巒重疊,山高坡陡,溝谷縱橫,溝深谷狹,地勢(shì)南北高、中間低;高原山地和河谷交錯(cuò)分布的地形地貌特征對(duì)區(qū)域的水熱、植被、土地等資源環(huán)境狀況具有明顯的再分配作用,同時(shí)也制約著人們對(duì)土地等自然資源的利用方式以及人類(lèi)的生產(chǎn)和居住方式。

2 數(shù)據(jù)準(zhǔn)備

本次研究選用三期遙感影像構(gòu)成時(shí)間序列來(lái)研究研究區(qū)的土地利用及其變化情況:1991年的Landsat-TM影像,2000年的Landsat-ETM+影像以及2005年的CBERS影像。三期遙感影像的處理流程如下:(1)預(yù)處理:包括波段組合、幾何校正、影像裁剪、色調(diào)匹配等;(2)土地利用現(xiàn)狀分類(lèi):在ERDAS IMAGINE 8.7平臺(tái)進(jìn)行監(jiān)督分類(lèi),分為:耕地、林地、草地、水域、建設(shè)用地和未利用地6類(lèi),并在ArcGIS 9.2中對(duì)分類(lèi)結(jié)果進(jìn)行人工修改完善,制成三期土地利用現(xiàn)狀分類(lèi)圖;(3)DEM分帶:將研究區(qū)DEM 分為3帶:＜4 000 m,4 000～5 000 m,＞5 000 m,并利用ArcGIS 9.2中柵格數(shù)據(jù)空間分析模塊中的統(tǒng)計(jì)工具,統(tǒng)計(jì)得到各高程帶上各土地利用類(lèi)型的面積(表1)。

表1 土地利用類(lèi)型在各高程帶上分布的面積 km2

3 土地利用的空間格局變化分析

3.1 土地利用的空間羅倫茲曲線

空間羅倫茲曲線是由羅倫茲于20世紀(jì)20年代發(fā)表用于研究工業(yè)集中化的一種頻率累計(jì)曲線[5],它可以用于刻畫(huà)空間單元的分布,也可用于對(duì)兩個(gè)空間單元的分布作對(duì)比,是研究離散區(qū)域數(shù)據(jù)集中化程度與均衡度的重要方法[7]。羅倫茲曲線被廣泛應(yīng)用于研究區(qū)域收入分配的差異。本文將不同高程帶的土地利用數(shù)據(jù)繪制羅倫茲曲線來(lái)進(jìn)行研究區(qū)土地利用類(lèi)型在不同高程帶上分布的集散程度及其變化研究。

下面以研究區(qū)1991年的草地為例,說(shuō)明羅倫茲曲線的繪制方法[8]。

(1)計(jì)算每一高程帶上的草地面積占研究區(qū)草地總面積的百分比及每一高程帶上的土地面積占研究區(qū)總面積的百分比(表3)。

表2 百分比計(jì)算值

(2)計(jì)算R值。R值為各高程帶的草地面積百分比與其土地面積百分比的比值(表4)。

(3)將高程帶按R值降序排列,逐次計(jì)算各高程帶草地面積的累積百分比和各高程帶土地面積的累積百分比(表5)。

表3 草地的R值

表4 累積百分比計(jì)算值

(4)以土地面積累積百分比為橫坐標(biāo),草地面積累積百分比為縱坐標(biāo)繪制草地的羅倫茲曲線。

按此方法分別繪制三期土地利用類(lèi)型的羅倫茲曲線圖(圖1)。

某土地利用類(lèi)型的羅倫茲曲線越接近均勻分布線(對(duì)角線)則其分布越均勻,越接近集中分布線其分布越集中。由圖1可以看出,研究區(qū)草地的羅倫茲曲線最接近均勻分布線,說(shuō)明其在各高程帶上分布最為均勻,受高程的影響最小;水域和未利用地次之;林地、建設(shè)用地的羅倫茲曲線最接近集中分布線,說(shuō)明其在高程分布上最為集中,即受高程的影響最大;耕地次之。在隨時(shí)間的變化方面,各土地利用類(lèi)型的羅倫茲曲線變化不大,其中變化最大的為2005年水域的羅倫茲曲線,這是由于2005年冰川積雪面積的增大所致。

圖1 不同年份土地利用羅倫茲曲線圖

3.2 土地利用的集中化指數(shù)

為了進(jìn)一步分析土地利用的空間格局變化,下面計(jì)算土地利用的集中化指數(shù)。集中化指數(shù)是基于羅倫茲曲線計(jì)算而來(lái)的用于描述地理數(shù)據(jù)分布的集中化程度的指數(shù)[7]。

式中:I——集中化指數(shù);A——羅倫茲曲線下方區(qū)域面積;M——數(shù)據(jù)集中于一點(diǎn)時(shí)洛倫茲曲線下方(整個(gè)正方形)的面積,此次計(jì)算中為10 000;R——數(shù)據(jù)均勻分布時(shí)洛倫茲曲線下方(對(duì)角線以下三角形)的面積,此次計(jì)算中為5 000。因此,公式(1)轉(zhuǎn)化為

可以看出,集中化指數(shù)的值域?yàn)閇0,1],其值越大,數(shù)據(jù)的分布越集中。

按照公式(2)計(jì)算得到各期各土地利用類(lèi)型的集中化指數(shù)(圖2)。

從圖2中的集中化指數(shù)印證了由羅倫茲曲線所得出的結(jié)論,同時(shí)還可以看出:1991-2005年,林地和建設(shè)用地的集中化指數(shù)最高且保持不變,說(shuō)明其分布的集中與分散程度不變,這是因?yàn)榱值睾徒ㄔO(shè)用地對(duì)高程的依賴(lài)性較強(qiáng),均分布在海拔低于4 000 m高程帶內(nèi);耕地和草地的集中化指數(shù)經(jīng)歷了先降后升的過(guò)程,分別由1991年的0.596和0.199降低到2000年的0.581和0.193,到2005年又分別增加到0.591和0.206,說(shuō)明二者的集中程度經(jīng)歷了先降后升的過(guò)程;水域的集中化指數(shù)在1991-2000年期間略有增加(從0.438增加到0.450),在2000-2005年期間又呈現(xiàn)出大幅的下降(從0.450降低到了0.413),即其集中程度先升后降,這主要是受不同時(shí)期遙感圖像上冰川積雪面積的增減所致;未利用地的集中化指數(shù)不斷增加,且增幅較大,表明其在高程帶上的分布逐漸趨于集中,受高程的制約程度逐漸下降。

圖2 土地利用集中化指數(shù)

4 結(jié)論

實(shí)驗(yàn)證明:將空間羅倫茲曲線運(yùn)用在土地利用在各地貌類(lèi)型的空間格局分析中,無(wú)需查看數(shù)據(jù)便能直觀明了地反映各土地利用類(lèi)型分布的集中與分散程度;此外,在此基礎(chǔ)上計(jì)算的集中化指數(shù)能更進(jìn)一步精確、定量地反映土地利用的集散分布及其細(xì)微的變化情況。本文采用地貌因素作為分析單元來(lái)研究土地利用及其變化,而地貌因素僅為影響土地利用類(lèi)型的重要因素之一,因此也為研究影響土地利用分布與變化的其它因素提供了方法上的借鑒。如,若運(yùn)用此方法將分析單元進(jìn)一步細(xì)化或選擇行政區(qū)域作為分析單元,配以農(nóng)業(yè)、經(jīng)濟(jì)、人口等方面的指標(biāo)進(jìn)行相關(guān)分析,或可對(duì)今后制定合理的土地利用規(guī)劃和土地資源的有效利用具有重要參考價(jià)值。

[1] 李秀彬.土地利用變化的解釋[J].地理科學(xué)進(jìn)展,2002,21(3):195-201.

[2] 李秀彬.全球環(huán)境變化研究的核心領(lǐng)域:土地利用/土地覆被變化的國(guó)際研究動(dòng)向[J].地理學(xué)報(bào),1996,51(6):553-554.

[3] William E R,William B M,TurnerI I B L.Modeling land use and cover as part of global environmental change[J].Climate Change,1994,28(1/2):45-64.

[4] 譚少華,倪紹祥.20世紀(jì)以來(lái)土地利用研究綜述[J].地域研究與開(kāi)發(fā),2006,25(5):84-88.

[5] 除多.山地土地利用/土地覆蓋變化研究-以西藏拉薩地區(qū)為例[M].北京:氣象出版社,2007:56.

[6] 張超,楊秉賡.計(jì)量地理學(xué)基礎(chǔ)[M].2版.北京:高等教育出版社,1991:50-51.

[7] 徐建華.計(jì)量地理學(xué)[M].北京:高等教育出版社,2006:35-39.

[8] 常捷,李慧典,黃裕婕,等.基于 RS和GIS的中尺度流域土地利用空間結(jié)構(gòu)變化:以伊洛河流域中部洛寧縣為例[J].測(cè)繪學(xué)院學(xué)報(bào),2004,21(1):44-45.