Industrial partition delimiting method and system

Abstract

The invention provides an industry partition delimiting method and system. The method comprises the following steps: acquiring spatial information data of an evaluation area; calculating and obtaining an ecological system value index of the evaluation region based on the spatial information data; calculating and obtaining a space development value index of the evaluation area based on the space information data; and according to the ecological system value index and the space development value index, performing partition delimitation of an evaluation region. According to the invention, quantitative evaluation can be carried out on the evaluation region, and the accuracy and reliability of partition delimiting are improved.

Description

Technical Field

[0001] The present invention relates to the technical field of regional assessment, and in particular to an industrial zoning demarcation method and system. Background Art

[0002] The implementation of ecological restoration projects will have an impact on the ecological pattern and affect the future industrial development of the entire region. How to optimize the future ecological space and industrial pattern of the region is extremely important for the ecological health and sustainable development of the region.

[0003] Current research on spatial development and demarcation basically only starts from a single perspective of industrial development or ecological protection. In addition, current research methods on spatial distribution generally analyze the distribution patterns of the importance of ecosystem services, without quantifying the evaluation indicators, and lack objectivity. Summary of the invention

[0004] The present invention provides an industrial zoning demarcation method and system to solve the defects of the prior art.

[0005] The present invention provides an industrial zoning demarcation method, comprising: S1: Collect spatial information data of the assessment area; S2: Calculate and obtain the ecosystem value index of the assessment area based on the spatial information data; S3: Calculate and obtain the spatial development value index of the assessment area based on the spatial information data; S4: Zoning the assessment area based on the ecosystem value index and the spatial development value index.

[0006] According to an industrial zoning demarcation method provided by the present invention, the spatial information data in step S1 specifically includes: first spatial information data, wherein the first spatial information data is remote sensing data; The second spatial information data includes topographic surveying and mapping data, hydrological monitoring data, land use data and vegetation index data.

[0007] According to an industrial zoning demarcation method provided by the present invention, step S1 further includes: Preprocessing the first spatial information data in the spatial information data specifically includes formatting the remote sensing data into a raster format and unifying the projection coordinate system and the spatial scale.

[0008] According to an industrial zoning demarcation method provided by the present invention, the ecosystem value index in step S2 specifically includes: A first value indicator, the first value indicator corresponding to the water conservation supply in the assessment area; a second value indicator, the second value indicator corresponding to the amount of soil erosion control in the assessment area; A third value indicator, wherein the third value indicator corresponds to the flood storage capacity of the assessment area.

[0009] According to an industrial zoning delineation method provided by the present invention, the expression of the first value index is:

[0010]

[0011] in, is the first evaluation indicator, For water conservation supply, is the operating cost per unit storage capacity of the reservoir, To evaluate the ecosystem type index value in the region, To assess the total number of ecosystem types in the region, For the The area of ​​the ecosystem, For the The amount of rainfall that produces runoff in this ecosystem, For the Surface runoff of ecosystems, For the Evapotranspiration of ecosystems.

[0012] According to an industrial zoning demarcation method provided by the present invention, the expression of the second value index is:

[0013]

[0014] in, is the second value indicator, is the sedimentation coefficient, is the soil erosion control amount, is the soil bulk density, is the unit reservoir desilting project cost, is the pollutant type index value, is the total number of pollutant types, For the The pure content of the pollutant, For the The unit treatment cost of each pollutant is is the rainfall erosivity factor, is the soil erodibility factor, is the slope length factor, is the slope factor, It is the vegetation cover management factor.

[0015] According to an industrial zoning demarcation method provided by the present invention, the expression of the third value index is:

[0016]

[0017] in, is the third value indicator, For flood storage capacity, is the engineering cost per unit storage capacity of the reservoir, For the Stormwater runoff from ecosystems is the unit area.

[0018] According to an industrial zoning demarcation method provided by the present invention, the spatial development value index in step S3 specifically includes: Distance from the preset main road, distance from the preset main river, terrain slope, altitude, precipitation, and vegetation index.

[0019] According to an industrial zoning demarcation method provided by the present invention, step S4 further comprises: S41: constructing an ecosystem value map according to the ecosystem value indicators; S42: constructing a space development value map according to the space development value index; S43: The ecosystem value map and the spatial development value map are coupled through an ArcGIS software raster calculator to complete zoning delineation.

