大型山地城市建設用地空間擴展及地形分異

呂志強，鄧 睿，卿珊珊

（重慶工商大學 旅游與國土資源學院，重慶400067）

大型山地城市建設用地空間擴展及地形分異

呂志強，鄧 睿，卿珊珊

（重慶工商大學 旅游與國土資源學院，重慶400067）

近年來，我國城市化進程已經進入加速發展階段，城市建設用地日益緊張，對山地城市來講，建設用地的供需矛盾尤為突出。以城市建設用地擴張的地形梯度變化為切入點，運用重慶市主城區2001年、2007年和2014年的遙感影像數據，通過混合像元分解提取建設用地。通過計算建設用地的分形維數、緊湊度指數和放射狀指數分析了主城區城市建設用地擴張水平方向的形態變化特征；結合區域的高程、地形起伏度，分析了建設用地的地形分異特征。結果表明：城市規模擴展迅速，建設用地數量從2001年的237.62 km2增加到2014年的818.85 km2；綜合分形維數指數、城市緊湊度指數、放射狀指數來看，城市擴展從2001—2007年的緊湊型擴展發展轉變成2007—2014年的外延擴展；主城區建設用地的擴展呈從低高程區的增長達到飽和狀態向高程較高的區域發展的態勢；從起伏度分析來看，建設用地無法向起伏度較大的區域大規模擴展。

建設用地擴展；混合像元分解；地形梯度；大型山地城市

城市建設用地擴張是城市發展在空間上最為顯著的特征，是城市地理學、土地科學、城市規劃學等多學科共同關注的領域[1]，作為 LUCC（Land-Use and Land-Cover Change）中變化最快、環境影響最劇烈的城市建設用地空間擴張直觀地體現了城市化和城市發展的空間過程與特征，開展這方面的研究對于合理配置土地資源、優化城市空間結構、實現城市可持續發展具有重要意義。

20世紀80年代以來，城市建設用地的急劇擴張成為中國城市化的突出特征和社會經濟加速發展的重要表現，針對此方面的研究主要集中在20世紀90年代。對城市建設用地擴張的研究主要集中在長三角、珠三角等經濟發達的沿海地區和內陸平原地區，例如上海市[2]、南京都市區[3-4]、無錫市[5]、常熟市[6]、北京市[7]、石家莊市[8]和長沙市[]等；黎夏等[10]在遙感多時相數據的基礎上通過熵的計算，定量地描述了珠江三角洲城市擴張的空間規律和擴張過程；周翔等[11]結合遙感影像與社會經濟數據分析了1990—2010年蘇錫常地區建設用地擴張的時空特征、空間模式及社會經濟驅動力。

隨著城市化進程的加快以及國家對建設用地需求的增長，城市建設逐漸由人口過度稠密的平原地區轉向人居環境更佳的山區。在中國，山地、丘陵和高原占了全國陸地面積的69%，分布在上述山地區域的城市稱為山地城市。因其地理條件的特殊性，山地城市的發展受自然條件限制，近幾年隨著經濟發展和城市化建設的不斷加強，更多的學者對山地城市的發展特征和規律進行研究。重慶作為西南地區的特大山地城市，近年來也有研究者開始了對重慶城市用地的研究。廖和平等[12]分析了重慶市直轄以后的城市空間擴展強度及其空間差異、城市空間擴展影響機制，定量揭示了城市空間擴張演變特征及其驅動機制的評價方法；牟鳳云等[13]運用多時相遙感數據對重慶市1978—2005年城市空間演化過程、動力機制、城鄉用地轉化和城市用地的內部分異進行了研究。此類研究多是針對建設用地水平方向的擴張特征進行研究，缺乏在豎直方向上的擴張特征的研究。

