A Method, Device, Equipment and Medium for Identifying the Service Scope of Urban Ecological Source Areas

Abstract

The present invention discloses a method, device, equipment and medium for identifying the service scope of urban ecological source areas. The method includes: obtaining land use data and ecosystem service value of the area to be identified, and constructing an ecological source area set; constructing a spatial analysis grid, and obtaining impedance evaluation data based on the resistance value, elevation data and slope data of grid cells; taking the impedance evaluation data as the cost, calculating the minimum cumulative cost distance from the grid cells to the ecological source areas, and obtaining a service impedance evaluation result; calculating the ecological service value of the ecological source areas to the grid cells by using a gravity model; dividing the maximum service scope of the ecological source areas according to the ecological source area corresponding to the maximum value in the ecological service value of the grid cells; identifying the edge areas served by the ecological source areas, and removing the edge areas from the maximum service scope of the ecological source areas to obtain the core area of the service of the ecological source areas. The present invention improves the accuracy and comprehensiveness of identifying the service scope of urban ecological source areas.

Description

Technical Field

[0001] The present invention relates to the technical field of urban ecological planning, and particularly to a method, device, equipment and medium for identifying the service range of urban ecological source areas. Background Art

[0002] Urban ecological source areas are the core components of urban ecological systems. Accurately identifying ecological source areas and reasonably dividing their service ranges can not only optimize land use, provide a scientific basis for urban planning, ensure the balanced distribution of ecosystem services, but also help promote cross-administrative regional ecological cooperation, provide guidance for the construction of ecological networks, thereby enhancing the connectivity and stability of the overall ecological system.

[0003] With the acceleration of the urbanization process and the increasing complexity of ecological environment problems, there are several limitations in the existing technologies for identifying urban ecological source areas and dividing their service ranges. The current methods for dividing the service range of ecological source areas are not accurate and comprehensive enough. The commonly used buffer analysis method is simple and intuitive, but it ignores the influence of factors such as terrain and land use on the spatial spread of ecological services, resulting in low accuracy of service range division; another method based on landscape connectivity, although considering spatial heterogeneity, mainly focuses on species migration and is difficult to comprehensively reflect the multi-dimensional service range of ecological source areas. In addition, most of the existing technologies only focus on one aspect of ecological source area identification or service range division, lacking a systematic method to organically combine the two. At the same time, the existing service range division methods are not fine enough to effectively distinguish the core area and the edge area within the service range of ecological source areas. Summary of the Invention

[0004] To solve the above technical problems, the present invention provides a method, device, equipment and medium for identifying the service range of urban ecological source areas, which combines ecological source area identification and service range division, improves the accuracy and comprehensiveness of the identification of the service range of ecological source areas, and can identify a more refined service range of ecological source areas.

[0005] An embodiment of the present invention provides a method for identifying the service range of urban ecological source areas, including:

[0006] Obtaining land use data and ecosystem service value of the area to be identified, and constructing an ecological source area set;

[0007] Constructing a spatial analysis grid in the area to be identified, where the spatial analysis grid includes a number of grid cells;

[0008] Based on the resistance value, elevation data and slope data corresponding to the grid cells, obtaining impedance evaluation data of the grid cells;

[0009] For each ecological source area in the ecological source area concentration, taking the impedance evaluation data as the cost, calculate the minimum cumulative cost distance from each grid cell to this ecological source area to obtain the service impedance evaluation result of the ecological source area;

[0010] According to the ecosystem service value and the service impedance evaluation result, use the gravity model to calculate the ecological service value of each ecological source area to the grid cell;

[0011] According to the ecological source area corresponding to the maximum value in the ecological service value of each grid cell, divide the grid cell to obtain the maximum service range of each ecological source area;

[0012] Based on the ecological service value, identify the marginal area of the ecological source area service, and exclude the marginal area from the maximum service range of the ecological source area to obtain the core area of the ecological source area service.

[0013] As an improvement of the above solution, the obtaining the land use data and ecosystem service value of the area to be identified and constructing the ecological source area set includes:

[0014] Obtain the land use data of the area to be identified, and based on the land use data, obtain the first area set corresponding to several natural land cover types in the area to be identified; delete the areas within the buffer range of the main traffic arteries in the first area set to obtain the second area set; merge the areas with the same natural land cover type in the second area set, and select the areas with an area greater than the first preset value in the merged area set to form the first data set;

[0015] Obtain the ecosystem service value of each area in the area to be identified, and according to the ecosystem service value, select the areas with an ecosystem service value greater than the second preset value in the area to be identified to form the second data set;

[0016] Take the intersection of the first data set and the second data set as the ecological source area set.

[0017] As an improvement of the above solution, the obtaining the impedance evaluation data of the grid cell based on the resistance value, elevation data, and slope data corresponding to the grid cell includes:

[0018] According to the resistance values of different land use areas preset, obtain the resistance value of the grid cell;

[0019] Obtain the elevation data and slope data corresponding to the grid cell;

[0020] Perform a weighted superposition operation on the resistance value, elevation data, and slope data of the grid cell to obtain the impedance evaluation data of the grid cell.

[0021] As an improvement to the above solution, the resistance values of the different land use areas are related to the land use types, urban layouts, and transportation networks of the corresponding land use areas; the land use types include at least one of forest land, water bodies, grassland, cultivated land, urban construction land, and other unused land; the higher the resistance value of the grid cell, the greater the obstruction of the corresponding area of the grid cell to the spread of ecological services.

