一种城市三维空间扩张模拟方法及装置

一种城市三维空间扩张模拟方法及装置
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.所述发展概率预测模块，用于将所述土地利用数据以及所述驱动力因子数据输入至已训练的发展概率预测模型中，以使所述发展概率预测模型预测出各土地利用类型在所
述模拟区域的每一像元中的发展概率；
32.所述像元迭代计算模块，用于对所述模拟区域内的所有像元进行迭代计算，预测每一所述像元经过模拟后的土地利用类型并在模拟后的土地类型为城市用地时，将所述空间因素数据输入到已训练的建筑物高度预测模型中，以使所述建筑物高度预测模型对像元上的建筑物的高度进行预测；其中，预测每一所述像元经过模拟后的土地利用类型具体包括：对所述模拟区域进行邻域划分，根据每一像元所在的邻域内各土地利用类型的占比、各土地利用类型在每一像元上的发展概率、各土地利用类型间的转换矩阵以及各土地类型的惯性系数，计算各土地利用类型在每一像元上的总分布概率，继而根据各土地利用类型在每一像元上的总分布概率确定每一像元模拟后的土地利用类型；
33.所述模拟结果生成模块，用于在所述模拟区域内的所有像元迭代完毕后，获得所述模拟区域的城市三维空间扩张模拟结果。
34.通过实施本发明实施例具有如下有益效果：
35.本发明实施例提供了一种城市三维空间扩张模拟方法及装置，所述方法基于已训练的发展概率预测模型对模拟区域内各土地利用类型在所有像元上的发展概率进行预测，然后根据各土地利用类型在每一像元上的发展概率、各土地利用类型间的转换矩阵以及各土地类型的惯性系数，计算各土地利用类型在每一像元上的总分布概率，根据各土地利用类型在每一像元上的总分布概率确定每一像元模拟后的土地利用类型，继而在像元经过模拟后的土地类型为城市用地时，通过已训练的建筑物高度模型对像元上建筑物的建筑物高度进行预测。与现有技术相比，在水平方向对每个像元的土地利用类型进行模拟预测之后，还对归属于城市用地的像元上的建筑物的高度进行了预测，实现了对城市扩张在水平方向和竖直方向上的模拟预测，使得对城市扩张的模拟更加的准确。
附图说明
36.图1是本发明一实施例提供的一种城市三维空间扩张模拟方法的流程示意图。
37.图2是本发明一实施例提供的一种城市三维空间扩张模拟装置的结构示意图。
具体实施方式
38.下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本邻域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
39.如图1所示，本发明一实施例提供了一种城市三维空间扩张模拟方法，至少包括如下步骤：
40.步骤s101:获取模拟区域的土地利用数据、模拟区域中每一像元所对应影响土地利用变化的驱动力因子数据以及模拟区域中每一像元所对应影响建筑物高度变化的空间因素数据。
41.步骤s102：将所述土地利用数据以及所述驱动力因子数据输入至已训练的发展概率预测模型中，以使所述发展概率预测模型预测出各土地利用类型在所述模拟区域的每一像元中的发展概率。
42.步骤s103：对所述模拟区域内的所有像元进行迭代计算，预测每一所述像元经过模拟后的土地利用类型并在模拟后的土地类型为城市用地时，将所述空间因素数据输入到已训练的建筑物高度预测模型中，以使所述建筑物高度预测模型对像元上的建筑物的高度进行预测；其中，预测每一所述像元经过模拟后的土地利用类型具体包括：对所述模拟区域进行邻域划分，根据每一像元所在的邻域内各土地利用类型的占比、各土地利用类型在每一像元上的发展概率、各土地利用类型间的转换矩阵以及各土地类型的惯性系数，计算各类土地利用类型在每一像元上的总分布概率，继而根据各类土地利用类型在每一像元上的总分布概率确定每一像元模拟后的土地利用类型。
43.步骤s104：在所述模拟区域内的所有像元迭代完毕后，获得所述模拟区域的城市三维空间扩张模拟结果。
44.对于步骤s101:通过对模拟区域的高分遥感影像进行预处理解译，获得模拟区域的土地利用数据；可以理解的土地利用数据包括土地利用类型以及各土地利用类型的面积，在本发明中土地利用类型，非建设用地、公共管理服务用地、商业用地、居住用地以及工业用地，其中，公共管理服务用地、商业用地、居住用地以及工业用地均归为本发明所定义的城市用地：当然具体城市用地的分类可以根据实际情况进行确定。
45.对于每一像元所对应的影响土地利用变化的驱动力因子数据，优选的在本发明中选取如下数据作为上述驱动力因子数据：模拟区域的地形数据(dem、坡度)、像元到城市中心的距离、像元到区县中心的距离、像元到铁路的距离、像元到高速公路的距离、像元到国道的距离、像元到省道的距离、像元到火车站的距离、像元到地铁站的距离、像元到海洋的距离、像元到湖泊的距离、像元到河流的距离、模拟区域的城市路网密度、模拟区域的人口密度、模拟区域的pm2.5和模拟区域的gdp。当然具体的驱动力因子数据，可以根据实际情况进行选定。
46.对于每一像元所对应的影响建筑物高度变化的空间因素数据，优选的在本发明中选取如下数据作为上述空间因素数据：模拟区域的地形数据(dem、坡度)、像元到公园绿地的距离、像元至水体的距离、模拟区域的房价、模拟区域的娱乐设施分布密度、模拟区域的超市分布密度、模拟区域的餐饮分布密度、模拟区域的工厂分布密度、模拟区域的商场分布密度、模拟区域的人口密度、模拟区域的医院分布密度、、模拟区域的公交站点分布密度、模拟区域的夜间灯光强度、像元到城市中心的距离、像元到区县中心的距离以及像元到城市道路的距离。当然具体的空间因素数据，可以根据实际情况进行选定。
47.对于步骤s102，首先对上述发展概率预测模型进行说明，在一个优选的实施例中，所述发展概率预测模型的构建方法，具体包括：获取所述模拟区域的历史土地利用数据以及所述驱动力因子数据；对所述历史土地利用数据以及所述驱动因子进行采样，获得采样数据；将所述采样数据输入到神经网络模型中进行训练，获得所述发展概率预测模型。优选的，通过均匀采样或分层采样对所述历史土地利用数据以及所述驱动因子进行采样，获得所述采样数据。优选的在将所述采样数据输入到神经网络模型中进行训练之前，还包括：对所述采样数据进行归一化处理。
48.为更好的说明发展概率预测模型的具体构建过程，以下列举实际示例进行进一步阐释：
49.以广东省深圳市所在的区域作为模拟区域，通过landsat影像解译生成，深圳市
2009年和2014年土地利用数据，并划分为非建设用地、公共管理服务用地、商业用地、居住用地和工业用地；其中，公共管理服务用地、商业用地、居住用地和工业用地均归为城市用地，2009年土地利用数据作为本发明上述历史土地利用数据，2014年土地利用数据作为验证数据。通过在地图软件上进行爬虫获取深圳市2009年和2014年的建筑物高度数据，将深圳市2009年的建筑物高度数据作为本发明的历史建筑物高度数据，将2014年的建筑物高度数据作为验证数据(这两个数据的具体用处，将在后文进行阐释)。所有的区域数据文件都要被统一成同尺度的栅格影像以方便进行像元尺度上的同质分析。紧接着选取上文所述的驱动力因子数据作为对深圳市进行城市扩张模拟时的驱动力因子数据；
50.在数据获取完毕后，在深圳市2009年土地利用数据和驱动力因子数据上进行随机点采样。根据深圳市2009年土地利用数据的特点选择采样方式，方法提供两种采样方式：1、均匀采样策略；2、分层采样策略。均匀采样策略使得各类土地利用类型的采样点一致，适用于不同土地利用类型面积差异较大的区域；分层采样法使得各类土地利用类型的采样点占总采样点的比例与各类初始土地利用类型与区域面积的比例一致，计算量较小，适合较大尺度区域采样。
51.示意性的，在本发明实施例中使用分层采样策略，使得各类土地利用类型的采样点与其土地利用类型面积所占比例一致，用采样点提取深圳市2009年土地利用数据和上述驱动力因子数据，从数据集中抽取10％作为采样点，通过采样点抽取驱动力因子以及其采样点对应的土地利用类型，生成上述采样数据。采样数据公式表示为：
52.x(l)＝[lucc(l)，x1(l)，x2(l)，x3(l)，...，x
w
(l)，...，x
m
(l)]
t
；
[0053]
其中lucc(l)表示第l个采样点的土地利用类型；x
w
(l)表示第l个采样点上的第w个驱动力因子，t为转置。
[0054]
紧接着对采样数据进行归一化处理，统一各个变量的量纲。归一化公式表示为：
[0055][0056]
其中max
w
和min
w
分别是第w个驱动力因子的最大和最小值，m为驱动因子的个数。
[0057]
采样数据归一化完毕之后，采用神经网络模型对归一后的采样数据进行训练，使得模型能达到更快的训练速度和训练精度。神经网络模型的公式如下：
[0058][0059]
其中η(n)是第n次迭代的学习率，e(n)和e(n
‑
1)是相邻两次迭代的神经网络输出的均方根误差，a，b，c是常数，取值范围分别为(1，2)、(0，1)、[1，1.1]。
[0060]
在本发明中，上述神经网络的输入层接收到驱动力因子数据后，将它们输出到隐藏层。本发明实施例中设隐藏层数为h＝12，隐藏层第v个神经元所收到的信号公式为：
[0061][0062]
其中net
v
(l)为隐藏层第v个神经元所收到的信号；ω
(w，v)
为输入层和隐藏层之间一一对应的参数，也就是两个层级间的权重值；隐藏层net
v
(l)对接收到的信号进行处理，然后将处理结果sig mod
v
输出至输出层。处理数据的函数通常为sig mod函数：
[0063][0064]
隐藏层输出函数与输出层输出函数间也有权值一一对应，输出层数数目u与本发明模拟的土地利用类别相等，土地利用类型共有5类即u＝5，假设在第q个输出层中，那么输出层的数据处理函数为：
[0065][0066]
net
q
(l)为输出层第q个神经元所收到的信号；θ
(v，q)
是隐藏层和输出层之间的权重参数，输出层接收到信号之后，会对这些信号产生响应，该响应值通过下式生成分布概率p(l，q)：
[0067][0068]
在本发明中p(l，q)为采样点训练神经网络的输出。网络训练完毕后，将所有土地利用变化驱动因子图层叠加并全图输入训练好的神经网络，将p(l，q)中的l转化成为像元的坐标(i，j)，得到如下公式：
[0069][0070]
其中p(i，j，q)为神经网络输出的第(i，j)个像元上第q种用地类型的发展概率。
[0071]
神经网络训练完成后，获得本发明上述发展概率预测模型，然后将深圳市当前的土地利用数据(即本发明上述模拟区域的土地利用数据)以及驱动力因子数据输入至训练好的发展概率预测模型中，预测出各土地利用类型在上述模拟区域的每一像元上的发展概率。
[0072]
优选的在模型训练完毕后，还可采用深圳市2014年的土地利用数据对模型进行验证，若验证不通过，则调整网络参数重新进行训练。
[0073]
对于步骤s103；首先对建筑物高度预测模型进行说明，在一个优选的实施例中，所述建筑物高度预测模型的构建方法，具体包括：获取所述模拟区域的历史建筑物高度数据以及所述空间因素数据；对所述历史建筑物高度数据以及所述空间因素数据进行采样，获得第二采样数据；根据所述第二采样数据通过随机森林回归算法进行模型训练，获得所述建筑物高度预测模型。优选的，通过随机采样或分层采样对所述历史建筑物高度数据以及所述空间因素数据进行采样，获得第二采样数据。优选的，所述根据所述第二采样数据通过随机森林回归算法进行模型训练，获得所述建筑物高度预测模型，具体包括：根据所述第二采样数据生成随机向量序列；根据所述第二采样数据中的每一子样本以及所述随机向量序列，分别构建每一子样本的回归树模型；根据每一子样本的回归树模型组合成初始随机森林回归模型，继而对所述初始随机森林回归模型进行训练，获得所述建筑物高度预测模型。
[0074]
依旧以深圳市为模拟区域为例，通过上文所描述的方法获取深圳市2009年的建筑物高度数据作为本发明的历史建筑物高度数据并结合上文所获取的空间因素数据，作为训
练建筑物高度预测模型的基础数据。
[0075]
在深圳市2009年的建筑物高度数据与上述空间因素数据上进行随机点采样，获得建筑物高度的采样数据(即本发明上述第二采样数据)；根据建筑物高度分布数据的特点选择采样方式，方法提供两种采样方式：1、随机采样策略；2、分层采样策略。
[0076]
示意性的，在本发明中使用随机采样策略，用采样点提取深圳市2009年的建筑物高度数据与上述空间因素数据，从数据集中抽取20％作为采样点，采样后的样本公式表示为：
[0077]
s(l)＝[bh(l)，s1(l)，s2(l)，s3(l)，...，s
w
(l)，...，s
m
(l)]
t
；
[0078]
其中bh(l)表示第l个采样点的建筑物高度；s
w
(l)表示第l个采样点上的第w个空间因素，t为转置。
[0079]
紧接着根据样本数据集s(l)生成随机向量序列θ
i
(i＝1，2，...，k)；
[0080]
然后采用bootstrap抽样方法从数据集s(l)中有放回地随机抽取k个子样本集，所有子样本集表示为s
i
(i＝1，2，...，k)；之后，假定子样本集相互独立分布，对每一个子样本集s
i
分别构建回归树模型{h(x，θ
i
)}；
[0081]
最后，由多个回归树组合{h(x，θ1)，h(x，θ2)，...，h(x，θ
k
)，}构成随机森林回归模型，模型的结果是所有决策树，h(x，θ2)回归结果的平均值，公式表示为：
[0082][0083]
其中；h(x)表示随机森林回归模型的预测值，{h(x，θ
i
)}表示第i个回归树模型的回归结果值。
[0084]
在本发明实施例中使用常用的精度评价指标检验随机森林回归对建筑物高度预测的精度，主要的评价指标为：拟合优度(决定系数，r2)、均方根误差(rmse)、平均绝对误差(mae)、解释方差得分(ev)。公式表示为：
[0085][0086][0087][0088][0089]
其中far
t，i
和far
p，i
分别为第i个建筑单元的实际容积率和预测容积率；为
所有建筑单元实际容积率的平均值，等于n为所有建筑单元的数量；var{far
t，i
‑
far
p，i
}和var{far
t，i
}分别为far
t，i
‑
far
p，i
和far
t，i
的方差。r2的取值范围为[0，1]，值越大，回归拟合效果越好。rmse和mae的取值范围都为[0，+∞]，值越小，回归拟合效果更好。ev用于解释因变量的方差变化，值越大说明回归拟合效果越好。
[0090]
在本发明实施例中，r2为0.8628，rmse为2.0180，mae为1.3707，ev为0.8632。结果表明随机森林回归的拟合效果较好。
[0091]
在高度预测模型构建完毕后，开始对模拟区域内所有像元进行迭代计算，预测每一像元经过模拟后的土地利用类型；在一个优选的实施例中，所述根据每一像元所在的邻域内各土地利用类型的占比、各土地利用类型在每一像元中的发展概率、各土地利用类型间的转换矩阵以及各土地类型的惯性系数，计算各类土地利用类型在每一像元上的总分布概率，具体包括：将所述每一像元所在的邻域内各土地利用类型的占比、各土地利用类型在每一像元中的发展概率、各土地利用类型间的转换矩阵以及各土地类型的惯性系数相乘，获得各类土地利用类型在每一像元上的总分布概率。优选的，所述根据各类土地利用类型在每一像元上的总分布概率，确定每一像元模拟后的土地利用类型，具体包括：根据各类土地利用类型在每一像元上的总分布概率构成轮盘，通过轮盘赌选择算法确定每一像元经过模拟后的土地利用类型。
[0092]
具体的，在迭代前提前预设好邻域大小、各土地利用类型间的转换矩阵以及各土地类型的惯性系数。
[0093]
示意性的在本发明实施例中的深圳市影像像元为30m，因而邻域范围选择7
×
7摩尔邻域。本发明实施例模拟的时长为：5年。迭代次数设为100次。
[0094]
迭代模拟区域当前土地利用数据的所有像元，计算每一像元所在邻域内包含的各土地利用类型在邻域内所占的比例(即上述每一像元所在的邻域内各土地利用类型的占比)，并与步骤s102中所得出的各土地利用类型在每一像元上的发展概率、各土地利用类型间的转换矩阵和各土地类型的惯性系数共同组成该像元上各类土地利用类型的总分布概率。
[0095]
本发明实施例中，通过邻域大小构建一个扫描窗口，然后统计扫描窗口内的各类像元的数量，进而来衡量每一像元所在的邻域内各土地利用类型的占比，具体公式如下：
[0096][0097]
其中为t第次迭代时，在第i行、第j列的像元为中心的n*n邻域内第q类土地利用类型的占比，con为条件函数，s
k
表示邻域内的当前被扫描的像元，该像元属于第k种土地利用类型；con(s
k
＝q)表示检测邻域内的当前土地利用像元类型是否为第q类；ω
k
为不同土地利用类型的预设的权重；则构成轮盘赌的各类土地利用类型在每一像元上的总分布概率的公式为：
[0098][0099]
其中表示第t第次迭代时，第q类土地类型在第(i，j)个像元上的总分布概率，
p(i，j，q)为由发展概率预测模型所输出的第q类土地类型在第(i，j)个像元中的发展概率，c(i’，j’)表示转换限制矩阵，i’，j’表示转换限制矩阵的坐标，表示在第t次迭代时，第q类土地类型的惯性系数。
[0100]
然后，该像元上的各类土地利用类型的总分布概率构成轮盘，通过轮盘赌选择，使区域内各种土地利用类型在像元上竞争，竞争获胜的土地利用类型占据该像元，作为该像元经过模拟后的土地利用类型；若像元上竞争获胜的土地利用类型为城市用地类型，使用上述的建筑物高度预测模型对该像元上的建筑物高度进行预测；否则，不需要进行预测。
[0101]
在本发明实施例中，建筑物高度预测的公式可以表示为：
[0102][0103]
其中bh
i，j
表示当第(i，j)个像元经过模拟后的土地利用类型为城市用地时，第(i，j)个像元上建筑物的预测高度，rfregression()为上述建筑物高度预测模型，为第(i，j)个像元上第w个空间因素。
[0104]
对于步骤s104，在迭代完影像中的全部有效像元，后获得城市三维空间扩张模拟结果并输出。
[0105]
通过本发明上述各实施例所述的城市三维空间扩张模拟方法，能够对模拟区域的水平方向和竖直方向的城市扩张进行模拟，提高了模拟结果的准确性。
[0106]
在上述方法项实施例的基础上，本发明对应提供了装置项实施例；
[0107]
如图2所示，本发明一实施例提供了一种城市三维空间扩张模拟装置，包括数据获取模块，发展概率预测模块、像元迭代计算模块以及模拟结果生成模块；
[0108]
所述数据获取模块，用于获取模拟区域的土地利用数据、模拟区域中每一像元所对应影响土地利用变化的驱动力因子数据以及模拟区域中每一像元所对应影响建筑物高度变化的空间因素数据；
[0109]
所述发展概率预测模块，用于将所述土地利用数据以及所述驱动力因子数据输入至已训练的发展概率预测模型中，以使所述发展概率预测模型预测出各土地利用类型在所述模拟区域的每一像元中的发展概率；
[0110]
所述像元迭代计算模块，用于对所述模拟区域内的所有像元进行迭代计算，预测每一所述像元经过模拟后的土地利用类型并在模拟后的土地类型为城市用地时，将所述空间因素数据输入到已训练的建筑物高度预测模型中，以使所述建筑物高度预测模型对像元上的建筑物的高度进行预测；其中，预测每一所述像元经过模拟后的土地利用类型具体包括：对所述模拟区域进行邻域划分，根据每一像元所在的邻域内各土地利用类型的占比、各土地利用类型在每一像元上的发展概率、各土地利用类型间的转换矩阵以及各土地类型的惯性系数，计算各土地利用类型在每一像元上的总分布概率，继而根据各土地利用类型在每一像元上的总分布概率确定每一像元模拟后的土地利用类型；
[0111]
所述模拟结果生成模块，用于在所述模拟区域内的所有像元迭代完毕后，获得所述模拟区域的城市三维空间扩张模拟结果。
[0112]
需说明的是，本发明上述装置项实施例是与本发明方法项实施例对应的，其能够实现本发明上述任意一项所述的城市三维空间扩张模拟方法。另，以上所描述的装置实施例仅仅是示意性的，其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开
的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。另外，本发明提供的装置实施例附图中，模块之间的连接关系表示它们之间具有通信连接，具体可以实现为一条或多条通信总线或信号线。本邻域普通技术人员在不付出创造性劳动的情况下，即可以理解并实施。
[0113]
以上所述是本发明的优选实施方式，应当指出，对于本技术邻域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也视为本发明的保护范围。