CN106251020A - Land_use change layout method for optimizing based on resosurces environment loading capacity - Google Patents

Land_use change layout method for optimizing based on resosurces environment loading capacity Download PDF

Info

Publication number
CN106251020A
CN106251020A CN201610634963.6A CN201610634963A CN106251020A CN 106251020 A CN106251020 A CN 106251020A CN 201610634963 A CN201610634963 A CN 201610634963A CN 106251020 A CN106251020 A CN 106251020A
Authority
CN
China
Prior art keywords
land
construction land
construction
use change
layout
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
Application number
CN201610634963.6A
Other languages
Chinese (zh)
Inventor
辜寄蓉
刘寅
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Information Center Of Sichuan Provincial Department Of Land And Resources
Original Assignee
Information Center Of Sichuan Provincial Department Of Land And Resources
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Information Center Of Sichuan Provincial Department Of Land And Resources filed Critical Information Center Of Sichuan Provincial Department Of Land And Resources
Priority to CN201610634963.6A priority Critical patent/CN106251020A/en
Publication of CN106251020A publication Critical patent/CN106251020A/en
Pending legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/29—Geographical information databases
    • 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
    • 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
    • G06Q10/00—Administration; Management
    • G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G06Q10/043—Optimisation of two dimensional placement, e.g. cutting of clothes or wood
    • 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—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/10—Services
    • G06Q50/16—Real estate
    • G06Q50/165—Land development

Abstract

The invention discloses a kind of Land_use change layout method for optimizing based on resosurces environment loading capacity, including: obtain survey region Land_use change topology data；Land_use change topology data is carried out spatialization process, forms Land_use change topology data spatialization evaluation index；According to Land_use change topology data spatialization evaluation index, set up resosurces environment loading capacity assessment indicator system；Use Field Using Fuzzy Comprehensive Assessment that resosurces environment loading capacity assessment indicator system carries out multidomain treat-ment, generate resosurces environment loading capacity three regions；Carry out construction land expansion simulation based on spatial database, form construction land suitability mode of extension；Construction land suitability mode of extension and resosurces environment loading capacity partitioned area are carried out spatial overlay analysis, according to analysis result, it is achieved preferred to Land_use change layout.

Description