關于城市建設用地空間擴展的研究都是基于城市建設用地范圍的提取，即土地利用分類，主要的土地利用分類方法有目視解譯和包括神經網絡分類方法、模糊分類方法以及決策樹分類方法等方法在內的遙感分類方法。隨著遙感技術的迅速發展，Ridd[14]提出將城市影像中的每個像元看成植被—不透水面—土壤這3種代表類型的線性組合，即V-I-S模型（Vegetation-impervious surface-soil）。此模型成為提取非滲透表面豐度的重要方法，Carlson等[15]利用植被覆蓋度與不透水面之間的關系，研究了適合于城市建成區不透水面信息提取的方法；Wu等[16]利用光譜混合分析法（SMA）對Landsat 7ETM+影像進行分解，提取出不透水面；錢樂祥等[17]用線性光譜模型有效地解決了亮度變化與陰影關聯的問題，展現了該方法的應用前景；李明杰等[18]在采用V-I-S模型與歸一化線性混合模型相結合的方法提取的非滲透表面豐度圖的基礎上，以地表溫度數據進行糾正，進一步提高豐度提取精度。

本研究嘗試運用混合像元分解的方法較高精度提取城市建設用地范圍，結合城市空間形態指數研究城市范圍的擴張過程和特征；根據提取端元的豐度值劃分建設用地密度等級，結合地形因素研究建設用地及建設用地等級在地形因素影響下的擴張水平和擴張特征。

1 研究區概況與數據來源

1.1 研究區概況

重慶市地跨東經105°17′—110°11′，北緯28°10′—32°13′，位于四川盆地、中國西南內陸青藏高原與長江中下游平原的過渡地帶，是全國著名的“山城”，是中國重要的中心城市之一，也是西南地區和長江上游最大的經濟中心城市和科技、文化、教育事業的中心。重慶市主城區位于川東平行嶺谷地區，明月山、銅鑼山、中梁山、縉云山等多條華鎣山余脈從東北向西南延伸，長江和嘉陵江自西向東切割山脈而過，是典型的“山水城市”。主城區面積5 473 km2，2013年常住人口為808.5萬人，人均GDP為6.67萬元，城鎮化率達到87.87%。

由于早期我國西南地區經濟較為落后，城市建設發展速度緩慢，重慶的城市布局形態呈現一種原始自然的點狀。隨著我國經濟的發展以及國家推出重點發展西南地區經濟的政策，山地城市建設進程加快，城市建設規模不斷加大。尤其是重慶成為中央直轄市后，城市逐漸由集中向分散發展，并向著與區域城市相融合、與外圍都市緊密聯合的方向發展。

重慶市主城區自改革開放，尤其是1997年設立為直轄市以來，城市化速度加快，區內建成區面積不斷擴大。2010年我國內陸地區第一個國家級開發開放新區兩江新區的成立使得主城區北部成為建設的熱點地區。本文因數據限制，采用處理過后的邊界截取大部分主城區作為研究區域，截掉部分土地覆被類型多為植被，對城市建設用地的擴張研究基本沒有影響。

1.2 數據來源與處理

重慶市主城區主要覆蓋LANDSAT TM/ETM/OLI軌道號為p128/r39和p128/r40的遙感影像中，本文為了體現重慶市直轄后的快速城市化過程，選擇2001年、2007年和2014年遙感影像圖作為數據源，其中2001年為ETM+數據，2007年為TM數據，2014年為OLI數據。首先在ENVI中對遙感影像進行幾何校正、配準、鑲嵌等處理，然后采用混合像元分解方法提取不透水面密度，并進一步提取建設用地。

圖1 研究區位置

為體現地形因素對山地城市空間演化的影響，本文將分析城市隨海拔高度和地形起伏度的變化呈現出來的空間拓展特點，故選擇DEM數據為輔助數據對山地城市的空間拓展進行分析。本文采用空間分辨率為85.6m的STRM（Shuttle Radar Topography Mission）DEM數據，在ArcGIS中計算地形起伏度并與建設用地密度相結合進行分類疊加、統計。

2 研究方法

2.1 混合像元分解

遙感影像是以像元為基本單位來檢測和獲取地物信息，混合像元是指一個像元中包括了兩種或兩種以上地物的現象。將每一像元分解，求得其覆蓋類型組分占像元面積的百分比的過程被稱為混合像元分解。本文采用歸一化線性光譜混合模型（NSMA）進行混合像元分解，首先對各波段的光譜值進行歸一化處理，然后依據線性光譜混合模型（LSMM）原理對混合像元進行高精度分解的方法。其中，光譜歸一化算法公式為：

線性光譜混合模型（LSMM）算法公式為：

