本发明专利技术属于地理信息系统的虚拟地理环境及仿真研究领域，涉及到城市空间发展的仿真及预测技术。本发明专利技术利用支持向量机技术设计了一种城市元胞用地性质变化决策函数，在此基础上提出了一种基于支持向量机的地理元胞模型并给出了基于该模型的城市空间发展模拟及预测的方法。该方法能够克服现有地理元胞模型的线性性缺陷，具有较快的仿真速度和较高的仿真预测精度。利用支持向量机技术建立城市元胞模型，用于解决地理元胞状态预测这一类典型非线性问题，有效的解决了传统的地理元胞模型线性性导致仿真预测结果误差过大问题，同现有的方法相比，对样本数量要求较少、计算效率更高、精度也更高。
【技术实现步骤摘要】
本专利技术属于地理信息系统的虚拟地理环境及仿真研究领域，涉及到城市空间发展的仿真及预测技术。
技术介绍
随着全球经济的持续发展，城市化和城市的可持续发展已成为全球城市发展的热点话题。城市空间作为城市各项功能的载体，其形态的合理性将对城市居民的生存质量产生决定性的影响，并进而影响到整个人类的生存和发展。在全球范围内的城市快速扩张过程中，城市空间形态急剧变化导致了诸如土地资源利用不当、耕地流失、生态环境破坏、污染严重、城市功能布局失衡等许多尖锐的现实问题的出现。开展对城市空间发展模型及模拟预测研究，有利于确保城市化与社会经济发展的协调和同步，优化土地利用，降低城市化进程中的盲目性，有效保证城市空间和社会经济的长期可持续发展。然而，由于城市空间变化具有显著的动态特征，不同阶段有不同特点。与此同时，城市空间形态演变又是一个典型的动态空间复杂系统，具有开放、自组织、动态、非平衡等特征，会受到政治、经济、文化、社会、自然、政策等多种因素的影响，其演变过程表现出高度的非线性性、动态性和复杂性。因此，现有的宏观、静态和确定的城市模型难以满足现代城市空间发展的要求。上世纪80至90年代，元胞自动机(Cellular Automata，CA)技术被引入到城市发展研究当中，形成了一个新的技术分支：地理元胞信息系统。CA以一种动态的视角来研究城市发展，并充分利用了计算机的并行计算能力来加速建模和仿真预测过程，有效的推动了城市空间研究的进程。但传统的CA模型转换规则具有典型的线性特征，在描述城市发展这一类复杂非线性系统时存在着较大的局限性，其仿真预测的结果存在着较大的误差。针对地理元胞模型的线性性的缺陷，逻辑回归、神经网络、决策树和粗糙集等新方法被引入到地理元胞模型转换规则的确定中，在一定程度上推动了地理元胞模型的发展。然而，逻辑回归模型要求转换规则的特征变量之间线性无关，这限制了它的应用范围；神经网络可以很好的解决非线性分类问题，但其自身属于黑盒方法，内部机理难以解释，也难以反映模型表述对象的内部机制；决策树和粗糙集等其余方法则对样本数量有较高的要求，在样本数量不足的情况下，其模型精度较低。支持向量机(Support Vector Machines，SVM)是一种建立在统计学习理论、结构风险最小化及VC维(Vapnik-Chervonenkis Dimension)理论基础上的机器学习方法。SVM具有理论基础坚实、数学模型简单明确的特点，在模式识别、函数估计、时序预测和回归分析等领域得到了深入的发展，其成功应用方向包括手写识别、人像识别、基因分类和二值分类等。SVM实际上是建立在统计学习理论基础上的一种数据挖掘及分类技术，其机理在于在给定的分类条件下寻找一个最优分类超平面，在确保分类精度的前提下，使得超平面两侧的空白区域最大化。地理元胞自动机是元胞自动机在城市空间发展的成功应用，城市空间发展可以在微观上看作是城市地理元胞状态的变化，地理元胞的状态可以简化为建设用地和非建设用地，因此地理元胞的状态预测可以看做一个典型的非线性分类问题。支持向量机在地理元胞自动机转换规则制定中的应用研究及基于支持向量机的元胞自动机模型的设计实现对于克服现有技术线性性缺陷及提高模型预测精度具有重要的理论和应用价值。
技术实现思路
本专利技术的目的在于克服传统地理元胞模型在具有典型非线性特征的城市地理现象领域应用存在的不足，利用支持向量机技术设计了一种城市元胞用地性质变化决策函数，在此基础上提出了一种基于支持向量机的地理元胞模型并给出了基于该模型的城市空间发展模拟及预测的方法。该方法能够克服现有地理元胞模型的线性性缺陷，具有较快的仿真速度和较高的仿真预测精度。本专利技术的技术方案是：S1、建立基于支持向量机的城市元胞自动机模型(Support Vector Machines based Geographical Cellular Automata，SVM-GCA)，具体为：S11、定义城市元胞状态：城市元胞状态为城市用地性质描述即城市元胞模型的因变量y，所述城市元胞状态包括两种状态：城市用地和非城市用地，分别用1和-1进行编码，1表示该元胞为城市用地，-1表示该元胞为非城市用地；S12、确定城市元胞邻域关系：采用Moore型邻域，元胞的上、下、左、右、左上、左下、右上、右下相邻的8个元胞称为该元胞的邻域，邻域半径为1；S13、确定城市元胞用地转换因子：城市元胞用地性质的转化取决于每个元胞当前的地理环境因子x1、经济增长因子x2、人口密度因子x3、土地价格因子x4、交通状况因子x5以及政策因子x6，其中，所述经济增长因子x2、人口密度因子x3、土地价格因子x4、交通状况因子x5和政策因子x6从城市统计数据中抽取；S14、设计城市元胞状态转换规则：城市元胞在t+1时刻转换为建设用地的概率为其中，F为决策函数，R为约束因子，Nt为t时刻元胞邻居影响因子；S2、进行城市遥感图像分类处理和用地转换因子确定及优化，具体为：利用Envi 4.0软件，通过非监督分类法按照四种用地类别对遥感图像进行分类处理，其中，所述四种用地类别为：水域及山体(包括河流湖泊等，编号为1)、农林用地(包括园林和植被、耕地和空地等，编号为2)、城乡建设用地(包括居住用地、商业、工业、公共设施等，编号为3)和交通用地(包括公路、铁路、机场等，编号为4)，并由此构成S13所述地理环境因子x1；S3、基于SVM-GCA的城市发展仿真预测，具体为：S31、基于GIS平台对城市遥感图像进行分类处理，获得城市用地基础地理信息，根据城市统计数据构造用地转换因子，在GIS平台中用6个图层分别描述6个影响因子；S32、构建城市元胞模型：元胞形状选择为正方形，利用GIS系统中栅格图像代表元胞空间，每个栅格对应一个元胞，元胞大小为L×L，元胞邻域类型确定为Moore型，即每个元胞有8个邻居元胞，其中，L为不为零的自然数；S33、根据S14所述转换规则，计算城市元胞在t+1时刻转换为建设用地的概率S34、得出城市元胞在t+1时刻用地性质，所述用地性质按如下方法确定：设城市元胞用地性质转变概率阈值为Pthreshold，当城市元胞在t+1时刻转换为建设用地的概率Pt+1≥Pthreshold时，则判定t+1时刻该元胞将变为城市用地，当Pt+1＜Pthreshold时，则判定t+1时刻该元胞用地性质不变，其中，所述Pthreshold通过城市发展历史数据拟合得到，为经验值；S4、根据S1-S3对城市元胞用地性质变化进行K次迭代，最终获得未来城市空间发展模拟结果，其中，所述K为不为零的自然数，是经验值。进一步地，S14所述t时刻邻居影响因子Nt的确定方法具体如下：t时刻的邻居影响因子Nt用元胞的邻居中城市用地元胞所占比例来表示，若第j个邻居元胞Cj为城市用地元胞，则S(Cj)＝1，否则，S(Cj)＝0，其中，S(Cj)为元胞状态函数。进一步地，S14所述约束因子R的确定方法如下：如果地理元胞用地性质无法转换为建设用地，则R＝0，否则R＝1。进一步地，S14所述决策函数F的确定方法如下：决策函数设训练样本集为{(xi，yi)，i＝1,2，…，l本文档来自技高网...
【技术保护点】
一种基于SVM‑GCA的城市空间发展模拟及预测方法，其特征在于，具体步骤如下：S1、建立基于支持向量机的城市元胞自动机模型(Support Vector Machines based Geographical Cellular Automata，SVM‑GCA)，具体为：S11、定义城市元胞状态：城市元胞状态为城市用地性质描述即城市元胞模型的因变量y，所述城市元胞状态包括两种状态：城市用地和非城市用地，分别用1和‑1进行编码，1表示该元胞为城市用地，‑1表示该元胞为非城市用地；S12、确定城市元胞邻域关系：采用Moore型邻域，元胞的上、下、左、右、左上、左下、右上、右下相邻的8个元胞称为该元胞的邻域，邻域半径为1；S13、确定城市元胞用地转换因子：城市元胞用地性质的转化取决于每个元胞当前的地理环境因子x1、经济增长因子x2、人口密度因子x3、土地价格因子x4、交通状况因子x5以及政策因子x6，其中，所述经济增长因子x2、人口密度因子x3、土地价格因子x4、交通状况因子x5和政策因子x6从城市统计数据中抽取；S14、设计城市元胞状态转换规则：城市元胞在t+1时刻转换为建设用地的概率为其中，F为决策函数，R为约束因子，Nt为t时刻元胞邻居影响因子；S2、进行城市遥感图像分类处理和用地转换因子确定及优化，具体为：利用Envi 4.0软件，通过非监督分类法按照四种用地类别对遥感图像进行分类处理，其中，所述四种用地类别为：水域及山体(包括河流湖泊等，编号为1)、农林用地(包括园林和植被、耕地和空地等，编号为2)、城乡建设用地(包括居住用地、商业、工业、公共设施等，编号为3)和交通用地(包括公路、铁路、机场等，编号为4)，并由此构成S13所述地理环境因子x1；S3、基于SVM‑GCA的城市发展仿真预测，具体为：S31、基于GIS平台对城市遥感图像进行分类处理，获得城市用地基础地理信息，根据城市统计数据构造用地转换因子，在GIS平台中用6个图层分别描述6个影响因子；S32、构建城市元胞模型：元胞形状选择为正方形，利用GIS系统中栅格图像代表元胞空间，每个栅格对应一个元胞，元胞大小为L×L，元胞邻域类型确定为Moore型，即每个元胞有8个邻居元胞，其中，L为不为零的自然数；S33、根据S14所述转换规则，计算城市元胞在t+1时刻转换为建设用地的概率S34、得出城市元胞在t+1时刻用地性质，所述用地性质按如下方法确定：设城市元胞用地性质转变概率阈值为Pthreshold，当城市元胞在t+1时刻转换为建设用地的概率Pt+1≥Pthreshold时，则判定t+1时刻该元胞将变为城市用地，当Pt+1＜Pthreshold时，则判定t+1时刻该元胞用地性质不变，其中，所述Pthreshold通过城市发展历史数据拟合得到，为经验值；S4、根据S1‑S3对城市元胞用地性质变化进行K次迭代，最终获得未来城市空间发展模拟结果，其中，所述K为不为零的自然数，是经验值。...
【技术特征摘要】
1.一种基于SVM-GCA的城市空间发展模拟及预测方法，其特征在于，具体步骤如下：S1、建立基于支持向量机的城市元胞自动机模型(Support Vector Machines based Geographical Cellular Automata，SVM-GCA)，具体为：S11、定义城市元胞状态：城市元胞状态为城市用地性质描述即城市元胞模型的因变量y，所述城市元胞状态包括两种状态：城市用地和非城市用地，分别用1和-1进行编码，1表示该元胞为城市用地，-1表示该元胞为非城市用地；S12、确定城市元胞邻域关系：采用Moore型邻域，元胞的上、下、左、右、左上、左下、右上、右下相邻的8个元胞称为该元胞的邻域，邻域半径为1；S13、确定城市元胞用地转换因子：城市元胞用地性质的转化取决于每个元胞当前的地理环境因子x1、经济增长因子x2、人口密度因子x3、土地价格因子x4、交通状况因子x5以及政策因子x6，其中，所述经济增长因子x2、人口密度因子x3、土地价格因子x4、交通状况因子x5和政策因子x6从城市统计数据中抽取；S14、设计城市元胞状态转换规则：城市元胞在t+1时刻转换为建设用地的概率为其中，F为决策函数，R为约束因子，Nt为t时刻元胞邻居影响因子；S2、进行城市遥感图像分类处理和用地转换因子确定及优化，具体为：利用Envi 4.0软件，通过非监督分类法按照四种用地类别对遥感图像进行分类处理，其中，所述四种用地类别为：水域及山体(包括河流湖泊等，编号为1)、农林用地(包括园林和植被、耕地和空地等，编号为2)、城乡建设用地(包括居住用地、商业、工业、公共设施等，编号为3)和交通用地(包括公路、铁路、机场等，编号为4)，并由此构成S13所述地理环境因子x1；S3、基于SVM-GCA的城市发展仿真预测，具体为：S31、基于GIS平台对城市遥感图像进行分类处理，获...
【专利技术属性】
技术研发人员：韩效，刘民岷，
申请(专利权)人：西南交通大学，
类型：发明
国别省市：四川;51
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。