Land_use change layout method for optimizing based on resosurces environment loading capacity
Technical field
Present invention relates particularly to a kind of Land_use change layout method for optimizing based on resosurces environment loading capacity.
Background technology
Land use planning is under the guidance of the idea of sustainable development, utilizes with Rational Land and protects, with full Foot human economic society growth requirement is the planning of basic goal, is to the land development in certain area coverage, utilizes, administers The general arrangement made on room and time with protection and layout, its open space planning relates to resources protection and development, resource cloth Office and the aspect such as structural adjustment, major project delimitation, and focus on the pass that is optimized for of land use structure adjustment and space layout Point, i.e. based on the visual angle producing, living, the ecological three lives coordinates, rationally determines the ratio of farming land, construction land, ecological land Relation and space layout；In general, the core of land use planning is to alleviate day by day fierce people-ground contradiction, meets warp Ji is social fast-developing and population increases the demand to land used, and construction land is carrying socio-economic development achievement and mankind's life Produce the main carriers lived, therefore coordinating production, on the basis of ecological land, it is achieved the optimization layout of construction land, be The important content of overall plan for land use, is also to solve the important channel of all kinds of problems in urbanization process.
Along with the propelling of urbanization process, construction land layout optimization is increasingly becoming the much-talked-about topic of research；1989, suddenly China's moral proposes Garden City theory, advocates Urban-rural Coordinated Development, serves important to the construction land layout in city ecotope Directive function, nineteen thirty-five relies special wide mu City Theory and " Athens Charter " within 1978, promulgated then to build in urban development If Land arrangement is expanded on further with the mutual relation of ecological environment；The nineties in 20th century, for solve city low-density without Sequence such as spreads at the urban issues, and the theory such as " compact development, astute increase " proposes in succession, delimit urban growth boundary, builds many fingers Target urban development Benefit evaluation system becomes the important means optimizing construction land layout；One-tenth along with computer technology Ripe, combine with GIS technology at the computer model such as cellular automata, neutral net and carry out the basis of construction land layout prediction On carry out construction land optimizing research and become new focus, rationalization and science for construction land layout optimization provide Technical support.
The Oregon of U.S.'s Portland proposes a set of land used preferred version, the Oregon area of Portland vast but / 3rd of local not enough Quanzhou that suitable people occupies, almost all population all concentrates in Williamette river valley, and population Nesting Zone is also agricultural development area, the lance under quick urbanization process, between urban construction, environmental conservation and agricultural development Shield is increasingly sharp-pointed, and formulating rational urban growth boundary becomes important content and the governmental duties of national planning；Ripple in 2004 The special blue delimitation result to urban growth boundary in 1979 is corrected, i.e. on the basis of existing urban growth boundary, if Four kinds of different development sights are determined: basic sight was i.e. increased by former pattern；Sight A advocates suitably expansion limit, existing city Boundary, and carry out greater compactness of exploitation at light rail and public bus network periphery；Sight B advocates to carry out on the basis of existing urban border Greater compactness of exploitation and the more effective utilization of soil traffic system；Sight C is Pattern of satellite town, advocates surrounding cities, Metropolitan Area The potential region increased as city；On the basis of carrying out situation setting and city's spatial structure prediction, for selecting optimum Urban development pattern, Portland proposes 46 concrete evaluation indexes in terms of eight, traffic, land used, air quality, water etc., and 46 indexs under different situations have been carried out detailed measuring and calculating and relative analysis, the distinct quality of various developmental patterns, thus Investigated by the public, the public select the most acceptable scheme, mix above multiple scheme on this basis and determine final city Growth border, city.The defect of this kind of method is: 1, selection indicators is more complicated, is difficult to value；2, subjective guideline is more, root According to the Characteristic Design of the U.S., it is not suitable for China's actual conditions.
Long Ying etc. use restrictive CA method to be simulated Beijing's urban morphology, set under four kinds of different sights Urban morphology: planning sight (control centre city, put forth effort develop new city), spread sight (the urban development mould of " stand large flat bread " Formula), grape cluster sight (maintaining on the basis of existing development scale, promote samll cities and towns and the exploitation of Along Traffic Trunk Lines), ecological protect Protect district's sight (avoiding on stream, to forbidding that construction district and basic farmland take, promoting the sustainable development in city)； On the basis of Scene Simulation, Long Ying etc. uses scenario analysis, by CA, Logistic and MonoLoop method, by often Secondary circulation increases the policy parameters such as cellular quantity, field policy coefficient, Kappa regression coefficient newly to the urban morphology under different sights Feasibility carry out the relative analysis of transverse and longitudinal, so that it is determined that the urban morphology of optimum；The defect of this kind of method is: only consider Urban Expansion, simulates multiple urban diffusion model, but is not extended the preferred of scheme.
Summary of the invention
For above-mentioned deficiency of the prior art, the invention provides a kind of Land_use change based on resosurces environment loading capacity Layout method for optimizing, the method can carry out the optimized placement scheme of Land_use change under the constraint of resosurces environment loading capacity Select, to determine optimal Land utilization's mode.
In order to reach foregoing invention purpose, the technical solution used in the present invention is:
A kind of Land_use change layout method for optimizing based on resosurces environment loading capacity is provided, including:
Obtain survey region Land_use change topology data；Described Land_use change topology data includes statistical data and space number According to；
Described Land_use change topology data is carried out spatialization process, forms Land_use change topology data spatialization evaluation and refer to Mark, sets up spatial database；
According to Land_use change topology data spatialization evaluation index, according to the intension criterion of resosurces environment loading capacity, set up Resosurces environment loading capacity assessment indicator system；
Use Field Using Fuzzy Comprehensive Assessment that described resosurces environment loading capacity assessment indicator system carries out multidomain treat-ment, generate money Environmental carrying capacity three regions, source, i.e. resource environment high-bearing capacity region, middle bearing capacity region and low bearing capacity region；
Carry out construction land expansion simulation based on described spatial database, form construction land suitability mode of extension；
Described construction land suitability mode of extension and resosurces environment loading capacity partitioned area are carried out spatial overlay analysis, According to analysis result, it is achieved preferred to Land_use change layout.
Described Land_use change topology data is carried out spatialization process, forms Land_use change topology data spatialization evaluation and refer to Mark, the concrete steps setting up spatial database include:
S11, the original fishing net of generation covering survey region；
S12, the original fishing net and the Land_use change topology data that generate are carried out space overlapping, generate in original fishing net each The spatial object value of individual grid, obtains population density, economic density, land resource enabling capabilities, Water resources supporting capacity, friendship Logical distribution of facilities, water conservancy projects distribution and afforestation rate evaluation space index；
S13, according in step S12 generate evaluation space index, according to windy and sandy soil data figure spot area and The area accounting of corresponding survey region similar cover type, is calculated the spatial density of each index, and computing formula is as follows:
SL=ZB × ZL
Wherein, SL is the index space density of each figure speckle carrying, and ZL is certain space-like object space at this grid Object value, ZB is the spatial area that certain space-like object is occupied in original grid；
S14, the index space density obtaining step S13 are standardized processing, and form Land_use change topology data empty Between change evaluation index.
According to Land_use change topology data spatialization evaluation index, according to the intension criterion of resosurces environment loading capacity, set up Concretely comprising the following steps of resosurces environment loading capacity assessment indicator system:
Intension criterion according to resosurces environment loading capacity, is classified as social economy's bearing capacity, Resources Carrying Capacity and environment Bearing capacity；
By Land_use change topology data spatialization evaluation index respectively according to social economy's bearing capacity, Resources Carrying Capacity and ring The intension criterion of border bearing capacity, sets up resosurces environment loading capacity assessment indicator system:
Economic density and population density index are divided into social economy's load bearing principle layer；
Land resource enabling capabilities, Water resources supporting capacity, means of transportation distributing and water conservancy projects distribution index are divided For Resources Carrying Capacity rule layer；
Index of forest coverage is divided into environmental carrying capacity rule layer.
Carry out construction land expansion simulation based on described spatial database, form the tool of construction land suitability mode of extension Body step is:
S21, determine construction land expansion change variable: land use pattern be changed to construction land and non-constructive land Between conversion,For the Land_use change state of t ij position, by the Land_use change state in t+1 momentTo t Land_use change stateChange be Land_use change state variable, use grid algebraically subtraction by Land_use change state Variable is converted to dichotomic variable,
Being converted into construction land by non-constructive land, variable is designated as 1；
Non-constructive land is converted into non-constructive land, and variable is designated as 0；
S22, determine construction land expansion variation: spatiality based on construction land layout, opening and adjustability Principle, build construction land influence of arrangement factor key element collection, including natural conditions, urban system, the main line of communication, transport node, Administrative obstruct, industrial factors, restriction element and program element；
S23, calculating variation power of influence:
For range attenuation class key element, use ArcGIS spatial analytical method to calculate the space length of variation, then use The space attraction force value of equation below calculating variation:
f k = e - β × dist k
Wherein, fkForce value is attracted for certain class variation space at kth space lattice；distk: kth space networks Lattice are to the space length value sum of all variations of certain class；β: take 0.0001；
For absolute restricted key element, grid algebraically subtraction is used to be converted into dichotomic variable, it is possible to build If the variation region of Land arrangement, attracting force value is 1；The variation region of construction land layout can not be carried out, inhale Gravitation value is 0；
Step S24, set up construction land expansion change with variation binary Logistic Spatial Regression Model:
L (P)=In (P/1-P)=β0+β1f1+β2f2+…+βnfn
Wherein, dependent variableOnly have two states, be expressed as construction land transition probability, have 0 and 1 two value；From change Amount i.e. variation is f1, f2, f3 ..., fn, and fn is the space captivation of all kinds of variations affecting construction land layout Value；β0The constant term obtained, β is returned for binary LogisticnThe each variation obtained is returned to construction for binary Logistic The weighted value of Land arrangement；
Step S25, construction land transition probability are predicted: binary Logistic based on construction land layout Yu variation Spatial Regression Model, structure construction land expansion probabilistic model:
P = 1 1 + e - z ,
Z=β0+β1f1+β2f2+…+βnfn；
Wherein, p is construction land transition probability；
Step S26, according to the relation between construction land layout and variation, carry out the construction land under different mode Layout simulation, calculates the construction land transition probability under different mode, and turns each spatial point construction land under each pattern Change probability averaged and standard deviation respectively, according to construction land transition probability meansigma methods and standard deviation, generate construction land Suitability mode of extension, specific as follows:
Construction land transition probability p >=construction land transition probability meansigma methods+2 times construction land transition probability standard deviation, Then it is set as construction land and controls extended area；
Construction land transition probability p >=construction land transition probability meansigma methods+1 times construction land transition probability standard deviation, Then it is set as construction land appropriateness extended area；
Construction land transition probability p >=construction land transition probability meansigma methods+0.5 times construction land transition probability standard Difference, then be set as the abundant extended area of construction land.
Described construction land suitability mode of extension and resosurces environment loading capacity partitioned area are carried out spatial overlay analysis, According to analysis result, it is achieved to preferably concretely comprising the following steps of Land_use change layout:
In the original fishing net of survey region, by construction land suitability mode of extension and resosurces environment loading capacity partition area Territory is overlapped, and the region that stack result is identical is added up and calculated its area, calculate the identical region of stack result with The area accounting of the original fishing net of survey region, according to area accounting, is carried out preferably Land_use change layout.
The method of the present invention establishes the resosurces environment loading capacity evaluation of a set of applicable land use planning in a creative way and refers to Mark, and set up the Land_use change placement scheme analogy method based on binary logistic homing method, thus provide one Kind under resosurces environment loading capacity retrains, carry out the preferred method of Land_use change placement scheme, it is achieved excellent to Land_use change layout Change design.
To retrain construction land layout, building optimum development space is target, constructs resource based on space scale Environment Carrying Capacity Assessment index system, spatial distribution the size based on its bearing capacity of measuring and calculating regional carrying capacity of resources and environments are drawn Determine high-bearing capacity district, middle bearing capacity district and low bearing capacity district three class bearing space, preferably establish bar for construction land layout Part.
On the basis of summing up the construction land influence of arrangement factor, to the pass between construction land layout and influent factor thereof System has carried out quantifying to judge, constructs construction land expansion suitability equation, with different developmental patterns, simulates different construction Land arrangement scheme.
With the evaluation result of resosurces environment loading capacity for constraint, use spatial choice, different development programs is commented Valency, based on producing, live and the principle of ecological three lives reconciling spatial scale development, finally determines the optimum side of town site layout Case.
Accompanying drawing explanation
Fig. 1 is the original fishing net figure of the covering survey region of one embodiment of the invention；
Fig. 2 is the original fishing net space overlapping figure with land use data of one embodiment of the invention；
Fig. 3 is the mode of extension of the construction land control extended area of the continuity formula layout pattern of one embodiment of the invention Figure；
Fig. 4 is the extension of the construction land appropriateness extended area of the continuity formula layout pattern of another embodiment of the present invention Ideograph；
Fig. 5 is the extension of the abundant extended area of construction land of the continuity formula layout pattern of another embodiment of the present invention Ideograph；
Fig. 6 is the extension of the construction land control extended area of the cities and towns guide type layout pattern of one embodiment of the invention Ideograph；
Fig. 7 is the construction land appropriateness extended area of the cities and towns guide type layout pattern of another embodiment of the present invention Mode of extension figure；
Fig. 8 is the abundant extended area of construction land of the cities and towns guide type layout pattern of another embodiment of the present invention Mode of extension figure；
Fig. 9 is the extension of the construction land control extended area of the traffic guide formula layout pattern of one embodiment of the invention Ideograph；
Figure 10 is the construction land appropriateness extended area of the traffic guide formula layout pattern of another embodiment of the present invention Mode of extension figure；
Figure 11 is the abundant extended area of construction land of the traffic guide formula layout pattern of another embodiment of the present invention Mode of extension figure.
Detailed description of the invention
Below the detailed description of the invention of the present invention is described, in order to those skilled in the art understand this Bright, it should be apparent that the invention is not restricted to the scope of detailed description of the invention, from the point of view of those skilled in the art, As long as various changes limit and in the spirit and scope of the present invention that determine, these changes are aobvious and easy in appended claim Seeing, all utilize the innovation and creation of present inventive concept all at the row of protection.
According to one embodiment of present invention, Land_use change layout method for optimizing based on resosurces environment loading capacity, specifically Step is as follows:
1, survey region Land_use change topology data is obtained；Land_use change topology data includes statistical data and space number According to；
An embodiment according to the application, the society in collection research region, economy, industry and resource data, the most such as Shown in table 1 below:
Table 1 Land_use change topology data
2, Land_use change topology data is carried out spatialization process, form Land_use change topology data spatialization evaluation and refer to Mark, sets up spatial database；
An embodiment according to the application, detailed process is:
2.1, with survey region as scope, the fishing net of the 500 × 500m covering whole survey region is created, as shown in Figure 1；
2.2, the fishing net generated in step 2.1 is carried out space overlapping with the Land_use change topology data of acquisition in step 1, Generate the spatial object value of each grid, obtain GDP value, the size of population, land resource distribution, distribution of water resources, traffic set Execute distribution, water conservancy projects distribution and forest cover distribution space evaluation index；
The fishing net of survey region and GDP data are carried out space overlapping, obtain each sky by an embodiment according to the application The GDP spatial data that layout is online, as shown in Figure 2；
2.3, by the evaluation space index of generation in step 2.2 according to the figure spot area of windy and sandy soil data and accordingly The area accounting of survey region similar cover type carries out being calculated the spatial density of each index, is economic density, people Mouth density, land resource enabling capabilities, Water resources supporting capacity, means of transportation distributing, water conservancy projects distribution and afforestation rate； Computing formula is as follows:
SL=ZB × ZL
Wherein, SL is the index space density of each figure speckle carrying in survey region, and ZL is that certain space-like object is at this net Spatial object value at lattice, ZB is the spatial area that certain space-like object is occupied in original grid；
2.4, it is standardized the index space density obtained in step 2.3 processing, its value is transformed between 0-1, Make each index be in the same order of magnitude, and then Comprehensive Correlation evaluation can be carried out, thus it is empty to define Land_use change topology data Between change evaluation index；
Use ArcGIS software, set up spatial database, the spatialization evaluation index data after above-mentioned process are put into State in data base.
3, according to the Land_use change topology data spatialization evaluation index formed in step 2.4, carry according to resource environment The intension criterion of power, sets up resosurces environment loading capacity assessment indicator system；Particularly as follows:
On the basis of metrics evaluation yardstick based on survey region statistical data, incorporate evaluation space yardstick, meanwhile, root According to the connotative definition to resosurces environment loading capacity, structure in terms of social economy's bearing capacity, Resources Carrying Capacity and environmental carrying capacity three Build resosurces environment loading capacity index system；
Wherein, social economy's bearing capacity, with real economy, population spatial distribution for characterizing, is resosurces environment loading capacity Core, for reflecting the diversity of socio-economic activity spatial distribution；
Resources Carrying Capacity be may utilize by region from natural resources and social resources two aspect, natural resources aspect Water, the soil resource reflection rigid support to future economy social development；Social resources aspect, by traffic and water resource infrastructure Spatial distribution characterize soft support to future economy social development；
Environmental carrying capacity is then characterized by forest cover: afforestation rate is the highest, shows Regional Eco-environmental Quality The best, preferably carry out ecological construction, and the relatively low probability that regioinvertions is construction land of afforestation rate is the biggest, The most suitable layout socio-economic activity in the case of other conditions are identical；
According to the intension criterion of above-mentioned social economy bearing capacity, Resources Carrying Capacity and environmental carrying capacity, will
Economic density and population density index are divided into social economy's load bearing principle layer；
Land resource enabling capabilities, Water resources supporting capacity, means of transportation distributing and water conservancy projects distribution index are divided For Resources Carrying Capacity rule layer；
Index of forest coverage is divided into environmental carrying capacity rule layer, thus sets up resosurces environment loading capacity evaluation index System, the most as shown in table 2:
Table 2 resosurces environment loading capacity assessment indicator system
4, use Field Using Fuzzy Comprehensive Assessment that the resosurces environment loading capacity assessment indicator system formed in step 3 is carried out subregion Process, generate resosurces environment loading capacity three regions, i.e. resource environment high-bearing capacity region, middle bearing capacity region and low bearing capacity Region.
5, carry out construction land expansion simulation based on the data in the spatial database formed in step 2, formed to build and use Ground suitability mode of extension, particularly as follows:
5.1, determine construction land expansion change variable:
For a certain concrete plot, in a certain concrete moment, can only there be two kinds of forms of expression, i.e. construction lands in this plot Or non-constructive land；The change of land use pattern i.e. shows as the transforming relationship between construction land and non-constructive land；
For the Land_use change state of t ij position, by the Land_use change state in t+1 momentSoil to t Ground utilization stateChange be Land_use change state variable, use grid algebraically subtraction Land_use change state to be become Amount is converted to dichotomic variable,
Being converted into construction land by non-constructive land, variable is designated as 1；
Non-constructive land is converted into non-constructive land, and variable is designated as 0；
5.2, construction land expansion variation is determined: during construction land influent factor collection builds, based on building If the spatiality of Land arrangement, opening and adjustability principle, build and comprise physical features, traffic, cities and towns grade, administrative resistance Hinder, industry development, the influent factor collection of many key elements such as planning, the identification for construction land influence of arrangement key element provides comparatively perfect Alternative elements, these alternative elements are chosen from the Land_use change topology data that step 1 is collected, shown in table 3 specific as follows:
Table 3 construction land influence of arrangement variable factors collection
5.3, variation power of influence size is calculated
Feature according to each variation and the mutual relation with construction land layout thereof, can be by construction land influence of arrangement Each key element that key element is concentrated is divided into range attenuation class key element and absolute restricted key element.
Range attenuation class key element mainly shows as the Space Elements such as cities and towns, traffic, administrative obstruction, and it is to construction land layout Impact be rendered as range attenuation trend, the nearest apart from these Space Elements, then condition of its exploitation is the best, and development cost is more Low, the probability producing construction land conversion is the biggest, otherwise the least；
Absolute restricted key element is also for forbidding building district's ecology key element, and such key element has pressure about to construction land layout Bundle, is the region that can not produce construction land expansion；
For range attenuation class key element, use ArcGIS spatial analytical method to calculate the space length of each variation, then adopt The space calculating variation by equation below attracts force value, and its value, between 0-1, is worth the biggest, is then converted into construction land Probability is the biggest, and vice versa, and its computing formula is:
f k = e - β × dist k
Wherein, fkForce value is attracted for certain class variation space at kth space lattice；Kth space lattice arrives The space length value sum of all variations of certain class；β: take 0.0001；
For absolute restricted key element, grid algebraically subtraction is used to be converted into dichotomic variable, it is possible to build If the variation region of Land arrangement, attracting force value is 1；The variation region of construction land layout can not be carried out, inhale Gravitation value is 0；
An embodiment according to the application, adopts for the variable in the construction land influence of arrangement variables set in step 3 Calculate its space in aforementioned manners and attract force value, as shown in table 4 below:
Table 4 construction land expansion variation captivation result table
5.4, construction land expansion change and variation binary Logistic Spatial Regression Model are set up
L (P)=In (P/1-P)=β0+β1f1+β2f2+…+βnfn
Binary Logistic Spatial Regression Model is provided with dependent variable and independent variable two class variable type, wherein, dependent variableHaving and only two states, be expressed as construction land transition probability P, represent by 0 and 1 two value, independent variable i.e. affects change Amount is f1, f2, f3 ..., fn, and fn is that the space of all kinds of variations affecting construction land layout attracts force value；β0For binary Logistic returns the constant term obtained, βnThe each variation obtained is returned to construction land layout for binary Logistic Weighted value；
5.5, construction land transition probability prediction
Spatial Regression Model based on the construction land layout set up in step 5.4 and variation relation is on each impact Variable and the judgement of weighted value thereof, structure construction land expansion probabilistic model:
P = 1 1 + e - z ,
Wherein, Z=β0+β1f1+β2f2+…+βnfn；P is construction land transition probability；
5.6, construction land placement scheme prediction
According to the relation between construction land layout and variation, carry out the construction land layout mould under different mode Intend, calculate the construction land transition probability under different mode, and to each spatial point construction land transition probability under each pattern Averaged and standard deviation respectively, according to construction land transition probability meansigma methods and standard deviation, generates construction land suitability Mode of extension, specific as follows:
Based on the driving relationship between construction land layout and variation, according to the basic mould of current construction Land arrangement Formula, on the basis of judging region practical development situation, can adjust power of influence or some variation of add drop of variation Construction land layout pattern under simulation different mode, such as continuity formula developmental pattern, satellite city developmental pattern and traffic guide formula Developmental pattern, carries out the construction land layout Path selection under different mode；
(1) continuity developmental pattern
Amount between construction land expansion and variation that continuity formula layout i.e. represents based on SPSS model regression result Change relation, determines that each point position in space produces the probability of new construction land expansion under relying on existing variation；
(2) cities and towns guide type layout
Cities and towns guide type layout is on the basis of regression result, based on construction land Present Layout Situation and town development feelings Condition, under the background of city Rapid Expansion, with satellite city, organic dispersion theory be rely on, by giving priority to center-town, emphasis Town and general town, promote the orderly balanced development in each grade cities and towns；
(3) traffic guide formula layout
Traffic guide formula layout is the Typical Representative of shaft orientation layout's pattern, based on construction land expansion and the main line of communication aobvious Write dependency, based on regression result, simulate, at enhancing road, construction land expansion effect is reduced the feelings of other key element effect Under condition, each point position in space produces the probability of construction land expansion；
Under various developmental patterns, the program prediction under different development scale can be carried out, i.e. calculating different layout mould Under formula on the basis of each point position in space transition probability, calculate meansigma methods and the standard of construction land transition probability under each pattern Difference, and determine that construction land controls extended area, construction land appropriateness extended area and construction land based on standard chosen below Fully extended area, and by result space, the concrete foundation of its scale selection is as follows:
Construction land transition probability p >=construction land transition probability meansigma methods+2 times construction land transition probability standard deviation, Then it is set as construction land and controls extended area；
Construction land transition probability p >=construction land transition probability meansigma methods+1 times construction land transition probability standard deviation, Then it is set as construction land appropriateness extended area；
Construction land transition probability p >=construction land transition probability meansigma methods+0.5 times construction land transition probability standard Difference, then be set as the abundant extended area of construction land.
An embodiment according to the application, obtains 9 kinds of layout result under three developmental patterns:
1. continuity formula layout result
As shown in Figures 3 to 5, under continuity formula layout pattern, construction land will present climing based on inner city Prolonging the external extension pattern of formula, its construction land distribution trend shows as that many south, north is few, the situation that Xi Duodong is few, regional coordination Loudness is poor.
2. cities and towns guide type layout result
As shown in Figure 6 to 8, under the guide type layout pattern of cities and towns, the state of development of stand, inner city large flat bread formula obtains Effective suppression, and the small and medium-sized cities such as center-town and emphasis obtain developing the most rapidly, cities and towns hierarchical system structure obtains Certain improvement.
3. traffic guide formula layout result
As shown in Figures 9 to 11, under traffic guide layout pattern, construction land layout presents in inner city being The heart, along national highway and provincial highway in axially spreading the trend of development, zonal distribution is obvious, away from the development in main line of communication region then by sternly Heavily suppress, present polarization trend.
6, construction land suitability mode of extension and resosurces environment loading capacity partitioned area are carried out spatial overlay analysis, root According to analysis result, it is achieved preferred to Land_use change layout, specific as follows:
In the original fishing net of survey region, by construction land suitability mode of extension and resosurces environment loading capacity partition area Territory is overlapped, and the region that stack result is identical is added up and calculated its area, calculate the identical region of stack result with The area accounting of the original fishing net of survey region, according to area accounting, is carried out preferably Land_use change layout；
An embodiment according to the application, according to " high-bearing capacity region emphasis layout, middle bearing capacity region appropriateness layout With a small amount of layout in low bearing capacity region " principle, construction land suitability mode of extension and resosurces environment loading capacity subregion are tied Fruit has carried out overlay analysis, it is thus achieved that in the different lower three kinds of placement schemes of layout pattern, construction land accounts for the carrying of different resource environment The area accounting of power functional areas, is carried out preferably placement scheme, and concrete preferred version is as shown in the table:
By the above-mentioned feature to three kinds of lower nine kinds of placement schemes of different mode and building in the range of different bearer power thereof If the comprehensive of land used accounting is analyzed and compare, cities and towns guide type minimum scheme and continuity formula minimum scheme are made the most at last For two big alternatives of construction land layout optimization, both schemes all have followed the basic principle of construction land layout, with In time, also reacted under the constraint of resosurces environment loading capacity, and construction land layout should pursue intensivization development, can appropriateness expansion build If land area.