式中：Riλ為第λ波段第i像元的光譜反射率（已知）；fki為對應于i像元組分第k個端元所占的分量值（待求）；Ckλ為第k個基本組分在第λ波段的光譜反射率；εiλ為殘余誤差值（即光譜的非模型化部分）；n為基本組分的數目。評價模型用殘差或均方根誤差RMS表示：

本研究在ENVI中，選擇植被、水體、建設用地為端元，對遙感圖片進行混合像元分解，采用選擇限制性條件的混合像元分解方法。提取出建設用地端元之后，根據端元的豐度值將建設用地分為高密度、中密度、低密度3個等級。

2.2 城市形態定量描述

在GIS平臺中，對得到的城市建設用地結果進行統計與分析，利用城市空間形態的定量模型來獲取城市形態的相關參數。

（1）城市空間形狀分形維數。

空間的分形維數可以描述城市邊界形狀的復雜性，反映出土地利用形狀的變化及土地利用受干擾的程度，是一個面積與周長的關系，計算公式為：

式中：St為t時期城市斑塊的分形維數；Pt，At分別為t時期城市斑塊的面積和周長；St的理論范圍在1～2，St值越大表示圖形形狀越復雜。

（2）城市空間緊湊度指數。

式中：K為城市的緊湊度；A為城市面積；P為城市輪廓周長。這一公式以圓形區域作為標準度量單位，圓形地物的緊湊度為1，其他任何形狀地物的緊湊度均小于1，如正方形的緊湊度為，地物離散程度越大，其緊湊度越低。緊湊度大的圖斑受外界干擾小，更容易保持內部資源的穩定性。

（3）Boyce-Clark形狀指數。

其基本思想是將研究的形狀與標準圓形形狀進行比較，得出一個相對指數的方法。這種方法是以面狀

式中：SBC為Boyce-Clark形狀指數；ri為某個圖形的優勢點（Vantage point）到圖形周界的半徑長度；n為具有相等角度差的輻射半徑的數量。n可以取不同的數量，數量越大，形狀指數值精度越高。對于該研究，城市形狀優勢點是圖形質心（Centroid），半徑的取值方向采用32方向，基本上可以避免一對多或多對一現象的出現。

2.3 起伏度計算

在計算地形起伏度之前需要確定起伏度的計算尺度，張偉等[18]在全國范圍內選取試驗區，在漸變尺度下對中國地形起伏度適宜計算尺度進行研究，認為基于SRTM和ASTER DEM這兩種不同的數據源，全國范圍內的地形起伏度適宜計算尺度分別為4.72，3.20 km2。本研究采用的是STRM（Shuttle Radar Topography Mission）DEM數據，結合重慶市主城區的典型性，研究中以3×3的計算窗口進行地形起伏度計算。

采用Arc Map中的空間分析（Spatial Analyst）模塊中的鄰域分析（Neighborhood Statistic）工具，先計算出3×3像元內的格網最大值，然后計算出其領域最小值，再利用模塊中的柵格計算工具（Raster Calculator）計算最大值與最小值高程差，就得到了該窗口的地形起伏度結果值。可表示為如下公式：

式中：RF是分析窗口內的地勢起伏度；Hmax是分析窗口內的最大高程值；Hmin是分析窗口內的最小高程值。

2.4 結合海拔高度和起伏度分級分區統計

參考李明杰等[19]的方法在GIS中將提取的建設用地分為低密度、中密度和高密度3個等級，本文研究的是建設用地相對密度在地形維度上的空間擴張，而不是具體建設用地圖斑密度的絕對數量的變化；采用平均面積加權法（AWAE）將高程分為8個區間，用1—8標記；采用1/2 Standard Deviation方法將起伏度分為10個區間，用1—10標記。結合海拔高度和起伏度分區，統計每一區間內3個建設用地等級所占數量和比例。地物的重心（或圖形的質心、行政區范圍內的人口中心、經濟中心等）為核心，向周邊延伸若干條半徑（半徑之間的夾角相等），計算各條半徑的變化情況，所以也被稱為半徑形狀指數，用放射狀指數SBC來表示。

3 結果與分析

3.1 建設用地擴展的總體特征