[0022] As an improvement to the above solution, for each ecological source area in the ecological source areas, taking the impedance evaluation data as the cost, calculating the minimum cumulative cost distance from each grid cell to this ecological source area, and obtaining the service impedance evaluation result of this ecological source area, including:

[0023] For each ecological source area in the ecological source areas, taking the impedance evaluation data of the grid cell as the cost raster, calculating the minimum value of the sum of the impedance evaluation data on the path from each grid cell passing through the cost surface to reach this ecological source area, obtaining the raster data of the minimum cumulative cost distance, and forming a raster layer, where the raster layer is the service impedance evaluation result of this ecological source area;

[0024] Integrating the raster layers of all ecological source areas in the ecological source areas to obtain the service impedance evaluation results of all ecological source areas.

[0025] As an improvement to the above solution, the gravity model is specifically:

[0026]

[0027] where F eco is the ecological service value of an ecological source area to a grid cell, V total is the total ecological system service value of the corresponding area of this grid cell, V eco is the ecological system service value of this ecological source area, R cost_dist is the service impedance evaluation result of this ecological source area within the area to be identified.

[0028] As an improvement to the above solution, identifying the edge area of the ecological source area service based on the ecological service value, including:

[0029] Based on the ecological service value, if the grid cell meets the preset threshold condition, it is determined that the grid cell belongs to the edge area of the ecological source area service; the formula for the threshold condition is:

[0030]

[0031] where F eco-max and F eco-minare respectively the maximum and minimum values of the ecological service values of the ecological source areas for the grid cell, and F eco-all is the sum of the ecological service values of all ecological source areas for the grid cell, and α is a preset relative difference threshold.

[0032] An embodiment of the present invention also provides an urban ecological source area service scope identification device, including:

[0033] An ecological source area construction module, configured to obtain land use data and ecosystem service value of the area to be identified, and construct an ecological source area set;

[0034] A grid construction module, configured to construct a spatial analysis grid in the area to be identified, and the spatial analysis grid includes a number of grid cells;

[0035] An impedance evaluation module, configured to obtain impedance evaluation data of the grid cell based on the resistance value, elevation data, and slope data corresponding to the grid cell;

[0036] A comprehensive impedance module, for each ecological source area in the ecological source area set, using the impedance evaluation data as the cost, calculating the minimum cumulative cost distance from each grid cell to the ecological source area, and obtaining the service impedance evaluation result of the ecological source area;

[0037] A service value calculation module, configured to calculate the ecological service value of each ecological source area for the grid cell by using a gravity model according to the ecosystem service value and the service impedance evaluation result;

[0038] A maximum service scope module, configured to divide the grid cell according to the ecological source area corresponding to the maximum value in the ecological service value of each grid cell, and obtain the maximum service scope of each ecological source area;

[0039] A fine service scope module, configured to identify the edge area of the ecological source area service based on the ecological service value, and remove the edge area from the maximum service scope of the ecological source area to obtain the core area of the ecological source area service.

[0040] An embodiment of the present invention also provides a computer device, including a processor and a memory, where a computer program is stored in the memory, and the computer program is configured to be executed by the processor, and when the processor executes the computer program, the urban ecological source area service scope identification method described in any one of the above is implemented.

[0041] An embodiment of the present invention also provides a computer-readable storage medium, where the computer-readable storage medium stores a computer program, and when the computer program runs, it controls the device where the computer-readable storage medium is located to execute the urban ecological source area service scope identification method described in any one of the above.

[0042] Compared with the prior art, the beneficial effects of a method, device, equipment and medium for identifying the service range of urban ecological source areas provided by the embodiments of the present invention are as follows: By combining the ecological patch method and the ecological importance evaluation method for ecological source area identification, high-quality ecological source areas can be identified, improving the accuracy and flexibility of ecological source area identification; By comprehensively evaluating multiple factors for the service impedance of ecological source areas, the comprehensiveness and accuracy of the identification of the service range of ecological source areas are improved; By constructing a gravity model to calculate the ecological service value, the accuracy of the identification of the service range of ecological source areas is further improved; By setting a relative difference threshold to divide the edge area and the core area of the ecological source area, a refined division of the service range of the ecological source area is achieved, providing a basis for the precise allocation of urban ecological protection resources and helping to improve resource utilization efficiency. The embodiments of the present invention can not only accurately identify high-value ecological source areas, but also scientifically divide their service ranges, improving the accuracy and comprehensiveness of the identification of the service range of ecological source areas, and providing strong technical support for urban ecological planning and management. BRIEF DESCRIPTION OF THE DRAWINGS

[0043] Figure 1 is a schematic flowchart of a method for identifying the service range of urban ecological source areas provided by the embodiments of the present invention;

[0044] Figure 2 is a schematic diagram of an ecological source area set in a method for identifying the service range of urban ecological source areas provided by the embodiments of the present invention;

[0045] Figure 3 is a schematic diagram of an impedance factor evaluation data set in a method for identifying the service range of urban ecological source areas provided by the embodiments of the present invention;

[0046] Figure 4 is a schematic diagram of the ecological service values of three certain ecological source areas in the area to be identified provided by the embodiments of the present invention;

[0047] Figure 5 is a schematic diagram of the maximum service range of an ecological source area in a method for identifying the service range of urban ecological source areas provided by the embodiments of the present invention;

