Urban ecological source service range identification method, device, equipment and medium

Abstract

The invention discloses an urban ecological source and land service range identification method, device and equipment and a medium, and the method comprises the steps: obtaining land utilization data and ecological system service value of a to-be-identified region, and constructing an ecological source and land set; constructing a spatial analysis grid, and obtaining impedance evaluation data based on the resistance value, the elevation data and the slope data of a grid unit; the impedance evaluation data is used as cost, the minimum cumulative cost distance from the grid units to the ecological source land is calculated, and a service impedance evaluation result is obtained; calculating an ecological service value of the ecological source to the grid unit by adopting a gravity model; dividing the maximum service range of the ecological source land according to the ecological source land corresponding to the maximum value in the ecological service values of the grid units; and identifying an edge region of the ecological source ground service, and removing the edge region from the maximum service range of the ecological source ground to obtain a core region of the ecological source ground service. According to the invention, the accuracy and comprehensiveness of urban ecological source service range identification are improved.

Description

Technical Field

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

[0002] Urban ecological sources are the core components of urban ecosystems. Accurately identifying ecological sources and reasonably dividing their service areas can not only optimize land use, provide a scientific basis for urban planning, and ensure the balanced distribution of ecosystem services, but also help promote ecological collaboration across administrative regions and provide guidance for the construction of ecological networks, thereby enhancing the connectivity and stability of the overall ecosystem.

[0003] With the acceleration of urbanization and the increasing complexity of ecological and environmental problems, existing technologies have several limitations in the identification of urban ecological sources and the division of service scopes. The current methods for dividing the service scope of ecological sources are not accurate and comprehensive enough. Although the commonly used buffer zone analysis method is simple and intuitive, it ignores the impact of factors such as terrain and land use on the spatial propagation of ecological services, resulting in low accuracy in the division of service scopes. Another method based on landscape connectivity, although it takes into account spatial heterogeneity, mainly focuses on species migration and is difficult to fully reflect the multi-dimensional service scope of ecological sources. In addition, most of the existing technologies only focus on one aspect of ecological source identification or service scope division, and lack a systematic method to organically combine the two. At the same time, the existing service scope division methods are not sophisticated enough and fail to effectively distinguish the core area and edge area within the service scope of the ecological source. Summary of the invention

[0004] In order to solve the above technical problems, the present invention provides a method, device, equipment and medium for identifying the service scope of an urban ecological source, which combines ecological source identification with service scope division, improves the accuracy and comprehensiveness of ecological source service scope identification, and can identify a more detailed ecological source service scope.

[0005] The embodiment of the present invention provides a method for identifying the service scope of an urban ecological source, including:

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

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

[0008] Obtaining impedance evaluation data of the grid unit based on the resistance value, elevation data, and slope data corresponding to the grid unit;

[0009] 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;

[0010] 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;

[0011] 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;

[0012] Based on the ecological service value, the edge area of ​​the ecological source service is identified, and the edge area is eliminated from the maximum service range of the ecological source to obtain the core area of ​​the ecological source service.

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

[0014] 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 buffer zone of the main traffic artery 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;

[0015] 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;

[0016] The intersection of the first data set and the second data set is used as an ecological source set.

[0017] As an improvement of the above solution, the impedance evaluation data of the grid unit is obtained based on the resistance value, elevation data and slope data corresponding to the grid unit, including:

[0018] Obtaining the resistance value of the grid unit according to the preset resistance values ​​of different land use areas;

[0019] Obtaining elevation data and slope data corresponding to the grid unit;

[0020] 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.

[0021] As an improvement of the above scheme, 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: forest land, water body, grassland, cultivated land, urban construction land and at least one of 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.

[0022] As an improvement of the above scheme, 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:

[0023] 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;

[0024] 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.

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

[0026]

[0027] 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.

[0028] As an improvement of the above solution, the step of identifying the edge area of ​​ecological source services based on the ecological service value includes:

[0029] 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:

[0030]

[0031] Among them, F eco-max With F eco-minare 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.

[0032] The embodiment of the present invention further provides a device for identifying the service range of an urban ecological source, including:

[0033] 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;

[0034] 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;

[0035] 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;

[0036] 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;

[0037] 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;

[0038] 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;