重慶市主城區位于長江和嘉陵江交匯處，城區被河流分割為3個主體部分，區域大小和地形條件的差異使各區域的擴展速度和方向差異較大。自改革開放以來，隨著經濟的發展和生產力布局的優化調整，重慶市主城區城市空間進入了前所未有的快速擴展階段。尤其是重慶市直轄以來，主城區建設用地面積在2001年約為237.62 km2，增加到2007年的328.99 km2，2014年的818.85 km2，其中2007—2014年建設用地擴張最為劇烈，增速為69.98 km2/a，和2001—2007年的15.23 km2/a相比，增速明顯，達到了2001—2007年的4.6倍。

從表1可以看出，2001—2007年主城區建設用地輪廓變化不大，高密度等級建設用地分布較廣；2007—2014年建設用地跨越了縉云山、中梁山、銅鑼山和明月山4條山脈大規模擴張，高密度等級比例嚴重下降。

表1 建設用地擴張總體特征 km2

2001—2007年主城區建設用地重心向東北遷移，由于山體阻隔，建設用地向正西和正東方向的拓展空間越來越小，主城向北擴張速度加快；2007—2014年主城區建設用地重心繼續向東北偏移，特別是2010年兩江新區的成立，使得城市向北擴張的趨勢繼續增強。

3.2 建設用地擴展的水平特征

在ArcGIS中提取出城市建設用地的面積和邊界周長，并計算出城市空間形態指數，計算結果見表2。

表2 主城區城市空間形態演化指數

從表2可以看出，2001—2014年，主城區城市形態分形維數為1.548～1.580，分形維數理論范圍為1～2，說明重慶市主城區的城市邊界形狀較復雜；但城市分形維數變化不大，說明城市形態穩定性較高，結合重慶的實際情況，說明城市一直在建設，邊界形態沒有穩定下來。2001—2014年，主城區緊湊度指數為0.003～0.004，遠遠低于1，由緊湊度指數的定義來看，緊湊度指數介于0～1，當緊湊度指數接近于1時，其形狀接近于圓形，城市內斑塊受外界干擾小，更容易保持內部資源的穩定性；而主城區的緊湊度遠遠小于1，說明城市形態離散程度大。同期城市形態放射狀指數變化較為明顯，并呈現出一直減小的趨勢。根據放射狀指數的定義來看，放射狀指數越趨于0，表明城市內各圖斑分布越緊湊，指數越大說明各圖斑分布越離散，較好地反映了區域內部的聯系強度。

重慶市主城區地貌多為丘陵山地，平坦地面較少，主要是沿江斷續階地和緩丘平壩。這就導致了城區的城市擴張不像平原城市均衡地向各個方向擴展，而主要是沿交通干線向外擴張，在不同方向上，擴張的速度和規模都有較大差異，整個城市呈跳躍式擴張。

（1）2001—2007年處于直轄后的加速發展階段。這一時期城市逐漸進入了快速擴張期，增加的建設用地面積約為91.37 km2。這個階段內分形維數下降，緊湊度指數增加，放射狀指數下降較為明顯，表明這一時期內空間擴張以填充型為主。結合實際情況來看，城市輪廓變化不大，江北機場建設逐漸完善，靠近機場的兩路片區得到迅速發展，渝北區介于唐家沱和觀音橋—人和之間的腹地、機場高速沿線地帶成為了城市擴張的主要區域。

（2）2007—2014年處于城市的快速擴張期。這一時期城市擴張進入了跳躍式空間擴張階段，增加的建設用地面積約為489.86 km2。階段內分形維數上升，緊湊度指數和放射狀指數下降，表明這一時期內空間擴張以外延發展為主。隨著兩江新區的成立，北部新區和兩路片區成為了主城區建設的熱點地區。此外，隨著微電子產業園、物流園及大學城的持續建設，主城區的西部得到快速發展；城區東南部由于茶園新區的建設也進行了快速擴張；城區東北部江北區復盛魚嘴片區異軍突起。城區跳躍性向城市外圍大規模擴張，城市外圍地區的城鎮開始承擔起中心城區外溢功能和新增職能，與中心城區協同發展。

3.3 建設用地擴展的地形分異

