基于RS和GIS的城市扩张研究.doc

基于RS和GIS的城市扩张研究 ——以福州市为例 摘要 采用1989年、2006年和2010不同历史阶段的遥感影像数据，借助于ArcGIS 软件和ERDAS软件，计算福州市建筑指数（NDBI）和植被指数（NDVI），通过两个指数相减运算，较精确的提取城镇用地信息。采用空间形态分析、地统计分析、等扇分析、景观测度及分形理论等研究方法,对福州市建成区及其周边城乡建设用地规模与形态的发展演变过程、特征及其原因进行研究，结果表明：福州市城乡建设用地由快速的、蔓延式的无序扩张向城市规划引导下的适度、有序扩张转变，城市规划对建设用地布局的引导作用增强。道路交通尤其是高速公路建设对城乡建设用地布局的影响较为明显。目前仍面临的突出问题是，城乡建设用地加速扩展的同时，耕地被各项建设用地占用，近年来缩减幅度大幅增加。 关键词: 城市扩张；演变特征；福州市 近年来，遥感与GIS 技术越来越被广泛地应用于土地利用、土地覆被变化、城市用地动态变化、城市规划等研究中，提高了城市用地发展研究的客观性、实时性和准确性，为城市规划的编制、监测和管理提供了有效途径和手段。采用所获取的1989年、2006 年和2010 年三期遥感影像数据，利用ArcGIS 软件，对福州市区及其周边城乡建设用地的发展演变过程进行对比分析，研究城乡建设用地在较长时间尺度内的演变特征及动力，以期能为科学合理地进行城乡规划、控制城乡建设用地规模提供一定的技术支撑。 1 研究区概况与研究方法 1.1研究区概况 福州市位于福建省东部的闽江下游沿海地区，东经118°08′～120°31′ ，北纬25°15′～26°29′之间。东濒东海，与台湾隔海相望，西邻三明市，和南平市，北接宁德市区，南连莆田市。 福州市面积为12152.5平方公里，海岸线曲折、漫长，大陆海岸线长度为920千米，居全省第二，海岛海岸线长度为390千米，居全省第一，福州海域宽阔，多港湾，岛屿。 福州市属亚热带海洋性季风气候，常年温暖湿润，雨量充沛，四季常青，全年无霜期达32天，年均气温在16～20°C之间，年均降水量为900～2100mm。 福州市是中国市场化程度和对外开放度较高的地区之一，改革开放以来特别是近年来经济社会持续快速协调健康发展，已初具经济繁荣、科技发达，设施完善、环境优美的现代都市风貌。 1.2 数据来源 选取福州市1989年Landsat TM与2006、2010年Landsat ETM遥感影像为研究数据。 1.3 研究方法 城市用地信息提取 城市建成区的提取是研究城市扩展的重要方面，城市建设用地的遥感信息提取通常采用人机交互的手工方式，自动化程度不高，许多方法如果不做后续处理，精度无法满足要求。近年来，基于指数的建筑用地提取技术相继提出，归一化的差值建筑覆盖指数NDBI、基于指数的-建筑用地指数IBI（Xu，2007）等，本文用NDBI-NDVI来表示城市的用地情况，公式 NDBI=(B5-B4)/(B5+B4) 式中，NDBI代表归一化的建筑覆盖指数，B4、B5分别为Landsat TM或ETM的第4、5波段的反射率。 NDVI=(B4-B3)/(B4+B3) 式中，NDVI代表归一化的植被指数，B4、B3分别为Landsat TM或ETM的第4、3波段的反射率。 一般情况下只要用NDBI就能反映出城市的建筑用地信息，但为了除掉植被的干扰，提高精度，将NDVI减去，得到研究区的建筑信息（图1），最后分别对各相关年份的地图进行矢量化，并结合影像判读来获得各年份城市建成区的矢量边界图，然后用其对所提取的建筑用地信息进行掩膜运算，获得了各年份的城市建成区范围（图2）。 图1 研究区不同年份的NDBI-NDVI 图2 研究区不同年份的建设用地分布 图3 1989-2010年城市扩张图 采用GIS 的空间分析方法对怀化城市空间形态进行定量测度，包括扩展强度指数、紧凑度指数、分维数、形状指数、等扇分析法。 （1） 扩展强度指数 计算城市扩展强度，比较不同时期扩展强弱和快慢，研究怀化城市扩展规模特征本文用城市年均增加面积与基年面积之比来表示扩展强度，公式为 AGR=((UA n+i - UA i )/n UA i) ?100% 式中：AGR 为城市年均扩展强度指数；UA n+i和UA i ，分别代表n+ i年和i年的城市面积； n 是以年为单位的时间。 （2）紧凑度指数 引入紧凑度公式BCI =2/ P，研究怀化城市扩展的空间形态变化特征。式中：BCI 表示紧凑度指数；A 表示建成区面积；P 指城市建成区轮廓的周长。 （3）分维数 分形维数表示一个集合在空间上的占有程度，维数越大，空间区域内含有分形单元的机会也就越大。常用的分维有：Hausdorff 维数、盒维数、相似维数、信息维数、关联维