一种“流空间”视角下的城镇开发边界模拟预测方法与流程

1.本发明涉及一种融合城市间“流空间”和遥感数据的城镇开发边界模拟和预测技术方法。


背景技术：

2.长久以来以规划手段在城市周围划定限制开发边界一直是城市管理者用以遏制城市蔓延的用地管理政策之一，2019年6月，自然资源部发布《城镇开发边界划定指南(试行)》，并将城镇开发边界定义为“在一定时期内因城镇发展需要，可以集中进行城镇开发建设，重点完善城镇功能的区域边界。城镇开发边界内可分为城镇集中建设区、城镇弹性发展区和特别用途区。”城镇开发边界既是一条“技术线”，又是一条“政策线”，是上级政府指导和约束下级政府国土空间开发行为的政策工具。以往的城镇开发边界划定面临着众多挑战，比如不同层级事权管理体制下的划定深度差异和传导机制效率问题。
3.从城镇开发边界的技术性而言，城镇开发边界划定的主要方法可分为质性法、指标体系法和模型法三种方法。其中质性法主要基于人口规模、设施成本、土地利用现状和环境敏感区等预测未来城镇用地面积，比如弗雷划定法；在指标法方面，不同组织和部门为了进一步科学划定城镇开发边界制定了不同的指标体系系统。而在模型法方面，近年来，随着遥感(remote sensing)和地理信息系统(geographic information system)的广泛应用，lbs数据、元胞自动机、多智能体、蚁群优化算法(ant colony optimization)等场景分析模拟技术均被用于进行城市增长模拟，此类城市增长模拟方法大多是基于城市发展的历史遥感影像图获取城市用地增长速度、强度和发展方向等，通过机器学习算法模拟不同情境下城市未来扩张状况。
4.元胞自动机模型(cellular automaton,ca)是从技术上模拟城镇开发边界的主要工具模型，ca模型能够适应不同管控逻辑，从城市“正向增长”角度和“底线控制”角度出发划定城镇开发边界。ca模型可以将自然因素、社会经济因素等多个方面城市空间扩张的驱动因素和限制因素纳入考虑。其中自然因素衡量指标主要包括到河流水系距离、地形地势(坡度、坡向)、生态资源、土地适宜性；社会经济因素包括到城市重要节点距离(包括城镇中心、商业中心、)、到城市主要道路距离、到铁路距离、邻域开发强度、土地利用政策。在传统的ca模型的基础上也不断衍生出新的技术方法，如sleuth模型、clue
‑
s模型、ann神经网络模型和flus等。
5.但在实际应用上，由于城市增长具有高度的复杂性，传统的ca模拟所碰到的最大问题是如何确定模型的结构和参数。城市空间扩张形态除了受到内生因素驱动之外，也同时会受到周边城市的外生因素影响，包括城镇体系结构。“外生性”的区域和城镇间因素对城市空间扩张的影响作用逐渐加强，但传统元胞自动机模型聚焦于单一研究对象的“内生性”影响因素。城镇间人口、生产和交通等要素的流动使得将某一特定区域视为一个研究整体时，容易忽略区域中不同空间单元的空间异质性差异，比如不同城市的城市扩张驱动力或限制因素权重存在一定差异。


技术实现要素：