[0048] Figure 6 is a schematic diagram of the stable core hinterland and potential vulnerable hinterland of an ecological source area in a method for identifying the service range of urban ecological source areas provided by the embodiments of the present invention;

[0049] Figure 7 is a schematic structural diagram of a device for identifying the service range of urban ecological source areas provided by the embodiments of the present invention;

[0050] Figure 8 is a schematic structural diagram of a computer device provided by the embodiments of the present invention. Detailed implementation manners

[0051] The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. Apparently, the described embodiments are only a part of the embodiments of the present invention, rather than all the embodiments. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present invention without creative efforts shall fall within the protection scope of the present invention.

[0052] Please refer to Figure 1 , Figure 1 which is a schematic flowchart of a method for identifying the service range of urban ecological source areas provided by an embodiment of the present invention. The method for identifying the service range of urban ecological source areas includes:

[0053] S1: Obtain the land use data and ecosystem service value of the area to be identified, and construct an ecological source area set;

[0054] As an optional embodiment, the obtaining the land use data and ecosystem service value of the area to be identified and constructing an ecological source area set includes:

[0055] Obtain the land use data of the area to be identified, and based on the land use data, obtain a first area set corresponding to several natural land cover types in the area to be identified; delete the areas within the buffer range of the main traffic arteries in the first area set to obtain a second area set; merge the areas with the same natural land cover type in the second area set, and select the areas with an area greater than a first preset value in the merged area set to form a first data set;

[0056] Obtain the ecosystem service value of each area in the area to be identified, and select the areas with an ecosystem service value greater than a second preset value in the area to be identified according to the ecosystem service value to form a second data set;

[0057] Take the intersection of the first data set and the second data set as the ecological source area set.

[0058] Preferably, the first preset value is the area value of the top 20% in terms of total area ratio; the second preset value is the ecosystem service value of the top 20% in terms of ecosystem service value ratio.

[0059] Specifically, as the area providing high-value ecological services, the ecological source area is a key element of the urban ecosystem. The embodiments of the present invention adopt the ecological patch method and the ecological importance evaluation method to identify potential ecological source areas, and then obtain the final ecological source area set through screening.

[0060] First, using the land use data of the city to be identified, select the areas corresponding to natural land cover types such as cultivated land, grassland, and water area to form the first area set; to exclude human interference, delete the areas within the 500-meter buffer zone of the main transportation arteries, where the main transportation arteries include national highways, provincial highways, county roads, expressways, railways, and urban first- and second-class roads; then, perform spatial geographical merging on the land with the same land use nature to form a preliminary ecological source data base plate, and then extract the areas with the top 20% total area ratio from it to form the first data set.

[0061] Secondly, evaluate the area to be identified through the ecosystem service value index, and extract the areas with the top 20% ecosystem service value ratio from it to form the second data set. Among them, the ecosystem service value can be calculated through the InVEST model (Integrated Valuation of Ecosystem Services and Tradeoffs model) or obtained from the Resource and Environment Science Data Platform of the Chinese Academy of Sciences. The ecosystem service value index aims to quantify the contribution degree of each ecological patch to the urban ecosystem.

[0062] Finally, perform a spatial overlay analysis on the first data set and the second data set, take the intersection of the two, and finally form the ecological source identification result, that is, the ecological source set, as Figure 2 shown. The embodiment of the present invention ensures that the identified ecological sources not only have sufficient area scale but also have high ecological value.

[0063] It should be noted that the embodiment of the present invention uses 20% as the relative threshold for identifying high-value ecological sources. This threshold is based on the Pareto principle, which shows that about 20% of the key areas provide 80% of the ecosystem services. However, 20% is not a fixed standard. The embodiment of the present invention designs an adjustable threshold mechanism: in cities with rich ecological resources, it can be increased to 25% to 30%; in cities with scarce resources, it can be reduced to 15% to 20%. In actual application, it can also be fine-tuned according to expert opinions and local ecological protection goals to ensure the flexibility and applicability of the method.

[0064] In the prior art, the ecological source identification methods adopted mostly rely on single or simple combined indicators, such as directly defining a specific land use type as an ecological source, or only evaluating the importance of ecological patches through landscape indices, without fully considering the functional value and service ability of the ecosystem, resulting in an incomplete evaluation of the importance of ecological sources. However, the embodiment of the present invention not only considers the spatial characteristics of ecological patches, such as area and connectivity, but also incorporates the value evaluation of ecosystem services, thus providing a more comprehensive and accurate identification of high-ecological-value source areas and overcoming the limitations of a single perspective.

[0065] S2: Construct a spatial analysis grid in the area to be recognized, where the spatial analysis grid includes a number of grid cells;

[0066] Specifically, establish a uniform fishing net of 1 km × 1 km in the area to be recognized as the basic spatial unit for subsequent analysis.

[0067] S3: Obtain the impedance evaluation data of the grid cell based on the resistance value, elevation data, and slope data corresponding to the grid cell;

[0068] As one optional embodiment, the obtaining the impedance evaluation data of the grid cell based on the resistance value, elevation data, and slope data corresponding to the grid cell includes:

[0069] Obtain the resistance value of the grid cell according to the resistance values of different land use areas preset;

[0070] Obtain the elevation data and slope data corresponding to the grid cell;

[0071] Perform a weighted superposition operation on the resistance value, elevation data, and slope data of the grid cell to obtain the impedance evaluation data of the grid cell.