Claims (5)

1. Land_use change layout method for optimizing based on resosurces environment loading capacity, is characterized in that, including:
Obtain survey region Land_use change topology data；Described Land_use change topology data includes statistical data and spatial data；
Described Land_use change topology data is carried out spatialization process, forms Land_use change topology data spatialization evaluation index, Set up spatial database；
According to Land_use change topology data spatialization evaluation index, according to the intension criterion of resosurces environment loading capacity, set up resource Environment Carrying Capacity Assessment index system；
Use Field Using Fuzzy Comprehensive Assessment that described resosurces environment loading capacity assessment indicator system carries out multidomain treat-ment, generate resource ring Bearing capacity three regions, border, i.e. resource environment high-bearing capacity region, middle bearing capacity region and low bearing capacity region；
Carry out construction land expansion simulation based on described spatial database, form construction land suitability mode of extension；
Described construction land suitability mode of extension and resosurces environment loading capacity partitioned area are carried out spatial overlay analysis, according to Analysis result, it is achieved preferred to Land_use change layout.
Land_use change layout method for optimizing based on resosurces environment loading capacity the most according to claim 1, is characterized in that, right Described Land_use change topology data carries out spatialization process, forms Land_use change topology data spatialization evaluation index, sets up sky The concrete steps of spatial database include:
S11, the original fishing net of generation covering survey region；
S12, the original fishing net and the Land_use change topology data that generate are carried out space overlapping, generate each net in original fishing net The spatial object value of lattice, obtains GDP value, the size of population, land resource distribution, distribution of water resources, means of transportation distributing, water conservancy Distribution of facilities and forest cover distribution space evaluation index；
S13, according in step S12 generate evaluation space index, according to the figure spot area of windy and sandy soil data and accordingly The area accounting of survey region similar cover type, is calculated the spatial density of each index, is economic density, population close Degree, the distribution of land resource enabling capabilities, Water resources supporting capacity, means of transportation distributing, water conservancy projects and afforestation rate；Calculate Formula is as follows:
SL=ZB × ZL
Wherein, SL is the index space density of each figure speckle carrying, and ZL is certain space-like object spatial object at this grid Value, ZB is the spatial area that certain space-like object is occupied in original grid；
S14, the index space density obtaining step S13 are standardized processing, and form Land_use change topology data spatialization Evaluation index.
Land_use change layout method for optimizing based on resosurces environment loading capacity the most according to claim 2, is characterized in that, root According to Land_use change topology data spatialization evaluation index, according to the intension criterion of resosurces environment loading capacity, set up resource environment and hold Concretely comprising the following steps of load power assessment indicator system:
Intension criterion according to resosurces environment loading capacity, is classified as social economy's bearing capacity, Resources Carrying Capacity and environment carrying Power；
Land_use change topology data spatialization evaluation index is held according to social economy's bearing capacity, Resources Carrying Capacity and environment respectively The intension criterion of load power, sets up resosurces environment loading capacity assessment indicator system:
Economic density and population density index are divided into social economy's load bearing principle layer；
Land resource enabling capabilities, Water resources supporting capacity, means of transportation distributing and water conservancy projects distribution index are divided into money Source load bearing principle layer；
Index of forest coverage is divided into environmental carrying capacity rule layer.
4. according to the Land_use change layout method for optimizing based on resosurces environment loading capacity described in claim 1 or 3, its feature It is to carry out construction land expansion simulation based on described spatial database, forms the concrete step of construction land suitability mode of extension Suddenly it is:
S21, determine construction land expansion change variable: being changed between construction land and non-constructive land of land use pattern Conversion,For the Land_use change state of t ij position, by the Land_use change state in t+1 momentSoil to t Utilization stateChange be Land_use change state variable, use grid algebraically subtraction Land_use change state variable to be turned It is changed to dichotomic variable,
Being converted into construction land by non-constructive land, variable is designated as 1；
Non-constructive land is converted into non-constructive land, and variable is designated as 0；
S22, determine construction land expansion variation: spatiality based on construction land layout, opening and adjustability principle, Build construction land influence of arrangement factor key element collection, including natural conditions, urban system, the main line of communication, transport node, administrative resistance Every, industrial factors, restriction element and program element；
S23, calculating variation power of influence:
For range attenuation class key element, use ArcGIS spatial analytical method to calculate the space length of variation, then use as follows The space attraction force value of formula calculating variation:
f k = e - β × dist k
Wherein, fkForce value is attracted for certain class variation space at kth space lattice；distk: kth space lattice arrives The space length value sum of all variations of certain class；β: control variable, takes 0.0001；
For absolute restricted key element, grid algebraically subtraction is used to be converted into dichotomic variable, it is possible to carry out building and use The variation region of ground layout, attracting force value is 1；The variation region of construction land layout, captivation can not be carried out Value is 0；
Step S24, set up construction land expansion change with variation binary Logistic Spatial Regression Model:
L (P)=In (P/1-P)=β0+β1f1+β2f2+…+βnfn
Wherein, dependent variableOnly have two states, be expressed as construction land transition probability P, have 0 and 1 two value；Independent variable is i.e. Variation is f1, f2, f3 ..., fn, and fn is that the space of all kinds of variations affecting construction land layout attracts force value；β0 The constant term obtained, β is returned for binary LogisticnThe each variation obtained is returned to construction land for binary Logistic The weighted value of layout；
Step S25, construction land transition probability are predicted: binary Logistic space based on construction land layout Yu variation Regression model, structure construction land expansion probabilistic model:
P = 1 1 + e - z ,
Wherein, Z=β0+β1f1+β2f2+…+βnfn；
Wherein, p is construction land transition probability；
Step S26, according to the relation between construction land layout and variation, carry out the construction land layout under different mode Simulation, calculates the construction land transition probability under different mode, and converts each spatial point construction land under each pattern general Rate averaged and standard deviation respectively, according to construction land transition probability meansigma methods and standard deviation, generate construction land suitable Property mode of extension, specific as follows:
Construction land transition probability p >=construction land transition probability meansigma methods+2 times construction land transition probability standard deviation, then by It is set to construction land and controls extended area；
Construction land transition probability p >=construction land transition probability meansigma methods+1 times construction land transition probability standard deviation, then by It is set to construction land appropriateness extended area；
Construction land transition probability p >=construction land transition probability meansigma methods+0.5 times construction land transition probability standard deviation, then It is set as the abundant extended area of construction land.
Land_use change layout method for optimizing based on resosurces environment loading capacity the most according to claim 4, is characterized in that, will Described construction land suitability mode of extension and resosurces environment loading capacity partitioned area carry out spatial overlay analysis, according to analyzing knot Really, it is achieved to preferably concretely comprising the following steps of Land_use change layout:
In the original fishing net of survey region, construction land suitability mode of extension is entered with resosurces environment loading capacity partitioned area Row superposition, adds up to the region that stack result is identical and calculates its area, calculates the identical region of stack result and research The area accounting of the original fishing net in region, according to area accounting, is carried out preferably Land_use change layout.
CN201610634963.6A 2016-08-05 2016-08-05 Land_use change layout method for optimizing based on resosurces environment loading capacity Pending CN106251020A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610634963.6A CN106251020A (en) 2016-08-05 2016-08-05 Land_use change layout method for optimizing based on resosurces environment loading capacity

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610634963.6A CN106251020A (en) 2016-08-05 2016-08-05 Land_use change layout method for optimizing based on resosurces environment loading capacity