[0020] The present invention also provides an industrial zoning demarcation system, comprising: Acquisition module: used to collect spatial information data of the assessment area; A first calculation module: used to calculate and obtain the ecosystem value index of the assessment area based on the spatial information data; A second calculation module: used for calculating and obtaining the spatial development value index of the assessment area based on the spatial information data; Coupling module: used to zoning the assessment area according to the ecosystem value index and the spatial development value index.

[0021] The present invention provides an industrial zoning demarcation method and system, which combines spatial development suitability and ecosystem service value to perform industrial development zoning demarcation. The present invention comprehensively considers the quantitative indicators of the ecosystem value and spatial development value of the assessment area, accurately describes and quantitatively evaluates the assessment area, and then constructs an ecosystem value map and a spatial development value map, and couples the two to achieve a visual presentation of the zoning demarcation results, thereby improving the accuracy and reliability of the zoning demarcation. BRIEF DESCRIPTION OF THE DRAWINGS

[0022] In order to more clearly illustrate the technical solutions in the present invention or the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0023] Figure 1 A schematic diagram of a process flow of an industrial zoning demarcation method provided by the present invention; Figure 2 A schematic diagram of the structure of an industrial zoning demarcation system provided by the present invention. DETAILED DESCRIPTION

[0024] In order to make the purpose, technical solution and advantages of the present invention clearer, the technical solution of the present invention will be clearly and completely described below in conjunction with the drawings in the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments, and they should not be understood as limitations on the present invention. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention. In the description of the present invention, it should be understood that the terms used are only for descriptive purposes and cannot be understood as indicating or implying relative importance.

[0025] Combine the following Figure 1 to Figure 2 Embodiments of the present invention are described.

[0026] like Figure 1 As shown, the present invention provides an industrial zoning demarcation method, comprising: S1: Collect spatial information data of the assessment area.

[0027] The spatial information data in step S1 specifically includes: first spatial information data, wherein the first spatial information data is remote sensing data; The second spatial information data includes topographic surveying and mapping data, hydrological monitoring data, land use data and vegetation index data.

[0028] In the process of data collection, topographic mapping data is collected through digital elevation models (DEMs) obtained by sensors carried by satellites or drones, because remote sensing technology can obtain topographic data on a large scale and efficiently; hydrological monitoring data is mainly obtained through hydrological monitoring stations, hydrological sensors and other equipment; and land use data and vegetation index data use GIS software to overlay and analyze remote sensing images, topographic mapping data, etc., to generate thematic maps such as land use maps, and at the same time, through field visits and investigations, to understand the current land use status, land ownership and other information; vegetation index data is obtained through multispectral remote sensing, using multispectral sensors carried by satellites or drones to obtain surface reflectance data to reflect the growth status and coverage of vegetation.

[0029] Wherein, step S1 also includes: Preprocessing the first spatial information data in the spatial information data specifically includes formatting the remote sensing data into a raster format and unifying the projection coordinate system and the spatial scale.

[0030] After collecting and acquiring remote sensing data, it is first necessary to unify it into a raster format. Raster data is a data model that divides geographic space into regular grid cells and assigns attribute values ​​to each cell to represent spatial distribution. In addition, since remote sensing data may come from different sensors and platforms, they may use different projection coordinate systems. Therefore, the present invention also unifies the coordinate system and spatial scale. The projection coordinate system defines how to convert points on the earth's surface into a plane coordinate system, which is crucial for subsequent data analysis and comparison.

[0031] S2: Calculate and obtain the ecosystem value index of the assessment area based on the spatial information data.

[0032] Furthermore, step S2 quantitatively evaluates the supply of each ecosystem service in the designated area and calculates the ecosystem service value through ecosystem service indicators including water conservation, soil conservation, and flood regulation.

[0033] The ecosystem value indicators in step S2 specifically include: The first value indicator corresponds to the water conservation supply in the assessment area.

[0034] The expression of the first value indicator is:

[0035]

[0036] in, is the first evaluation indicator, For water conservation supply, is the operating cost per unit storage capacity of the reservoir, To evaluate the ecosystem type index value in the region, To assess the total number of ecosystem types in the region, For the The area of ​​the ecosystem, For the The amount of rainfall that produces runoff in this ecosystem, For the Surface runoff of ecosystems, For the Evapotranspiration of ecosystems.

[0037] For the first value indicator, this embodiment calculates the supply of water conservation services by using the water balance equation in the assessment area, and for the value calculation, the shadow engineering method is used to calculate the water conservation value by simulating the cost required to build water conservancy facilities with a water storage capacity equivalent to the water conservation capacity.