[0039] The refined service range module is used to identify the edge area of ​​the ecological source service based on the ecological service value, and remove the edge area from the maximum service range of the ecological source to obtain the core area of ​​the ecological source service.

[0040] An embodiment of the present invention also provides a computer device, including 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, the method for identifying the service scope of an urban ecological source as described in any one of the above items is implemented.

[0041] An embodiment of the present invention also provides a computer-readable storage medium, which stores a computer program, wherein when the computer program is running, the device where the computer-readable storage medium is located is controlled to execute any of the above-mentioned urban ecological source service scope identification methods.

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

[0043] Figure 1 It is a flow chart of a method for identifying the service scope of an urban ecological source provided by an embodiment of the present invention;

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

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

[0046] Figure 4 It is a schematic diagram of ecological service values ​​of three ecological sources in the area to be identified provided by an embodiment of the present invention;

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

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

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

[0050] Figure 8 It is a structural schematic diagram of a computer device provided by an embodiment of the present invention. DETAILED DESCRIPTION

[0051] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. 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.

[0052] See also Figure 1 , Figure 1 The following is a flow chart of a method for identifying the service scope of an urban ecological source provided by an embodiment of the present invention. The method for identifying the service scope of an urban ecological source includes:

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

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

[0055] 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 buffer zone of the main traffic artery 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;

[0056] 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;

[0057] The intersection of the first data set and the second data set is used as an ecological source set.

[0058] Preferably, the first preset value is the area value that accounts for the top 20% of the total area; the second preset value is the ecosystem service value that accounts for the top 20% of the ecosystem service value.

[0059] Specifically, ecological sources, as areas that provide high-value ecological services, are key elements of urban ecosystems. The embodiment of the present invention uses the ecological patch method and the ecological importance assessment method to identify potential ecological sources, and then obtains the final ecological source set through screening.

[0060] First, the land use data of the city to be identified are used to screen out areas corresponding to natural land cover types such as cultivated land, grassland and water to form the first area set; in order to eliminate human interference, the areas within the 500-meter buffer zone of the main traffic arteries are deleted, where the main traffic arteries include: national highways, provincial highways, county roads, expressways, railways, and urban primary and secondary roads; then, the land uses with the same land use nature are spatially and geographically merged to form a preliminary ecological source data base, and then the areas with the top 20% of the total area are extracted from it to form the first data set.

[0061] Secondly, the identified areas are evaluated by the ecosystem service value index, and the areas with the top 20% of ecosystem service value are extracted to form the second data set. Among them, the ecosystem service value can be calculated by the InVEST model (an integrated evaluation model of ecosystem services and trade-offs) or obtained from the Resource and Environmental Science Data Platform of the Chinese Academy of Sciences. The ecosystem service value index aims to quantify the contribution of each ecological patch to the urban ecosystem.

[0062] Finally, the first and second data sets are spatially overlaid and analyzed, and the intersection of the two is taken to finally form the ecological source identification result, namely the ecological source set, such as Figure 2 The embodiment of the present invention ensures that the identified ecological source area not only has a sufficient area scale, but also has a high ecological value.

[0063] It should be noted that the embodiment of the present invention adopts 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 key areas provide 80% of 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 used mostly rely on single or simple combination of 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, which fail to fully consider the functional value and service capacity of the ecosystem, resulting in an incomplete assessment of the importance of ecological sources. The embodiment of the present invention not only considers the spatial characteristics of ecological patches, such as area, connectivity, etc., but also incorporates the value assessment of ecosystem services, thereby providing a more comprehensive and accurate identification of high ecological value sources and overcoming the limitations of a single perspective.

[0065] S2: constructing a spatial analysis grid in the area to be identified, wherein the spatial analysis grid includes a plurality of grid cells;

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

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

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

[0069] Obtaining the resistance value of the grid unit according to the preset resistance values ​​of different land use areas;

[0070] Obtaining elevation data and slope data corresponding to the grid unit;

[0071] 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.

[0072] As one of the optional embodiments, 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.