[0072] As one optional embodiment, the resistance values of different land use areas are related to the land use types, urban layouts, and transportation networks of the corresponding land use areas; the land use types include at least one of forest land, water body, grassland, cultivated land, urban construction land, and other unused land; the higher the resistance value of the grid cell, the greater the obstacle to the spread of ecological services in the area corresponding to the grid cell.

[0073] Specifically, the impact of the ecological source area on the surrounding environment is a complex process restricted by multiple factors, which mainly include: natural factors (such as terrain elevation and slope), land use types (such as forest land, water area, agricultural land, and construction land), urban planning and layout (such as the accessibility of transportation networks), and distance attenuation, etc. In order to comprehensively evaluate the impact of these impedance factors on the services of the ecological source area, the embodiments of the present invention obtain the impedance evaluation data of each grid cell based on the resistance value, elevation data, and slope data, and form an impedance factor evaluation data set, as Figure 3 shown.

[0074] First, rasterize and reclassify the land use data in the area to be recognized, where the raster size is the same as the size of the above grid cell; comprehensively consider factors such as land use types, urban layouts, and transportation networks, and assign resistance values to different types of land use. Specifically, as shown in Table 1, these resistance values reflect the impact degree of different land types on ecological processes. Thus, the resistance value corresponding to each grid cell can be obtained.

[0075] Table 1 Resistance values of different land uses

[0076]

[0077] It should be noted that the specific resistance values can be adjusted according to the actual situations and key concerns of different cities to ensure the accuracy and practicality of the evaluation.

[0078] After obtaining the resistance values, the elevation and slope factors are also incorporated into the resistance evaluation. Specifically, the elevation data and slope data of each grid cell are obtained from the land use data of the area to be identified, and weighted superposition operations are performed by assigning weights to the resistance values, elevation data, and slope data to calculate the impedance evaluation data, forming an impedance factor evaluation data set corresponding to the area to be identified.

[0079] S4: For each ecological source area in the ecological source area set, using the impedance evaluation data as the cost, calculate the minimum cumulative cost distance from each grid cell to this ecological source area to obtain the service impedance evaluation result of this ecological source area;

[0080] As one optional embodiment, for each ecological source area in the ecological source area set, using the impedance evaluation data as the cost, calculate the minimum cumulative cost distance from each grid cell to this ecological source area to obtain the service impedance evaluation result of this ecological source area, including:

[0081] For each ecological source area in the ecological source area set, using the impedance evaluation data of the grid cell as the cost raster, calculate the minimum value of the sum of the impedance evaluation data on the path from each grid cell passing through the cost surface to reach this ecological source area to obtain the raster data of the minimum cumulative cost distance, forming a raster layer, and this raster layer is the service impedance evaluation result of this ecological source area;

[0082] Integrate the raster layers of all ecological source areas in the ecological source area set to obtain the service impedance evaluation results of all ecological source areas.

[0083] Specifically, in order to quantify the difficulty of service acquisition in each area to the ecological source area and thus evaluate the impact of impedance factors on the spread of ecological services, the embodiment of the present invention uses the impedance evaluation data as the cost raster and combines the cost distance analysis function of GIS (Geographic Information System) to calculate the minimum value of the sum of the impedance evaluation data on the path when each grid cell in the area to be identified passes through the cost surface to reach the ecological source area, that is, the minimum cumulative cost distance. Specifically, the following operations are performed one by one for each ecological source area in the ecological source area set:

[0084] Input the ecological source area set as the source data into GIS;

[0085] Use the impedance factor evaluation data set corresponding to the area to be recognized as the data set of the cost raster;

[0086] Utilize the cost distance analysis function of GIS to generate raster data representing the minimum cumulative cost distance from grid cells to the nearest ecological source area.

[0087] After iteration, multiple raster layers are generated. Each raster layer represents the service impedance distribution of a specific ecological source area. Finally, these layers are integrated to form the service impedance evaluation result of all ecological source areas for the area to be recognized. This service impedance evaluation result is a raster data set. If there are n ecological source areas, the service impedance evaluation result contains n raster files. The embodiments of the present invention can reflect the service impedance characteristics of each ecological source area in detail and comprehensively consider the influence of all source areas, which helps to comprehensively understand the spatial distribution and acquisition difficulty of urban ecological services.

[0088] S5: According to the ecosystem service value and the service impedance evaluation result, use the gravity model to calculate the ecological service value of each ecological source area for the grid cell;

[0089] Specifically, in order to determine the service scope of each ecological source area, the embodiments of the present invention innovatively apply the gravity model to the division of the ecological service scope. This model was originally widely used in the field of economic geography to describe the economic interaction between cities. The embodiments of the present invention construct a model that can accurately describe the spatial distribution of ecological services by analogizing the ecosystem service value to "mass" and the spatial resistance (i.e., the service impedance evaluation result) to "distance".

[0090] As an optional embodiment, the gravity model is specifically:

[0091]

[0092] Among them, F eco is the ecological service value of a certain ecological source area for a certain grid cell, V total is the total ecosystem service value of the area corresponding to the grid cell, V eco is the ecosystem service value of this ecological source area, R cost_dist is the service impedance evaluation result of this ecological source area within the area to be recognized.

[0093] Furthermore, through the calculation of the gravity model, the ecological service value of each ecological source area within the area to be recognized can be obtained. As Figure 4 shown, it is a schematic diagram of the ecological service values of three ecological source areas provided by the embodiments of the present invention within the area to be recognized.