[0038] In a specific embodiment, the grid data of land use, rainfall, and evapotranspiration are obtained through public data, and the surface runoff data are calculated by rainfall (mm) and the average surface runoff coefficient. The calculation process is carried out in the ArcGIS 10.2 raster calculator tool, and the obtained surface runoff coefficient is shown in Table 1.

[0039]

[0040] Table 1 Surface runoff coefficient A second value indicator corresponds to the amount of soil erosion control in the assessment area. The expression of the second value indicator is:

[0041]

[0042] in, is the second value indicator, is the sedimentation coefficient, is the soil erosion control amount, is the soil bulk density, is the unit reservoir desilting project cost, is the pollutant type index value, is the total number of pollutant types, For the The pure content of the pollutant, For the The unit treatment cost of each pollutant is is the rainfall erosivity factor, is the soil erodibility factor, is the slope length factor, is the slope factor, It is the vegetation cover management factor.

[0043] For the second evaluation index, the embodiment of the present invention uses the replacement cost method, i.e., the cost accounting of pollutant treatment to calculate the value of reducing non-point source pollution based on the soil conservation amount and the nitrogen and phosphorus content in the soil, where To reduce the value of silt deposition, To reduce the value of non-point source pollution.

[0044] For the above-mentioned soil erosion control amount, the RUSLE soil loss equation is used in the assessment area to calculate the soil conservation service supply. After standardization, the ecological degradation space of specific ecosystem services is formed. The specific calculation formula is shown above.

[0045] For the calculation of soil conservation service supply, the soil conservation amount is first selected, that is, the amount of soil erosion reduced by the ecosystem, and the difference between the potential soil erosion and the actual soil erosion is used as the evaluation index of the soil conservation function of the ecosystem. Among them, the actual soil erosion refers to the amount of soil erosion under the current surface vegetation coverage, and the potential soil erosion refers to the amount of soil erosion that may occur without surface vegetation coverage.

[0046] In the RUSLE soil and water loss equation, the calculation principles of each factor are as follows.

[0047] Rainfall erosivity factor The monthly scale calculation formula is used for quantification, and the specific expression is:

[0048] In the above formula, is the average annual rainfall, in mm; is the average monthly rainfall in month i, in mm.

[0049] Soil erodibility factor , the calculation formula is as follows:

[0050]

[0051] In the formula, represents the soil erodibility factor before correction, represents the corrected soil erodibility factor, , , They are the percentage contents of clay particles <0.002 mm, powder particles 0.002 mm to 0.05 mm, and sand particles 0.05 mm to 2 mm. is the percentage content of organic carbon; Vegetation cover factor : Its value is between 0 and 1. The smaller the value, the stronger the soil and water conservation capacity of vegetation and the smaller the soil loss; Forest, bush, grassland The value is calculated based on the vegetation coverage of the study area in that year:

[0052] In the formula, is the vegetation coverage.

[0053] A third value indicator, wherein the third value indicator corresponds to the flood storage capacity of the assessment area.

[0054] Wherein, the expression of the third value indicator is:

[0055]

[0056] in, is the third value indicator, For flood storage capacity, is the engineering cost per unit storage capacity of the reservoir, For the Stormwater runoff from ecosystems is the unit area.

[0057] For the third value indicator, in the assessment area, the flood storage equation is used to measure the supply of flood storage services through the flood storage capacity of vegetation and water bodies, and a specific service degradation space is formed after standardization. The calculation method of flood storage value is to calculate the economic value of flood storage to reduce the threat of floods, and use the replacement cost method, that is, the construction cost of the reservoir to calculate the flood storage value of the ecosystem. The value of flood storage is calculated using the replacement engineering method and the construction cost of the reservoir, and then the third value indicator and the corresponding calculation formula of flood storage can be obtained.

[0058] In a specific embodiment, a daily rainfall exceeding 50 mm is defined as a rainstorm event, the interpolation process is completed in the ArcGIS 10.2 Inverse Distance Weighted (IDW) tool, the rainstorm surface runoff raster data is the product of the rainstorm rainfall and the runoff coefficient, and the calculation process is completed in the ArcGIS 10.2 Raster Calculator tool. For the calculation of the supply of flood storage services, the calculation process is also performed in the ArcGIS 10.2 Raster Calculator tool.

[0059] After completing the calculation of the above multiple value indicators, it is necessary to combine the factors to obtain the overall evaluation result of the ecosystem service value. The combination of factors adopts the following formula:

[0060] In the formula, It is the comprehensive evaluation result of regional ecological suitability; For the The evaluation value of each factor; For the The weight of a factor.