[0073] Specifically, the impact of ecological sources on the surrounding environment is a complex process, which is restricted by many factors, including 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 accessibility of transportation network) and distance attenuation. In order to comprehensively evaluate the impact of these impedance factors on ecological source services, the embodiment of the present invention obtains the impedance evaluation data of each grid unit based on the resistance value, elevation data and slope data to form an impedance element evaluation data set, such as Figure 3 shown.

[0074] First, the land use data in the identified area is rasterized and reclassified, where the grid size is consistent with the size of the above grid unit; considering factors such as land use type, urban layout and transportation network, resistance values ​​are assigned to different types of land use, as shown in Table 1. These resistance values ​​reflect the degree of influence of different land types on ecological processes. In this way, the resistance value corresponding to each grid unit 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 conditions and focus of different cities to ensure the accuracy and practicality of the assessment.

[0078] After obtaining the resistance value, the elevation and slope elements are also included in the resistance evaluation. Specifically, the elevation data and slope data of each grid unit are obtained from the land use data of the area to be identified, and the weights corresponding to the resistance value, elevation data, and slope data are assigned to perform weighted superposition operations to calculate the impedance evaluation data and form the impedance element evaluation data set corresponding to the area to be identified.

[0079] S4: for each ecological source in the ecological source set, taking the impedance evaluation data as the cost, calculating the minimum cumulative cost distance from each grid unit to the ecological source, and obtaining the service impedance evaluation result of the ecological source;

[0080] As one of the optional embodiments, for each ecological source in the ecological source set, 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:

[0081] 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;

[0082] 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.

[0083] Specifically, in order to quantify the difficulty of obtaining services from each region to the ecological source, and thus evaluate the impact of impedance factors on the dissemination of ecological services, the embodiment of the present invention uses impedance evaluation data as a cost grid, combined with the cost distance analysis function of GIS (geographic information system), to calculate the minimum value of the sum of impedance evaluation data on the path passed by each grid unit in the identified area from the cost surface to the ecological source, that is, the minimum cumulative cost distance. Specifically, the following operations are performed one by one for each ecological source in the ecological source concentration:

[0084] The ecological source site set is input into GIS as source data;

[0085] The impedance element evaluation dataset corresponding to the area to be identified is used as the dataset of the cost raster;

[0086] The cost distance analysis function of GIS is used to generate raster data representing the minimum cumulative cost distance from the grid cell to the nearest ecological source.

[0087] After iteration, multiple raster layers are generated, each raster layer represents the service impedance distribution of a specific ecological source, and finally these layers are integrated to form the service impedance evaluation results of all ecological sources for services within the identified area. The service impedance evaluation result is a raster dataset. If there are n ecological sources, the service impedance evaluation result contains n raster files. The embodiment of the present invention can carefully reflect the service impedance characteristics of each ecological source, and comprehensively consider the impact of all sources, which helps to fully understand the spatial distribution and acquisition difficulty of urban ecological services.

[0088] S5: 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;

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

[0090] As one of the optional embodiments, the gravity model is specifically:

[0091]

[0092] Among them, F eco is the ecological service value of a certain ecological source to a certain 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.

[0093] Furthermore, through the calculation of the gravity model, the ecological service value of each ecological source in the area to be identified can be obtained, such as Figure 4 As shown, it is a schematic diagram of the ecological service values ​​of three ecological sources in the area to be identified provided by an embodiment of the present invention.

[0094] The embodiment of the present invention takes into account the service capacity of the ecological source itself and incorporates the spatial resistance factor, ultimately quantifying the relative impact of different ecological sources on the same area, which can more accurately reflect the spatial attenuation characteristics of ecological services, and thus achieve more accurate identification of the service scope of the ecological source.

[0095] S6: dividing 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;

[0096] Specifically, based on the calculation results of the gravity model, the principle of "selecting the largest to belong" is adopted for each grid unit, and it is divided into the service range of the ecological source that has the greatest impact on it. Among them, for a certain grid unit, the calculation results of the gravity model can obtain the ecological service value of each ecological source to the grid unit, and the maximum value of these ecological service values ​​is screened out, and the grid unit is divided into the service range of the ecological source corresponding to the maximum value. After all grid units are divided into the corresponding ecological source range one by one, the maximum service range of each ecological source can be obtained, such as Figure 5 shown.

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

