基于GIS和RS的城市空間形態變化研究
——以西昌市為例

裴勇軍，周光紅，田學軍

(1.衢州市測繪院，浙江衢州324000;2.成都理工大學，四川 成都610059;3.中國冶金地質總局地球物理勘查院，河北保定071000)

一、前 言

城市空間形態變化研究，是城市發展趨勢分析的基礎，對城市遠景規劃有十分重要的指導作用。它也是城市人類活動演變以及區域城鎮化分布研究的重點內容和基礎內容。目前，基于GIS和RS的城市空間形態變化是關于城市發展趨勢研究的熱點，是對城鎮化和人類活動趨勢研究的有效手段。研究利用GIS和RS首次在西昌市基于不同時段影像、解譯數據和地形數據，完成了圖像數據、圖形數據和屬性數據的疊合，并在此基礎上，進行了西昌市空間形態變化和城鎮化的分析與評價。

二、研究方法及數據處理

一個城市的空間形態變化除了受到自然因素的影響，更多地受到人類活動和經濟等因子的影響。論文在研究西昌市空間形態發展模式時把影響西昌市空間形態變化的影響因子概括為自然因子和人類活動因子。人類活動因子主要是通過土地利用變化率來表示，作為研究的一個驅動力變化因子，通過分析1989年TM數據和1999年ETM數據求得土地利用變化率，土地利用變化率越高的地方，說明人類活動越活躍。自然因子主要是從西昌市的DEM數據提取坡度、坡向、地面粗糙度、相對高程作為研究的驅動力因子。坡度大小是制約生產力空間布局的重要因子，對交通、經濟及人類聚居影響很大;坡向影響人類居住地的選擇，也是直接造成局部地區氣候特征的差異，影響農業的發展;地表粗糙度是反映地形的起伏變化和侵蝕程度的宏觀地形因子;相對高程是地勢起伏大小的指標，地勢起伏對人類的居住習慣及城市空間形態變化有著重要的影響。最后通過AHP按不同權重疊加5個因子得到西昌市空間形態變化圖。研究方法思路圖見圖1。

1.人類活動因子——土地利用變化率的提取

(1)遙感圖像的預處理及監督分類

研究所采用的數據為于1989年05月11日拍攝的西昌市Landsat TM數據，以及于1999年11月23日拍攝的西昌市Landsat ETM數據。在Envi軟件中通過圖像幾何校正、圖像裁剪和拼接、多波段數據融合及波段選擇、圖像增強完成對TM和ETM圖像的預處理。在實地考察、參考其他學者的成果和反復解譯的基礎上建立各地物的解譯標志，利用Envi軟件采用最大似然(Likelihood Classification)分類器進行監督分類，再采用人工目視解譯的方法進行局部修改得到分類圖。

(2)土地利用變化率的提取

研究中土地利用變化的速度主要是通過土地利用類型的時空變化進行度量，它既可表征單一土地利用類型的時序變化，也可對區域土地利用動態的總體狀況及其區域分異進行分析。表達式如下

式中，Ka、Kb分別為區域某一特定土地利用類型研究期初及研究期末的面積;T為研究時段長;R為單一土地利用類型變化率。該指數可以刻畫區域土地利用變化程度。分類統計得出表1的結果。

表1 研究區域土地利用類型面積及其變化

利用ArcGIS軟件對分類圖進行分析即可得出土地利用變化率。土地利用變化率的值越高表示該地的土地利用變化越快，也表明該地人類活動越活躍，空間形態變化越明顯。為了便于研究，研究在ArcGIS中采用自然裂點法將土地利用變化率分為5個級別，即無變化、變化較低、變化一般、變化較快、變化顯著(如圖2所示)。

2.自然因子——坡度、坡向、相對高程、地面粗糙度的提取

(1)坡度因子的提取

研究采用的DEM數據利用經過幾何校正的地形圖提取的高程數據生成的。基于DEM的坡度提取通常在3像素×3像素的DEM柵格分析窗口中，采用幾何平面來擬合或差分計算的方法進行。研究在ArcGIS軟件中提取的坡度如圖3所示。

圖2 土地利用變化率分類圖

圖3 西昌市坡度圖

(2)坡向因子的提取

研究在ArcGIS軟件中基于DEM數據提取坡向。為了便于統一計算，需將坡向根據四川省坡向分級標準進行重分類合并為4類，即平緩地、陽坡、半陽坡、陰坡。最后按照此類分法提取出坡向分布如圖4所示。

圖4 西昌市坡向分級圖

(3)地面粗糙度的提取

地面粗糙度一般定義為地表單元的曲面面積與其在水平面上的投影面積之比。研究在ArcGIS軟件中基于DEM提取的地面粗糙度如圖5所示。

圖5 西昌市地面粗糙度

(4)相對高程的提取

本文采用的是領域統計提取方法［1］在ArcGIS軟件中提取西昌市的相對高程，整個研究區域的相對高程數據可以利用ArcGIS的柵格計算器求取，這種方法提取出來的相對高程數據由于采用了領域統計分析方法，因而提取出來的相對高程數據沒有被分級，數值具有連續性(如圖6所示)。

3.利用層次分析法(AHP)獲得空間形態變化圖

研究基于層次分析法求的以上5個因子的權重，解決了評價要素的不同屬性、不同度量標準、不同定性與定量標準等的統一化和規范化問題。

(1)建立層次結構模型

根據前面提取出的因子，將分類指標體系分為目標層、指標層(因素層)兩個層次。目標層為空間形態變化圖。指標層為:坡度、坡向、地面粗糙度、相對高程、土地利用變化率。

(2)構造判斷矩陣

