遥感见证城市扩张

“经纬石旁话遥测”是《武汉大学学报·信息科学版》的官方订阅号

“经纬”石是武汉大学星湖园旁由宁津生院士亲笔题词的一块石头

本文改编自学术论文《2000-2016年中国地级以上城市高精度城区边界遥感提取及时空扩展分析》

发表于《武汉大学学报·信息科学版》2018年12期 “测绘遥感高层论坛”专刊

宁晓刚1,3 王浩1 张翰超1,2,3  刘娅菲1 庞博1  郝铭辉1

1. 中国测绘科学研究院

2. 武汉大学遥感信息工程学院

3. 城市空间信息工程北京市重点实验室

城区边界和城区面积是城镇化的重要表征和扩展分析的基础。然而城区边界存在概念和提取标准不统一、精度较低、可比性较差等问题。

为此，提出了基于高分辨率影像和地理信息资料辅助的城区半自动化提取方法，充分利用高分辨率影像上的城市景观特征、先验地理信息知识和一系列标准规则，以得到精度高、一致性强的数据。

以中国337个地级以上城市为研究区，采用该方法得到了2000年、2005年、2010年、2016年4期城区边界成果，并开展了城区时空扩展及用地效率等相关分析。

结果表明：①16年间城区扩展迅速，城区主要集中分布在东部和中部，东西部地区差异大；②城市用地效率与城镇化发展水平显著相关，城区扩展以外延型为主；③大多城市城区扩展超前于人口增长，少量城市城区扩展滞后于人口增长；④城区扩展以占用耕地为主。

========以下为全文，建议收藏后阅读，大约需要15 min=======

随着城镇化发展，城镇人口不断增多，伴随而来的是城市区域的迅速扩张。

特别是2000年以来，中国的城镇化发展进入了新的阶段。中国“十三五”规划、国家新型城镇化规划(2014—2020)等提出要合理控制城镇开发边界、保护基本农田，通过优化城市内部空间结构，促进城市紧凑集约、高效绿色发展，以解决城市盲目扩张，土地城镇化快于人口城镇化，建设用地粗放低效，城镇化与工业化不匹配，城镇空间分布和规模结构不合理，资源环境承载力不足，生态环境问题突出等诸多问题。

合理控制城镇增长边界，优化城市内部空间结构，促进城市紧凑发展已经成为城市可持续发展的内在要求。开展城镇化监测与分析，对于准确把握城镇化进程，科学开展城市规划和建设，促进城市可持续发展有着极为重要的作用。

城镇化可从人口城镇化和土地城镇化两方面进行度量。人口城镇化可采用城市化率(城市人口与总人口的比率)表示，土地城镇化可采用非农化的城区面积表示。相比于人口城镇化产生的人口转移这一难以度量的指标，土地城镇化即城市空间扩展是城镇化作用于地理空间的直接结果，更容易被客观度量。

遥感技术具有真实、客观、现势性强、成本低等优势，已成为土地城镇化监测的重要技术手段。

城镇化监测首要任务即确定城区面积或城区边界，然而城区或城区边界相关概念众多，这些定义有着不同的出发点和特点。许多研究将土地利用分类结果中的建设用地或不透水面作为城区。还有一些研究人员综合考虑城市用地的功能定义城区。

城镇化监测采用的影像数据源众多。早期，MODIS影像是城市化监测的主要数据源。随后，Landsat系列等中分辨率卫星影像被广泛使用。随着越来越多的高分辨率卫星发射升空，许多研究人员尝试利用高分辨率遥感影像开展城区提取，但由于其自动化提取的能力限制，多用于单个城市尺度。另外，在获取城区范围的基础上，也有不少学者进行了城市扩展模式和特征的分析。

上述研究为城镇化监测与分析所需的城区边界提取提供了丰富的方法，但仍存在不足。主要体现在：城区边界提取的方法和标准不统一，导致城市和区域间城区提取精度和面积比较困难；中低分辨率影像得到的城区边界提取精度较低；高分辨率影像的城区边界提取算法难以适应于全国尺度。

