CN113177354A - 基于改进ann耦合异步元胞自动机的城市动态扩张模拟方法 - Google Patents
基于改进ann耦合异步元胞自动机的城市动态扩张模拟方法 Download PDFInfo
- Publication number
- CN113177354A CN113177354A CN202110433496.1A CN202110433496A CN113177354A CN 113177354 A CN113177354 A CN 113177354A CN 202110433496 A CN202110433496 A CN 202110433496A CN 113177354 A CN113177354 A CN 113177354A
- Authority
- CN
- China
- Prior art keywords
- ann
- particles
- cell
- hpso
- particle
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 230000001413 cellular effect Effects 0.000 title claims abstract description 40
- 238000000034 method Methods 0.000 title claims abstract description 24
- 238000004088 simulation Methods 0.000 title claims abstract description 24
- 230000008878 coupling Effects 0.000 title claims abstract description 5
- 238000010168 coupling process Methods 0.000 title claims abstract description 5
- 238000005859 coupling reaction Methods 0.000 title claims abstract description 5
- 239000002245 particle Substances 0.000 claims abstract description 87
- 238000006243 chemical reaction Methods 0.000 claims abstract description 22
- 238000005457 optimization Methods 0.000 claims abstract description 14
- 238000005192 partition Methods 0.000 claims abstract description 11
- 238000011160 research Methods 0.000 claims abstract description 10
- 238000000638 solvent extraction Methods 0.000 claims abstract description 4
- 238000013528 artificial neural network Methods 0.000 claims description 24
- 230000008859 change Effects 0.000 claims description 18
- 230000006399 behavior Effects 0.000 claims description 13
- 230000000694 effects Effects 0.000 claims description 13
- 238000011161 development Methods 0.000 claims description 9
- 230000006870 function Effects 0.000 claims description 9
- 230000007704 transition Effects 0.000 claims description 6
- 230000009286 beneficial effect Effects 0.000 claims description 4
- 238000004422 calculation algorithm Methods 0.000 claims description 4
- 230000003993 interaction Effects 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 claims description 3
- 230000002093 peripheral effect Effects 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 238000005065 mining Methods 0.000 claims description 2
- 230000001537 neural effect Effects 0.000 claims description 2
- 239000000284 extract Substances 0.000 abstract description 2
- 230000008569 process Effects 0.000 description 3
- 238000007477 logistic regression Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000001364 causal effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 238000007637 random forest analysis Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
- G06F18/232—Non-hierarchical techniques
- G06F18/2321—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
- G06F18/23213—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions with fixed number of clusters, e.g. K-means clustering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Artificial Intelligence (AREA)
- General Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Business, Economics & Management (AREA)
- Software Systems (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Mathematical Physics (AREA)
- Computational Linguistics (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- Tourism & Hospitality (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Medical Informatics (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Probability & Statistics with Applications (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- General Business, Economics & Management (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种基于改进ANN耦合异步元胞自动机的城市动态扩张模拟方法,属于地理信息技术领域,S1:通过K‑means方法对研究区域进行分区;S2:在每个分区内构建ANN‑HPSO‑CA模型得到每个分区的转换规则;S3:将ANN‑HPSO‑CA模型与异步元胞自动机模型耦合。本发明将一种与改进ANN与异步元胞自动机CA模型相结合,通过使用粒子群结合CA优化之后的ANN提取转换规则,并且考虑了因地理位置和经济差异所导致的空间异质性与时间异质性,更加符合真实的城市扩张模拟,能够提升城市扩张模拟精度高,为相关领域决策者提供了数据支撑和参考。
Description
技术领域
本发明属于地理信息技术领域,涉及一种基于改进ANN耦合异步元胞自动机的城市动态扩张模拟方法。
背景技术
元胞自动机(Cellular Automata,CA)是一种时间、空间、状态都离散,空间相互作用和时间因果关系都为局部的网格动力学模型,具有模拟复杂***时空演化过程的能力。CA的核心问题是获取有效的转换规则。如今有许多城市扩张模拟方法和模型,例如逻辑回归、人工神经网络和随机森林算法等与CA模型耦合。逻辑回归能解释驱动因子对城市扩张模拟的重要性,但处理非线性、复杂的扩张模拟变化时模拟精度较低。所以采用使用粒子群结合CA优化后的人工神经网络ANN提取转换规则来进行城市扩张模拟,另外由于每一个城市栅格由于经济,地理位置的差异,存在空间异质性的情况,并且由于一些元胞在演化的过程中由于演化速率的不同,也同时存在时间差异性,传统的元胞自动机与现实模拟有较大的的差异,导致模拟精度较低。
发明内容
有鉴于此,本发明的目的在于提供一种基于改进ANN与分区异步元胞自动机模型的城市动态扩张模拟方法。
为达到上述目的,本发明提供如下技术方案:
一种基于改进ANN耦合异步元胞自动机的城市动态扩张模拟方法,包括以下步骤:
S1:通过K-means方法对研究区域进行分区;
S2:通过基于ANN-HPSO-CA模型得到每个分区的转换规则;
S3:将ANN-HPSO-CA模型与异步元胞自动机模型耦合。
进一步,在所述步骤S1中,将基期年与末期年两年的研究区域的数据进行叠加,得到土地利用变化数据,对所述土地利用变化数据与其他影响土地利用变化的驱动因子数据进行K-means空间聚类分析,聚类结果中的每一个类(簇)对应元胞空间的一个分区。
进一步,所述步骤S2具体包括以下步骤:
S21:基于人行为粒子群优化和元胞自动机的改进人工神经网络ANN-HPSO-CA,在原粒子群优化算法PSO中引入全局最坏粒子,即在每次迭代中,整个种群中适应性最差的粒子,表示为Gworst:
Gworst(t)=argmax{f(Pbest1),f(Pbest2),…,f(PbestN)}
f(·)代表对应粒子的适应度值;
引入服从标准正态分布的学习系数r3,即r3∈(0,1),若r3>0,则认为是一个有推力的学习系数,有助于提高粒子的飞行速度;若r3<0,则认为是惩罚的学习系数,降低粒子的飞行速度;若r3=0,表示这些坏习惯或行为对粒子没有影响;将两个加速度系数r1和r2分别替换为随机学习系数r1和r2:
r1和r2是两个值在[0,1]之间的随机数,r1+r2=1;对于一个粒子,在每一个d=1,2,…,D中,r1、r2、r3的值是相同的,不同粒子的r1、r2、r3的值是不同的;如果溢出边界,就将当前的值设置为边界值;
S22:将S1得到的聚类的结果输入到ANN-HPSO-CA中进行处理,得到每一个分区的转换规则。
进一步,所述ANN-HPSO-CA为:在HPSO中引入CA,通过CA提高搜索能力来挖掘更多潜在搜索空间,使得粒子有更多的天赋来摆脱原来的轨迹;
在HPSO的第一层循环中,所有搜索空间中的粒子被放在一个给定的二维网格中,对于每一个粒子,CA计算找出当前粒子邻域区域内的局部最好粒子和局部最坏粒子并设定为对应的初始值,然后计算当前粒子的适应度值,如果当前粒子的适应度值比局部最好粒子的适应度值更高,就根据CA转换规则对当前粒子进行更新;
HPSO的第二层循环将ANN的权重和偏差参数初始化,赋值给粒子,在HPSO-CA中找出最优粒子;在计算粒子的适应度值时,将神经网络中的代价函数应用于HPSO中的适应度函数,最后,全局最优粒子即为神经网络中的最优权值和偏差。
进一步,所述S3具体包括以下步骤:
S31:首先通过基期年的土地利用数据与末期年的土地利用数据得到土地利用变化数据,然后通过标准格网划分的方法将研究区域划分为很多标准格网,得到一个速率格网;
S32:通过结合土地利用变化数据按照以下公式得到元胞转换速率:
其中Vi为速率格网i中的元胞演化速率,Ci为速度格网i中发生变化的栅格(元胞)数目,Si为速率格网i中非空栅格(元胞)的数目;
S33:通过下列公式计算元胞的演化间隔:
其中Intervali即为速率格网Vi中的元胞演化间隔;
S34:定义最大元胞演化间隔MaxI,然后把所有元胞的演化间隔都控制在[1,MaxI]区间中,用下列的公式展示:
即,如果元胞的演化间隔大于MaxI,则该元胞的演化间隔就为MaxI,否则就为Intervali;
S35:异步元胞自动机CA模型中每个元胞的总体转换概率由五个部分组成:城市发展适宜性Ps、邻域效应Ω、约束条件Rc、随机因子RA以及通过元胞演化间隔得到的所述城市发展适宜性通过ANN-HPSO-CA模型计算得到;所述邻域效应代表n×n邻域范围内的中心元胞与周边元胞的相互作用;具体表现为中心元胞n×n邻域范围内城市元胞的密度;元胞i在时刻t的邻域效应表示为:
限制因素约束了特殊土地利用类型向城市的转换;限制因素表示为:
其中当元胞i在时刻t满足扩张约束条件时,con()取值为1,否则取值为0;RA为0-1之间的随机数;
元胞i在时刻t的总体转换概率表示为:
在每次迭代中计算每个元胞的总体转换概率,并根据概率值选择最大的前N个非城市元胞进行转换;根据基期年和末期年的土地利用变化数据得到城市扩张量M,模型迭代次数则为即达到实际城市扩张量时模型停止模拟,输出模拟结果。
本发明的有益效果在于:本发明将一种与改进ANN与异步元胞自动机CA模型相结合,通过使用粒子群结合CA优化之后的ANN提取转换规则,并且考虑了因地理位置和经济差异所导致的空间异质性与时间异质性,更加符合真实的城市扩张模拟,能够提升城市扩张模拟精度高,为相关领域决策者提供了数据支撑和参考。
本发明的其他优点、目标和特征在某种程度上将在随后的说明书中进行阐述,并且在某种程度上,基于对下文的考察研究对本领域技术人员而言将是显而易见的,或者可以从本发明的实践中得到教导。本发明的目标和其他优点可以通过下面的说明书来实现和获得。
附图说明
为了使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明作优选的详细描述,其中:
图1为基于改进ANN耦合异步元胞自动机的诚实动态扩张模拟***结构示意图。
具体实施方式
以下通过特定的具体实例说明本发明的实施方式,本领域技术人员可由本说明书所揭露的内容轻易地了解本发明的其他优点与功效。本发明还可以通过另外不同的具体实施方式加以实施或应用,本说明书中的各项细节也可以基于不同观点与应用,在没有背离本发明的精神下进行各种修饰或改变。需要说明的是,以下实施例中所提供的图示仅以示意方式说明本发明的基本构想,在不冲突的情况下,以下实施例及实施例中的特征可以相互组合。
其中,附图仅用于示例性说明,表示的仅是示意图,而非实物图,不能理解为对本发明的限制;为了更好地说明本发明的实施例,附图某些部件会有省略、放大或缩小,并不代表实际产品的尺寸;对本领域技术人员来说,附图中某些公知结构及其说明可能省略是可以理解的。
本发明实施例的附图中相同或相似的标号对应相同或相似的部件;在本发明的描述中,需要理解的是,若有术语“上”、“下”、“左”、“右”、“前”、“后”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此附图中描述位置关系的用语仅用于示例性说明,不能理解为对本发明的限制,对于本领域的普通技术人员而言,可以根据具体情况理解上述术语的具体含义。
请参阅图1,为基于改进ANN与分区异步元胞自动机模型的城市动态扩张模拟方法,该方法包括以下步骤:
S1:通过K-means方法对研究区域进行分区;
S2:通过基于ANN-HPSO-CA模型的到每个分区的转换规则;
S3:将ANN-HPSO-CA模型与异步元胞自动机模型耦合。
可选的,所述S1具体为:
将基期年与末期年两年的研究区域的数据进行叠加,得到土地利用变化数据,对土地利用变化数据与其他影响土地利用变化的驱动因子数据进行K-means空间聚类分析,聚类结果中的每一个类(簇)就对应元胞空间的一个分区;
可选的,所述S2具体为:
S21:基于人行为粒子群与元胞自动机模型原理:基于人行为粒子群优化和元胞自动机的改进人工神经网络是一种全新的改进人工神经网络,通过在传统粒子群优化算法的速度方程中引入了全局最坏粒子,并且赋予其随机权重,使其服从标准正态分布,该优化方案有助于平衡PSO的勘测开发能力。
考虑到我们现实生活中存在着一些有坏习惯或者行为的人,同时,我们都只到这些坏习惯或行为会给周围的人带来一些影响。如果我们都知道这些坏习惯或行为会给周围的人带来一些影响。如果我们从这些坏习惯或行为中得到警告,这对我们是有益的。相反,如果我们从这些坏习惯或行为中学习,它对我们是有害的。因此,我们必须对这些坏习惯或行为给出一个客观和理性的看法。所以作者在原PSO中引入了全局最坏粒子,在每次迭代中,它是整个种群中适应性最差的粒子,表示为Gworst:
Gworst(t)=argmax{f(Pbest1),f(Pbest2),…,f(PbestN)}
f(·)代表对应粒子的适应度值。为了模拟人的行为并充分利用最坏情况,引入学习系数r3,它服从标准正态分布;即r3属于(0,1),若r3>为0,则认为是一个有推力的学习系数,有助于提高粒子的飞行速度;因此,可以提高粒子的探测能力。反之,若r3<0,则认为是惩罚的学习系数,会降低粒子的飞行速度;因此,加大开发力度是有益的。如果r 3=0,表示这些坏习惯或行为对粒子没有影响。同时,为了降低所求解问题的参数敏感性,我们将两个加速度系数r1和r2分别替换为随机学习系数r1和r2。
r1和r2是两个值在[0,1]之间的随机数,r1+r2=1。对于一个粒子,在每一个d=1,2,…,D中,r1、r2、r3的值是相同的。但是不同粒子r1、r2、r3的值是不同的,因为每次迭代中都会重新生成。如果溢出边界,我们就将其设置为边界值。
ANN-HPSO-CA模型原理:因为粒子在搜索过程中往往会陷入局部最优,为了克服上述缺点,更好地平衡勘测开发,在HPSO中引入了CA,通过CA提高搜索能力来挖掘更多潜在搜索空间,使得粒子有更多的天赋来摆脱原来的轨迹。在HPSO的第一层循环中,所有搜索空间中的粒子被放在一个给定的二维网格中,对于每一个粒子,CA会计算找出当前粒子邻域区域内的局部最好粒子和局部最坏粒子并设定为对应的初始值,然后计算当前粒子的适应度值,如果当前粒子的适应度值比局部最好粒子的适应度值更高,那么就根据CA转换规则对当前粒子进行更新。
HPSO的第二层循环将ANN的权重和偏差参数初始化,赋值给粒子,在HPSO-CA中找出最优粒子。在计算粒子的适应度值时,将神经网络中的代价函数应用于HPSO中的适应度函数,最后,全局最优粒子也就是神经网络中的最优权值和偏差。
S22:将S1得到的聚类的结果输入到ANN-HPSO-CA中进行处理,得到每一个分区的转换规则。
可选的,所述S3具体为:
S31:首先通过基期年的土地利用数据与末期年的土地利用数据得到土地利用变化数据然后通过标准格网划分的方法将研究区域划分为很多标准格网,得到一个速率格网;
S32:通过结合土地利用变化数据按照以下公式得到元胞转换速率:
Vi为速率格网i中的元胞演化速率,Ci为速度格网i中发生变化的栅格(元胞)数目,Si为速率格网i中非空栅格(元胞)的数目。
S33:通过下列公式计算元胞的演化间隔:
其中Intervali即为速率格网Vi中的元胞演化间隔。
S34:定义了最大元胞演化间隔MaxI,然后把所有元胞的演化间隔都控制在[1,MaxI]区间中。可用下列的公式展示:
即是如果元胞的演化间隔大于MaxI,则该元胞的演化间隔就为MaxI,否则就为Intervali。
S35:异步元胞自动机CA模型中每个元胞的总体转换概率由五个部分组成:城市发展适宜性Ps、邻域效应Ω与约束条件Pc和随机因子RA以及通过元胞演化间隔得到的城市发展适宜性通过ANN-HPSO-CA模型计算得到;邻域效应代表n×n邻域范围内的中心元胞与周边元胞的相互作用;具体表现为中心元胞n×n邻域范围内城市元胞的密度;元胞i在时刻t的邻域效应表示为:
限制因素约束了特殊土地利用类型向城市的转换;限制因素表示为:
其中当元胞i在时刻t满足扩张约束条件时,con()取值为1,否则取值为0;RA为0-1之间的随机数;
元胞i在时刻t的总体转换概率表示为:
在每次迭代中计算每个元胞的总体转换概率,并根据概率值选择最大的前N个非城市元胞进行转换。根据基期年和末期年的土地利用变化数据得到城市扩张量M,模型迭代次数则为即达到实际城市扩张量时模型停止模拟,输出模拟结果。
最后说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本技术方案的宗旨和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (5)
1.一种基于改进ANN耦合异步元胞自动机的城市动态扩张模拟方法,其特征在于:包括以下步骤:
S1:通过K-means方法对研究区域进行分区;
S2:通过基于ANN-HPSO-CA模型得到每个分区的转换规则;
S3:将ANN-HPSO-CA模型与异步元胞自动机模型耦合。
2.根据权利要求1所述的基于改进ANN耦合异步元胞自动机的城市动态扩张模拟方法,其特征在于:在所述步骤S1中,将基期年与末期年两年的研究区域的数据进行叠加,得到土地利用变化数据,对所述土地利用变化数据与其他影响土地利用变化的驱动因子数据进行K-means空间聚类分析,聚类结果中的每一个类对应元胞空间的一个分区。
3.根据权利要求1所述的基于改进ANN耦合异步元胞自动机的城市动态扩张模拟方法,其特征在于:所述步骤S2具体包括以下步骤:
S21:基于人行为粒子群优化和元胞自动机的改进人工神经网络ANN-HPSO-CA,在原粒子群优化算法PSO中引入全局最坏粒子,即在每次迭代中,整个种群中适应性最差的粒子,表示为Gworst:
Gworst(t)=argmax{f(Pbest1),f(Pbest2),…,f(PbestN)}
f(·)代表对应粒子的适应度值;
引入服从标准正态分布的学习系数r3,即r3∈(0,1),若r3>0,则认为是一个有推力的学习系数,有助于提高粒子的飞行速度;若r3<0,则认为是惩罚的学习系数,降低粒子的飞行速度;若r3=0,表示这些坏习惯或行为对粒子没有影响;将两个加速度系数r1和r2分别替换为随机学习系数r1和r2:
r1和r2是两个值在[0,1]之间的随机数,r1+r2=1;对于一个粒子,在每一个d=1,2,…,D中,r1、r2、r3的值是相同的,不同粒子的r1、r2、r3的值是不同的;如果溢出边界,就将当前的值设置为边界值;
S22:将S1得到的聚类的结果输入到ANN-HPSO-CA中进行处理,得到每一个分区的转换规则。
4.根据权利要求1-3任一所述的基于改进ANN耦合异步元胞自动机的城市动态扩张模拟方法,其特征在于:所述ANN-HPSO-CA为:在HPSO中引入CA,通过CA提高搜索能力来挖掘更多潜在搜索空间,使得粒子有更多的天赋来摆脱原来的轨迹;
在HPSO的第一层循环中,所有搜索空间中的粒子被放在一个给定的二维网格中,对于每一个粒子,CA计算找出当前粒子邻域区域内的局部最好粒子和局部最坏粒子并设定为对应的初始值,然后计算当前粒子的适应度值,如果当前粒子的适应度值比局部最好粒子的适应度值更高,就根据CA转换规则对当前粒子进行更新;
HPSO的第二层循环将ANN的权重和偏差参数初始化,赋值给粒子,在HPSO-CA中找出最优粒子;在计算粒子的适应度值时,将神经网络中的代价函数应用于HPSO中的适应度函数,最后,全局最优粒子即为神经网络中的最优权值和偏差。
5.根据权利要求1所述的基于改进ANN耦合异步元胞自动机的城市动态扩张模拟方法,其特征在于:所述S3具体包括以下步骤:
S31:首先通过基期年的土地利用数据与末期年的土地利用数据得到土地利用变化数据,然后通过标准格网划分的方法将研究区域划分为很多标准格网,得到一个速率格网;
S32:通过结合土地利用变化数据按照以下公式得到元胞转换速率:
其中Vi为速率格网i中的元胞演化速率,Ci为速度格网i中发生变化的栅格(元胞)数目,Si为速率格网i中非空栅格(元胞)的数目;
S33:通过下列公式计算元胞的演化间隔:
其中Intervali即为速率格网Vi中的元胞演化间隔;
S34:定义最大元胞演化间隔MaxI,然后把所有元胞的演化间隔都控制在[1,MaxI]区间中,用下列的公式展示:
即,如果元胞的演化间隔大于MaxI,则该元胞的演化间隔就为MaxI,否则就为Intervali;
S35:异步元胞自动机CA模型中每个元胞的总体转换概率由五个部分组成:城市发展适宜性Ps、邻域效应Ω、约束条件Pc、随机因子RA以及通过元胞演化间隔得到的所述城市发展适宜性通过ANN-HPSO-CA模型计算得到;所述邻域效应代表n×n邻域范围内的中心元胞与周边元胞的相互作用;具体表现为中心元胞n×n邻域范围内城市元胞的密度;元胞i在时刻t的邻域效应表示为:
限制因素约束了特殊土地利用类型向城市的转换;限制因素表示为:
其中当元胞i在时刻t满足扩张约束条件时,con()取值为1,否则取值为0;RA为0-1之间的随机数;
元胞i在时刻t的总体转换概率表示为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110433496.1A CN113177354A (zh) | 2021-04-19 | 2021-04-19 | 基于改进ann耦合异步元胞自动机的城市动态扩张模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110433496.1A CN113177354A (zh) | 2021-04-19 | 2021-04-19 | 基于改进ann耦合异步元胞自动机的城市动态扩张模拟方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113177354A true CN113177354A (zh) | 2021-07-27 |
Family
ID=76924297
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110433496.1A Pending CN113177354A (zh) | 2021-04-19 | 2021-04-19 | 基于改进ann耦合异步元胞自动机的城市动态扩张模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113177354A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114442998A (zh) * | 2022-01-27 | 2022-05-06 | 东南大学 | 一种开源软件项目的演化建模方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112163367A (zh) * | 2020-08-20 | 2021-01-01 | 同济大学 | 融合萤火虫算法和元胞自动机的城市扩张模拟预测方法 |
CN112487724A (zh) * | 2020-12-11 | 2021-03-12 | 重庆邮电大学 | 基于分区和改进cnn-ca模型的城市动态扩张模拟方法 |
-
2021
- 2021-04-19 CN CN202110433496.1A patent/CN113177354A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112163367A (zh) * | 2020-08-20 | 2021-01-01 | 同济大学 | 融合萤火虫算法和元胞自动机的城市扩张模拟预测方法 |
CN112487724A (zh) * | 2020-12-11 | 2021-03-12 | 重庆邮电大学 | 基于分区和改进cnn-ca模型的城市动态扩张模拟方法 |
Non-Patent Citations (6)
Title |
---|
YUE WANG: "An improved artificial neural network based on human-behaviour particle swarm optimization and cellular automata", 《HTTPS://WWW.SCIENCEDIRECT.COM/SCIENCE/ARTICLE/PII/S095741741930572X》 * |
井长青: "耦合神经网络与元胞自动机的城市土地利用动态演化模型", 《干旱区域研究 》 * |
刘明皓: "ANN-CA模型改进及城市用地动态模拟", 《重庆邮电大学学报(自然科学版)》 * |
李沁等: "基于改进粒子群算法和元胞自动机的城市扩张模拟――以南京为例", 《长江流域资源与环境》 * |
柯新利等: "基于空间数据挖掘的分区异步元胞自动机模型研究", 《中国图象图形学报》 * |
王月: "改进的粒子群优化算法及其应用研究", 《中国优秀硕士学位论文全文数据库 (信息科技辑)》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114442998A (zh) * | 2022-01-27 | 2022-05-06 | 东南大学 | 一种开源软件项目的演化建模方法 |
CN114442998B (zh) * | 2022-01-27 | 2024-07-19 | 东南大学 | 一种开源软件项目的演化建模方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Gharehchopogh | An improved Harris Hawks optimization algorithm with multi-strategy for community detection in social network | |
Guan et al. | An artificial-neural-network-based, constrained CA model for simulating urban growth | |
CN112884226A (zh) | 基于多智能体算法的国土空间格局模拟规划方法及*** | |
CN111260118A (zh) | 一种基于量子粒子群优化策略的车联网交通流量预测方法 | |
CN110473592A (zh) | 基于图卷积网络的有监督的多视角人类协同致死基因预测方法 | |
CN110164129B (zh) | 基于gernn的单交叉口多车道交通流量预测方法 | |
CN108090624B (zh) | 一种改进元胞自动机的城市生态安全模拟与预测的方法 | |
CN110826244B (zh) | 模拟轨道交通影响城市生长的共轭梯度元胞自动机方法 | |
CN109492796A (zh) | 一种城市空间形态自动分区方法与*** | |
CN117313931B (zh) | 一种基于火灾场景下地铁站台乘客疏散时间的预测方法 | |
CN116401941B (zh) | 一种地铁站闸机疏散能力的预测方法 | |
Li et al. | Genetic algorithms for determining the parameters of cellular automata in urban simulation | |
US20230020112A1 (en) | Relating complex data | |
CN116756828A (zh) | 一种城市空间规划方法、介质及*** | |
CN109840551A (zh) | 一种用于机器学习模型训练的优化随机森林参数的方法 | |
Hu et al. | Dynamically Optimized Unstructured Grid (DOUG) for Analog Ensemble of numerical weather predictions using evolutionary algorithms | |
CN112487724A (zh) | 基于分区和改进cnn-ca模型的城市动态扩张模拟方法 | |
CN117556979B (zh) | 基于群体智能搜索的无人机平台与载荷一体化设计方法 | |
CN113177354A (zh) | 基于改进ann耦合异步元胞自动机的城市动态扩张模拟方法 | |
Cai et al. | The mining of urban hotspots based on multi-source location data fusion | |
CN112633591B (zh) | 一种基于深度强化学习的空间搜索方法及设备 | |
CN113743659A (zh) | 一种基于成分法和马尔可夫元胞自动机的城市布局预测方法及应用 | |
CN115691140B (zh) | 一种汽车充电需求时空分布的分析与预测方法 | |
CN116167254A (zh) | 基于城市大数据的多维城市仿真推演方法和*** | |
CN106649521A (zh) | 一种基于任务负载的数据迁移方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210727 |