[0098] Specifically, after obtaining the maximum service range of the ecological source, the situation in which each grid unit is simultaneously affected by multiple ecological sources is analyzed to identify the potential vulnerable hinterland, i.e., the marginal area, in the ecological source.

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

[0100] 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:

[0101]

[0102] 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.

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

[0104] Specifically, when a grid cell satisfies When , it is considered that the grid unit is similarly affected by multiple ecological sources and is not clearly attributed to any ecological source. It is defined as a potential vulnerable hinterland, that is, a 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 cities with more fragile ecosystems, the threshold can be lowered to more strictly define the stable core area.

[0106] Furthermore, the potentially vulnerable hinterland is removed from the maximum service range of each ecological source, and the remaining area is the stable core hinterland of the ecological source, that is, the core area. The core area is dominated by a specific ecological source, and the ecological services are relatively stable and sufficient. Figure 6 As shown, it is a schematic diagram of the stable core hinterland and the potentially vulnerable hinterland of the ecological source provided by the embodiment of the present invention. The stable core hinterland represents the main beneficiary area of ​​the ecological source service and should be protected and maintained; while the potentially vulnerable hinterland is an area where the ecological service is relatively weak 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 service by quantitatively analyzing the service impact range of the ecological source, which not only clarifies the core impact area of ​​the ecological source, but also identifies the transition zone where the service is relatively weak, providing more targeted guidance for urban ecological management.

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

[0108] In terms of urban master planning and green space system optimization, this invention can be used to identify and protect key ecological sources, assist in determining ecological protection red lines, and rationally layout urban construction land and ecological land. By analyzing the service scope of ecological sources, planners can identify areas with insufficient ecological service coverage and plan new parks and green spaces in a targeted manner, thereby optimizing green space layout, improving the spatial balance of ecological services, and creating a more livable environment for citizens.

[0109] In the planning and implementation of ecological restoration projects, this invention can help determine priority restoration areas. By identifying potential vulnerable hinterlands, city managers can develop targeted restoration plans, such as establishing ecological corridors in vulnerable areas between ecological sources, or increasing green space in marginal areas, which will help 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 impact of large-scale construction projects on surrounding ecological sources and their service areas, which helps to formulate more accurate environmental protection measures and minimize the negative impact of projects on important ecological areas. At the same time, this analysis can provide a scientific basis for project site selection and scale decision-making, helping to balance development needs and ecological protection.

[0111] In terms of cross-administrative ecological management and coordination, this invention can help identify key ecological sources and their service areas across multiple administrative regions, and provide a scientific basis for ecological collaborative management. This will help break down administrative barriers, promote the overall coordination of regional ecological governance, and provide support for the establishment of a more equitable and effective cross-regional ecological compensation mechanism.

[0112] In urban resilience assessment and adaptive planning, this invention can be used to analyze the distribution and function of key ecological sources and help identify ecologically vulnerable areas. This information provides an 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 sources, urban development strategies can be adjusted in a timely manner, the adaptability and resilience of urban ecosystems can be enhanced, and the overall resilience of the city can be improved.

[0113] Correspondingly, the present invention also provides a device for identifying the service scope of an urban ecological source, which can implement all the processes of the method for identifying the service scope of an urban ecological source in the above-mentioned embodiment.

[0114] See also Figure 7 , Figure 7 1 is a schematic diagram of a structure of a device for identifying a service range of an urban ecological source provided by an embodiment of the present invention. The device for identifying a service range of an urban ecological source includes:

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

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

[0117] An impedance evaluation module 703, 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;

[0118] The comprehensive impedance module 704 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;

[0119] A service value calculation module 705 is 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;

[0120] The maximum service range module 706 is 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;

[0121] The refined service range module 707 is used to identify the edge area of ​​the ecological source service based on the ecological service value, and remove the edge area from the maximum service range of the ecological source to obtain the core area of ​​the ecological source service.

[0122] Preferably, the ecological source construction module 701 is specifically used for:

[0123] 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 buffer zone of the main traffic artery 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;

[0124] 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;

[0125] The intersection of the first data set and the second data set is used as an ecological source set.

[0126] Preferably, the impedance evaluation module 703 is specifically used for:

[0127] Obtaining the resistance value of the grid unit according to the preset resistance values ​​of different land use areas;