[0094] The embodiments of the present invention consider the service capabilities of ecological source areas themselves, incorporate spatial resistance factors, and finally quantify the relative impacts of different ecological source areas on the same region, which can more accurately reflect the spatial attenuation characteristics of ecological services, and further achieve a more accurate identification of the service scope of ecological source areas.

[0095] S6: Divide each grid cell according to the ecological source area corresponding to the maximum value in the ecological service value of each grid cell to obtain the maximum service scope of each ecological source area;

[0096] Specifically, based on the calculation results of the gravity model, for each grid cell, the principle of "selecting the largest and attributing" is adopted, and it is divided into the service scope of the ecological source area that has the greatest impact on it. Among them, for a certain grid cell, the ecological service values of each ecological source area on this grid cell can be obtained from the calculation results of the gravity model, and the maximum value among these ecological service values is selected. Then, this grid cell is divided into the service scope of the ecological source area corresponding to this maximum value. After all grid cells are divided into the corresponding ecological source area ranges one by one, the maximum service scope of each ecological source area can be obtained, as Figure 5 shown.

[0097] S7: Based on the ecological service value, identify the edge area of the ecological source area service, and exclude the edge area from the maximum service scope of the ecological source area to obtain the core area of the ecological source area service.

[0098] Specifically, after obtaining the maximum service scope of the ecological source area, analyze the situation where each grid cell is affected by multiple ecological source areas at the same time to identify the potential vulnerable hinterland in the ecological source area, that is, the edge area.

[0099] As one of the optional embodiments, the identifying the edge area of the ecological source area service based on the ecological service value includes:

[0100] Based on the ecological service value, if the grid cell meets the preset threshold condition, it is determined that the grid cell belongs to the edge area of the ecological source area service; the formula of the threshold condition is:

[0101]

[0102] Among them, F eco-max and F eco-min are respectively the maximum value and the minimum value in the ecological service values of the ecological source area on the grid cell, F eco-all is the sum of the ecological service values of all ecological source areas on the grid cell, and α is a preset relative difference threshold.

[0103] Preferably, the relative difference threshold α is set to 25%.

[0104] Specifically, when a certain grid cell satisfies it is considered that the multiple ecological source areas affecting this grid cell have similar impacts, and it is not clear which ecological source area it belongs to. It is defined as a potential vulnerable hinterland, that is, the marginal area.

[0105] It should be noted that the relative difference threshold α can be adjusted according to the specific characteristics of the city and the research objectives. For example, for a city with a relatively fragile ecosystem, this threshold can be lowered to more strictly define the stable core area.

[0106] Furthermore, the potential vulnerable hinterland is excluded from the maximum service range of each ecological source area, and the remaining area is the stable core hinterland of this ecological source area, that is, the core area. The core area is mainly affected by a specific ecological source area, and the ecological services are relatively stable and sufficient. As Figure 6 shown, it is a schematic diagram of the stable core hinterland and the potential vulnerable hinterland of the ecological source area provided by the embodiment of the present invention. The stable core hinterland represents the main beneficiary area of the ecological source area services and should be protected and maintained with emphasis; while the potential vulnerable hinterland is the area with relatively weak ecological services, and targeted ecological restoration and strengthening measures need to be taken. The embodiment of the present invention effectively distinguishes the core area and the marginal area of the ecological services by quantitatively analyzing the service influence range of the ecological source area, not only clarifies the core influence area of the ecological source area, but also identifies the transitional zone with relatively weak services, providing more targeted guidance for urban ecological management.

[0107] The method for identifying the service range of the urban ecological source area provided by the embodiment of the present invention has broad application prospects and beneficial effects in multiple fields:

[0108] In terms of urban overall planning and green space system optimization, the present invention can be used to identify and protect key ecological source areas, assist in determining the ecological protection red line, and rationally layout urban construction land and ecological land. By analyzing the service range of the ecological source area, planners can identify areas with insufficient ecological service coverage, and plan new park green spaces in a targeted manner, so as to optimize the green space layout, improve the spatial balance of ecological services, and create a more livable environment for citizens.

[0109] In the planning and implementation of ecological restoration projects, the present invention can help determine the priority restoration areas. By identifying the potential vulnerable hinterland, urban managers can formulate targeted restoration plans, such as establishing ecological corridors in the vulnerable areas between ecological source areas, or increasing the green space area in the marginal areas, which helps to improve the connectivity and stability of the overall ecosystem and maximize the benefits of limited restoration resources.

[0110] For environmental impact assessment and project decision-making, the present invention can be used to analyze the potential impacts of large-scale construction projects on surrounding ecological source areas and their service ranges, helping to formulate more accurate environmental protection measures and minimize the negative impacts of projects on important ecological regions. At the same time, such analysis can provide a scientific basis for project site selection and scale decision-making, assisting in balancing development needs and ecological protection.

[0111] In terms of cross-administrative region ecological management and coordination, the present invention can help identify key ecological source areas and their service ranges spanning multiple administrative regions, providing a scientific basis for ecological collaborative management. This helps to break down administrative barriers, promote overall coordination of regional ecological governance, and can support the establishment of a more fair and effective cross-regional ecological compensation mechanism.