为解决上述问题，本文提出了一种计算机自动识别和人工解译相结合的标准化城区边界半自动化提取方法，以中国337个地级以上城市为监测区，利用高分辨率影像和一致的提取标准，获得了2000年、2005年、2010年、2016年4期的城区边界；城区面积被用于土地城镇化表征，并结合其他专题数据对全国地级以上城市的时空扩展特征进行分析。

1

城区数据和提取方法

选择中国337个地级以上城市(包括地级市、自治州、地区、盟，不含港澳台)作为研究区。选取其市辖区(没有市辖区的选择政府所在地的县)作为城区提取范围，存在行政区划调整的城市以2016年市辖区为准。

根据中国经济区域划分为东部、中部、西部和东北4大地区。根据《第一财经周刊》提出的城市分级结果，将地级以上城市分为1线城市、新1线城市、2线城市、3线城市、4线城市、5线城市(见图 1)。

地级以上城市是全国经济、政治和人口集聚中心，其城镇化能够代表本地区的最高水平。21世纪初期，地级以上城市经济社会快速发展，城区扩展迅速，在中国城镇化空间扩展监测与分析中具有代表性。

图  1  全国337个地级以上城市及地域分布

数据准备和处理

收集高分辨率遥感影像、人口数据、国内生产总值(gross domestic product, GDP)数据、1：10万比例尺土地利用数据、全国第一次地理国情普查成果数据、2016年基础性地理国情监测成果数据，用于全国地级以上城市城区边界提取和城市空间扩展分析。

其中，高分辨率遥感影像主要采用优于2m分辨率影像，覆盖整个研究区域，具体情况见表 1，表 1中数据类型按照使用的优先顺序进行排列，影像时相多为当年度6月份左右，来源于原国家测绘地理信息局和各省级测绘地理信息部门。坐标系均转化到CGCS2000大地坐标系。

人口和GDP数据来源于各年份的中国统计年鉴；1：10万比例尺土地利用数据来源于中国科学院资源环境科学数据中心，土地利用合并成耕地、建设用地、林草地、水域和未利用地；全国第一次地理国情普查成果数据和2016年基础性地理国情监测数据来源于原国家测绘地理信息局，包括地表覆盖、城市道路要素、地名地址、乡镇和县级行政区划界线，为城区边界提取提供了很好的参考。

高分辨率遥感影像经过正射纠正、影像融合、影像镶嵌、裁切等处理步骤，所有专题数据经过几何纠正重投影，以保障与影像的匹配。

表１ 本文用到的主要遥感影像数据

城区边界提取

城区边界采用计算机自动提取和人工解译修正相结合的提取方法。

采用融合直角点和直角边特征对高分辨率遥感影像自动提取居民点多边形；将居民点多边形与区政府所在地叠加，区政府所在地落到的集中连片居民点多边形作为初始城区边界；参照地理国情监测的分类体系和标准，按照城区边界提取原则和规则对初始城区边界人工目视解译修正和编辑，以获取最终的城区边界(见图 2)。最终的城区边界包括了中心城区和飞地型城区。

初始城区边界的修正原则和规则如下。

图  2  城区边界提取流程图

1) 行政区划边界限定原则。城区必须位于城市行政区划范围内，不能超过行政区划边界。

2) 城区边界线勾绘原则。城区边界线优先沿道路、河流等线性地物的边界勾绘，不能跨越房屋建筑(区)和构筑物。城区边界线内侧不能有耕地存在。当城区边界落到河岸时，如果河岸周边为城市建设用地绿化景观，则城区边界线沿河流的高水位线勾绘，连同绿化部分归入城区；否则沿河岸附近的道路勾绘。线性地物界线不明显时，城区边界优先沿具有明显界线的完整地块勾绘，如院落或围墙。若地块边界不明显或整个地块内未建成面积比例较大，可按实际建设用地边界勾绘。

3) 中心城区判定原则。中心城区即中心集中连片区，通过道路一直向外延伸，直到具有城市景观的建设用地与中心集中连片区边缘的距离超过50 m或不再具有城市景观。

4) 城市景观判定原则。城市景观主要包括城市房屋建筑区、城市道路、城市绿地、城市广场、停车场、体育场等景观，其中由城市道路围成的街区是最主要的城市景观特征。

