CN102254106A - City extension simulating method based on cellular automata - Google Patents
City extension simulating method based on cellular automata Download PDFInfo
- Publication number
- CN102254106A CN102254106A CN2011102274233A CN201110227423A CN102254106A CN 102254106 A CN102254106 A CN 102254106A CN 2011102274233 A CN2011102274233 A CN 2011102274233A CN 201110227423 A CN201110227423 A CN 201110227423A CN 102254106 A CN102254106 A CN 102254106A
- Authority
- CN
- China
- Prior art keywords
- cellular
- land
- neighbours
- road
- unused
- 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
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a city extension simulating method based on a cellular automaton. In the city extension simulating method, the iterative process of the tradition cellular automaton is improved, the cell neighbor space is layered, and neighbors in all directions and other factors influencing cell transformation are showed in a matrix form. The city extension predicting method disclosed by the invention has remarkable advantages: firstly, the method is superior to the traditional point-by-point cell judgment method, and the prediction efficiency is increased remarkably; in addition, the matrix form is introduced to the iterative process, thus the distance factor of a city extension system is simplified, and the attribute judgment of a central cell is simplified, and the simulation result is totally consistent to that of the traditional traversing method.
Description
Technical field
The present invention relates to city expansion simulation field, relate in particular to the city expansion analogy method of a kind of employing based on cellular automaton.
Background technology
Cellular automaton (cellular automata is called for short CA) is that a kind of time, space and state all disperse, and has powerful spatial modeling ability and arithmetic capability, can simulate the dynamic system with space-time characteristic
[1], and extensively applied to the simulation that the city is expanded.But cellular automaton causes the standard cellular automaton to adopt the mode of " traversal " cellular unit all the time in iterative process with characteristic and " from top to bottom " operating mechanism thereof of simple local rule Simulation of Complex phenomenon
[2]This kind alternative manner simple possible, be convenient to move on computers, but when cellular quantity is big, adopt the mode inefficiency of traversal, and be not easy to consider distance factor, when taking geographical phenomenon into account and be subjected to many-sided factor affecting, also acquire a certain degree of difficulty, influenced the popularization that the ingot automaton model is used.And city system dynamic evolution process has the complicacy of height, and be subjected to multifactorial influences such as nature, society, economy, culture, politics, the introducing of a large amount of factor factors and distance factor makes and utilizes the inefficiency problem of cellular automaton simcity expansion more obvious.
The iterative process of tradition cellular automaton is specific as follows:
If cellular space
, wherein,
Be cellular space size,
Expression the
iRow, the
jThe cellular of row,
Cellular next moment state representation in center is
, wherein
,
Be the discrete time,
,
For cellular exists
,
State constantly,
Be neighbours
[1]Then be in size
The cellular space in, the judgement that utilizes traditional cellular Automation Model to finish a step-length need be carried out
Inferior
Functional operation is to judge next moment state of each cellular, and therefore when cellular quantity was big, iteration efficient was lower.
The list of references that relates in the literary composition is as follows:
[1] Zhou Chenghu, Sun Zhanli, Xie Yichun. geographical cellular automaton research [M]. Beijing: Science Press, 1999.
[2] Li Xia, Ye Jiaan, Liu Xiaoping etc. geographical simulation system: cellular automaton and multiple agent [M]. Beijing: Science Press, 2007.
Summary of the invention
At the deficiency that has inefficiency in the prior art based on the city expansion analogy method of traditional cellular automaton, the iterative process that the present invention is directed to traditional cellular automaton is improved, layering is carried out in cellular neighbours space, the neighbours that represent each orientation with the form of matrix influence the factor factor of cellular conversion with other, and provide a kind of based on expansion analogy method in the city improved cellular automaton of above-mentioned matrix operation, high efficiency.
In order to solve the problems of the technologies described above, the present invention adopts following technical scheme:
A kind of city expansion analogy method based on the improved cellular automaton of matrix may further comprise the steps:
S1 is according to the city simulation drawing, urban land is divided into industrial land, merchant takes six types of ground, residential land, water body, road and unused lands, marks off according to city planning and plans the key area;
The definition of S2 land used cellular attribute layer:
Calculate the distance of each unused land cellular respectively apart from one-level road and secondary road
,
, and definition land used cellular is apart from the road distance attribute
, wherein,
Be sign of operation, can determine according to different cities Characteristics of Development and condition;
Be cellular space size; Right simultaneously
Carry out normalization and obtain the land used cellular apart from road distance attribute layer
Wherein,
rRadius for neighbours' scope;
,
Represent the in neighbours' scope respectively
Row, the
The cellular of row,
,
Expression
Constantly
Place's cellular
Individual neighbours' cellular,
,
Expression
Constantly
Place's cellular
Individual neighbours' cellular, from top to bottom
Increase successively, from left to right
Increase successively;
Figure 1 shows that neighbours' position view, Figure 2 shows that to be in the center cellular
Neighbours' matrix of position, the upper left corner
S4 is right
Moment neighbours' matrix
Carry out the iteration of a step-length, obtain the unused land cellular and exist
State constantly
, wherein,
The representing matrix computing;
S5 is according to neighbours' matrix
, calculate the transition probability that the unused land cellular is converted into urban land
:
Wherein,
Represent the distance of some neighbours' matrixes apart from the center cellular,
Expression
tThe set of all neighbours' matrixes constantly,
lExpression the
lIndividual neighbours' cellular,
,
,
Being sign of operation, is any computing in addition subtraction multiplication and division or other complex calculation, can determine according to different cities Characteristics of Development and condition;
S6 works as
The time,
, unused land is converted into construction land; Otherwise,
, unused land does not transform; Wherein,
Decide by urban development speed.
Adopt the conversion of map algebraically distance of obstacle among the above-mentioned steps S2-1, calculate the distance of each land used cellular respectively apart from one-level road and secondary road
,
Compared with prior art, the present invention has following advantage and beneficial effect:
The inventive method is improved at the iterative process of traditional cellular automaton, and layering is carried out in cellular neighbours space, represents that with the form of matrix the neighbours in each orientation influence the factor factor of cellular conversion with other.Expansion Forecasting Methodology in city of the present invention has significant advantage, and at first, this method is better than the mode that cellular is judged in traditional pointwise, has significantly improved forecasting efficiency; In addition, introduced matrix in iterative process, simplified the difficulty of considering the distance factor in the expanding system of city and having simplified center cellular property determine, analog result and traditional traversal method are in full accord.
Description of drawings
Fig. 1 is mole type neighbours position views;
Fig. 2 is neighbours' matrix synoptic diagram;
Fig. 3 is the simcity synoptic diagram in certain city;
Fig. 4 be embodiment apart from the attribute synoptic diagram, wherein, figure (a) be among the embodiment each cellular apart from the one-level road apart from the attribute synoptic diagram, scheme (b) be among the embodiment each cellular apart from secondary road apart from the attribute synoptic diagram;
Fig. 5 is consuming time-neighbours' radius size fitted figure of embodiment.
Embodiment
Figure 3 shows that the simcity synoptic diagram in certain city, will be example with simcity synoptic diagram shown in Figure 3 below, adopt the inventive method respectively and come the simcity expansion based on the method for traditional cellular automaton.Comprise in the simcity that industrial land, merchant take ground, residential land, water body, road and unused land several types, mark off the planning key area according to city planning simultaneously; Choose the land used cellular apart from road distance and planning factors as the cellular attribute; Road changes the bigger and road influence degrees at different levels of the influence of construction land into to unused land and varies in size, and is transformed to the basis with map algebraically distance of obstacle, calculates each cellular respectively apart from the one-level road
And the distance of secondary road
, as shown in Figure 4; Definition land used cellular is apart from the road distance attribute
, represent the combined influence of roads at different levels to the land used cellular, right simultaneously
Carry out normalization, obtain
The attribute layer, wherein,
Be sign of operation, can be any computing in addition subtraction multiplication and division or other complex calculation, can determine according to different cities Characteristics of Development and condition.The artificial interference is one of factor of influence very important in the urban development, the limitation attribute of definition planning here layer
As follows:
In the method simcity expansion of employing based on traditional cellular automaton, consider planning factors and apart from the influence of road distance, judge that simultaneously neighbours' cellular is the construction land number, and adopt random number to judge whether the center cellular is converted into construction land, so that the result meets the randomness of real world development.Specific algorithm is that to last row of last column, each iteration judges that all whether the center cellular is unused land, if then carry out next step judgement, otherwise jumps into another circulation from first row, the first row iteration; If for unused land calculate its neighbours' cellular and other factor factor pairs its be converted to the effect degree of urban land
S I, j T+1 'If this probability is greater than a certain random number, then the center element dysuria with lower abdominal colic turns to urban land.
The inventive method is an example with mole type (Moore type) neighbours and adiabatic border, at first defines neighbours' matrix and is
, then the definable transformation rule is:
Wherein
,
,
Be sign of operation, can be any computing in addition subtraction multiplication and division or other complex calculation, can determine according to different cities Characteristics of Development and condition.Try to achieve according to following formula
Be concrete numeric type, need corresponding with cellular state.Here adopt stochastic matrix to determine cellular state:
When the matrix operation result more than or equal to
The time, unused land is converted into construction land, otherwise does not change.Wherein,
Decide by urban development speed.Table 1 is depicted as the comparison consuming time of adopting the inventive method and the expansion of traditional cellular automaton simcity.
The comparison consuming time of table 1 the inventive method and traditional cellular automation method
Adopt least square method to obtain the inventive method and traditional cellular automation method quadratic fit function consuming time, as follows:
Wherein,
Be the consuming time of traditional cellular automation method;
Be the consuming time of the inventive method;
It is neighbours' radius.
Table 2 is consuming time-neighbours' radius size match table
? | Method based on traditional cellular automaton | The inventive method |
Constant term | 2.267 | -3.67×10 -2 |
Once | 8.79×10 -2 | 7.32×10 -2 |
Quadratic term | — | 7.50×10 -2 |
Residual | 0.0643 | 0.0022 |
From table 1 and fitting result shown in Figure 5, at radius
The time, adopt the inventive method iteration once take time less than method based on traditional cellular automaton.But along with the increase of radius, the once required time of the inventive method iteration significantly increases.Main cause is, though need to adopt the method for traversal based on the method for traditional cellular automaton, causes its iteration time once longer, and the increase of radius is very little to the influence of traversal; And the inventive method needs to redefine each neighbours' matrix in each step-length, and along with the increase of radius, needs neighbours' matrix of definition to be quadratic power and increases, and is as shown in table 2, so trend along with the radius rising consuming time is obvious.From fitting result, based on being once linear relationship between the consuming time and radius size of the method for traditional cellular automaton, and be quadratic function relation between the consuming time and radius size of the inventive method.Because in actual applications, choose r=2 or 3 usually,
Situation very rare, so in the general applied research, if
The time, the inventive method still is better than the method based on traditional cellular automaton.
Claims (2)
1. the city expansion analogy method based on cellular automaton is characterized in that, may further comprise the steps:
S1 is according to the city simulation drawing, urban land is divided into industrial land, merchant takes six types of ground, residential land, water body, road and unused lands, marks off according to city planning and plans the key area;
The definition of S2 land used cellular attribute layer:
Calculate the distance of each unused land cellular respectively apart from one-level road and secondary road
,
, and definition land used cellular is apart from the road distance attribute
, wherein,
Be sign of operation, can determine according to different cities Characteristics of Development and condition;
Be cellular space size; Right simultaneously
Carry out normalization and obtain the land used cellular apart from road distance attribute layer
Wherein,
rRadius for neighbours' scope;
,
Represent the in neighbours' scope respectively
Row, the
The cellular of row,
,
Expression
Constantly
Place's cellular
Individual neighbours' cellular,
,
Expression
Constantly
Place's cellular
Individual neighbours' cellular, from top to bottom
Increase successively, from left to right
Increase successively;
S4 is right
Moment neighbours' matrix
Carry out the iteration of a step-length, obtain the unused land cellular and exist
State constantly
, wherein,
The representing matrix computing;
S5 is according to neighbours' matrix
, calculate the transition probability that the unused land cellular is converted into urban land
:
Wherein,
Represent the distance of some neighbours' matrixes apart from the center cellular,
Expression
tThe set of all neighbours' matrixes constantly,
lExpression the
lIndividual neighbours' cellular,
,
,
Be sign of operation, can determine according to different cities Characteristics of Development and condition;
2. the city expansion analogy method based on cellular automaton according to claim 1 is characterized in that:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011102274233A CN102254106A (en) | 2011-08-10 | 2011-08-10 | City extension simulating method based on cellular automata |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011102274233A CN102254106A (en) | 2011-08-10 | 2011-08-10 | City extension simulating method based on cellular automata |
Publications (1)
Publication Number | Publication Date |
---|---|
CN102254106A true CN102254106A (en) | 2011-11-23 |
Family
ID=44981367
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2011102274233A Pending CN102254106A (en) | 2011-08-10 | 2011-08-10 | City extension simulating method based on cellular automata |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102254106A (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104156537A (en) * | 2014-08-19 | 2014-11-19 | 中山大学 | Cellular automaton urban growth simulating method based on random forest |
CN104992041A (en) * | 2015-08-06 | 2015-10-21 | 武汉大学 | City expansion boundary prediction method based on space syntax |
CN106251014A (en) * | 2016-07-29 | 2016-12-21 | 西南交通大学 | Development of urban space based on SVM GCA simulation and Forecasting Methodology |
CN110009257A (en) * | 2019-04-17 | 2019-07-12 | 青岛大学 | Multiple dimensioned variable window cellular Automation Model based on urban traffic blocking sprawling analysis |
CN110688726A (en) * | 2019-08-19 | 2020-01-14 | 华南师范大学 | Spatial orientation self-adaptive city expansion simulation method, system and storage medium |
CN111797523A (en) * | 2020-06-30 | 2020-10-20 | 南京图申图信息科技有限公司 | Land use change-oriented multi-level vector cellular automaton modeling method |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060058968A1 (en) * | 2004-09-14 | 2006-03-16 | Gangadhar Deepak K | Protein structure alignment using cellular automata |
CN1987906A (en) * | 2005-12-22 | 2007-06-27 | 中国煤炭地质总局航测遥感局 | Method for dynamicaly predicting land use change |
CN101436204A (en) * | 2008-12-12 | 2009-05-20 | 中山大学 | City evolvement simulation implementing method based on paralleling elementary cell automatic machine |
-
2011
- 2011-08-10 CN CN2011102274233A patent/CN102254106A/en active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060058968A1 (en) * | 2004-09-14 | 2006-03-16 | Gangadhar Deepak K | Protein structure alignment using cellular automata |
CN1987906A (en) * | 2005-12-22 | 2007-06-27 | 中国煤炭地质总局航测遥感局 | Method for dynamicaly predicting land use change |
CN101436204A (en) * | 2008-12-12 | 2009-05-20 | 中山大学 | City evolvement simulation implementing method based on paralleling elementary cell automatic machine |
Non-Patent Citations (2)
Title |
---|
王海军 等: "利用地图代数和数据场拓展元胞自动机理论", 《武汉大学学报·信息科学版》, vol. 35, no. 12, 31 December 2010 (2010-12-31), pages 1474 - 1477 * |
罗平 等: "地理特征元胞自动机及城市土地利用演化研究", 《武汉大学学报·信息科学版》, vol. 29, no. 6, 30 June 2004 (2004-06-30) * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104156537A (en) * | 2014-08-19 | 2014-11-19 | 中山大学 | Cellular automaton urban growth simulating method based on random forest |
CN104992041A (en) * | 2015-08-06 | 2015-10-21 | 武汉大学 | City expansion boundary prediction method based on space syntax |
CN104992041B (en) * | 2015-08-06 | 2017-11-28 | 武汉大学 | A kind of urban sprawl Boundary Prediction method based on Space Syntax |
CN106251014A (en) * | 2016-07-29 | 2016-12-21 | 西南交通大学 | Development of urban space based on SVM GCA simulation and Forecasting Methodology |
CN110009257A (en) * | 2019-04-17 | 2019-07-12 | 青岛大学 | Multiple dimensioned variable window cellular Automation Model based on urban traffic blocking sprawling analysis |
CN110009257B (en) * | 2019-04-17 | 2023-09-08 | 青岛大学 | Multi-scale variable window cellular automaton model based on urban traffic congestion spreading analysis |
CN110688726A (en) * | 2019-08-19 | 2020-01-14 | 华南师范大学 | Spatial orientation self-adaptive city expansion simulation method, system and storage medium |
CN110688726B (en) * | 2019-08-19 | 2023-01-10 | 华南师范大学 | Spatial orientation self-adaptive city expansion simulation method, system and storage medium |
CN111797523A (en) * | 2020-06-30 | 2020-10-20 | 南京图申图信息科技有限公司 | Land use change-oriented multi-level vector cellular automaton modeling method |
CN111797523B (en) * | 2020-06-30 | 2024-03-22 | 南京图申图信息科技有限公司 | Land utilization change-oriented multi-level vector cellular automaton modeling method |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102254106A (en) | City extension simulating method based on cellular automata | |
CN101435708B (en) | Real time roaming and scaling method of vector map in mobile navigation equipment | |
Wang et al. | A novel approach for the layout of electric vehicle charging station | |
Hong et al. | A decision support model for reducing electric energy consumption in elementary school facilities | |
Amado et al. | Solar energy integration in urban planning: GUUD model | |
Liu et al. | Impact of urban form on building energy consumption and solar energy potential: A case study of residential blocks in Jianhu, China | |
Liu et al. | [Retracted] Digital Twin‐Based Ecogreen Building Design | |
Wang et al. | Positivity and difference of influence of built environment around urban park on building energy consumption | |
Liu et al. | A multi-objective optimization framework for designing urban block forms considering daylight, energy consumption, and photovoltaic energy potential | |
CN107730027A (en) | A kind of load forecasting method and device based on region building property | |
CN102254105A (en) | Urban sprawl forecasting method based on cloud model cellular automata | |
CN110288133A (en) | Planning substation automatic addressing method based on distant view year saturation loading distribution map | |
CN106355271A (en) | Optimal spatial allocation method of land resources applied to land use planning | |
AU2021102983A4 (en) | Method and system for measuring accessibility of ecological effect of green space on building energy consumption carbon emission reduction | |
Zhang | RETRACTED: Energy-saving design of civil engineering buildings based on FPGA and embedded system | |
Shuai et al. | Does ecological network affect the construction land spatial development intensity?-An empirical study based on the main urban area of Jinan, China | |
CN109934513B (en) | Irregular airport facing industrial area layout system and method based on multi-agent evolution algorithm | |
Ascione et al. | Energy assessment in town planning: urban energy maps | |
CN103195161A (en) | Construction method of basic module unit of building structure body | |
Javanroodi et al. | Evaluating the impacts of urban form on the microclimate in the dense areas | |
Ghasemi Sangi et al. | Investigation of the Height Distribution Effect of Residential Complex Blocks on Optimization of Cooling and Heating Loads (Tehran, District 9) | |
Zheng et al. | Extension and application of space syntax a case study of urban traffic network optimizing in Beijing | |
Yang et al. | Progress on geographic cellular automata model | |
CN112365072B (en) | Integrated design scheme generation method and system for park integrated energy system | |
Tabkhi et al. | Investigating the feasibility of using nanomaterials in electricity and energy systems for urban use |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20111123 |