[0112] In urban resilience assessment and adaptive planning, the present invention can be used to analyze the distribution and functions of key ecological source areas, helping to identify ecologically vulnerable areas. This information provides important reference for formulating urban adaptive strategies, such as strengthening ecological construction in vulnerable areas and improving the city's ability to cope with extreme weather events. By dynamically monitoring changes in ecological source areas, urban development strategies can be adjusted in a timely manner to enhance the adaptability and resilience of the urban ecosystem and improve the overall resilience of the city.

[0113] Correspondingly, the present invention also provides a device for identifying the service range of urban ecological source areas, which can implement all processes of the method for identifying the service range of urban ecological source areas in the above embodiments.

[0114] Please refer to Figure 7 , Figure 7 which is a schematic structural diagram of a device for identifying the service range of urban ecological source areas provided by an embodiment of the present invention. The device for identifying the service range of urban ecological source areas includes:

[0115] An ecological source area construction module 701, configured to obtain land use data and ecosystem service value of the area to be identified, and construct an ecological source area set;

[0116] A grid construction module 702, configured to construct a spatial analysis grid in the area to be identified, and the spatial analysis grid includes a plurality of grid cells;

[0117] An impedance evaluation module 703, configured to obtain impedance evaluation data of the grid cell based on the resistance value, elevation data, and slope data corresponding to the grid cell;

[0118] A comprehensive impedance module 704, configured to, for each ecological source area in the ecological source area set, use the impedance evaluation data as the cost, calculate the minimum cumulative cost distance from each grid cell to the ecological source area, and obtain the service impedance evaluation result of the ecological source area;

[0119] A service value calculation module 705, configured to calculate the ecological service value of each ecological source area for the grid cell by using a gravity model according to the ecological system service value and the service impedance evaluation result;

[0120] A maximum service range module 706, configured to divide the grid cell according to the ecological source area corresponding to the maximum value in the ecological service value of each grid cell, so as to obtain the maximum service range of each ecological source area;

[0121] A fine service range module 707, configured to identify the edge area of the ecological source area service based on the ecological service value, and remove the edge area from the maximum service range of the ecological source area to obtain the core area of the ecological source area service.

[0122] Preferably, the ecological source area construction module 701 is specifically configured to:

[0123] Obtain the land use data of the area to be identified, and based on the land use data, obtain a first area set corresponding to several natural land cover types in the area to be identified; delete the areas within the buffer range of the main traffic artery in the first area set to obtain a second area set; merge the areas with the same natural land cover type in the second area set, and select the areas with an area greater than a first preset value in the merged area set to form a first data set;

[0124] Obtain the ecological system service value of each area in the area to be identified, and select the areas with an ecological system service value greater than a second preset value in the area to be identified according to the ecological system service value to form a second data set;

[0125] Take the intersection of the first data set and the second data set as the ecological source area set.

[0126] Preferably, the impedance evaluation module 703 is specifically configured to:

[0127] Obtain the resistance value of the grid cell according to the resistance values of different land use areas preset;

[0128] Obtain the elevation data and slope data corresponding to the grid cell;

[0129] Perform a weighted superposition operation on the resistance value, the elevation data, and the slope data of the grid cell to obtain the impedance evaluation data of the grid cell.

[0130] Preferably, the resistance values of the different land use areas are related to the land use types, urban layouts, and transportation networks of the corresponding land use areas; the land use types include at least one of forest land, water body, grassland, cultivated land, urban construction land, and other unused land; the higher the resistance value of the grid cell, the greater the obstacle to the spread of ecological services in the area corresponding to the grid cell.

[0131] Preferably, the comprehensive impedance module 704 is specifically configured to:

[0132] For each ecological source area where the ecological source areas are concentrated, taking the impedance evaluation data of the grid cells as cost rasters, calculating the minimum value of the sum of the impedance evaluation data on the path from each grid cell passing through the cost surface to reach this ecological source area, obtaining raster data of the minimum cumulative cost distance, and forming a raster layer, where the raster layer is the service impedance evaluation result of this ecological source area;

[0133] Integrating the raster layers of all ecological source areas where the ecological source areas are concentrated to obtain the service impedance evaluation results of all ecological source areas.

[0134] Preferably, the gravity model is specifically:

[0135]

[0136] where F eco is the ecological service value of an ecological source area to a grid cell, V total is the total ecological system service value of the area corresponding to this grid cell, V eco is the ecological system service value of this ecological source area, R cost_dist is the service impedance evaluation result of this ecological source area within the area to be identified.

[0137] Preferably, identifying the marginal area of the ecological source area service based on the ecological service value includes:

[0138] Based on the ecological service value, if the grid cell meets the preset threshold condition, it is determined that the grid cell belongs to the marginal area of the ecological source area service; the formula of the threshold condition is:

[0139]

[0140] where F eco-max and F eco-min are respectively the maximum value and the minimum value in the ecological service values of the ecological source area to the grid cell, F eco-all is the sum of the ecological service values of all ecological source areas to the grid cell, and α is a preset relative difference threshold.

[0141] In specific implementation, the working principle, control process, and achieved technical effects of the urban ecological source area service scope recognition device provided in the embodiments of the present invention are correspondingly the same as those of the urban ecological source area service scope recognition method in the above embodiments, and will not be elaborated herein.

[0142] See Figure 8 , Figure 8 is a structural block diagram of a computer device provided in an embodiment of the present application. The computer device includes: a processor 801, a memory 802, and a computer program stored in the memory 802 and executable on the processor 801. When the processor 801 executes the computer program, it implements the steps in the embodiments of the above urban ecological source area service scope recognition method. Alternatively, when the processor 801 executes the computer program, it implements the functions of each module / unit in the above device embodiments.

