裴勇軍，周光紅，田學(xué)軍

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

一、前 言

城市空間形態(tài)變化研究，是城市發(fā)展趨勢(shì)分析的基礎(chǔ)，對(duì)城市遠(yuǎn)景規(guī)劃有十分重要的指導(dǎo)作用。它也是城市人類活動(dòng)演變以及區(qū)域城鎮(zhèn)化分布研究的重點(diǎn)內(nèi)容和基礎(chǔ)內(nèi)容。目前，基于GIS和RS的城市空間形態(tài)變化是關(guān)于城市發(fā)展趨勢(shì)研究的熱點(diǎn)，是對(duì)城鎮(zhèn)化和人類活動(dòng)趨勢(shì)研究的有效手段。研究利用GIS和RS首次在西昌市基于不同時(shí)段影像、解譯數(shù)據(jù)和地形數(shù)據(jù)，完成了圖像數(shù)據(jù)、圖形數(shù)據(jù)和屬性數(shù)據(jù)的疊合，并在此基礎(chǔ)上，進(jìn)行了西昌市空間形態(tài)變化和城鎮(zhèn)化的分析與評(píng)價(jià)。

二、研究方法及數(shù)據(jù)處理

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

1.人類活動(dòng)因子——土地利用變化率的提取

(1)遙感圖像的預(yù)處理及監(jiān)督分類

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

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

研究中土地利用變化的速度主要是通過土地利用類型的時(shí)空變化進(jìn)行度量，它既可表征單一土地利用類型的時(shí)序變化，也可對(duì)區(qū)域土地利用動(dòng)態(tài)的總體狀況及其區(qū)域分異進(jìn)行分析。表達(dá)式如下

式中，Ka、Kb分別為區(qū)域某一特定土地利用類型研究期初及研究期末的面積;T為研究時(shí)段長(zhǎng);R為單一土地利用類型變化率。該指數(shù)可以刻畫區(qū)域土地利用變化程度。分類統(tǒng)計(jì)得出表1的結(jié)果。

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

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

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

(1)坡度因子的提取

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

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

圖3 西昌市坡度圖

(2)坡向因子的提取

研究在ArcGIS軟件中基于DEM數(shù)據(jù)提取坡向。為了便于統(tǒng)一計(jì)算，需將坡向根據(jù)四川省坡向分級(jí)標(biāo)準(zhǔn)進(jìn)行重分類合并為4類，即平緩地、陽坡、半陽坡、陰坡。最后按照此類分法提取出坡向分布如圖4所示。

圖4 西昌市坡向分級(jí)圖

(3)地面粗糙度的提取

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

圖5 西昌市地面粗糙度

(4)相對(duì)高程的提取

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

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

研究基于層次分析法求的以上5個(gè)因子的權(quán)重，解決了評(píng)價(jià)要素的不同屬性、不同度量標(biāo)準(zhǔn)、不同定性與定量標(biāo)準(zhǔn)等的統(tǒng)一化和規(guī)范化問題。

(1)建立層次結(jié)構(gòu)模型

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

(2)構(gòu)造判斷矩陣

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

圖6 西昌市相對(duì)高程

表2 層次分析法的判斷矩陣標(biāo)度及其含義

確定評(píng)價(jià)因子權(quán)重判斷矩陣

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

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

(3)一致性檢驗(yàn)

由于客觀事物的復(fù)雜性或?qū)κ挛镎J(rèn)識(shí)的片面性，通過判斷矩陣求出特征向量，需要對(duì)判斷矩陣進(jìn)行一致性和隨機(jī)性檢驗(yàn)，檢驗(yàn)公式

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

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

CR小于0.10時(shí)認(rèn)為判斷矩陣具有一致性;否則，不具有一致性，需要重新構(gòu)造判斷矩陣。對(duì)于本次研究中構(gòu)建的判斷矩陣，CR為0.013 3小于0.1，表明判斷矩陣具有滿意的一致性。

判斷矩陣構(gòu)建后，通過計(jì)算各矩陣的特征值和特征向量，并進(jìn)行一致性檢驗(yàn)，得出各指標(biāo)的權(quán)重。指標(biāo)體系和權(quán)重的具體內(nèi)容見表3。

表3 因子指標(biāo)體系和權(quán)重

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

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

圖7 西昌市空間形態(tài)變化分類圖

三、結(jié) 論

通過對(duì)得到的西昌市空間形態(tài)變化圖的分析可以看出:

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

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

表4 坡度、坡向、地面粗糙度和相對(duì)高程對(duì)城市空間形態(tài)變化關(guān)系表

可以看出人類居住地逐漸往低坡度的平原地區(qū)發(fā)展，坡度在45°下，空間形態(tài)變化較大。坡向?yàn)槠骄彽睾桶腙柶驴臻g形態(tài)變化較大。相對(duì)高程在200～450 m的地方變化較大，2000 m以上的地方基本上無變化。

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

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

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

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

［3］祁燕.基于遙感與GIS的北京市城鎮(zhèn)化進(jìn)程的動(dòng)態(tài)研究［D］.北京:北京林業(yè)大學(xué)，2009.

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

［5］李立娜.西昌市土地資源可持續(xù)利用評(píng)價(jià)研究［J］.安徽農(nóng)業(yè)科學(xué).2009，37(36):18133-18135.

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

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

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