5) 飞地型城区判定原则。飞地型城区即飞地型集中连片区，是城区的重要组成部分。其具有以下特征：①通过道路与中心城区联系紧密；②具有城市景观特征；③面积大于60 ha且至少包含以下地理单元：城市级的街道办事处、行政部门、工矿区、开发区、科研单位、大专院校等。

6) 城区边界逆向提取原则。首先，完成最新时相的城区边界提取；其次，按现势性从新到旧的顺序，采用逆增量更新技术获得城区边界变化部分；最后，将后时相边界与城区边界变化部分融合得到前时相的城区边界，以保证多期边界之间的拓扑准确。

计算机自动提取和人工解译修正相结合的提取方法大大节省了人工目视和手工勾绘获得初始城区的效率。经测试，自动识别获取初始城区边界相比人工目视手工解译效率提升近2倍。连同初始城区边界解译修正这一环节，综合提升效率30%。图 3为本研究的城区内部典型地物示例。

图３ 城区地物解译标志示例

城市空间扩展分析指标

采用城区扩展速度、扩展强度、地均GDP、城市用地增长弹性系数4个指标，从城区时空扩展过程、城区土地利用效率变化、城区扩展协调性、城区扩展区域土地利用变化4个方面对城区时空扩展进行分析。

1）城区扩展速度。城区扩展速度V_i为某一 时间段内某城市的城区面积年增长速度，反映了单位时间内不同城区扩展快慢的绝对差异。

式中，ΔU_i,j 为 j 时段 i 单元的城区扩展面积；Δt_i 为 j 时段的时间跨度。

2）城区扩展强度。扩展强度 N_i为某一时间段内某城市城区面积相对于基期城区的年扩展比例，反映了单位时间内不同城区扩展快慢的相对差异。

式中，ΔUi,j 为就时段 i 单元的城区扩展面积；Δti 为 j 时段的时间跨度；M_i为基期 i 单元的城区面积。

3）地均GDP。地均GDP反映单位城区面积上产出的GDP效益。 

4）城市用地增长弹性系数。城市用地增长弹性系数 K 是某时期内城区面积增长率与城区人口增长率之比，是测定城区扩展是否协调的重要指标。国内学者通常用城市用地增长弹性系数来衡量城市扩张与人口增长之间的协调关系，认为城市用地增长弹性系数为1.12比较合理。

式中，R_G 为城区面积年均增长率；R_P 为城市人口 年均增长率；A_t 表示某时期末的城区面积；A₀ 表 示某时期初的城区面积；Δt表示时间跨度。

2

城区边界提取精度分析

本研究提取了全国337个地级以上城市的城区边界矢量成果，其可靠性和高精度主要体现在以下3个方面。

1) 技术方法保障提取结果的高精度。本研究提取城区边界首先采用融合直角点和直角边特征的居民点自动提取方法，保证了提取结果的一致性。另外，提取原则和标准严格、规范，可操作性强。解译人员需严格按照6项规范化的提取原则和具体技术规定操作，能够保证城区边界提取结果的可靠性。

2) 高分辨率遥感影像作为主体数据保障提取结果的高精度。采用的遥感影像分辨率高，能够有效避免误分、漏分现象。已有研究结果表明，分辨率对于城区边界提取精度具有重要的影响，高分辨率影像相比MODIS和Landsat影像，能够有效避免城区边界的误分、漏分，且能够将城区与城乡结合部、乡镇建设用地等区分开来。

3) 高质量的辅助数据保障提取结果的高精度。难点区域有真实可靠的辅助数据支持，保障提取成果的准确性。城区边界提取最大的难点在于集中连片、城市景观和飞地型城区的判断。本研究通过参考全国第一次地理国情普查成果和2016年基础性地理国情监测成果数据中的城镇功能单元数据、开发区数据、保税区数据、城市道路数据等，很好地解决了上述城区3大特征的判断。

3

城区时空扩展分析

城区时空扩展过程

2016年全国地级以上城市的城区总面积达到了35637.2

km²，占全国行政辖区总面积的2.0%，相比于2000年增长了122.6%。从城市、省份、城市群、经济区4个尺度对全国地级以上城市的城区面积时空扩展过程进行分析。