Publications (1)

Publication Number Publication Date
CN106251020A true CN106251020A (en) 2016-12-21

Family

ID=58077772

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610634963.6A Pending CN106251020A (en) 2016-08-05 2016-08-05 Land_use change layout method for optimizing based on resosurces environment loading capacity

Country Status (1)

Country Link
CN (1) CN106251020A (en)

Cited By (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107220768A (en) * 2017-05-26 2017-09-29 北京师范大学 A kind of construction land Space Expanding method and system
CN107423566A (en) * 2017-07-27 2017-12-01 北京师范大学 A kind of land-use optimization method for taking into account ecological water saving and contamination control
CN107730099A (en) * 2017-09-30 2018-02-23 四川师范大学 A kind of space planning method for establishing model
CN107886216A (en) * 2017-10-12 2018-04-06 天津大学 A kind of land carrying capacity analysis method based on Evaluation of Land Use Suitability
CN108537441A (en) * 2018-04-09 2018-09-14 中国地质大学（武汉） Land Use Transition economic society Effect Evaluation and monitoring information integrated system
CN108805471A (en) * 2018-08-21 2018-11-13 北京师范大学 Evaluation method for water resources carrying capacity based on the analysis of hybrid system interactively
CN109858180A (en) * 2019-02-26 2019-06-07 中国农业大学 A kind of area crops water consumption spatial framework optimum design method
CN110674975A (en) * 2019-08-30 2020-01-10 华北水利水电大学 Spatial layout method and device for reducing carbon emission in urban planning
CN111160794A (en) * 2019-12-31 2020-05-15 中国地质科学院探矿工艺研究所 Geological environment bearing capacity analysis method based on AHP-GIS coupling analysis
CN111310340A (en) * 2020-02-19 2020-06-19 中南大学 Urban area interaction abnormal relation identification method and equipment based on human movement
CN112487724A (en) * 2020-12-11 2021-03-12 重庆邮电大学 Urban dynamic expansion simulation method based on partition and improved CNN-CA model
WO2021103322A1 (en) * 2019-11-29 2021-06-03 青岛理工大学 Auxiliary design system for new district planning of famous historical city
CN112966908A (en) * 2021-02-23 2021-06-15 深圳市城市规划设计研究院有限公司 Function hybrid evaluation method and device for construction land and storage medium
CN113095574A (en) * 2021-04-15 2021-07-09 南京经纬地诚土地规划咨询有限公司 Territorial space universe renovation method
US20210294935A1 (en) * 2017-11-10 2021-09-23 Autodesk, Inc. Techniques for automatically generating designs having characteristic topologies for urban design projects
CN113837913A (en) * 2021-09-09 2021-12-24 胡月明 Method and device for determining resource bearing capacity key threshold of cultivated land of villages and small towns
CN115375203A (en) * 2022-10-25 2022-11-22 山东省地质矿产勘查开发局八〇一水文地质工程地质大队(山东省地矿工程勘察院) Comprehensive analysis system for multi-element water resource optimization configuration
CN116467788A (en) * 2023-05-10 2023-07-21 悉地（苏州）勘察设计顾问有限公司 Barrier-free environment construction management method and system
CN117113514A (en) * 2023-10-17 2023-11-24 中天昊建设管理集团股份有限公司 Building design layout optimizing system based on BIM technology

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103824226A (en) * 2014-03-05 2014-05-28 广州大学 Ploughing function evaluating method based on GIS (geographic information system)

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103824226A (en) * 2014-03-05 2014-05-28 广州大学 Ploughing function evaluating method based on GIS (geographic information system)

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
黄中杰等: "泸州市区域发展强度空间格局", 《经济研究导刊》 *

Cited By (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107220768A (en) * 2017-05-26 2017-09-29 北京师范大学 A kind of construction land Space Expanding method and system
CN107423566A (en) * 2017-07-27 2017-12-01 北京师范大学 A kind of land-use optimization method for taking into account ecological water saving and contamination control
CN107730099A (en) * 2017-09-30 2018-02-23 四川师范大学 A kind of space planning method for establishing model
CN107886216A (en) * 2017-10-12 2018-04-06 天津大学 A kind of land carrying capacity analysis method based on Evaluation of Land Use Suitability
US20210294935A1 (en) * 2017-11-10 2021-09-23 Autodesk, Inc. Techniques for automatically generating designs having characteristic topologies for urban design projects
US11922099B2 (en) * 2017-11-10 2024-03-05 Autodesk, Inc. Techniques for automatically generating designs having characteristic topologies for urban design projects
CN108537441A (en) * 2018-04-09 2018-09-14 中国地质大学（武汉） Land Use Transition economic society Effect Evaluation and monitoring information integrated system
CN108805471A (en) * 2018-08-21 2018-11-13 北京师范大学 Evaluation method for water resources carrying capacity based on the analysis of hybrid system interactively
CN109858180A (en) * 2019-02-26 2019-06-07 中国农业大学 A kind of area crops water consumption spatial framework optimum design method
CN110674975B (en) * 2019-08-30 2023-04-07 华北水利水电大学 Spatial layout method and device for reducing carbon emission in urban planning
CN110674975A (en) * 2019-08-30 2020-01-10 华北水利水电大学 Spatial layout method and device for reducing carbon emission in urban planning
WO2021103322A1 (en) * 2019-11-29 2021-06-03 青岛理工大学 Auxiliary design system for new district planning of famous historical city
CN111160794A (en) * 2019-12-31 2020-05-15 中国地质科学院探矿工艺研究所 Geological environment bearing capacity analysis method based on AHP-GIS coupling analysis
CN111310340A (en) * 2020-02-19 2020-06-19 中南大学 Urban area interaction abnormal relation identification method and equipment based on human movement
CN111310340B (en) * 2020-02-19 2022-08-16 中南大学 Urban area interaction abnormal relation identification method and equipment based on human movement
CN112487724A (en) * 2020-12-11 2021-03-12 重庆邮电大学 Urban dynamic expansion simulation method based on partition and improved CNN-CA model
CN112966908A (en) * 2021-02-23 2021-06-15 深圳市城市规划设计研究院有限公司 Function hybrid evaluation method and device for construction land and storage medium
CN113095574A (en) * 2021-04-15 2021-07-09 南京经纬地诚土地规划咨询有限公司 Territorial space universe renovation method
CN113837913B (en) * 2021-09-09 2023-11-17 胡月明 Method and device for determining key threshold value of rural cultivated land resource bearing capacity
CN113837913A (en) * 2021-09-09 2021-12-24 胡月明 Method and device for determining resource bearing capacity key threshold of cultivated land of villages and small towns
CN115375203A (en) * 2022-10-25 2022-11-22 山东省地质矿产勘查开发局八〇一水文地质工程地质大队(山东省地矿工程勘察院) Comprehensive analysis system for multi-element water resource optimization configuration
CN116467788A (en) * 2023-05-10 2023-07-21 悉地（苏州）勘察设计顾问有限公司 Barrier-free environment construction management method and system
CN116467788B (en) * 2023-05-10 2023-11-21 悉地（苏州）勘察设计顾问有限公司 Barrier-free environment construction management method and system
CN117113514A (en) * 2023-10-17 2023-11-24 中天昊建设管理集团股份有限公司 Building design layout optimizing system based on BIM technology
CN117113514B (en) * 2023-10-17 2024-01-09 中天昊建设管理集团股份有限公司 Building design layout optimizing system based on BIM technology

Similar Documents

Publication Publication Date Title
CN106251020A (en) Land_use change layout method for optimizing based on resosurces environment loading capacity
Grădinaru et al. Green infrastructure in strategic spatial plans: Evidence from European urban regions
Jin et al. Designing a spatial pattern to rebalance the orientation of development and protection in Wuhan
Zhang et al. Planning an ecological network of Xiamen Island (China) using landscape metrics and network analysis
QuanLi et al. Agent-based modeling and simulations of land-use and land-cover change according to ant colony optimization: a case study of the Erhai Lake Basin, China
CN103136603A (en) Intelligent land utilization layout optimal allocation method
CN106845623A (en) A kind of electric power wireless private network base station planning method based on artificial fish-swarm algorithm
Tang et al. Dynamic evolution and scenario simulation of habitat quality under the impact of land-use change in the Huaihe River Economic Belt, China
Yuan et al. Regional land-use allocation using a coupled MAS and GA model: From local simulation to global optimization, a case study in Caidian District, Wuhan, China
Rondinelli et al. Integrated regional development planning: Linking urban centres and rural areas in Bolivia
Yang et al. Layout optimization and multi-scenarios for land use: An empirical study of production-living-ecological space in the Lanzhou-Xining City Cluster, China
Choudhury et al. GIS augmented computational intelligence technique for rural cluster electrification through prioritized site selection of micro-hydro power generation system
CN101419623A (en) Geographical simulation optimizing system
Dai et al. Agent-based model of land system: Theory, application and modelling framework
Bazan-Krzywoszańska et al. Economic conditions for the development of energy efficient civil engineering using RES in the policy of cohesion of the European Union (2014–2020). Case study: The town of Zielona Gora
Li et al. Multi-objective location-scale optimization model and solution methods for large-scale emergency rescue resources
CN103577899B (en) A kind of service combining method combined with QoS based on creditability forceast
Yang et al. Construction of ecological security pattern based on ecosystem service evaluation and minimal cumulative resistance model: a case study of Hefei City, China
Cao Development potential evaluation for land resources of forest tourism based on fuzzy ahp method
Zhai et al. Characteristic analyses, simulations and predictions of land use in poor mountainous cities: a case study in the central area of Chengde County, China
Kononova Geoecological marketing of tourist-recreational zones of cities territories
Zhou et al. An interval fuzzy land-use allocation model (IFLAM) for Beijing in association with environmental and ecological consideration under uncertainty
Zambrano-Asanza et al. Multicriteria decision analysis in geographic information systems for identifying ideal locations for new substations
Zhang et al. Eco-city planning in China: A review of policies and cases 2009–2015
CN109635912B (en) Pipeline layout method and device based on bidirectional ant colony and dynamic programming

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into 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: 20161221