[0061] S3: Calculate and obtain the spatial development value index of the assessment area based on the spatial information data.

[0062] The space development value index in step S3 specifically includes: Distance from the preset main road, distance from the preset main river, terrain slope, altitude, precipitation, and vegetation index.

[0063] In this embodiment, the selected spatial development value index is as above. Specifically, first, in the evaluation area, the distance to the main road is calculated, because the distance to the main road determines the land development cost and location advantage. Generally, the closer the land development is to the road, the easier it is to form a spatial agglomeration effect. For the distance to the river, the construction and development of land have a great impact on the ecological environment of the nearby river. In principle, the development land should be as far away from the river as possible to avoid damage to the river ecosystem. For the calculation of slope, the terrain is one of the limiting factors affecting land construction and development. The larger the slope, the lower the feasibility of development and construction. The slope is derived from DEM data. For altitude: the impact of altitude on vegetation coverage, production and life, the higher the altitude, the less conducive to development. For rainfall, rainfall affects the distribution of water resources in the basin and is a very important factor for plant growth. For the calculation of vegetation coverage, the height reflects the vegetation coverage, which is one of the important indicators for measuring the living environment and one of the ecological conditions for human production and life. The vegetation index is often expressed by the normalized difference vegetation index (NDVI), and the calculation formula is:

[0064] In the formula, is the reflectivity in the near-infrared band, It is the reflectance of the red light band. After the NDVI value is calculated, it needs to be standardized. The standardized value is 0~100.

[0065] After obtaining the relevant data of the above indicators, it is also necessary to combine the factors to obtain the overall evaluation results. The combination of factors uses the following formula:

[0066] In the formula, The comprehensive evaluation results of the suitability of national land space development are For the The evaluation value of the factor, For the In a specific embodiment, the specific weight settings of the above indicators are shown in Table 2.

[0067]

[0068] Table 2 Index weight settings in spatial development value indicators S4: Zoning the assessment area based on the ecosystem value index and the spatial development value index.

[0069] Wherein, step S4 further comprises: S41: constructing an ecosystem value map according to the ecosystem value indicators; S42: constructing a space development value map according to the space development value index; S43: The ecosystem value map and the spatial development value map are coupled through an ArcGIS software raster calculator to complete zoning delineation.

[0070] The actual meaning of step S2 and step S3 is to perform quantitative calculations on multiple data indicators from the two aspects of protection and development, and finally to integrate the evaluation results of the two aspects of data indicators in step S4 to perform regional industrial zoning calculations.

[0071] In a specific embodiment, the present invention sorts the ecosystem value indicators and spatial development value indicators in steps S2 and S3, and divides the indicators of the two aspects into three levels: high, medium and low according to the division threshold of 1 / 3, and then demarcates the regions according to the divided indicator levels.

[0072] The present invention sets the area with high ecosystem service value and high development suitability as an ecotourism area; sets the area with high ecosystem service value and moderate development suitability as an industrial scientific research base; sets the area with high ecosystem service value and low development suitability as a nature reserve; sets the area with low ecosystem service value and high development suitability as an industrial area; sets the area with low ecosystem service value and moderate development suitability as an agricultural product processing area; sets the area with low ecosystem service value and low development suitability as an agricultural area; sets the area with moderate ecosystem service value and high national land space development suitability as a modern agricultural park; sets the area with moderate ecosystem service value and moderate national land space development suitability as an ecological industrial park; and sets the area with moderate ecosystem service value and low national land space development suitability as an ecological agricultural area.

[0073] For the intuitiveness of the data, in steps S41 to S43, corresponding visualization and data processing can be performed. The specific process is to first process and analyze the corresponding spatial data based on the ecosystem value index calculated in step S2, and use ArcGIS to interpolate or classify the processed spatial data, so as to convert the discrete ecosystem value index into a continuous or classified spatial distribution map, and use ArcGIS's mapping function to display the calculated ecosystem value in the form of a map, so as to intuitively reflect the value of the ecosystem in different regions. Then for the spatial development value index mentioned in step S3, it is necessary to generate a spatial development value map based on the comprehensive evaluation results using ArcGIS mapping function to display the spatial development potential or suitability of each region in the evaluation area. Finally, according to the planning objectives and actual needs, design a suitable calculation formula or logic rule, couple the raster data of the two layers, and the area with high ecosystem value and suitable spatial development value can be demarcated as a specific industrial zone by weighted summation, logical judgment, etc.