（1）地形重心轉移。在GIS中分別統計各年份建設用地所在的海拔高度和起伏度的最小值、最大值、平均值和標準差，結果如表3、表4所示。

從建設用地所在海拔高度的最大值來看（表3），2001—2007年建設用地向高度較高的區域擴展。此后，因為針對區內的各類自然保護區、風景名勝區、森林公園區和“四山”禁建區等高海拔區域的生態保護政策，到了2014年建設用地所在的海拔高度下降到低于2001年，其中的建設用地向其他用地轉化。結合建設用地所在海拔高度的平均值和標準差來看，2001—2014年，建設用地整體持續向海拔高度較高的區域擴展，所在高度也相對集中。

從建設用地所在地形起伏度的最大值來看（表4），隨著經濟發展和技術水平的提高，所在起伏度也經歷了增大爾后又降低的過程。結合建設用地所在起伏度的平均值和標準差來看，和海拔高度不同，起伏度的平均值基本保持不變；標準差相差不大，但保持逐漸減小，說明建設用地所在起伏度相對而言更加集中。

表3 建設用地所在海拔高度特征 m

表4 建設用地所在地形起伏度特征

（2）建設用地擴展與海拔高度。在GIS中對研究區的高程分區，分別統計2001年、2007年和2014年的建設用地面積。

從每一高程區間建設用地的分布情況來看（圖2），主城區建設用地在每一高程區間內逐年增長，尤其是2007—2014年增幅較大。從每一年建設用地的分布情況看，主城區建設用地呈現偏態分布，主要分布在1，2，3區，分別占了全部建設用地面積的75.3%，74%和73.4%。可以看出，3個年份建設用地面積在1—3區所占比重相差不大，說明主城區建設用地按照一定的速度和規律穩定地進行擴展。2001年和2007年建設用地面積分布區間都是從1區增長到2區達到最高，然后逐漸下降。直到2014年建設用地面積分布區間從1區逐漸增長到3區達到最高，然后逐漸下降。說明隨著年份的增加，主城區建設用地在低高程區的增長達到飽和狀態，逐漸向高程較高的區域發展。

在高程分區的基礎上統計2001年和2014年主城區建設用地密度等級，并進行比較分析，詳見圖3。

從圖3來看，2001年主城區建設用地主要分布在高程1區、2區和3區。建設用地在每一個高程分區內密度分布較為均勻，在高程1—4區均是中密度等級建設用地最多，低密度等級次之，高密度等級最少；隨著高程的增加，擴展難度加大，建設用地密度等級分布呈現不同的趨勢，在高程5—8區低密度等級建設用地最多，中密度等級次之，高密度等級最少。

從2014年來看，主城區建設用地仍然主要分布在高程1—3區，但高程4區的比重有增加的趨勢。建設用地整體分布和在每一高程區間內的分布出現極為明顯的不均勻，高密度等級建設用地比重只有4.5%。在高程2—5區間內中密度等級建設用地最多，低密度等級次之，高密度等級最少；而在高程最低的1區和高程達到一定高度的6—8區內低密度等級的建設用地最多，中密度次之，高密度等級最少，高程8區從圖3可以看出幾乎沒有，說明在研究區內高程的最高點仍然有建設用地但排列極為稀疏。

結合重慶主城區的實際情況，2001—2014年重慶直轄以來，隨著經濟的飛速發展，建設用地急速擴展，經歷了“擴展—飽和—擴展”的過程。建設用地呈現出由適宜性最強的低高程區慢慢向適宜性較強的較高高程區擴展的趨勢。且隨著建設用地的擴展，新建設用地的開發建設，出現了大面積的待建和新建建設用地，使得原有的高密度等級建設用地轉化為中密度等級的建設用地；或是由于整個建設用地面積的增加，基數增大，使得原有的高密度等級建設用地比重嚴重下降。

（3）建設用地擴展與地形起伏度。在GIS中對研究區的起伏度分區，分別統計2001年、2007年和2014年的建設用地面積，得到如下結果。

圖3 2001年和2014年高程分區統計建設用地密度等級比例及標準差