按城市级别分析，城区平均面积由高到低依次为1线城市、新1线城市、2线城市、3线城市、4线城市、5线城市。扩展面积最大的3个城市依次是上海、深圳、天津。从城区面积占比分析，上海市的城区面积占比最大，其次是天津市。2000—2016年，新1线城市扩展面积最大，达4820.0 km2, 2线城市与3线城市次之。1线城市16年间扩展速度为119.9 km2/a，5线城市扩展速度为140.0 km2/a(见图 4)。

图４ 全国地级以上城市2000-2016年城区面积空间分布图

按省份分析，2016年城市城区面积最大的前三甲省份分别是江苏、山东和广东，城区面积最小的3个省份是西藏、青海和海南。海南省城区面积小，但其城区面积占比高。

2000—2016年时间段中，城市扩展面积和扩展速度最大的前三甲省份分别是江苏、山东、广东，扩展面积和扩展速度最小的3个省份分别为西藏、青海、海南。扩展面积最大的江苏和扩展面积最小的西藏两者扩展速度相差50余倍。

扩展强度最大的前三甲省份依次为宁夏、云南、江西，扩展强度分别为16.0%、14.8%、14.0%。扩展强度较小的后3个省市为北京、黑龙江、西藏，扩展强度分别为2.1%、3.1%、3.6%。

按不同城市群分析，长三角、珠三角、京津冀、哈长、成渝5大城市群地级以上城市的城区面积呈阶段递增。

长三角城市群城区面积在5大城市群中增长最为迅速，且城区面积占比涨幅较大。珠三角城市群与长三角城市群类似，城区面积增长迅速且城区面积占比涨幅大。长三角和珠三角城市群城区面积占比远高于其他城市群，表明该两个城市群的城镇化程度远高于其他地区。京津冀地区的城区面积增长最为缓慢，城区面积占比几乎保持不变，且在5大城市群中最低，表明京津冀城市群土地城镇化综合水平低，且城镇化发展不足，亟需京津冀协同发展战略来带动。哈长城市群的城区面积增长速度和城区面积占比仅高于京津冀城市群，其城区总面积在5个城市群中最少。

按不同经济区分析，东北、东部、中部、西部4个地区的城区面积均呈阶段递增。

2016年，东部地区的城市城区面积最大，东北地区城区面积最小。城区面积占比方面，2016年东部地区最大为3.83%，西部地区最小为0.94%，表明西部地区的土地城镇化水平在全国4大区域仍处于最低。2000—2016年，各地区扩展面积与扩展速率从高到低依次为东部地区、西部地区、中部地区、东北地区(见图 5)。东部地区扩展面积为10224.5 km2，东北地区扩展面积为1617.7 km2。2000—2016年，从不同经济区分析，各地区扩展强度由高到低依次为西部地区、东部地区、东北地区、中部地区。西部地区扩展强度达9.01%，中部地区扩展强度仅为1.21%，面积占比增长最为缓慢。

图  5  全国地级以上城市2000—2016年城区扩展面积空间分布图

城区土地利用效率变化

对城区土地利用的经济效率进行分析。2000—2016年，全国各省份地均GDP，即单位城区面积GDP呈持续增长趋势。2016年, 地均GDP排名前三甲的省市依次为北京、广东、吉林，排名后三位的省份依次是青海、西藏、海南。排名第一的北京和排名最后的青海地均GDP相差近200倍。

按经济区分析，2000—2016年各地区地均GDP呈递增趋势。中部地区的地均GDP由2000年4大地区中的最后一位攀升至2016年的第二位。东部地区在各阶段地均GDP均处于全国领先。

按城市级别分析，1线城市的地均GDP在各阶段均领先全国其他城市，5线城市的地均GDP在各阶段均落后于全国其他城市，这与其实际城镇化发展水平一致。