6.本发明针对现有通过元胞自动机模型模拟和预测城镇开发边界的不足，提供一种结合城镇“流空间”的城镇开发边界模拟和预测方法，包括以下步骤：
7.(1)运用研究区域范围内的遥感影像和城镇建设用地空间扩散驱动因素数据，通过logit模型计算各空间因素影响权重，用于提供步骤(5)，具体计算公式如下：
[0008][0009]
式中p
logit
为土地利用类型变更状态，其中p
logit
取值1或0，1为非城镇建设用地转换为城镇建设用地，0则表示未变更，x为城镇建设用地空间扩散驱动因素，β为空间因素影响系数，k为城镇建设用地空间扩散驱动因素数量
[0010]
(2)计算研究区域范围内限制土地利用类型转变因素，用于提供步骤(5)，如河流水域和双评价因素等，创建变量p
restriction
，取值为0(禁止转变为城镇建设用地)或1(允许转变为城镇建设用地)。
[0011]
(3)运用重力模型对城市流进行计算，具体公式如下：
[0012][0013][0014]
式中i
i
和i
j
分别表示城镇i和城镇j的流强度，d
ij
为城镇i,j的距离，α为参数,i
flow(x,y)
为城镇i中元胞xy的流强度，d
i(xy)
为元胞xy到城镇i中心距离，d
j(xy)
为元胞xy到城镇j中心距离；
[0015]
(4)根据步骤(3)计算获得的城镇间流强度，构建城镇等级体系，形成不同等级城镇建设用地分配指标体系，用于提供步骤(5)，具体公式如下：
[0016]
p
threshold
＝p*k
i
[0017]
式中p
threshold
为“流空间”视角下城镇i的城镇建设用地转换概率，取值为0
‑
1，p为研究范围内基准城镇建设用地转换概率，k
i
为城镇体系中城镇i的权重。
[0018]
(5)制定约束性元胞自动机模型转换规则，模拟和预测未来城镇开发边界，具体公式如下：
[0019][0020][0021]
元胞自动机约束规则如下：
[0022]
if any cell is not already developedto urban
[0023]
&if
[0024]
then cell i is developed with some other probability
[0025]
为在t时刻，栅格ij周围n*n的邻域对其影响作用，为栅格ij在t时刻是城镇建设用地的转换概率，n为每次迭代中允许用地类型发生转换的栅格数量。
[0026]
本发明的优点是:
[0027]
(1)通过将“流”强度纳入元胞自动机模型中，相较于传统logit
‑
ca模型，可将城镇网络纳入城镇开发边界模拟和预测中；
[0028]
(2)解决以往元胞自动机基于历史数据预测未来情景的问题，通过拆解用地转换概率公式，规划师可将城镇结构等主观规划意图纳入城镇空间扩张预测中；
[0029]
(3)方法适用性、延展性好，算法基于主流软件可应用于元胞自动机模型的众多领域和具体问题，亦可嵌入现有城镇开发边界模拟和预测软件中，如flus和clue
‑
s等。
附图说明
[0030]
图1：本发明的流程示意图
[0031]
图2：西安地表覆盖遥感数据(左2010，右2020)
[0032]
图3：西安城镇空间扩张驱动力因素
[0033]
图4：西安各区县间人流强度
[0034]
图5：西安限建区
[0035]
图6：西安市2020土地覆盖模拟图
[0036]
图7：西安市2030土地覆盖预测图
具体实施方式
[0037]
下面结合实例和附图对本发明做进一步详细说明：
[0038]
一种在遥感数据基础上运用统计模型结合城市流数据，通过python实现对市县城镇开发边界模拟和预测技术的方法.
[0039]
案例选取西安市为研究对象，如附图1所示，算法步骤如下：
[0040]
(1)获取30米分辨率地表覆盖数据globeland30中2010和2020年西安市地表覆盖遥感数据，见附图2。
[0041]
(2)选取影响西安市城镇空间扩张的驱动力因素，本实例中驱动力因素包括“到铁路距离”、“到高速路距离”、“到道路距离”、“坡度”、“坡向”和“到政府机构距离”，在arc gis 10.7中处理为30米栅格，见附图3。
[0042]
(3)运用logit模型计算各空间驱动力因素影响权重和用地转换概率，具体计算公式如下：
[0043][0044]
式中p
logit
为土地利用类型变更状态，比如1为非城镇建设用地转换为城镇建设用地，0则表示未变更，x为城镇建设用地空间扩散驱动因素，包括“到铁路距离”、“到高速路距离”、“到道路距离”、“坡度”、“坡向”和“到政府机构距离”，β为空间因素影响系数。
[0045]
通过计算得到各类土地利用类型转变概率，结果见表1。
[0046]
表1各类用地利用类型转变概率统计
[0047]
(4)获取西安手机信令数据，计算各个区县间人流强度，见附图4。
[0048]
(5)运用重力模型对城市流进行计算，具体公式如下：
[0049][0050][0051]
式中i
i
和i
j
分别表示城镇i和城镇j的流强度，d
ij
为城镇i,j的距离，α为参数(此实例中取值1),i
flow(x,y)
为城镇i中元胞xy的流强度，d
i(xy)
为元胞xy到城镇i中心距离，d
j(xy)
为元胞xy到城镇j中心距离。
[0052]
通过计算得到人流强度对各个栅格的影响作用，结果见表2：
[0053] 平均值标准差最小值最大值流强度0.00000006070.00000054500.00000000000.0000178000
[0054]
表2人流强度对各个栅格的影响作用结果
[0055]
(6)选取西安城镇扩张限制因素，本实例中选取“西安城镇建设适宜性评价”、“农业生产适宜性评价”、“生态保护重要性等级”、“河流水系”和“国家、省级文物保护单位范围因素”构建西安市限建区，见附图5。
[0056] 耕地林地草地湿地水体人造地表平均值0.1309180.0043180.004021//0.08707标准差0.0128930.0248880.006261//0.03738
[0057]
(7)根据城镇间流强度，构建城镇等级体系，形成不同等级城镇建设用地分配指标体系，具体公式如下：
[0058]
p
threshold
＝p*k
i
[0059]
式中p
threshold
为“流空间”视角下城镇i的城镇建设用地转换概率，取值为0
‑
1，p为研究范围内基准城镇建设用地转换概率，k
i
为城镇体系中城镇i的权重，在本实例中，西安市不同区县用地分配权重见表3：
[0060]
表3不同区县权重表
[0061][0062]
(8)制定约束性元胞自动机模型转换规则，模拟和预测未来城镇开发边界，具体公式如下：
[0063][0064][0065]
元胞自动机约束规则如下：
[0066]
if any cell is not already developed to urban
[0067]
&if
[0068]
then cell i is developed with some other probability
[0069]
为在t时刻，栅格ij周围n*n的邻域对其影响作用，为栅格ij在t时刻是城镇建设用地的转换概率，n为每次迭代中允许用地类型发生转换的栅格数量，p
logit
和i
flow(x,y)
均采用标准化处理。在本实例中，选取n为3，用地类型转换栅格数量表见表4，最终生成西安市2020土地覆盖模拟图(见附图6)，并根据2010
‑
2020年间西安市土地覆盖类型变化情况，同比预测2030年西安市各类土地覆盖类型变化情况，得到2030土地覆盖预测图(见附图7)。
[0070]
表4用地类型转换栅格数量表
[0071][0072]