[0074] like Figure 2 As shown, the present invention also provides an industrial zoning demarcation system, comprising: Collection module 100: used to collect spatial information data of the evaluation area; The first calculation module 200 is used to calculate and obtain the ecosystem value index of the assessment area based on the spatial information data; The second calculation module 300 is used to calculate and obtain the spatial development value index of the assessment area based on the spatial information data; Coupling module 400: used to zoning the assessment area according to the ecosystem value index and the spatial development value index.

[0075] The device embodiments described above are merely illustrative, wherein the units described as separate components may or may not be physically separated, and the components displayed as units may or may not be physical units, that is, they may be located in one place, or they may be distributed on multiple network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the scheme of this embodiment. Ordinary technicians in this field can understand and implement it without paying creative labor.

[0076] Through the description of the above implementation methods, those skilled in the art can clearly understand that each implementation method can be implemented by means of software plus a necessary general hardware platform, and of course, can also be implemented by hardware. Based on this understanding, the above technical solution is essentially or the part that contributes to the prior art can be embodied in the form of a software product, and the computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, a disk, an optical disk, etc., including a number of instructions for a computer device (which can be a personal computer, a server, or a network device, etc.) to execute the methods described in each embodiment or some parts of the embodiments.

[0077] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for delineating industrial zones, characterized in that: include: S1: Collect spatial information data of the assessment area; S2: Calculate and obtain the ecosystem value index of the assessment area based on the spatial information data; S3: Calculate and obtain the spatial development value index of the assessment area based on the spatial information data; S4: Zoning the assessment area based on the ecosystem value index and the spatial development value index.

2. The method for delineating industrial zones according to claim 1, characterized in that: The spatial information data in step S1 specifically includes: first spatial information data, wherein the first spatial information data is remote sensing data; The second spatial information data includes topographic surveying and mapping data, hydrological monitoring data, land use data and vegetation index data.

3. The method for delineating industrial zones according to claim 2, characterized in that: Step S1 also includes: Preprocessing the first spatial information data in the spatial information data specifically includes formatting the remote sensing data into a raster format and unifying the projection coordinate system and the spatial scale.

4. The method for delineating industrial zoning according to claim 1, characterized in that: The ecosystem value indicators in step S2 specifically include: A first value indicator, the first value indicator corresponding to the water conservation supply in the assessment area; a second value indicator, the second value indicator corresponding to the amount of soil erosion control in the assessment area; A third value indicator, wherein the third value indicator corresponds to the flood storage capacity of the assessment area.

5. The method for delineating industrial zoning according to claim 4, characterized in that: The expression of the first value indicator is: in, is the first evaluation index, For water conservation supply, is the operating cost per unit storage capacity of the reservoir, To evaluate the ecosystem type index value in the region, To assess the total number of ecosystem types in the region, For the The area of ​​the ecosystem, For the The amount of rainfall that produces runoff in this ecosystem, For the Surface runoff of ecosystems, For the Evapotranspiration of ecosystems.

6. The method for delineating industrial zones according to claim 4, characterized in that: The expression of the second value indicator is: in, is the second value indicator, is the sedimentation coefficient, is the soil erosion control amount, is the soil bulk density, is the unit reservoir desilting project cost, is the pollutant type index value, is the total number of pollutant types, For the The pure content of the pollutant, For the The unit treatment cost of each pollutant is is the rainfall erosivity factor, is the soil erodibility factor, is the slope length factor, is the slope factor, It is the vegetation cover management factor.

7. The method for delineating industrial zones according to claim 5, characterized in that: The expression of the third value indicator is: in, is the third value indicator, For flood storage capacity, is the engineering cost per unit storage capacity of the reservoir, For the Stormwater runoff from ecosystems is the unit area.

8. The method for delineating industrial zones according to claim 1, characterized in that: The space development value index in step S3 specifically includes: Distance from the preset main road, distance from the preset main river, terrain slope, altitude, precipitation, and vegetation index.

9. The method for delineating industrial zones according to claim 1, characterized in that: Step S4 further comprises: S41: constructing an ecosystem value map according to the ecosystem value indicators; S42: constructing a space development value map according to the space development value index; S43: The ecosystem value map and the spatial development value map are coupled through an ArcGIS software raster calculator to complete zoning delineation.

10. An industrial zoning system, characterized in that: include: Acquisition module: used to collect spatial information data of the assessment area; A first calculation module: used to calculate and obtain the ecosystem value index of the assessment area based on the spatial information data; A second calculation module: used for calculating and obtaining the spatial development value index of the assessment area based on the spatial information data; Coupling module: used to zoning the assessment area according to the ecosystem value index and the spatial development value index.

Images