采用凸壳原理对337个地级以上城市的扩展类型进行判断并分析。如图 6所示, 2000—2016年，填充型城市数量占比为26.7%，包括太原市、佳木斯市、宜昌市、白山市、盘锦市、延边朝鲜族自治州、十堰市、辽阳市、兰州市、枣庄市、大同市、锦州市、宝鸡市、洛阳市等90个城市。外延型城市数量占比为73.3%，包括吐鲁番地区、阜阳市、博尔塔拉蒙古自治州、北京市、天津市、惠州市、防城港市、吴忠市、六安市、鹤壁市、清远市、厦门市、中卫市等247个城市。2000—2016年, 全国地级以上城市空间扩展以外延型为主，城区形态趋于分散，城区土地利用效率有所下降(见图 7，其中新疆的城市用地增长弹性系数真实值是112.77，因显示原因表示为16.00)。

图 6 2000-2016年全国地级以上城市扩展类型空间分布图

图7 全国各省市2000-2016年地级以上城市用地增长弹性系数图

城区扩展协调性

2000—2016年，新疆、重庆、浙江、江苏、湖北、上海、辽宁、天津、安徽、江西、云南、宁夏、福建、内蒙、广西、甘肃、山西、贵州、湖南、四川、陕西、黑龙江、河南、吉林共24个省份的城市用地增长弹性系数远高于最优值1.12，表明这些省市城区扩展超前于人口增长。海南、青海、河北、西藏、北京共5个省市的城市用地增长弹性系数远低于最优值1.12，表明这些省市城区扩展滞后于人口增长。仅广东、山东两个省市的城市用地增长弹性系数接近最优值1.12。新疆、重庆的城市用地增长弹性系数分别为112.77和12.33，表明这两个省市16年间城区扩展远远高于人口增长速度。

城区扩展区域土地利用变化

2000—2016年，全国地级以上城市空间扩展占用最多的是耕地，占扩展总面积的56.4%。其次是建设用地，占扩展总面积的23.9%；水域及未利用地占用面积最少，占比为9.7%。

全国城区扩展以占用耕地为主的省份达28个。剩余的3个省市中，北京占用建设用地比例为50.8%，比占用耕地面积略高6.6%，表明北京近16年的扩展以合并农村居民点等建设用地和占用耕地并存；广东以占用水域及未利用地为主，占比为44.1%；西藏以占用林草地为主，占比达56.1%。

2000—2016年，城区扩展占用耕地面积排前10的省市为浙江、重庆、安徽、河南、吉林、江苏、上海、宁夏、山东、山西；占用林草地面积排名前10的省市为西藏、湖南、海南、新疆、内蒙、贵州、江西、云南、广西、福建；占用建设用地面积排名前10的省市为北京、陕西、四川、广西、河北、青海、海南、辽宁、福建、山西；占用水域及未利用地面积排名前10的省市有广东、甘肃、天津、湖北、贵州、新疆、黑龙江、福建、宁夏、内蒙(见图 8)。

图8 2000-2016年全国各省市空间扩展占用土地比例图

4

结语

本文提出了一种基于高分辨率遥感影像和地理信息资料辅助的计算机自动识别和人工解译相结合的标准化城区边界半自动化提取方法，获得了中国337个地级以上城市的2000年、2005年、2010年、2016年4期的高精度的城区边界。利用本研究得到的城区边界与专题数据相结合，开展了中国地级以上城市的城市时空扩展分析，结论如下：

①中国城镇化进程在2000-2016年间发展迅速，城区主要集中分布在东部和中部，东西部地区差异大。

②中国地级以上城市用地效率与城镇化发展水平显著相关，城市用地效率逐步上升。

③大部分城市城区扩展超前于人口增长，少量城市城区扩展滞后于人口增长。

④新增城区侵占耕地问题严重，占用建设用地也有较大比重。

进一步的研究工作包括：

①优化城区边界提取算法，将人工解译规则转化为机器学习算法，以提升城区边界提取效率。

②探讨中国地级以上城市的城区边界数据应用于政府部门和民众的途径和方法。

作者简介 

宁晓刚   博士, 研究员, 主要从事城市地理国情监测和遥感应用研究。[e-mail]

引用格式 

宁晓刚, 王浩, 张翰超, 刘娅菲, 庞博, 郝铭辉. 2000-2016年中国地级以上城市高精度城区边界遥感提取及时空扩展分析[J].武汉大学学报 · 信息科学版, 2018，43(12): 1916-1926