[0143] Exemplarily, the computer program can be divided into one or more modules / units. The one or more modules / units are stored in the memory 802 and executed by the processor 801 to complete the present application. The one or more modules / units can be a series of computer program instruction segments capable of completing specific functions, and these instruction segments are used to describe the execution process of the computer program in the computer device.

[0144] The computer device may include, but is not limited to, a processor 801 and a memory 802. Those skilled in the art can understand that the schematic diagram is only an example of the computer device and does not constitute a limitation on the computer device. It may include more or fewer components than shown, or combine certain components, or different components. For example, the computer device may further include input / output devices, network access devices, buses, etc.

[0145] The processor 801 may be a central processing unit (CPU), or may also be other general-purpose processors, digital signal processors (DSPs), application specific integrated circuits (ASICs), field programmable gate arrays (FPGAs), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. The general-purpose processor may be a microprocessor or the processor may also be any conventional processor, etc. The processor 801 is the control center of the computer device, connecting various parts of the entire computer device through various interfaces and lines.

[0146] The memory 802 can be used to store the computer programs and / or modules. By running or executing the computer programs and / or modules stored in the memory 802, and by invoking the data stored in the memory 802, the processor 801 realizes various functions of the computer device. The memory 802 mainly includes a program storage area and a data storage area. Among them, the program storage area can store an operating system, application programs required for at least one function (such as a sound playback function, an image playback function, etc.); the data storage area can store data created according to the use of the mobile phone (such as audio data, phone book, etc.). In addition, the memory 802 can include high-speed random access memory, and can also include non-volatile memory, such as a hard disk, a memory, a plug-in hard disk, a Smart Media Card (SMC), a Secure Digital (SD) card, a Flash Card, at least one magnetic disk storage device, a flash memory device, or other volatile solid-state storage devices.

[0147] Among them, if the modules / units integrated in the computer device are implemented in the form of software function units and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on such an understanding, to implement all or part of the processes in the above embodiment methods of the present application, it can also be completed by a computer program instructing relevant hardware. The computer program can be stored in a computer-readable storage medium. When the computer program is executed by the processor 801, the steps of the above various method embodiments can be implemented. Among them, the computer program includes computer program code, and the computer program code can be in the form of source code, object code, an executable file, or some intermediate form, etc. The computer-readable medium can include: any entity or device capable of carrying the computer program code, a recording medium, a USB flash drive, a mobile hard disk, a magnetic disk, an optical disc, a computer memory, a read-only memory (ROM), a random access memory (RAM), an electrical carrier signal, a telecommunication signal, and a software distribution medium, etc.

[0148] The embodiment of the present invention also provides a computer-readable storage medium. The computer-readable storage medium includes a stored computer program. Among them, when the computer program runs, it controls the device where the computer-readable storage medium is located to execute the urban ecological source area service range identification method described in any one of the above embodiments.

[0149] The embodiments of the present invention provide a method, device, equipment and medium for identifying the service scope of urban ecological source areas, and the beneficial effects are as follows: By combining the ecological patch method and the ecological importance evaluation method for ecological source area identification, high-quality ecological source areas can be identified, improving the accuracy and flexibility of ecological source area identification; By comprehensively considering various factors to evaluate the service impedance of ecological source areas, the comprehensiveness and accuracy of the identification of the service scope of ecological source areas are improved; By constructing a gravity model to calculate the ecological service value, the accuracy of the identification of the service scope of ecological source areas is further improved; By setting a relative difference threshold to divide the edge area and the core area of ecological source areas, the refined division of the service scope of ecological source areas is realized, providing a basis for the precise allocation of urban ecological protection resources and helping to improve resource utilization efficiency. The embodiments of the present invention can not only accurately identify high-value ecological source areas, but also scientifically divide their service scopes, improving the accuracy and comprehensiveness of the identification of the service scope of ecological source areas, and providing strong technical support for urban ecological planning and management.

[0150] The above is the preferred embodiment of the present invention. It should be noted that for those of ordinary skill in the art, without departing from the principle of the present invention, several improvements and refinements can be made, and these improvements and refinements are also regarded as the protection scope of the present invention.

Claims

1. A method for identifying the service scope of an urban ecological source, characterized in that: include: Obtain land use data and ecosystem service value of the area to be identified and construct an ecological source set; Constructing a spatial analysis grid in the area to be identified, wherein the spatial analysis grid includes a plurality of grid cells; Obtaining impedance evaluation data of the grid unit based on the resistance value, elevation data, and slope data corresponding to the grid unit; For each ecological source in the ecological source collection, the impedance evaluation data is used as the cost, and the minimum cumulative cost distance from each grid unit to the ecological source is calculated to obtain the service impedance evaluation result of the ecological source; According to the ecosystem service value and the service impedance evaluation results, the gravity model is used to calculate the ecological service value of each ecological source area to the grid unit; According to the ecological source area corresponding to the maximum value of the ecological service value of each grid unit, the grid unit is divided to obtain the maximum service range of each ecological source area; Based on the ecological service value, identifying the edge area of ​​the ecological source service, and removing the edge area from the maximum service range of the ecological source to obtain the core area of ​​the ecological source service; Wherein, the step of identifying the edge area of ​​ecological source services based on the ecological service value includes: Based on the ecological service value, if the grid unit meets a preset threshold condition, it is determined that the grid unit belongs to the edge area of ​​the ecological source service; the formula of the threshold condition is: Among them, F eco-max With F eco-min are the maximum and minimum values ​​of the ecological service value of the ecological source to the grid unit, respectively, eco-all is the sum of the ecological service values ​​of all ecological sources to the grid unit, and α is the preset relative difference threshold.