判斷矩陣是表示本層所有因素針對上一層某一個因素的相對重要性的比較。判斷矩陣的元素用Santy的1～9標度方法給出。層次分析法的判斷矩陣標度及其含義見表2。

圖6 西昌市相對高程

表2 層次分析法的判斷矩陣標度及其含義

確定評價因子權重判斷矩陣

利用方根法計算判斷矩陣aij的最大特征根及其對應特征向量，經計算，矩陣的特征向量為:W=(0.192 1，0.097 3，0.109 7，0.118 9，0.482 0)此即為坡度，坡向，相對高程，地面粗糙度，土地利用變化率分別對空間形態變化趨勢的權重。

計算判斷矩陣的最大特征根為5.059 5。

(3)一致性檢驗

由于客觀事物的復雜性或對事物認識的片面性，通過判斷矩陣求出特征向量，需要對判斷矩陣進行一致性和隨機性檢驗，檢驗公式

式中，CR為判斷矩陣的隨機一致性比率;RI為判斷矩陣的平均隨機一致性指標［1］，由大量試驗給出。CI判斷矩陣的一致性指示，由式(3)確定

式中，λmax為判斷矩陣的最大特征根;m為判斷矩陣的階數。

CR小于0.10時認為判斷矩陣具有一致性;否則，不具有一致性，需要重新構造判斷矩陣。對于本次研究中構建的判斷矩陣，CR為0.013 3小于0.1，表明判斷矩陣具有滿意的一致性。

判斷矩陣構建后，通過計算各矩陣的特征值和特征向量，并進行一致性檢驗，得出各指標的權重。指標體系和權重的具體內容見表3。

表3 因子指標體系和權重

(4)求取空間形態變化圖

根據上述分析，結合西昌市空間動態變化的5個影響因子的量化結果，根據層次分析法確定的各因子權重，在ArcGIS軟件中進行多因子分級加權運算，即可得到西昌市空間動態變化趨勢圖 。值越大，說明該地發展速度越快，為了便于分析研究，本文采用自然裂點法將西昌市空間形態變化趨勢圖分為5個級別，即無變化、變化較低、變化一般、變化較快、變化顯著(如圖7所示)。

圖7 西昌市空間形態變化分類圖

三、結 論

通過對得到的西昌市空間形態變化圖的分析可以看出:

1)空間形態變化較大的地方主要集中在城鎮邊緣，可以看出西昌市城鎮變化方向為城鎮不斷擴大，大量侵占耕地和草地。從時間變化特征來看，研究區1989—1999年土地利用類型面積有較大幅度的變化，其中耕地減少幅度較大，10年間共減少1 829.1 0 hm2。 林 地 面 積 顯 著 增 加，共 增 加2 497.116 hm2，其原因主要是1998年以后國家實施天然林保護工程和退耕還林政策，研究區大量的陡坡耕地生態退耕所致，同時研究區發揮光熱資源優勢，進行產業結構調整，發展經濟林木也是導致林地增加的重要原因［6］。草地面積明顯減少，共減少了2 287.42 hm2。建設用地呈現出快速增加的趨勢，10年面積增加2 042.24 hm2;水域和未利用地分別減少了 98.55 hm2和331.26 hm2。

2)將空間形態變化圖與坡度圖、坡向圖、地面粗糙度圖和相對高程圖進行分析統計，可以得到坡度、坡向、地面粗糙度和相對高程對城市空間形態變化關系表(見表4)。

表4 坡度、坡向、地面粗糙度和相對高程對城市空間形態變化關系表

可以看出人類居住地逐漸往低坡度的平原地區發展，坡度在45°下，空間形態變化較大。坡向為平緩地和半陽坡空間形態變化較大。相對高程在200～450 m的地方變化較大，2000 m以上的地方基本上無變化。

3)空間形態變化活躍的地方主要集中在瀘黃高速公路兩旁、邛海周圍及城鎮化邊緣。研究結果與實際情況相符合，表明研究方法能夠較準確地預測出城市空間形態變化趨勢。

通過此項研究表明:在西昌市首次建立了基于不同時段影像、解譯數據和地形數據的基礎上，完成了圖像數據、圖形數據和屬性數據的疊合，并在此基礎上，進行了西昌市空間形態變化和城鎮化的分析與評價。研究利用RS技術和GIS技術相結合的技術方法，可以方便快速地得出土地利用類型率，以及坡度、坡向、地面粗糙度和相對高程等自然因子對城市空間形態變化的影響，能夠較準確地預測出城市空間形態變化趨勢，對城鎮土地利用管理與城鎮遠景規劃有十分重要的指導作用。

［1］楊斌.“數字山地”框架下的山地本體及數字化分類研究——以四川省為例［D］.成都:成都理工大學，2009.

［2］鄧書斌.ENVI遙感圖像處理方法［M］.北京:科學出版社，2010.

［3］祁燕.基于遙感與GIS的北京市城鎮化進程的動態研究［D］.北京:北京林業大學，2009.

［4］胡玉福，鄧良基，張世熔.基于RS和GIS的西昌市土地利用及景觀格局變化［J］.農業工程學報，2011，27(10):322-327.

［5］李立娜.西昌市土地資源可持續利用評價研究［J］.安徽農業科學.2009，37(36):18133-18135.

［6］彭建，柳昆，鄭付濤等.基于AHP的地下空間開發利用適宜性評價［J］.地下空間與工程學報，2010，6(4):688-694.

［7］湯國安，楊昕.地理信息系統空間分析實驗教程［M］.北京:科學出版社，2006.

［8］周可法，吳世新.基于RS和GIS技術下城鎮空間變化分析及應用研究［J］.干旱區地理，2002，25(1):61-64.