[0128] Obtaining elevation data and slope data corresponding to the grid unit;

[0129] 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.

[0130] Preferably, 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.

[0131] Preferably, the comprehensive impedance module 704 is specifically used for:

[0132] 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;

[0133] 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.

[0134] Preferably, the gravity model is specifically:

[0135]

[0136] 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.

[0137] Preferably, the step of identifying the edge area of ​​ecological source services based on the ecological service value includes:

[0138] 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:

[0139]

[0140] 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.

[0141] In the specific implementation, the working principle, control process and technical effect achieved by the urban ecological source service scope identification device provided in the embodiment of the present invention are the same as those of the urban ecological source service scope identification method in the above-mentioned embodiment, and will not be repeated here.

[0142] See also Figure 8 , Figure 8 801 is a structural block diagram of a computer device provided in an embodiment of the present application, and 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, the steps in the above-mentioned urban ecological source service range identification method embodiment are implemented. Alternatively, when the processor 801 executes the computer program, the functions of each module / unit in the above-mentioned device embodiments are implemented.

[0143] Exemplarily, the computer program may be divided into one or more modules / units, which are stored in the memory 802 and executed by the processor 801 to complete the present application. The one or more modules / units may be a series of computer program instruction segments capable of completing specific functions, which 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 will appreciate that the schematic diagram is merely an example of a computer device and does not constitute a limitation on the computer device. The computer device may include more or fewer components than shown in the diagram, or may combine certain components, or different components. For example, the computer device may also include an input / output device, a network access device, a bus, etc.

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

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

[0147] Wherein, if the module / unit integrated in the computer device is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the present application implements all or part of the processes in the above-mentioned embodiment method, and can also be completed by instructing the relevant hardware through a computer program. The computer program can be stored in a computer-readable storage medium, and the computer program can implement the steps of the above-mentioned various method embodiments when executed by the processor 801. Wherein, the computer program includes computer program code, and the computer program code can be in source code form, object code form, executable file or some intermediate form. The computer-readable medium may include: any entity or device capable of carrying the computer program code, recording medium, U disk, mobile hard disk, disk, optical disk, computer memory, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), electric carrier signal, telecommunication signal and software distribution medium, etc.

[0148] An embodiment of the present invention also provides a computer-readable storage medium, which includes a stored computer program, wherein when the computer program is running, the device where the computer-readable storage medium is located is controlled to execute the urban ecological source service scope identification method described in any of the above embodiments.

[0149] The embodiment of the present invention provides a method, device, equipment and medium for identifying the service scope of an urban ecological source, and its beneficial effects are: by combining the ecological patch method and the ecological importance evaluation method to identify the ecological source, it is possible to identify high-quality ecological sources, thereby improving the accuracy and flexibility of ecological source identification; by comprehensively evaluating the service impedance of the ecological source based on multiple factors, the comprehensiveness and accuracy of the ecological source service scope identification are improved; by constructing a gravity model to calculate the ecological service value, the accuracy of the ecological source service scope identification is further improved; by setting a relative difference threshold to divide the edge area and the core area of ​​the ecological source, the refined division of the ecological source service scope is achieved, which provides a basis for the precise allocation of urban ecological protection resources and helps to improve resource utilization efficiency. The embodiment of the present invention can not only accurately identify high-value ecological sources, but also scientifically divide their service scopes, improve the accuracy and comprehensiveness of the ecological source service scope identification, and provide strong technical support for urban ecological planning and management.

[0150] The above is a preferred embodiment of the present invention. It should be pointed out that a person skilled in the art can make several improvements and modifications without departing from the principle of the present invention. These improvements and modifications are also considered to be within the scope of protection 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, the edge area of ​​the ecological source service is identified, and the edge area is removed from the maximum service range of the ecological source to obtain the core area of ​​the ecological source service.

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 buffer zone of the main traffic artery 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 whose area is 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. The method for identifying the service scope of an urban ecological source as claimed in claim 1, characterized in that: 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 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 to-be-identified area, 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; The refined service range module is used to identify the edge area of ​​the ecological source service based on the ecological service value, and remove the edge area from the maximum service range of the ecological source to obtain the core area of ​​the ecological source service.

9. 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 7.

10. 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 7 is implemented.

Images