從每一年份建設用地分布情況來看（圖4），年份增加建設用地隨之逐漸增長，分散到每一區間仍然可以看出2007—2014年建設用地急劇增長的情況。從每一起伏度區間建設用地分布情況來看，起伏度1區間建設用地面積變化不大，重慶為典型山地城市，地勢平坦的起伏度1區間開發較早，已經近乎飽和，總量變化不大。建設用地擴展程度最劇烈的是起伏度2區間、3區間和4區間，數量從2001年的205.36 km2增加2007年的285.69 km2以及2014年的717.77 km2，分別占總面積的86.4%，87.0%和87.7%，比重稍微有所增加但相差不大。這3個起伏度區間開發難度相對較低，在現有的技術水平和經濟投入下這部分的土地轉換為建設用地的可能性較大。其中絕對數量增加最多的是第2區間，主要分布渝北區和兩江新區，也是主城區開發的重點區域。建設用地面積變化最小的是起伏度8區間、9區間和10區間，建設用地面積僅從2001年的2.67 km2增加到2014年的5.10 km2。從難易程度和開發的可行性方面來講，這3個起伏度區間建設開發的可能性最小，所以變化最小。與高程分區不同，隨著起伏度的增加，建設用地擴展規模和速度迅速減小，建設用地無法向起伏度較大的區域大規模擴展。

圖4 起伏度分區統計建設用地面積

圖5是在起伏度分區的基礎上統計2001年和2014年主城區建設用地密度等級，并進行比較分析。從圖5可以看出，主城區建設用地高密度等級在各個起伏度分區比例下降明顯，說明隨著城市的擴展，使得城市建設用地總量這一基數快速增加，同時低密度和中密度等級建設用地的用地數量增加更為顯著；隨著城市的發展，城市生態環境的變化，例如城市綠地的增加，樹冠的生長以及舊建筑樓頂積水等，使得較早年份的高密度等級建設用地轉化為非高密度等級，導致高密度等級建設用地的比例降低。

圖5 2001年和2014年起伏度分區統計建設用地密度等級比例及標準差

4 結論

（1）2001—2014年，隨著經濟發展、城市化建設的不斷加強和技術水平的提高，在一定程度上打破了山地城市自然條件的限制，重慶市主城區建設用地跨越4大山脈快速擴展，經歷了“擴展—飽和—擴展”的過程，建設用地在已有規模內達到飽和后進行外延式擴展。

（2）建設用地在不同高程區間和起伏度區間下分布差異特征明顯，隨著開發技術水平的進步，建設用地有向高程較高的區域擴展的趨勢；但無論技術如何進步，起伏度較大的區域建設用地數量增加不明顯，說明建設用地無法大規模向起伏度較大的區域擴展。從不同年份建設用地在高程和起伏度各個區間所占比例來看，主城區建設用地按照一定的速度和規模穩定地進行擴展。

文章雖然采用改進的混合像元分解模型并且在此基礎上反復試驗提取城市建設用地范圍，以期提高提取精度，但也會把城市化水平較高的區域內公園、公共綠地等誤分為非建設用地，或是因遙感影像采集的季節不同，水體干涸后的灘涂等誤分為建設用地，希望在以后的研究中繼續改進。

[1]曹廣忠，郜曉雯，劉濤.都市區與非都市區的城鎮用地增長特征：以長三角地區為例[J].人文地理，2011，26（5）：65-70.

[2]張心怡，劉敏，孟飛.基于RS和GIS的上海城建用地擴展研究[J].長江流域資源與環境，2006，15（1）：29-33.

[3]姚士謀.南京大都市空間演化與地域結構發展策略[J].地理學與國土研究，2001，17（3）：7-11.

[4]陳江龍，高金龍，徐夢月，等.南京大都市區建設用地擴張特征與機理[J].地理研究，2014，33（3）：427-438.

[5]薛東前，王傳勝.無錫城市用地擴展的時空特征與趨勢分析[J].資源科學，2003，25（1）：9-14.

[6]馬榮華，陳雯，陳小卉，等.常熟市城鎮用地擴展分析[J].地理學報，2004，59（3）：418-426.

[7]顧朝林.北京土地利用/覆蓋變化機制研究[J].自然資源學報，1999，14（4）：307-312.

