CN117010751B - Urban ring land utilization toughness assessment method based on ring division - Google Patents
Urban ring land utilization toughness assessment method based on ring division Download PDFInfo
- Publication number
- CN117010751B CN117010751B CN202310932730.4A CN202310932730A CN117010751B CN 117010751 B CN117010751 B CN 117010751B CN 202310932730 A CN202310932730 A CN 202310932730A CN 117010751 B CN117010751 B CN 117010751B
- Authority
- CN
- China
- Prior art keywords
- toughness
- index
- land
- ring
- ring layer
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000009826 distribution Methods 0.000 claims abstract description 25
- 238000004458 analytical method Methods 0.000 claims abstract description 5
- CIWBSHSKHKDKBQ-JLAZNSOCSA-N Ascorbic acid Chemical compound OC[C@H](O)[C@H]1OC(=O)C(O)=C1O CIWBSHSKHKDKBQ-JLAZNSOCSA-N 0.000 claims description 30
- 238000010276 construction Methods 0.000 claims description 27
- 230000007704 transition Effects 0.000 claims description 26
- 238000004364 calculation method Methods 0.000 claims description 24
- 230000002093 peripheral effect Effects 0.000 claims description 20
- ZZZCUOFIHGPKAK-UHFFFAOYSA-N D-erythro-ascorbic acid Natural products OCC1OC(=O)C(O)=C1O ZZZCUOFIHGPKAK-UHFFFAOYSA-N 0.000 claims description 15
- 229930003268 Vitamin C Natural products 0.000 claims description 15
- 238000011161 development Methods 0.000 claims description 15
- 235000019154 vitamin C Nutrition 0.000 claims description 15
- 239000011718 vitamin C Substances 0.000 claims description 15
- 239000011159 matrix material Substances 0.000 claims description 4
- 239000002131 composite material Substances 0.000 claims description 3
- 238000010586 diagram Methods 0.000 claims description 2
- 238000011156 evaluation Methods 0.000 abstract description 9
- 238000007619 statistical method Methods 0.000 abstract description 7
- 239000010410 layer Substances 0.000 description 94
- 238000011160 research Methods 0.000 description 4
- 230000000877 morphologic effect Effects 0.000 description 3
- 230000002776 aggregation Effects 0.000 description 2
- 238000004220 aggregation Methods 0.000 description 2
- 238000012512 characterization method Methods 0.000 description 2
- 239000012792 core layer Substances 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000004445 quantitative analysis Methods 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013467 fragmentation Methods 0.000 description 1
- 238000006062 fragmentation reaction Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
Classifications
-
- 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/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- 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/24—Querying
- G06F16/245—Query processing
- G06F16/2458—Special types of queries, e.g. statistical queries, fuzzy queries or distributed queries
- G06F16/2462—Approximate or statistical queries
-
- 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
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- 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
-
- 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/26—Government or public services
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A30/00—Adapting or protecting infrastructure or their operation
- Y02A30/60—Planning or developing urban green infrastructure
Abstract
The invention belongs to the technical field related to regional space toughness analysis, and discloses a regional urban ring land utilization toughness assessment method. The method comprises the following steps: s1, dividing a city ring to be treated into three ring layers, and determining a toughness target corresponding to each ring layer and a toughness index for measuring each toughness target; s2, dividing each ring layer into a plurality of grid units, and calculating the toughness index of each grid unit; calculating the weight value of each toughness index, calculating the comprehensive toughness index corresponding to each grid unit by using each toughness index and the weight obtained by corresponding each toughness index, and visualizing the index value to obtain the spatial distribution of toughness; s3, calculating a variation coefficient, a skewness coefficient and a kurtosis coefficient of the comprehensive toughness index, and realizing statistical analysis of the toughness index. The method solves the problems that in the prior art, the land utilization toughness of the urban ring layer is not quantitatively evaluated, the evaluation result is difficult to reveal the difference of the internal toughness of the urban ring, and the space planning is difficult to be guided in a targeted manner.
Description
Technical Field
The invention belongs to the technical field related to regional space toughness analysis, and particularly relates to a regional land utilization toughness assessment method for a regional ring.
Background
With the continuous development of urbanization, the urban ring becomes an important space carrier for new urbanization, and meanwhile, the vulnerability problem faced by the urban ring with high urbanization is increasingly prominent. As an advanced form of urban regional spatial form evolution, urban rings have obvious ring layer spatial structures, the development stage of each ring layer, the spatial element aggregation characteristic and the spatial layout mode are different, the main risks are different, and the problems of knowing the risk difference inside the urban rings and improving the vulnerability of the urban rings become common attention.
The toughness emphasizing system has the key properties of diversity, redundancy, flexibility, self-organization and the like in the face of resistance, recovery and adaptability in risk disturbance. Land utilization toughness is a component of urban space toughness, focusing on the capability of toughness features such as the number of land utilization types, diversity of space structures and functions, connectivity, robustness and redundancy, etc., to maintain the basic functions of the system unchanged and to enable recovery and adaptation. The urban ring land utilization toughness focuses on how to buffer and eliminate risks through the coordinated configuration of land utilization system elements, and has important significance in promoting the urban ring safety toughness development.
At present, urban ring toughness assessment is still in a starting stage, and most of researches focus on toughness levels of subsystems such as society and economy, and most of the researches adopt statistical indexes to assess urban ring toughness, and focus on urban ring space toughness, especially land utilization toughness is insufficient. The evaluation scale is mainly based on urban ring global toughness evaluation, the toughness difference between different ring layers of the urban ring cannot be reflected, and recognition and quantitative analysis of urban ring land utilization toughness space characteristics and ring layer toughness difference are the basis for optimizing urban ring land utilization layout and realizing urban ring safety integrated development.
Therefore, how to provide a method for evaluating the toughness of urban ring land utilization of a layered ring is a problem that needs to be solved by those skilled in the art.
Disclosure of Invention
Aiming at the defects or improvement demands of the prior art, the invention provides a regional urban ring land utilization toughness assessment method, which solves the problems that the quantitative assessment of the regional urban ring land utilization toughness is insufficient, the toughness difference inside the urban ring is difficult to reveal as the assessment result, and the regional planning practice of the urban ring is difficult to guide in a targeted manner in the prior art.
In order to achieve the above object, according to the present invention, there is provided a method for evaluating toughness of urban land use of a split-level urban land, comprising the steps of:
s1, dividing a city ring to be treated into three ring layers, namely a core ring layer, a transition ring layer and a peripheral ring layer, and determining a toughness target corresponding to each ring layer and a toughness index for measuring each toughness target;
s2, dividing each ring layer into a plurality of grid units, and calculating a toughness index corresponding to each grid unit so as to obtain the toughness index corresponding to each grid unit in each ring layer; for each ring layer, calculating and obtaining a weight value of each toughness index by using the toughness index numerical distribution corresponding to each grid unit, calculating and obtaining a comprehensive toughness index corresponding to each grid unit by using each toughness index and the weight obtained by the corresponding toughness index, and visualizing the comprehensive toughness index to obtain a toughness space distribution diagram of each grid;
s3, calculating a variation coefficient, a skewness coefficient and a kurtosis coefficient corresponding to each ring layer by utilizing the comprehensive toughness index of the grid corresponding to each ring layer, and evaluating the toughness of the urban ring to be treated.
Further preferably, in step S1, the toughness targets corresponding to the core ring layer are diversity, connectivity and redundancy, the toughness targets corresponding to the transition ring layer are diversity, robustness and high efficiency, and the toughness targets corresponding to the peripheral ring layer are diversity, accessibility and stability.
Further preferably, the toughness index measuring the diversity is a land utilization type diversity index, the toughness index measuring the connectivity is a form connectivity, the toughness index measuring the redundancy is a land development strength and an infrastructure density, the toughness index measuring the robustness is a construction land crush degree, the toughness index measuring the high efficiency is an industrial average yield value, the toughness index measuring the reachability is a medical institution reachability and a road network density, and the toughness index measuring the stability is a land utilization ecological service function index.
Further preferably, in step S2, the dividing each ring layer into a plurality of grid units is performed using a fishing net tool of ArcGIS.
Further preferably, in step S2, the calculation of the toughness index of the core hoop layer is performed in the following manner:
R GM ＝1-∑f a 2 /(∑f a ) 2
D＝S/C
wherein R is GM Is a diversified index, f a Is the area of the a-th land use type, DI is the "source" plaque average distance index; d, d gh The distance from the grid g in the "source" plaque to the grid h in the "sink" plaque, p and q are the total number of grids of the "source" plaque and the "sink" plaque respectively, D is the regional land development intensity, S is the construction land area, and C is the regional total area; infrastructure density was calculated using ArcGIS nuclear density analysis tools.
Further preferably, in step S2, the calculation of the toughness index of the transition ring layer is performed in the following manner:
PD＝NP/S
A＝T/SI
wherein PD is the degree of breakage of the construction land, NP is the number of plaques of the construction land, S is the total area of the construction land, A is the industrial average yield value, T is the regular industrial total yield value, and SI is the total area of the regional industrial land; the calculation of the diversified indexes is the same as the calculation method of the diversified indexes of the core circle.
Further preferably, in step S2, the calculation of the toughness index of the outer circumferential layer is performed in the following manner:
R d ＝L/C
ESV k ＝Σ f A k ×VC k
ESV f ＝∑ f A k ×VC kf
ESV＝∑∑A k ×VC kf
wherein R is d The road network density is L is the total length of the roads of county roads and above, and C is the total area of the area; ESV (electronic service Voltage) k 、ESV f ESV is the ecological service value of the k-th type land, the ecological service value of the f-th service function and the total ecological system service value respectively; a is that k Is the k land use type area; VC (vitamin C) kf The ecological service value of the f service unit area of the k type land use type; diversity ofThe index is the same as the calculation method of the core ring layer, and the accessibility of the medical institution is obtained by calculating an OD cost matrix in the ArcGIS.
Further preferably, in step S2, the weight value is determined by an entropy method.
Further preferably, in step S2, the integrated toughness index is performed according to the following relation:
wherein F is the land utilization toughness level comprehensive index, ω of the grid j Is the weight of the j index, p ij Is a normalized value of the j-th index of grid i, i=1, 2,.. m, n is the total number of grids, m is the index number of each turn.
Further preferably, in step S3, the coefficient of variation, the skewness coefficient and the kurtosis coefficient are respectively performed according to the following relations:
c v ＝σ/μ
wherein c v Is the coefficient of variation, σ is the standard deviation, μ is the average value, SK is the skewness coefficient, μ is the average value, M 0 Is the mode, K is the kurtosis coefficient, i=1, 2, &.. i is the number of grids, n is the total number of grids, F i Is the land use toughness level composite index value for grid i.
In general, compared with the prior art, the above technical solution conceived by the present invention has the following beneficial effects:
1. the invention provides a regional urban ring land utilization toughness assessment method, which is characterized in that a regional ring land utilization toughness assessment framework is constructed, land utilization toughness index data of different ring layers are obtained through calculation based on land utilization data, public service facility POI data, statistical annual survey data and the like, and the spatial distribution characteristics of the land utilization toughness of a core ring layer, a transition ring layer and a peripheral ring layer are obtained through standardized processing and weighted superposition, so that statistical analysis is further carried out on the toughness distribution difference among the ring layers, and the method has important significance in quantifying urban ring space toughness and optimizing urban ring land utilization layout;
2. the regional space difference of toughness is further excavated by carrying out regional layer assessment on the toughness of the urban ring to conform to the structural characteristics of the development of the urban ring layer, and the regional space difference is an important foundation for implementing the layout optimization of the urban ring layer. Land utilization toughness is calculated and evaluated from grid scale, so that the space aggregation and difference characteristics of toughness can be reflected well, and the regional weak toughness area of the urban ring can be identified;
3. the method can evaluate the land utilization toughness level of the urban area relatively pertinently, and meanwhile, the land utilization data, the POI data of the public service facilities and the statistical yearbook data are all open source data, so that the acquisition difficulty is low, the evaluation method is suitable for the scale of the urban area, the operation cost is low, and the operation is simple and convenient.
Drawings
FIG. 1 is a flow chart of a method of urban ring land utilization toughness assessment for a split ring layer constructed in accordance with a preferred embodiment of the invention;
FIG. 2 is a spatial characterization of the core horizon land utilization toughness assessment index constructed in accordance with a preferred embodiment of the present invention;
FIG. 3 is a spatial characterization of the toughness assessment index for a transitional loop layer constructed in accordance with a preferred embodiment of the present invention;
FIG. 4 is a spatial signature of a peripheral containment land use assessment index constructed in accordance with a preferred embodiment of the present invention;
FIG. 5 is a graph of core course land utilization toughness assessment results constructed in accordance with a preferred embodiment of the present invention;
FIG. 6 is a plot of land utilization toughness assessment results for a transition ring constructed in accordance with a preferred embodiment of the present invention;
FIG. 7 is a plot of land utilization toughness assessment of a peripheral hoop constructed in accordance with a preferred embodiment of the present invention.
Detailed Description
The present invention will be described in further detail with reference to the drawings and examples, in order to make the objects, technical solutions and advantages of the present invention more apparent. It should be understood that the specific embodiments described herein are for purposes of illustration only and are not intended to limit the scope of the invention. In addition, the technical features of the embodiments of the present invention described below may be combined with each other as long as they do not collide with each other.
As shown in fig. 1, a method for evaluating toughness of urban ring land use of a layered ring includes:
s1, constructing a regional urban ring land utilization toughness assessment framework. The method comprises the following specific steps:
s11, urban ring layer division. The space structure research of the existing Wuhan urban ring is combined, the urban ring is divided into three parts of a core, a transition ring layer and a peripheral ring layer, and the space range of each ring layer is defined. Setting up
S12, constructing a toughness target layer. And analyzing the land utilization characteristics of each ring layer and providing a toughness target. The core ring layer toughness targets comprise diversity, connectivity and redundancy; the transition ring layer toughness targets comprise diversity robustness and high efficiency; the peripheral hoop toughness targets include diversity, accessibility, and stability. Wherein:
diversity is a systematic goal of urban ring land utilization. Diversity refers to the property of a system having or having different types of elemental composition, which is the basis for toughness. The diversity of land utilization types reflects the diversity of urban functions and can meet various demands of urban systems.
The urban ring core ring layer is a highly dense urban area, and the land layout presents the characteristics of concentrated connection of construction lands and broken plaque mosaic distribution of ecological lands. The highly concentrated economic factors and population make the space concentration of the core layer risk stronger, and the requirements on redundancy and quick response capability of the service land functions are higher; at the same time, the high density of the built space and the continuous watertight ground increase ecological risks. Connectivity and redundancy are therefore the primary toughness targets for the core hoop layer.
The transition ring layer does not form a highly connected built-up space, and mainly comprises a city border area of a core city, a new city and a built-up area of surrounding cities. The construction land and the non-construction land of the area are distributed in a staggered way, and the occupation of the industrial and mining land is relatively high; the construction land is inefficient and has significant fragmentation characteristics. Robustness and efficiency are thus the primary toughness targets for the transition ring layer.
The peripheral ring layer is mainly formed by continuously distributing non-construction lands such as cultivated lands and ecological lands, and the urban and rural construction lands are distributed in plaque and point shapes. The unbalanced distribution of land makes regional population more likely to face the problems of insufficient infrastructure, insufficient traffic accessibility, etc., and poor response capability in the face of impact. Meanwhile, the development and construction of the punctiform land affect the ecological service supply capability of the stable agricultural land to a certain extent. Thus, accessibility and stability are the primary toughness targets for the outer hoop layer.
S13, constructing a toughness index layer. And (5) screening land utilization toughness indexes of the core ring layer, the transition ring layer and the peripheral ring layer based on the toughness targets.
TABLE 1 Metropolis Ring land use toughness assessment index System with ring layers
The core circle layer diversity index comprises a land utilization type diversity index; the connectivity index includes a morphological connectivity index; redundancy indicators include land development strength and infrastructure density.
The transition circle layer diversity index comprises a land utilization type diversity index; the robustness index comprises the construction land breaking degree; the high efficiency index includes industrial average.
The peripheral circle layer diversity index comprises a land utilization type diversity index; the stability index comprises a land utilization ecological service function index; reachability metrics include medical facility reachability and road network density.
S2, dividing urban land utilization toughness evaluation basic units.
S21, dividing the range of the core ring layer, the transition ring layer and the peripheral ring layer into a plurality of grids by using a fishing net tool of the ArcGIS, and dividing the ranges into grids of 1km x 1km in the embodiment to form a basic unit for urban ring land utilization toughness evaluation.
S22, calculating the land utilization toughness indexes of the core ring layer, the transition ring layer and the peripheral ring layer.
A. And (5) calculating the toughness index of the core ring land. The calculated index value of the land utilization toughness of the core ring layer and the spatial distribution thereof are shown in fig. 2.
The core layer land utilization toughness index comprises a land utilization type diversity index, a form connectivity index, a land development strength and an infrastructure density. Wherein:
(1) Land use type diversification index. And obtaining each land use type and area thereof based on land use investigation data, and calculating and obtaining a land use type diversity index based on the land use type data. The calculation formula is as follows:
R GM ＝1-∑f a 2 /(∑f a ) 2
wherein R is GM For a diversified index, f a Is the area of the a-th land use type.
(2) The method for calculating the morphological connectivity index comprises the following steps:
extracting construction land, woodland, grassland and water based on land utilization coverage data; calculating the distance from the center of each grid in the 'source' landscape (construction land) to the nearest grid center in the 'sink' landscape (forest land, grassland and water body), adding up and dividing by the number of the 'source' landscape grids to obtain an 'source-sink' landscape average distance index, and taking the 'source-sink' landscape average distance index as a form connectivity index. The calculation formula is as follows:
where DI is the "source" plaque average distance index; d, d gh The distance from grid g in the "source" patch to grid h in the "sink" patch, and p and q are the total number of grids of the "source" patch and the "sink" patch, respectively.
(3) Land development intensity. The land development strength is expressed as a ratio of the area of the constructed land of the area to the total area of the area based on the land use coverage data. The calculation formula is as follows:
d=s/C, where D is the regional land development intensity, S is the construction land area, and C is the total regional area.
(4) Infrastructure density. Based on the distribution data of the POIs (points of interest Point of Interest) of the medical service facilities, on the basis of the 1km grid of the core circle layer, calculating the spatial distribution density of the infrastructure by using an ArcGIS nuclear density analysis tool, and obtaining the density of the infrastructure by using the grid value.
B. And (5) calculating toughness indexes of the land utilization of the transition ring layer. The calculated toughness index value of the land utilization of the transition ring layer and the spatial distribution thereof are shown in fig. 3.
The land utilization toughness indexes of the transition ring layer comprise land utilization type diversified indexes, construction land crushing degree and industrial land average yield values. Wherein:
(1) And (4) calculating the land use type diversification index as in S22.
(2) Degree of breakage of construction land. Based on land utilization coverage data, plaque density is calculated by means of Fragstats software, and the degree of land breakage for construction is reflected in a plaque density index. The formula is:
PD＝NP/S
where PD is the degree of breakage of the construction land, NP is the number of construction land patches, and S is the total area of the construction land.
(3) Industrially homogeneous production values. Based on the statistical yearbook data and the land survey data, the total industrial yield value is expressed in terms of a ratio of the total industrial yield value to the area of the industrial land. The formula is:
A＝T/SI
wherein A is the industrial average production value, T is the general industrial production value, and SI is the general area of the regional industrial land.
C. And calculating toughness indexes of the land of the peripheral ring layer. The calculated toughness index value of the land utilization of the transition ring layer and the spatial distribution thereof are shown in fig. 4.
The land utilization toughness indexes of the transition ring layer comprise land utilization type diversification indexes, medical institution accessibility, road network density and land utilization ecological service value indexes. Wherein:
(1) S22, land utilization type diversity index calculation is the same as that of S;
(2) Medical facility accessibility. And calculating the reachable range of the medical institution vehicle for 30min by means of an OD cost matrix in the ArcGIS based on the POI distribution data of the medical institution and the road network of the medical institution.
(3) Road network density. Based on the road vector data, the ratio of the road length to the total area of the area above county roads is calculated to represent the road network density. The calculation formula is as follows:
R d ＝L/C
wherein R is d The road network density is L is the total length of the roads of county and above, and C is the total area of the area.
(4) The land use ecological service value index is calculated based on the land use coverage data and the ecological service value data of each land use type unit area, and the ecological service values of the land use types are added to obtain the land use ecological service value index. The calculation formula is as follows:
ESV k ＝∑ f A k ×VC k
ESV k ＝∑ f A k ×VC kf
ESV＝∑∑A k ×VC kf
in ESV of k 、ESV f ESV is respectively the ecological service value of the kth land, the ecological service value of the f service function and the total ecological system service value; a is that k Area for the kth land use type; VC (vitamin C) kf And f, serving the ecological service value of the unit area for the k type land utilization type. Wherein, equivalent ecological system service value, eachThe ecological service value of the land utilization type in unit area is obtained by referring to different provinces of data.
S23, evaluating the land utilization toughness of the urban ring.
A. The entropy method determines the index weight. The method comprises the following specific steps:
(1) And (5) normalizing the index value. And uniformly collecting sample points by using a fishing net tool of ArcGIS software, and standardizing the collected index point set data based on a range method. The formula for normalizing the forward index is as follows:
the negative index normalization formula is:
wherein F is ij Is the value of the jth index of the ith grid, F jmax And F jmin The maximum value and the minimum value of the matrix array where the j index is located are respectively.
(2) And (5) calculating the index weight. And (3) inputting the standardized index point set data into a Python entropy method operation program to obtain the weights of all indexes. The entropy method comprises the following calculation steps:
calculating the specific gravity of the jth index of the ith grid:
calculating index information entropy:
calculating information redundancy:
d j ＝1-e j
and (5) calculating index weight:
wherein n is the grid number of the research area, and m is the evaluation index number of each circle of layer.
B. And (5) analyzing the comprehensive toughness indexes of each ring layer. Based on the index weight, calculating a comprehensive toughness index, and dividing toughness grades by adopting a natural breaking point method to obtain comprehensive toughness index values and spatial distribution of the comprehensive toughness index values of a core ring layer, a transition ring layer and a peripheral ring layer, wherein the comprehensive toughness index values are shown in fig. 5, 6 and 7. The calculation formula of the comprehensive toughness index is as follows:
wherein F is the land utilization toughness level comprehensive index of the grid, and the higher the toughness level is, the higher the numerical value is, omega j Weight of j index, p ij The standard value of the j index of the grid i is obtained, n is the number of grids, and m is the index number of each circle of layers.
S3, toughness between the ring layers is compared and analyzed. And carrying out statistical analysis on three groups of toughness evaluation results of the core ring layer, the transition ring layer and the peripheral ring layer, selecting a variation coefficient, a skewness coefficient and a kurtosis coefficient, and carrying out quantitative analysis on numerical distribution, concentrated trend and peak characteristics of the urban ring layer land utilization toughness level. The variation coefficient reflects the discrete degree of the toughness numerical distribution, the skewness coefficient describes the morphological characteristics of the normal distribution of the toughness numerical, and the steepness degree of the probability distribution of the kurtosis coefficient toughness numerical. Wherein, the calculation formula of the variation coefficient is:
c v ＝σ/μ
wherein, c v The coefficient of variation, σ is the standard deviation, and μ is the average value.
The calculation formula of the skewness coefficient is as follows:
wherein SK isSkewness coefficient, mu is average value, M 0 Mode, σ is standard deviation. SK (SK)>0 is positive (right) bias, SK<0 is negative (left) bias; the greater the absolute value, the greater the degree of deflection.
The calculation formula of kurtosis coefficient is:
where K is a kurtosis coefficient, σ is a standard deviation, μ is an average value, i=1, 2 i Is the land use toughness level composite index value for grid i. The range of kurtosis is [1, ++ ], the data kurtosis value completely obeying normal distribution is 3, the peak value>3 indicates a larger outlier. The following is the calculated results of statistical analysis of the toughness values of each layer using the Wuhan urban ring as an implementation case, as shown in Table 2 below:
table 2 statistical analysis of the toughness values of each turn
The evaluation result of the embodiment reflects the spatial distribution characteristics of land utilization toughness of the Wuhan urban ring, and identifies low-value areas of toughness of the core ring layer, the transition ring layer and the peripheral ring layer; meanwhile, the statistical analysis of the toughness values among the ring layers reflects the discrete and concentrated differences of the toughness values of each ring layer, wherein the numerical distribution of the toughness levels of the peripheral ring layers is most dispersed; the toughness levels of the three ring layers show right deviation characteristics, and the concentration trend of low toughness values is more remarkable; the difference of the toughness level inside each ring layer is larger, and the difference of the transition ring layers is the most obvious.
The invention provides a regional urban ring land utilization toughness assessment method, which is characterized in that land utilization toughness index data of different ring layers are obtained through calculation by constructing a regional ring land utilization toughness assessment framework, and land utilization toughness spatial distribution characteristics of a core ring layer, a transition ring layer and a peripheral ring layer are obtained through standardized treatment and weighted superposition, so that statistical analysis is further carried out on toughness distribution differences among ring layers, and the method has important significance in quantifying urban ring space toughness and optimizing urban ring land utilization layout.
It will be readily appreciated by those skilled in the art that the foregoing description is merely a preferred embodiment of the invention and is not intended to limit the invention, but any modifications, equivalents, improvements or alternatives falling within the spirit and principles of the invention are intended to be included within the scope of the invention.
Claims (5)
1. A method for evaluating the toughness of urban ring land utilization of a layered ring is characterized by comprising the following steps:
s1, dividing a city ring to be treated into three ring layers, namely a core ring layer, a transition ring layer and a peripheral ring layer, and determining a toughness target corresponding to each ring layer and a toughness index for measuring each toughness target;
s2, dividing each ring layer into a plurality of grid units, and calculating a toughness index corresponding to each grid unit so as to obtain the toughness index corresponding to each grid unit in each ring layer; for each ring layer, calculating and obtaining a weight value of each toughness index by using the toughness index numerical distribution corresponding to each grid unit, calculating and obtaining a comprehensive toughness index corresponding to each grid unit by using each toughness index and the weight obtained by the corresponding toughness index, and visualizing the comprehensive toughness index to obtain a toughness space distribution diagram of each grid;
s3, calculating a variation coefficient, a skewness coefficient and a kurtosis coefficient corresponding to each ring layer by using the comprehensive toughness indexes of the grids corresponding to each ring layer, so as to evaluate the toughness of the urban ring to be treated;
in step S1, the toughness targets corresponding to the core ring layer are diversity, connectivity and redundancy, the toughness targets corresponding to the transition ring layer are diversity, robustness and high efficiency, and the toughness targets corresponding to the peripheral ring layer are diversity, accessibility and stability;
the toughness indexes for measuring the diversity are land utilization type diversity indexes, the toughness indexes for measuring the connectivity are form connectivity, the toughness indexes for measuring the redundancy are land development strength and infrastructure density, the toughness indexes for measuring the robustness are construction land crushing degree, the toughness indexes for measuring the high efficiency are industrial average yield values, the toughness indexes for measuring the accessibility are medical institution accessibility and road network density, and the toughness indexes for measuring the stability are land utilization ecological service function indexes;
in step S2, the toughness index of the core ring layer is calculated as follows:
D＝S/C
wherein R is GM Is a diversified index, f a Is the area of the a-th land use type, DI is the "source" plaque average distance index; d, d gh The distance from the grid g in the "source" plaque to the grid h in the "sink" plaque, p and q are the total number of grids of the "source" plaque and the "sink" plaque respectively, D is the regional land development intensity, S is the construction land area, and C is the regional total area; calculating the density of the infrastructure by using an ArcGIS nuclear density analysis tool;
in step S2, the calculation of the toughness index of the transition ring layer is performed in the following manner:
PD＝NP/S
A＝T/SI
wherein PD is the degree of breakage of the construction land, NP is the number of plaques of the construction land, S is the total area of the construction land, A is the industrial average yield value, T is the regular industrial total yield value, and SI is the total area of the regional industrial land; the calculation of the diversified indexes is the same as the calculation method of the diversified indexes of the core circle;
in step S2, the calculation of the toughness index of the outer circumferential layer is performed in the following manner:
R d ＝L/C
ESV k ＝∑ f A k ×VC k
ESV f ＝∑ f A k ×VC kf
ESV＝∑∑A k ×VC kf
wherein R is d The road network density is L is the total length of the roads of county roads and above, and C is the total area of the area; ESV (electronic service Voltage) k 、ESV f ESV is the ecological service value of the k-th type land, the ecological service value of the f-th service function and the total ecological system service value respectively; a is that k Is the k land use type area; VC (vitamin C) kf The ecological service value of the f service unit area of the k type land use type; the diversity index is the same as the calculation method of the core ring layer, and the accessibility of the medical institution is obtained by calculating an OD cost matrix in the ArcGIS.
2. The method for evaluating toughness of urban area land use of a layered urban area according to claim 1, wherein in step S2, each layered area is divided into a plurality of grid units, and the division is performed by using a fishing net tool of ArcGIS.
3. A method for evaluating toughness of a regional land use in a metro area according to claim 1 or 2, wherein in step S2, the weight value is determined by an entropy method.
4. A method for evaluating toughness of a layered urban land use according to claim 1 or 2, wherein in step S2, the integrated toughness index is performed according to the following relation:
wherein F is the land utilization toughness level comprehensive index, ω of the grid j Is the weight of the j index, p ij Is a normalized value of the j index of grid i,i=1, 2, & gt..n, j=1, 2, & gt..m, i is the number of grids, n is the total number of grids, and m is the index number of each turn layer.
5. The method for evaluating toughness of urban land use of a layered urban area according to claim 1 or 2, wherein in step S3, the coefficient of variation, the coefficient of skewness and the coefficient of kurtosis are respectively performed according to the following relations:
c v ＝σ/μ
wherein c v Is the coefficient of variation, σ is the standard deviation, μ is the average value, SK is the skewness coefficient, μ is the average value, M 0 Is the mode, K is the kurtosis coefficient, i=1, 2, &.. i is the number of grids, n is the total number of grids, F i Is the land use toughness level composite index value for grid i.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310932730.4A CN117010751B (en) | 2023-07-27 | 2023-07-27 | Urban ring land utilization toughness assessment method based on ring division |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310932730.4A CN117010751B (en) | 2023-07-27 | 2023-07-27 | Urban ring land utilization toughness assessment method based on ring division |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117010751A CN117010751A (en) | 2023-11-07 |
| CN117010751B true CN117010751B (en) | 2024-03-22 |
Family
ID=88572179
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310932730.4A Active CN117010751B (en) | 2023-07-27 | 2023-07-27 | Urban ring land utilization toughness assessment method based on ring division |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117010751B (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107145724A (en) * | 2017-04-25 | 2017-09-08 | 武汉大学 | A kind of spacial analytical method of the green trace integration service ability assay in city |
| CN113487181A (en) * | 2021-07-05 | 2021-10-08 | 上海市园林科学规划研究院 | Method for evaluating ecological safety pattern of urbanized area |
| CN115689293A (en) * | 2022-11-15 | 2023-02-03 | 中国科学院地理科学与资源研究所 | Urban waterlogging toughness evaluation method based on pressure-state-response framework |
| CN115841042A (en) * | 2022-12-21 | 2023-03-24 | 福州大学 | Urban ecological toughness three-dimensional evaluation method constructed based on toughness theory |
| CN115860213A (en) * | 2022-11-30 | 2023-03-28 | 华南理工大学 | Multi-scenario city emergency sampling place site selection optimization method, system, equipment and medium |
| CN116128350A (en) * | 2023-01-03 | 2023-05-16 | 北京交通大学 | Evaluation method and evaluation device for biological diversity value of railway corridor |
-
2023
- 2023-07-27 CN CN202310932730.4A patent/CN117010751B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107145724A (en) * | 2017-04-25 | 2017-09-08 | 武汉大学 | A kind of spacial analytical method of the green trace integration service ability assay in city |
| CN113487181A (en) * | 2021-07-05 | 2021-10-08 | 上海市园林科学规划研究院 | Method for evaluating ecological safety pattern of urbanized area |
| CN115689293A (en) * | 2022-11-15 | 2023-02-03 | 中国科学院地理科学与资源研究所 | Urban waterlogging toughness evaluation method based on pressure-state-response framework |
| CN115860213A (en) * | 2022-11-30 | 2023-03-28 | 华南理工大学 | Multi-scenario city emergency sampling place site selection optimization method, system, equipment and medium |
| CN115841042A (en) * | 2022-12-21 | 2023-03-24 | 福州大学 | Urban ecological toughness three-dimensional evaluation method constructed based on toughness theory |
| CN116128350A (en) * | 2023-01-03 | 2023-05-16 | 北京交通大学 | Evaluation method and evaluation device for biological diversity value of railway corridor |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117010751A (en) | 2023-11-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Foroozesh et al. | Assessment of sustainable urban development based on a hybrid decision-making approach: Group fuzzy BWM, AHP, and TOPSIS–GIS | |
| CN110334864B (en) | GIS-based multi-specification-in-one urban and rural space partition method | |
| Gigović et al. | Application of the GIS-DANP-MABAC multi-criteria model for selecting the location of wind farms: A case study of Vojvodina, Serbia | |
| Zhu et al. | Flood disaster risk assessment based on random forest algorithm | |
| Liu et al. | Toward classification-based sustainable revitalization: Assessing the vitality of traditional villages | |
| CN110096563B (en) | Ecological network construction method based on ecological service logic coding | |
| CN110781267A (en) | Multi-scale space analysis and evaluation method and system based on geographical national conditions | |
| Lin et al. | Fine-scale mapping of urban ecosystem service demand in a metropolitan context: A population-income-environmental perspective | |
| CN110298520A (en) | A kind of farmland construction difficulty determines method | |
| Tian et al. | Suburban identification based on multi-source data and landscape analysis of its construction land: A case study of Jiangsu Province, China | |
| Jia et al. | Study on ecological evaluation of urban land based on GIS and RS technology | |
| CN115641718A (en) | Short-term traffic flow prediction method based on bayonet flow similarity and semantic association | |
| Lyu et al. | Linking regional sustainable development goals with ecosystem services to identify ecological security patterns | |
| Sun et al. | Response of ecosystem service value to land use/cover change in the northern slope economic belt of the Tianshan Mountains, Xinjiang, China | |
| Wang et al. | Multidimensional spatiotemporal evolution detection on China’s rural poverty alleviation | |
| Zhang et al. | Prioritizing sponge city sites in rapidly urbanizing watersheds using multi-criteria decision model | |
| CN113627654A (en) | Urban ecological corridor construction method and device based on suitability and connectivity | |
| Shen et al. | Improving landscape ecological network connectivity in urbanizing areas from dual dimensions of structure and function | |
| Ren et al. | Land use evaluation considering soil properties and agricultural infrastructure in black soil region | |
| CN117010751B (en) | Urban ring land utilization toughness assessment method based on ring division | |
| CN113627818A (en) | Park green space construction project comprehensive benefit evaluation method based on urban relocation | |
| CN113327021A (en) | TOPSIS entropy weight method-based comprehensive land evaluation method and device | |
| CN113112152A (en) | State-of-the-earth space planning current situation assessment method | |
| CN116911507A (en) | Urban block vitality evaluation method and storage medium based on space projection pursuit | |
| CN115146990B (en) | Urban activity quantitative evaluation method integrating multisource geographic big data |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |