1、引用格式：王克晓，周蕊近 年重庆主城区不透水面变化与驱动力分析 山地学报，（）：，（）：山 地 学 报 文章编号：（）：收稿日期（）：；改回日期（）：基金项目（）：重庆市技术创新与应用发展专项（）（）作者简介（）：王克晓（），男，河南平顶山人，硕士，助理研究员，主要研究方向：资源环境遥感。（），：通讯作者（）：周蕊（），女，硕士，研究员，主要研究方向：农业信息化。（），：近 年重庆主城区不透水面变化与驱动力分析王克晓，周 蕊（重庆市农业科学院 农业科技信息研究所，重庆 ）摘要：城市不透水面反映城市建设用地的利用强度和承载功能变化。不同城市的不透水面时空演化特征和空间分异机制因地貌不同而有所差异
2、。重庆主城区是典型的“山城”城区，以往相关研究多关注城市建设用地，缺乏对城市不透水面的深入研究。把握重庆主城区不透水面的扩展演化过程及驱动力机制有助于理解其城市化发展进程，指导城区发展规划。本研究选取重庆主城区 年不同时期的不透水面数据，基于扩展变化测度指标、象限分析、景观格局指标以及地理探测器等方法，研究了不透水面时空变化特点及其驱动力因素。结果表明：（）重庆主城区不透水面扩展速度和扩展强度持续增加，年面积从 增加到 ，年均扩展速度 ，整体扩展强度指数达到 ；年以后，其扩展速度和强度均有所减缓，扩展强度下降为 。（）重庆市主城区不透水面的扩展方向主要为向北和向西发展，具有明显扩展优势的主要为
3、渝北区和沙坪坝区；在 年、年和 年三个时间段内，其扩展方式经历了内部填充外向扩展内部填充的转变，以东北和北向的扩展辐射范围最大为 ，其次西部辐射范围达到 。（）随着城市化进程发展，重庆主城区的不透水面斑块不断产生与融合，扩展过程仍处于不稳定状态，形状复杂度不断提升，斑块聚集度逐步上升。（）主城区的不透水面主要分布在 以下的缓坡区。主城区不透水面的空间分异是多因素共同作用的结果，人口密度、等社会经济因素是重庆主城区不透水面分布分异的主要因素，高程、坡度等自然因素在一定程度上影响着城区不透水面的空间分异，但影响力相对较弱。本研究结果可为重庆主城区城市化建设和发展规划提供科学依据。关键词：不透水面；
4、时空特征；驱动力；重庆中图分类号：文献标志码：城市不透水面的时空格局变化反映了一个城市的发展水平 ，是预测城市发展方向的风向标。明确城市不透水面的演变规律对揭示城市发展轨迹、完善用地结构、指导城市规划与建设具有重要意义 。随着遥感对地观测技术的发展，城市不透水面的遥感提取、时空变化及驱动机制等方面的研究已经有了理论成果。例如，针对不透水面的遥感提取技术与方法，有代表性的包括影像范数结合差值指数与阈值分割算法 、优化影像光谱指数和阈值的不透水面提取指数 以及多尺度特征融合的深度学习方法 ；在时空变化特征分析方面，引入不同密度等级的转移矩阵动态度和中位数中心 ，采用扩展强度和景观指数 ，结合不透水
5、面平均中心、标准差椭圆与景观格局指数 ，综合空间形态分析、等扇距分析及分形分析 ；在土地覆盖类型变化驱动力分析方面，主要是考虑自然环境、人口、经济、政策和城市生态环境因素与不透水面变化的关系 ，借助灰色关联分析 、地理加权回归 、结构方程模型 、地理探测器 等方法开展驱动力分析，以确定不同驱动因子对不同土地覆盖类型的影响程度。在实际应用方面，基于上述方法，学者们在多个城市已进行有益尝试：基于时序 影像获取，利用扩张速率与强度、核密度、中位数中心点指标分析莆田市城区多年不透水面其时空变化特征，并结合灰色关联度模型探究其相关驱动因素 ；利用光谱混合分解法对多期 影像开展西安市 和哈尔滨市中心城区
6、不透水面提取、时空变化特征以及驱动因素分析。不同城市的不透水面时空演化特征和空间分异机制会因所依托地貌不同而有所差异。重庆主城区是典型的“山城”地貌，加上长江、嘉陵江两大江流交汇于此，使其城市发展形态和不透水面的扩展变化呈现出独特性 。自 年重庆直辖，区域经济快速发展，主城区城市化进程加快，城镇建设用地和城市不透水面持续扩张。但已有研究 关注点多集中在城市建设用地整体的时空变化，没有从不透水面时空格局变化的角度来认识城市发展进程。探讨重庆主城区不透水面的扩展特点、时空格局变化、驱动力因素及其分布的空间分异，有助于深入了解建设用地利用强度和承载功能变化，理解城市发展路径，科学把握城市扩展演化过程
7、及驱动机制。本研究利用 年、年、年、年、年和 年不同年份的不透水面数据，借助扩展变化测度指标、象限分析以及景观指数方法，分析重庆市直辖后主城区近 年来不透水面的时空变化特点，并基于地理探测器模型量化不同驱动因素之间的交互作用及其对不透水面空间分布的贡献率，探讨主城区不透水面的时空格局演变及其驱动关系，为重庆主城区城市化建设和发展规划提供科学依据。研究区与数据源 研究区概况重庆市位于长江上游地区（、），主城区面积 （图 ），以“山城”地貌为主，分布有缙云山脉、中梁山脉、铜锣山脉和明月山脉四大山脉，长江、嘉陵江两大江流交汇于此，地势起伏较大，地形切割破碎。重庆主城区作为国家中心城市，是我国西南地区
8、经济发展的重要支撑点之一。图 重庆主城区位置图 ，数据源获取土地利用类型数据源采用重庆市 年、年、年、年、年和 年六个时间节点的 分辨率的中国逐年土地覆盖数据集（，）（：）。该数据集是基于 （）上可获得的 数据，根据反射光谱、光谱指数、物候和地理位置等时空特征，结合随机森林分类器得到分类结果，主要包括农田、草地、不透水面、荒地、水体、灌木、冰雪、森林和湿地九种土地利用类型。其中，不透水面的平均精度得分超过 ，总体精度达到 ，且与现山地学报 卷 第 期有的土地覆被专题产品表现出良好的一致性 。本研究先利用重庆主城区行政矢量边界对六期 数据集进行裁剪，后通过地类编码提取获得城区不透水面空间分布。为
9、分析不透水面面积与经济发展和人口增长的相关性，本研究还采用了 年重庆主城区历年的 和人口规模数据，数据来源于对应年份的年鉴统计资料（：）。为进一步分析其空间分异性，选取研究区的高程、坡度、人口密度和 等因子进行驱动力分析。其中，人口密度和 为 年数据，来源于中国科学院资源环境科学与数据中心（：），空间分辨率均为 ；坡度数据基于 的 数据计算获取；对应的不透水面分布数据通过 年的 土地覆盖数据集提取。研究方法 不透水面扩展变化测度指标为反映研究区域内不透水面随时间的扩展变化程度，分别选取扩展速度指数 和扩展强度指数 表示研究时段内平均每年增加的不透水面面积和增长率 ，计算公式如下：（）（）式中，
10、、分别为研究初期和末期的城区不透水面增长面积；为年份间隔。象限方位分析法为比较重庆主城区在不同方向上的不透水面变化差异，以重庆渝中区解放碑为中心，东西方向、南北方向分别为横轴和纵轴，按照八个方向绘制主城区不同时期的不透水面面积变化雷达方位图，分析主城区不同方向的不透水面变化情况。景观格局指数为了解主城区不透水面在空间尺度上的景观分布格局变化，选取包括斑块数量（）、斑块密度（）、最大斑块指数（）、斑块平均大小（）、总边界长度（）、平均边界密度（）、景观形状指数（）、斑块平均邻接指数（）和聚集度指数（）九个景观指标来反映的重庆主城区不透水面的空间分布格局情况。其中，、和 指标用于衡量不透水面斑块的
11、扩张与融合情况，、用于比较不透水面斑块整体形状的复杂程度，和 两个指标来衡量不透水面斑块间的联系性和聚集力 。上述指标的计算通过景观分析软件 完成，各指标计算方法及其生态学意义可参考该软件帮助文档。地理探测器地理探测器包括因子探测、交互作用探测、风险区探测和生态探测四部分，其原理和计算公式见文献 ，可用于探测城区不透水面扩张影响因素的差异。本研究基于其“因子探测”和“交互作用探测”功能，以不透水面积分布为探测因变量，以驱动因子为自变量，进行单因子贡献率及因子交互作用定量分析，探讨重庆主城区不透水面分布的空间分异特点，从而明确各因子对其空间分异的相对重要性。本研究选用高程、坡度、以及人口密度等因
12、素进行驱动力分析。利用自然断裂法将高程、人口密度、等数值型数据进行五级离散化，坡度数据按照 第三次全国国土调查技术规程（）进行五级离散化；然后通过间隔 的格网点对上述因子图层和 土地覆盖图层进行多值提取，并对格网点中的土地覆盖类型属性进行二值处理（代表不透水面，表示其他）；最后基于不透水类型格网点和各量化因子开展地理探测分析。结果与分析 不透水面面积动态扩展变化分析如表 所示，重庆主城区 年不透水面 面 积 呈 现 逐 渐 增 长 趋 势，总 面 积 共 增 加 ，增长了 倍，在主城区中面积占比提升至 ，年均扩展速度 ，整体扩展强度指数达到 。其中，年，伴随着直辖后城市化的发展，主城区不透水面
13、面积快速增长，扩展速度 ，扩展强度 ；年期间，随着城市化进程的加快，不透水面积扩展速度随之稳定增加，年扩展强度 ，近 年重庆主城区不透水面变化与驱动力分析表 研究区不透水面面积及其变化 年份不透水面积 所占比例 时间段扩展速度（）扩展强度 年 年 年 年 年 保持在 以上；年以后，重庆主城区城市建成区扩展步伐放缓，不透水面扩展速度随之减小，扩展强度指标下降为 。不透水面象限分布及拓展方向分析如图 和图 所示，重庆市主城区在不同时间图 年研究区不透水面变化雷达图 段内不透水面的扩展侧重方向变化不大，整体呈现向北和向西发展的趋势。年，不透水面的扩展方向集中在西南西西北方向，主要分布在九龙坡区、沙坪
14、坝区和渝北区等距离象限中心 范 围 的 近 距 离 地 带 内；年 和 年，主城区不透水面的主要扩展方向均在北西北西方向，以渝北区、北碚区和沙坪坝区为主，辐射范围扩大到距离中心 范围；年，不透水面扩展变化主要在东北、北、西、西南、东等五个方位表现为外向式扩展，并以东北和北向的扩展幅度最大，对应分布于渝北区和江北区，而图 研究区不透水面时空变化图 在西北方位表现为内部填充式扩展。年，除南向为外向扩展外，其他七个方位均为内部填充式扩展，面积变化最大的仍是渝北和沙坪坝两区。总体上看，重庆市主城区不透水面面积最大的三个辖区为渝北区、九龙坡区和沙坪坝区，但具有明显扩展优势的分别是位于北东北方向的渝北区和
15、西向的沙坪坝区，这一扩展方向与已有研究 关于重庆城市用地重心的变化方向基本一致。不透水面要素景观格局分析如表 所示，不透水面斑块数量（）、斑块密度（）、最 大 斑 块 指 数（）、斑 块 平 均 大 小（）指标逐年增大，说明期间不透水面斑块大量出现，且与周边斑块有不断融合的迹象，印证山地学报 卷 第 期表 年研究区不透水面景观格局指标变化 年份 个 （个 ）（个 ）了在大量新兴建设区出现的同时，城市不透水面斑块不断融合扩大的现象。总边界长度（）、平均边界密度（）以及景观形状指数（）变化趋势基本一致，呈现增加的态势，说明重庆主城区不透水面的形状伴随着不透水面斑块的融合，复杂度在逐渐上升。尤其 年
16、以后，、和 三个指标的上升幅度更加显著，说明此阶段不透水面的扩张并未稳定，主城区的城市发展仍处于快速发展时期。平均邻接指数（）和聚集度指数（）一致增大，说明 年重庆主城区不透水面之间的联系越来越紧密，不透水面间的聚集度逐步上升。不透水面扩展变化驱动力分析 不透水面变化驱动因素分析 自然条件重庆主城区被缙云山脉、中梁山脉、铜锣山脉和明月山脉等四大山脉切割，加上长江与嘉陵江两江交汇，严重限制了主城区建成区的连片发展，也成为影响重庆主城区不透水面扩展方向的主要原因。铜梁山山脉将主城区东部和南部分割开来，使得重庆主城区的地势呈南北走势。地势相对平坦的北面，成为城市扩张的理想方向，也是不透水面扩展的主要
17、方向。另外，位于主城区西侧的处于缙云山山脉与中梁山山脉之间的南北条状地势平坦区也是不透水面的主要扩展区域。基于地形分布指数 的 年重庆主城区不透水面在坡度梯度区间上的统计信息（表 ）显示：不透水面中有 分布于 ，分布于 ，在 内的不透水面占 。可见，主城区不透水面的分布优势主要集中在 以下的缓坡区。不透水面的分布指数在坡度梯度上总体呈现减少趋势，分布指数均大于 的区域主要集中在 以内的缓坡区。因此，坡度是影响重庆主城区不透水面扩展规模的主要影响因素。表 研究区不透水面坡度梯度分布指数 坡度梯度（）不透水面分布面积 占不透水总面积比例 分布指数 以上 经济发展经济发展是城市化的主要动力，不透水面
18、的扩展是城市化发展的一种重要体现。重庆市是我国西部唯一的直辖市。年以来，在国家“西部大开发”战略和“一带一路”等一系列政策的支持下，主城区经济快速发展，城市化进程明显加快，不透水面不断扩张。主城区 年的 数据与对应年份不透水面面积的相关性决定系数 为 （图 ），二者具有极强的相关性，说明经济的发展是重庆主城区不透水面扩展的重要推动力量。人口增长人口增长是城市化的重要表现，也是推动城市发展的重要推动力量，对城区不透水面的变化有重要影响 。在城市经济发展的吸引和刺激下，城市 ，近 年重庆主城区不透水面变化与驱动力分析图 不透水面面积与城区 和人口数量的关系 人口持续增长，带动住房、道路等基础设施以
19、及教育、医疗等公共服务设施的建设需求，从而导致城区不透水面面积的增加。主城区 年人口数据与对应年份不透水面面积的相关性决定系数为 （图 ），人口增长与不透水面扩展态势高度同步。重庆主城区人口的逐年增加有力推动了城市不透水面的扩展。规划政策 年，横跨北碚区、渝北区和江北区的重庆北部新区设立，形成了四个工业园、四个商务区和三个科技园的产业功能布局，重庆城市发展重心开始向北偏移，相应辖区的不透水面面积明显增加。年，两江新区批准成立，加快了重庆向北发展的序幕。同年，两江新区东北部的龙兴工业园和鱼复工业园相继开始建设。与此同时，年位于缙云山和中梁山山脉槽谷地带的重庆大学城片区批准建设，年沙坪坝区西部新城
20、正式成立，年市级重大文化产业项目落地西部新城，伴随着这一系列城市发展规划实施，如今的重庆大学城片区已成为西部地区重要的科学研究中心和科技创新基地。年，随着西部（重庆）科学城和成渝地区双城经济圈国家战略的启动，重庆向西发展势头正盛。综上可见，重庆主城区向北和向西发展的政策是主城区不透水面扩张方向的重要推动因素。驱动因子对不透水面分布的空间分异分析 单驱动因子的解释力分析由于政策因素难以在空间上进行量化，本研究选用高程、坡度、以及人口密度等因素进行驱动力分析。如表 所示，各驱动因子按照 统计量由大到小排序为人口密度、高程、坡度，且各因子对主城区不透水面分布的影响均通过显著性检验（）。从驱动因子对不
21、透水面的分布影响来看，人口密度的解释力最高为 ，明显占主导作用；其次是 ，其解释力为 ，是影响不透水面空间分异的重要因素；高程和坡度等自然因素的解释力分别为 和 ，一定程度上影响着不透水面的空间分异。因子探测结果表明，在重庆主城区社会经济因素对不透水面的空间分异影响程度较大，而高程与坡度等自然因素的影响相对较弱，但其因子探测结果能够通过显著性检验，这也反映了重庆主城区典型的“山城”地貌特色。表 研究区不透水面空间分异驱动因子探测结果 因子探测指标驱动因子高程坡度 人口密度 统计量 值 交互因子的解释力分析交互探测可以量化驱动因子共同作用对因变量的解释力。如表 所示，人口密度与坡度交互作用对不透
22、水面空间分异的影响最强，因子交互探测的 值最高为 ；交互程度对不透水面空间分异作山地学报 卷 第 期用的影响程度较高的还有人口密度高程（值 ）、人口密度 （值 ）、高程（值 ）、坡度（值 ），而高程坡度两个自然因子交互作用的 值为 。上述因子交互作用表明，人口密度与 等社会经济因素是影响重庆主城区不透水面分布的主要因素，而自然因素对其分布的影响却有限。同时，任意两个因子的交互作用下均表现出 （）（），（），表明上述因子两两交互作用的类型主要为交互增强型，说明引起主城区不透水面空间分异结果的并非单一影响因子，而是多因素共同作用形成。表 研究区不透水面空间分异驱动因子交互探测结果 驱动因子各驱动因
23、子交互作用 值高程坡度 人口密度高程 坡度 人口密度 讨论 山地城市不透水面扩展演化过程的独特性城市不透水面的扩展强度、方向和时空结构特征可以体现出城市的发展水平及其社会经济发展状况。自重庆市直辖后，重庆市主城区不透水面积的扩展速度和扩展强度持续增加。年以后，主城区不透水面扩展速度和强度均有所放缓，但 产值仍然保持在高速增长。这是响应国家经济内涵式发展要求，经济发展方式由城市规模和速度效益向城市质量和品质效益转变取得的成效。近 年来，重庆市主城区在不同时间段内不透水面的扩展侧重方向变化不大，但不同时期不透水面扩展方式有所不同，扩展方式也先后经历了内部填充外向扩展内部填充的转变过程。这种在扩展方
24、向上的差异性和方式上蔓延类型是与之前学者关于重庆主城区建设用地扩展方向 和山地城市建设用地扩展蔓延方式 的研究结果是一致的，而又与平原城市建设用地的蔓延方式 有所不同。重庆市主城区不透水面联系愈加紧密，斑块聚集度不断上升，反映出主城区城市建设用地在多中心模式形成的同时，也在有组织地扩张。山地地貌对城市不透水面空间分布形态的重塑自然条件对于城市建设用地的扩张起着关键作用。随着城镇化水平的提升，经济、人口条件的作用逐渐突显，控制着建设用地的扩张规模和方向等 。本研究基于地理探测法分析了高程、坡度等自然因素与人口密度、等社会经济因素对主城区不透水面分布的影响程度，其驱动因子的重要性与通过灰色关联分析
25、得出的重庆市主城区建设用地扩展驱动因子 的结果一致。在人口增长与经济发展两大推动因素的驱使和相关城市开发政策的指引下，受限于“山城”地貌和“两江”两大主要自然环境因素，主城区不透水面空间分布呈现出“山水之城”的独特形态。不透水面变化驱动机制对城市规划政策的指导 年以后，重庆主城区不透水面的扩展速度和扩展强度有所减弱，扩展模式主要表现为内部填充式，景观格局上表现出不透水面斑块的扩张、融合和聚集度的上升，产值和所吸纳的人口规模依然处于连续高速增长状态。结果说明粗放式城市开发模式不再是城市经济发展和人口城市化的唯一途径。在国家可持续发展要求下，规划政策和产业政策是引领城市未来发展的关键环节。当前，西
26、部（重庆）科学城连片规划建设高新产业基地，加紧引进高科技人才发展高新科技产业，正是开发城市未来经济发展潜力，提升城市未来竞争力的重要决策。结论（）近 年 来，重 庆 主 城 区 不 透 水 面 从 增加到 ，增大 倍，在主城区中的面积占比提升至 ，年均扩展速度 ，扩展强度指数为 。年以后，重庆主城区城市建成区扩展步伐放缓，不透水面扩展速度随之减小，扩展强度有所下降。（）近 年来，重庆主城区不透水面扩展方向整体呈现向北和向西发展趋势，具有明显扩展优势 ，近 年重庆主城区不透水面变化与驱动力分析的主要为渝北区和沙坪坝区。在 年、年和 年的三个时间段内，主城区不透水面的扩展方式经历了内部填充外向扩展
27、内部填充的转变过程。（）近 年来，随着城市化的不断发展，重庆主城区不透水面斑块不断经历新生与融合的周期过程，扩展仍不稳定，复杂度不断提升，城市发展处于快速发展期，联系愈加紧密，斑块聚集度逐步上升，城市凝聚力不断提高。（）重庆的“山城”地貌和“两江”是影响重庆主城区不透水面扩展方向的主要原因，坡度是影响重庆主城区不透水面扩展规模的主要限制因素，分布在 以下缓坡区的不透水面占 ；经济发展规模与人口数量是推动主城区不透水面扩展的主要驱动因素，历年 、人口规模与对应不透水面面积的相关性决定系数 分别为 和 。社会经济因素是主城区不透水面空间分异的主要影响因素，高程、坡度等自然因素的影响相对较弱，主城区
28、不透水面的空间分异是多因素共同作用的结果。规划政策对重庆主城区不透水面的扩张方向具有巨大推动作用。参考文献（）李益敏，杨舒婷，吴博闻，等昆明市呈贡区不透水面时空变化及驱动力分析 自然资源遥感，（）：，（）：刘美，杜国明，于凤荣，等哈尔滨城乡梯度建设用地结构变化及不透水面遥感监测分析 遥感技术与应用，（）：，（）：段潘，张飞，刘长江基于 的新疆典型城市不透水面提取及空间差异分析 遥感学报，（）：，（）：崔雯卓，林晨，陈曦，等半城市化地区不透水面遥感提取优化方法研究 地理与地理信息科学，（）：，（）：王晶，高帅，郭亮，等基于多尺度特征融合的 网络高分影像不透水面提取研究 遥感技术与应用，（）：，（
29、）：宋小琪，林志垒 年厦门市不透水面时空演变特征 亚热带资源与环境学报，（）：，（）：王姗姗，蒲英霞，李升峰，等利用长时序流域不透水面数据集探索城市增长模式 以秦淮河流域为例 遥感技术与应用，（）：，（）：杨丽媛，陈洋波 年深圳市不透水面信息提取及其时空格局变化分析 水土保持通报，（）：，（）：沈非，袁甲，黄薇薇，等基于地学信息图谱的合肥市城市扩展时空特征及驱动力分析 长江流域资源与环境，（）：，（）：尹振南，周廷刚，李洪忠，等粤港澳大湾区不透水面变化及城市扩展研究 地理与地理信息科学，（）：，（）：周鹏，谢元礼，高志远，等西安市不透水面的变化及其驱动力 水土保持通报，（）：山地学报 卷 第
30、期 ，（）：吴玉刚华北平原农地利用变化及驱动力分析 中国农业资源与区划，（）：，（）：赵昊，刘平辉鄱阳湖生态经济区土地城镇化时空格局特征与驱动因素分析 江西农业学报，（）：，（）：许冬凡，张加龙，刘钱威，等 西南省会城市扩展时空变化及驱动因素分析 科学技术与工程，（）：，（）：赵小风，刘梦丽，王柏源基于地理探测器的城市建设用地规模时空差异及影响因素分析 中国土地科学，（）：，（）：许承权，白雅晴，田玉刚，等 基于 指数的 年莆田市城区不透水面时空演变分析 地域研究与开发，（）：，（）：，（）：国洪磊，周启刚基于 和 的城市用地扩展的时空特征研究 以重庆主城区为例 重庆工商大学学报（自然科学版）
31、，（）：，（），（）：李建辉，周启刚，陈倩，等基于 和 的建设用地变化研究 以重庆市主城九区为例 重庆工商大学学报（自然科学版），（）：，（），（）：嵇涛，杨华，何太蓉重庆主城区建设用地扩展的时空特征及驱动因子分析 长江流域资源与环境，（）：，（）：吕志强，代富强，周启刚，等 近 年来典型山地城市建设用地蔓延的时空特征分析 水土保持研究，（）：，（）：，（）：程结海，胡盼，袁占良基于夜光遥感数据的中原城市群建成区时空变化分析 测绘通报，（）：，（）：，（）：中华人民共和国自然资源部第三次全国国土调查技术规程：北京：地质出版社，：，李靖，廖和平，蔡进，等山地城市边缘土地利用格局与变化的地形梯度分布特征分析 以重庆市巴南区为例 长江流域资源与环境，（）：，近 年重庆主城区不透水面变化与驱动力分析 ，：，（）：蔡恩香，寇紫薇，孟珂，等 年郑州城市扩张时空特征分析 河南农业大学学报，（）：，（）：周棱楠，潘安，罗芳 年成都市建设用地扩张时空演变及驱动力分析 西南林业大学学报（自然科学），（）：，（），（）：，（，）：，；，；，（）；，（）；，（），山地学报 卷 第 期 ，（），；，：；（责任编辑钟雨倩），近 年重庆主城区不透水面变化与驱动力分析