2. The method for identifying the service scope of an urban ecological source as claimed in claim 1, characterized in that: The obtaining of land use data and ecosystem service value of the area to be identified and the construction of an ecological source set include: Acquire land use data of the area to be identified, and based on the land use data, acquire a first area set corresponding to a plurality of natural land cover types in the area to be identified; delete areas in the first area set that are within the traffic artery buffer zone to obtain a second area set; merge areas with the same natural land cover type in the second area set, and select areas in the merged area set that have an area greater than a first preset value to form a first data set; Acquire the ecosystem service value of each area in the to-be-identified area, and select, according to the ecosystem service value, areas in the to-be-identified area where the ecosystem service value is greater than a second preset value to form a second data set; The intersection of the first data set and the second data set is used as an ecological source set.

3. The method for identifying the service scope of an urban ecological source as claimed in claim 1, characterized in that: The step of obtaining impedance evaluation data of the grid unit based on the resistance value, elevation data, and slope data corresponding to the grid unit includes: Obtaining the resistance value of the grid unit according to the preset resistance values ​​of different land use areas; Obtaining elevation data and slope data corresponding to the grid unit; A weighted superposition operation is performed on the resistance value, the elevation data and the slope data of the grid unit to obtain impedance evaluation data of the grid unit.

4. The method for identifying the service scope of an urban ecological source as claimed in claim 3, characterized in that: The resistance values ​​of different land use areas are related to the land use type, urban layout and transportation network of the application area; The land use type includes: at least one of forest land, water body, grassland, cultivated land, urban construction land and other unused land; the higher the resistance value of the grid unit, the greater the obstacle to the dissemination of ecological services in the area corresponding to the grid unit.

5. The method for identifying the service scope of an urban ecological source as claimed in claim 1, characterized in that: For each ecological source in the ecological source collection, the impedance evaluation data is used as the cost, and the minimum cumulative cost distance from each grid unit to the ecological source is calculated to obtain the service impedance evaluation result of the ecological source, including: For each ecological source in the ecological source concentration, the impedance evaluation data of the grid unit is used as a cost grid, and the minimum value of the sum of the impedance evaluation data of each grid unit on the path from the cost surface to the ecological source is calculated to obtain the grid data of the minimum cumulative cost distance, and form a grid layer, which is the service impedance evaluation result of the ecological source; The grid layers of all ecological sources in the ecological source set are integrated to obtain the service impedance evaluation results of all ecological sources.

6. The method for identifying the service scope of an urban ecological source as claimed in claim 1, characterized in that: The gravity model is specifically: Among them, F eco is the ecological service value of an ecological source to a grid unit, V total is the sum of the ecosystem service values ​​of the area corresponding to the grid unit, V eco The ecosystem service value of the ecological source, R cost_dist It is the service impedance evaluation result of the ecological source in the area to be identified.

7. A device for identifying the service range of an urban ecological source, characterized in that: include: The ecological source construction module is used to obtain land use data and ecosystem service value of the area to be identified and construct an ecological source set; A grid construction module, used to construct a spatial analysis grid in the area to be identified, wherein the spatial analysis grid includes a plurality of grid units; An impedance evaluation module, used to obtain impedance evaluation data of the grid unit based on the resistance value, elevation data and slope data corresponding to the grid unit; A comprehensive impedance module is used to calculate the minimum cumulative cost distance from each grid unit to the ecological source for each ecological source in the ecological source set, using the impedance evaluation data as the cost, to obtain the service impedance evaluation result of the ecological source; A service value calculation module, used to calculate the ecological service value of each ecological source to the grid unit using a gravity model according to the ecosystem service value and the service impedance evaluation result; A maximum service range module, used to divide the grid unit according to the ecological source corresponding to the maximum value of the ecological service value of each grid unit, to obtain the maximum service range of each ecological source; A refined service range module, for identifying the edge area of ​​the ecological source service based on the ecological service value, and removing the edge area from the maximum service range of the ecological source to obtain the core area of ​​the ecological source service; Wherein, the step of identifying the edge area of ​​ecological source services based on the ecological service value includes: Based on the ecological service value, if the grid unit meets a preset threshold condition, it is determined that the grid unit belongs to the edge area of ​​the ecological source service; the formula of the threshold condition is: Among them, F eco-max With F eco-min are the maximum and minimum values ​​of the ecological service value of the ecological source to the grid unit, respectively, eco-all is the sum of the ecological service values ​​of all ecological sources to the grid unit, and α is the preset relative difference threshold.

8. A computer device, characterized in that: It includes a processor and a memory, wherein a computer program is stored in the memory, and the computer program is configured to be executed by the processor, and when the processor executes the computer program, it implements the urban ecological source service scope identification method as described in any one of claims 1 to 6.

9. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a computer program, wherein when the device where the computer-readable storage medium is located executes the computer program, the method for identifying the service scope of an urban ecological source as described in any one of claims 1 to 6 is implemented.

Images