[8]肖捷穎，葛京鳳，沈彥俊，等.基于GIS的石家莊市城市土地利用擴展分析[J].地理研究，2003，22（6）：789-798.

[9]林目軒，陳秧分，師迎春，等.大城市內部建設用地擴張差異及其原因：以長沙市區為例[J].經濟地理，2006，26（5）：836-841.

[10]黎夏，葉嘉安.利用遙感監測和分析珠江三角洲的城市擴張過程：以東莞市為例[J].地理研究，16（4）：56-62.

[11]周翔，陳亮，象偉寧.蘇錫常地區建設用地擴張過程的定量分析[J].應用生態學報，2014，25（5）：1422-1430.

[12]廖和平，彭征，洪惠坤，等.重慶市直轄以來的城市空間擴展與機制[J].地理研究，2007，26（6）：1137-1146.

[13]牟鳳云，張增祥，譚文彬.基于遙感和GIS的重慶市近30年城市形態演化特征分析[J].云南地理環境研究，2008，20（5）：1-5.

[14]Ridd M K.Exploring a VIS（vegetation-impervious surfacesoil）model for urban ecosystem analysis through remote sensing：comparative anatomy for cities[J].International Journal of Remote Sensing，1995，16（12）：2165-2185.

[15]Carlson T N，Arthur S T.The impact of land use-land cover changes due to urbanization on surface microclimate and hydrology：a satellite perspective[J].Global and Planetary Change，2000，25（1）：49-65.

[16]Wu C，Murray A T.Estimating impervious surface distribution by spectral mixture analysis[J].Remote Sensing of Environment，2003，84（4）：493-505.

[17]錢樂祥，崔海山.運用歸一化光譜混合模型分析城市地表組成[J].國土資源遙感，2006，18（2）：64-68.

[18]張偉，李愛農.基于DEM的中國地形起伏度適宜計算尺度研究[J].地理與地理信息科學，2012，28（4）：8-12.

[19]李明杰，錢樂祥，陳健飛.非滲透表面豐度提取方法運用：以廣州市海珠區為例[J].遙感信息，2011（2）：36-40.

Spatial Expansion and Topographic Differentiation of Urban Construction Land in Large Mountainous Cities

LYU Zhiqiang，DENG Rui，QING Shanshan
（Department of Land Resource Management，Chongqing Technology and Business University，Chongqing400067，China）

The urbanization process in China has entered the accelerating stage in recent years，the urban construction land is increasingly scarce，especially in the mountainous city，the contradiction between supply and demand of construction lands is obvious.Based on the terrain gradient change of urban construction land expansion，we used remote sensing image data of the main urban area in 2001，2007 and 2014 in Chongqing City to extract the urban construction lands by using the method of decomposition of mixed pixels.The fractal dimension of construction land，compact index and the radial index were computed to analyze the main urban construction land expansion form change characteristics of horizontal direction.Combining urban construction land with the regional elevation and relief，we analyzed the construction land differentiation characteristics of terrain.The results showed that the city scale expanded rapidly，and the area of construction land increased from 237.62 km2to 818.85 km2during the period from 2001 to 2014.With respect to comprehensive fractal dimension index，compact city index，the radial index in general，urban expansion is compact extension type，but during the period from 2007 to 2014，the urban expansion shifted toward the epitaxial extension type.Urban expansion appeared as compact extension type in the period from 2001 to 2007，but it turned to epitaxial extension type in the period from 2007 to 2014.The urban construction land in low elevation area reached to the saturated state，while high elevation area is a growing momentum，from the point of rolling degree analysis，construction lands are unable to expand to the larger area.

construction land expansion；mixed pixel decomposition；terrain gradient；large mountainous cities

F301.2;P208

A

1005-3409（2017）01-0232-07

2015-12-15

2016-02-01

重慶市社會科學規劃培育項目（2014PY63）；國家自然科學基金（41101155，41301351）；重慶市教委科技項目（KJ1400623）；教育部人文社科項目（12XJCZH005）；重慶市高等學校青年骨干教師資助計劃（2011）

呂志強（1978—），男，山東萊蕪人，博士，副教授，主要從事3S應用、環境遙感、資源環境經濟與農村可持續發展研究。E-mail：lv0726＠126.com