An ecological safety pattern key area identification method, device, medium and product

Abstract

The application discloses an ecological security pattern key area identification method, equipment, medium and product, relates to the field of ecological security pattern construction, and comprises the following steps: determining an ecological source according to geographical space data, adopting a morphological space pattern identification technology, an area threshold and connectivity analysis; generating a comprehensive ecological resistance surface of the ecological source; determining an ecological corridor according to the ecological source and the comprehensive ecological resistance surface, adopting a minimum cost path model; constructing an ecological security pattern network graph based on the ecological source and the ecological corridor; performing ecological security pattern network topology feature analysis according to the ecological security pattern network graph; determining a comprehensive ecological pressure index according to an ecological pressure index factor; and determining ecological security pattern partitioning according to the network topology structure and the comprehensive ecological pressure index. The application can fully consider the interaction and connection of the ecological security pattern at a local scale.

Description

Technical Field

[0001] This application relates to the field of ecological security pattern construction, and in particular to a method, device, medium and product for identifying key areas of ecological security pattern. Background Technology

[0002] Constructing an ecological security pattern is an effective measure to analyze the structure and function of regional ecosystems and maintain ecosystem services. Existing technologies mainly focus on the identification of source areas and the extraction of corridors to construct ecological security networks, paying more attention to the identification of geographical entities, but lacking analysis of the strength and direction of functional connections between elements. This results in limitations in the spatial pattern of the ecological security pattern and the coupling and feedback of various ecological elements in the network. At the same time, the ecological security patterns currently constructed are usually regarded as planar structures, ignoring their inherent hierarchical differentiation (such as the functional gradient between the core protected area and the human activity area and the buffer zone between the core protected area and the human activity area). Summary of the Invention

[0003] The purpose of this application is to provide a method, device, medium and product for identifying key areas of ecological security patterns, to realize spatial analysis of regional ecological security patterns, and to accurately identify core ecological sources and ecological corridors that play a key role in the stability and connectivity of ecological security patterns, while fully considering the interaction and connection of ecological security patterns at the local scale.

[0004] To achieve the above objectives, this application provides the following solution:

[0005] Firstly, this application provides a method for identifying key areas in an ecological security pattern, the method comprising:

[0006] Obtain geospatial data of the area to be identified; the geospatial data includes: land use data, slope, elevation, water network, and road network;

[0007] Based on geospatial data, morphological spatial pattern identification technology, area thresholds, and connectivity analysis were used to determine the ecological source areas;

[0008] Based on geospatial data, resistance factor indicators are determined; and based on these indicators, the weights of the resistance factor indicators are determined using the analytic hierarchy process (AHP); then, a weighted superposition is performed to generate a comprehensive ecological resistance surface for the ecological source area; the resistance factor indicators include: land use type, slope, distance from railways and highways, elevation, and distance from water bodies;

[0009] Based on the ecological source area and the comprehensive ecological resistance surface, the minimum cost path model is used to determine the ecological corridor;

[0010] An ecological security pattern network diagram is constructed based on ecological source areas and ecological corridors; the ecological security pattern network diagram uses ecological source areas as network nodes and ecological corridors as edges connecting network nodes.

[0011] The ecological security pattern network topology characteristics were analyzed based on the ecological security pattern network diagram, and the analysis results were obtained; the analysis results include the community structure and the core nodes and edge nodes within the community structure;

[0012] Ecological pressure index factors of ecological source areas are obtained, and a comprehensive ecological pressure index of ecological source areas is determined based on ecological pressure index factors; the ecological pressure index factors include: soil erosion, population concentration pressure factor and land use intensity factor.

[0013] Based on the network topology of the ecological security pattern network diagram and the comprehensive ecological pressure index, the ecological security pattern zoning of the areas to be identified is determined.

[0014] Optionally, the step of determining the ecological source area based on geospatial data using morphological spatial pattern recognition technology, area thresholds, and connectivity analysis specifically includes:

[0015] The types of ecological source areas are determined based on land use data; these types include: forest land, grassland, shrubland, water bodies, and wetlands.

[0016] Based on the type of ecological source area, morphological spatial pattern identification technology is used to determine the core area patches in the ecological space;

[0017] Core patches are obtained by screening core areas in the ecological space based on area thresholds;

[0018] Connectivity analysis was performed on the core patches to determine the patch importance index;

[0019] Ecological source areas are determined based on core patches with a patch importance index greater than 1.

[0020] Optionally, the connectivity analysis of the core plaques to determine the plaque importance index specifically includes:

[0021] Using formula Determine the plaque importance index dPC;

[0022] Where PC is the potential connectivity index. , where n is the total number of core plaques. and A represents the area of ​​the i-th and j-th core patches, respectively. L P represents the total area of ​​the region to be identified. ijPCremove represents the maximum product probability of all paths between the i-th and j-th core patches, and PCremove is the probability connectivity index of the remaining core patches after removing a single core patch.

[0023] Optionally, the step of determining ecological corridors based on the ecological source area and the comprehensive ecological resistance surface using a minimum cost path model specifically includes:

[0024] Using formula Determine the minimum cost path model;

[0025] Where LCP is the final minimum cost path generated from the j-th ecological source to the i-th ecological source, and D ij R represents the spatial distance traversed from the j-th ecological source area to the i-th ecological source area. i Let represent the ecological resistance value of the i-th ecological source site in a certain direction in space, where m and n refer to the number of rows and columns traversed by the i-th to j-th ecological source sites in the ecological resistance raster data, where m is the number of rows and n is the number of columns. This is a function that takes the minimum value.

[0026] Optionally, the step of performing network topology analysis of the ecological security pattern based on the ecological security pattern network diagram to obtain the analysis results specifically includes:

[0027] Based on the ecological security pattern network diagram, the GN algorithm is used to determine the community structure;

[0028] Core-periphery structure analysis is performed on each community structure in the ecological security pattern network diagram to obtain the core nodes and peripheral nodes within the community structure.

[0029] Optionally, the step of obtaining the ecological pressure index factors of the ecological source area and determining the comprehensive ecological pressure index of the ecological source area based on the ecological pressure index factors specifically includes:

[0030] Using formula Determine the amount of soil erosion; among which, For the soil erosion modulus of the grid, For vegetation coverage, The rainfall erosivity index, The soil erodibility index. Factors related to soil and water conservation measures;

[0031] Using formula Determine population agglomeration pressure factors Where pd is the population density;

[0032] Using formula Determine land use intensity factors ;in, Let A be the area of ​​the i-th land use type in the grid, and let A be the area of ​​a single grid cell. Assign a value to the land use intensity of the i-th land use type;

[0033] Using formula Determine the comprehensive ecological pressure index (EP) of the ecological source area; among which, This is the normalized value of soil erosion. This is the normalized value of the population agglomeration pressure factor. This is the normalized value of the land use intensity factor.

[0034] Optionally, determining the ecological security pattern zoning of the area to be identified based on the network topology of the ecological security pattern network map and the comprehensive ecological pressure index specifically includes:

[0035] Based on the comprehensive ecological pressure index, the natural discontinuity method is used to classify the ecological pressure levels of ecological source areas.

[0036] Based on the network topology of the ecological security pattern network diagram and the comprehensive ecological pressure index, the ecological source areas are divided into multiple ecological security clusters according to the ecological security pattern zoning principle.

[0037] Within each ecological security cluster, an ecological security core area and an ecological security control area are defined.

[0038] Based on the ecological pressure levels of the core ecological security area, the ecological security control area, and the ecological source area, ecologically fragile areas and ecologically sensitive areas are identified.

[0039] Secondly, this application provides a computer device, including: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the method for identifying key areas of ecological security patterns.

[0040] Thirdly, this application provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the aforementioned method for identifying key areas of ecological security patterns.

[0041] Fourthly, this application provides a computer program product, including a computer program that, when executed by a processor, implements the aforementioned method for identifying key areas of ecological security patterns.

[0042] According to the specific embodiments provided in this application, this application has the following technical effects:

[0043] This application provides a method, device, medium, and product for identifying key areas of ecological security patterns. By deeply analyzing the interactions and connections between ecological source areas and ecological corridors, it constructs an ecological security pattern network map that reflects the structure and interrelationships of regional ecological security patterns at different scales. Through topological features and ecological pressure analysis of the ecological security pattern network map, spatial analysis of the regional ecological security pattern can be achieved. This application simplifies the regional ecological security network into a binary structure by analyzing the community structure and the core and peripheral nodes within the community structure. It reveals the spatial self-organization characteristics of ecological processes and the hierarchical characteristics of elements within the network, helping to quantify the differences in functional importance among nodes (i.e., ecological source areas) from the perspective of maintaining the structural stability of the ecological security pattern network map, and determining the interactions and connections of ecological security patterns at local scales. Furthermore, it effectively tracks the impact of human activities on the regional ecological security pattern based on the comprehensive ecological pressure index of ecological source areas. By combining the impact of human activities with the characteristics of the network's own structure and its dynamic evolution over many years, this application can accurately identify core ecological sources and ecological corridors that play a key role in the stability and connectivity of the ecological security pattern, providing targets for differentiated ecological protection and restoration management measures, and effectively revealing the internal mechanisms of ecological network evolution or degradation. Attached Figure Description

[0044] To more clearly illustrate the technical solutions in the embodiments of this application or the prior art, the drawings used in the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of this application. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0045] Figure 1 This is a schematic diagram of a method for identifying key areas of ecological security patterns in one embodiment of this application;

[0046] Figure 2 A schematic diagram of the community structure determined by the GN (Girvan-Newman) algorithm (the circles at the bottom represent nodes in the ecological security pattern network). Detailed Implementation

[0047] The technical solutions of the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of this application, and not all embodiments. Based on the embodiments of this application, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of this application.

[0048] To make the above-mentioned objectives, features and advantages of this application more apparent and understandable, the application will be further described in detail below with reference to the accompanying drawings and specific embodiments.

[0049] In one exemplary embodiment, such as Figure 1 As shown, a method for identifying key areas in an ecological security pattern is provided, comprising the following steps S101 to S108. Wherein:

[0050] S101, Obtain geospatial data of the area to be identified; the geospatial data includes: land use data, slope, elevation, water network, and road network;

[0051] S102. Based on geospatial data, morphological spatial pattern identification technology, area threshold and connectivity analysis are used to determine the ecological source area;

[0052] S102 specifically includes:

[0053] S21, Determine the ecological source land type based on land use data; the ecological source land type includes: forest land, grassland, shrubland, water body and wetland;

[0054] S22. Based on the type of ecological source area, the morphological spatial pattern identification (MSPA) technique is used to identify core area patches in the ecological space; core area patches have the potential to become ecological source areas.

[0055] S23, core patches in the ecological space are screened based on the area threshold to obtain core patches; the area threshold can be 10 square kilometers.

[0056] S24, perform connectivity analysis on the core patches to determine the patch importance index;

[0057] S24 specifically includes:

[0058] Using formula Determine the plaque importance index dPC;

[0059] Where PC is the potential connectivity index. , where n is the total number of core plaques. and A represents the area of ​​the i-th and j-th core patches, respectively. L P represents the total area of ​​the region to be identified. ij PCremove represents the maximum product probability of all paths between the i-th and j-th core patches, and PCremove is the probability connectivity index of the remaining core patches after removing a single core patch.

[0060] S25, determine the ecological source area based on core patches with a patch importance index greater than 1.

[0061] S103, determine the resistance factor indicators based on geospatial data; and based on the resistance factor indicators, use the Analytic Hierarchy Process (AHP) to determine the weights of the resistance factor indicators; then perform weighted superposition to generate a comprehensive ecological resistance surface of the ecological source area; the resistance factor indicators include: land use type, slope, distance from railway and highway, elevation, and distance from water body;

[0062] That is, the comprehensive ecological resistance surface of the ecological source area is obtained by superimposing the various resistance factor indicators based on their weights.

[0063] S104. Based on the ecological source area and the comprehensive ecological resistance surface, the minimum cost path model is used to determine the ecological corridor;

[0064] S104 specifically includes:

[0065] Using formula Determine the minimum cost path model;

[0066] Where LCP is the final minimum cost path generated from the j-th ecological source to the i-th ecological source, and D ij R represents the spatial distance traversed from the j-th ecological source area to the i-th ecological source area. i Let represent the ecological resistance value of the i-th ecological source site in a certain direction in space, where m and n refer to the number of rows and columns traversed by the i-th to j-th ecological source sites in the ecological resistance raster data, where m is the number of rows and n is the number of columns. This is a function that takes the minimum value.

[0067] S105, Construct an ecological security pattern network diagram based on ecological source areas and ecological corridors; The ecological security pattern network diagram uses ecological source areas as network nodes and ecological corridors as edges connecting network nodes.

[0068] S106, Based on the ecological security pattern network diagram, perform network topology feature analysis of the ecological security pattern to obtain analysis results; the analysis results include the community structure and the core nodes and edge nodes within the community structure;

[0069] S106 specifically includes:

[0070] S61. Based on the ecological security pattern network diagram, the GN algorithm is used to determine the community structure;

[0071] The process of determining the structure of a community is as follows:

[0072] (1) Construct an undirected and unweighted ecological security pattern network graph G=(N,L) based on the ecological source areas and ecological corridors. Ecological source areas are abstracted as network nodes, and the ecological corridors connecting ecological source areas are abstracted as network edges. When there is an ecological corridor connecting any two ecological source areas, the edge value is set to 1; when there is no ecological corridor connecting them, the edge value is set to 0. N is the set of abstracted network nodes in the complex network model G. ,in, Let represent the i-th network node, and L represent the abstract set of network edges. The number of connections between each node and other network nodes satisfies ... .

[0073] (2) Detecting community structure in a graph based on the GN (Girvan-Newman) algorithm. Community structure reflects the clustering characteristics of a network. A community is a subgraph of a graph, where nodes within a community are tightly connected, while connections to external nodes are sparse. Select any node in the graph sequentially, and calculate the Edge Betweenness Centrality (EBC) of all edges based on the shortest path. Edge Betweenness Centrality is defined as the number of shortest paths through an edge in a network, and is an indicator of the importance of edges in the network. For each node in the graph, repeat the above steps, and sum the n EBC values ​​of each edge and divide by 2 to obtain the final EBC value of each edge. Remove the edge with the highest EBC value and calculate the modularity Q of the network community structure. Iterate the above process until no edge can be removed. The final result is represented as a hierarchical clustering tree, and as shown in... Figure 2 As shown. Each independent network structure corresponds to a horizontal division of the dendrogram, in order to... Figure 2 Taking the horizontal line as an example, each branch intersecting the horizontal line represents a community composed of a group of nodes, and all nodes at the bottom of the branch are members of that community. Specifically, the core structure is detected using the `cluster_edge_betweenness()` function in the R language.

[0074] (3) Using the maximum modularity Q-value as the clustering objective, determine the appropriate number of communities m ( Modularity is a method for measuring the strength of a network community structure; a higher value indicates a better community partitioning effect. Modularity is defined as the ratio of the total number of edges within a community to the total number of edges in the network, minus an expected value. This expected value represents the ratio of the total number of edges within a community to the total number of edges in the network when the network is set to random.

[0075] ;

[0076] In the formula, Points The degree. When the point When there are edges connecting them, ,otherwise Let m be the total number of edges in the network. Assume... Representing points respectively The two clubs they belonged to, when When in a club ,on the contrary, 0. The specific calculations can be performed using the `modularity` package in Python.

[0077] S62 performs a core-periphery structure analysis on each community structure in the ecological security pattern network diagram to obtain the core nodes and peripheral nodes within the community structure.

[0078] Core-edge structure analysis is performed on the community structures within the network to identify network nodes in a core position. A core-edge structure is a special structure consisting of a group of interconnected nodes, characterized by a tightly connected center and sparsely distributed edges, with all edge points connected to the core points. Based on the community structure partitioning results, a symmetric binary network containing nodes and edges is constructed for each community structure. . For community i-network A collection of network nodes , , This represents the set of edges within community structure i. Input this into the UCINET social network analysis software. The matrix performs core-periphery structure analysis on each community structure according to the menu path network (N) > core / periphery (Y) > absolute (C), and outputs the nodes contained in the core and periphery of each community structure.

[0079] S107, Obtain the ecological pressure index factors of the ecological source area, and determine the comprehensive ecological pressure index of the ecological source area based on the ecological pressure index factors; the ecological pressure index factors include: soil erosion, population concentration pressure factor and land use intensity factor; wherein, the ecological pressure index factors are determined based on raster data and station data with a spatial resolution of 1km for climate, population and land use.

[0080] S107 specifically includes:

[0081] Using formula Soil erosion was determined based on the fundamental principles of the general soil and water conservation equation, selecting indicators such as precipitation erosiveness, soil erodibility, slope and slope length, and surface vegetation cover. For the soil erosion modulus of the grid, For vegetation coverage, The rainfall erosivity index, The soil erodibility index. Factors related to soil and water conservation measures;

[0082] in, , It is the normalized vegetation index; Information contributed by a completely vegetated surface. Information contributed by surfaces with no vegetation cover and Take the maximum and minimum values ​​of the NDVI raster data respectively (remove negative values).

[0083] ;

[0084] ;

[0085] In the formula, R represents the annual average annual rainfall erosion force; Let be the erosivity of rainfall in the j-th half-month; j is the 24 half-months of a year; i is the number of days with erosive rainfall in the j-th half-month, i = 1, 2, ..., m; α is a parameter, α = 0.3937 in the warm season and α = 0.3101 in the cold season. The rainfall on the i-th erosive day of the j-th half-month was obtained by Kriging interpolation using daily rainfall observation data from meteorological stations. The specific calculation method was performed using ArcGIS software, and the spatialization of rainfall was completed according to the menu path Geostatistical Analyst > Kriging.

[0086] ;

[0087] ;

[0088] In the formula, K represents the soil erodibility factor before correction, and K represents the soil erodibility factor after correction. , , and These represent the percentage content of clay, silt, sand, and organic carbon, respectively. Soil texture percentages can be extracted from the soil data attribute items.

[0089] Assuming that human pressure on natural ecosystems increases logarithmically with increasing population density, and that pressure increases at a rate of 1000 people / km²... 2 When saturation is reached. Therefore, a population density greater than 1000 people / km² is considered suitable. 2 The grid population density pressure factor is assigned a value of 10, and the population density is less than 1000 people / km². 2 The grid population density pressure factor, using the formula To determine; where pd is the population density;

[0090] Land use intensity was divided into five levels and assigned values: unused land was assigned 1; forest and grassland 2; water area 3; cultivated land 4; and construction land 5. The formula was then used. Determine land use intensity factors ;in, Let A be the area of ​​the i-th land use type in the grid, and let A be the area of ​​a single grid cell. Assign a value to the land use intensity of the i-th land use type;

[0091] When the three types of factors contribute equally to the ecological pressure index, the three types of factors are normalized, and the comprehensive ecological pressure index score of each grid is the sum of the scores of the three types of indicators:

[0092] ;

[0093] in, This is the normalized value of soil erosion. This is the normalized value of the population agglomeration pressure factor. This is the normalized value of the land use intensity factor.

[0094] S108. Based on the network topology of the ecological security pattern network map and the comprehensive ecological pressure index, determine the ecological security pattern zoning of the area to be identified.

[0095] S108 specifically includes:

[0096] S81. Based on the comprehensive ecological pressure index, the natural discontinuity method is used to classify the ecological pressure level of the ecological source area, thereby obtaining the ecological pressure level of the ecological source area; the ecological pressure level of the ecological source area includes: low ecological pressure area, general ecological pressure area and high ecological pressure area.

[0097] S82. Based on the network topology and comprehensive ecological pressure index of the ecological security pattern network diagram, the ecological source areas are divided into multiple ecological security clusters, i.e., first-level zoning, using the ecological security pattern zoning principle. The ecological security pattern zoning principle includes: consistency and heterogeneity of ecological spatial pattern; high correlation and integrity of ecological network system.

[0098] The consistency and heterogeneity of ecological spatial patterns are that, spatially, each community structure, ecological processes, ecological functions, and ecological structures are relatively consistent, while there are significant differences in ecological structures and functions among different communities.

[0099] The high degree of interconnectedness and integrity of ecological networks are the core of dynamic evolution and mutual feedback mechanisms. Changes in the information state of a node in the network can be propagated and flow through the structure of interactions between nodes, thereby causing changes in the entire system.

[0100] S83 divides each ecological security cluster into an ecological security core area and an ecological security control area, i.e., a two-level zoning.

[0101] Based on the internal organizational patterns of the ecological security network, the primary and secondary partitioning identifies ecological security clusters with macroscopic spatial characteristics according to the local cohesion characteristics of network communities. Furthermore, based on core-periphery structural features, it identifies the ecological security core area and ecological security control area within each cluster. An ecological security cluster is a collaborative unit formed by multiple spatially adjacent and internally connected ecological elements, including ecological source areas and ecological corridors, and serves as a functional module of the regional ecological security pattern. The ecological security core area consists of nodes located in the core structure of the network community and the connections between these nodes. It is characterized by high connectivity and strong influence, and can play a global control role, serving as a key unit for maintaining the overall stability, connectivity, and function of the ecological security cluster. The ecological security control area consists of nodes located in the peripheral structure of the network community and the connections within the peripheral structure, as well as the connections between the core and peripheral structures. It is mostly a transition zone, structurally dependent on the core area, and its function is mainly buffering and filtering, making it susceptible to external disturbances.

[0102] S84, based on the ecological pressure levels of the core ecological security area, the ecological security control area, and the ecological source area, determines the ecologically fragile area and the ecologically sensitive area, i.e., the three-level zoning.

[0103] The third-level zoning, based on the second-level zoning, further divides ecologically vulnerable areas and ecologically sensitive areas by considering the degree of ecological pressure impact and spatial distribution, and using ecological source areas as the basis. When an ecological source area exists within a second-level zoning, and more than 30% of its area is located in an area of ​​high ecological pressure, that ecological source area is designated as an ecologically vulnerable area. When ecologically vulnerable areas are distributed within the core ecological security area, other ecological source areas within the core area, as well as ecological source areas within the ecological security control area of ​​that ecological security cluster, are all designated as ecologically sensitive areas. When ecologically vulnerable areas are distributed within the ecological security control area, all ecological source areas within the core ecological security area of ​​that ecological security cluster are designated as ecologically sensitive areas.

[0104] In one exemplary embodiment, a computer device is provided, which may be a server or a terminal. The computer device includes a processor, memory, input / output interfaces (I / O), and a communication interface. The processor, memory, and I / O interfaces are connected via a system bus, and the communication interface is connected to the system bus via the I / O interfaces. The processor of the computer device provides computing and control capabilities. The memory of the computer device includes non-volatile storage media and internal memory. The non-volatile storage media stores an operating system, computer programs, and a database. The internal memory provides an environment for the operation of the operating system and computer programs in the non-volatile storage media. The I / O interfaces of the computer device are used for exchanging information between the processor and external devices. The communication interface of the computer device is used for communication with external terminals via a network connection. When the computer program is executed by the processor, it implements a method for identifying key areas of an ecological security pattern.

[0105] In one exemplary embodiment, a computer device is provided, including a memory and a processor, wherein the memory stores a computer program, and the processor executes the computer program to implement the steps in the above-described method embodiments.

[0106] In one exemplary embodiment, a computer-readable storage medium is provided storing a computer program that, when executed by a processor, implements the steps in the above-described method embodiments.

[0107] In one exemplary embodiment, a computer program product is provided, including a computer program that, when executed by a processor, implements the steps in the above-described method embodiments.

[0108] It should be noted that the user information (including but not limited to user device information, user personal information, etc.) and data (including but not limited to data used for analysis, data stored, data displayed, etc.) involved in this application are all information and data authorized by the user or fully authorized by all parties, and the collection, use and processing of the relevant data must comply with relevant regulations.

[0109] Those skilled in the art will understand that all or part of the processes in the above embodiments can be implemented by a computer program instructing related hardware. The computer program can be stored in a non-volatile computer-readable storage medium. When executed, the computer program can include the processes of the embodiments described above. Any references to memory, databases, or other media used in the embodiments provided in this application can include at least one of non-volatile and volatile memory. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory, optical memory, high-density embedded non-volatile memory, resistive random access memory (ReRAM), magnetic random access memory (MRAM), ferroelectric random access memory (FRAM), phase change memory (PCM), graphene memory, etc. Volatile memory can include random access memory (RAM) or external cache memory, etc. By way of illustration and not limitation, RAM can take many forms, such as Static Random Access Memory (SRAM) or Dynamic Random Access Memory (DRAM).

[0110] The databases involved in the embodiments provided in this application may include at least one type of relational database and non-relational database. Non-relational databases may include, but are not limited to, blockchain-based distributed databases. The processors involved in the embodiments provided in this application may be general-purpose processors, central processing units, graphics processing units, digital signal processors, programmable logic devices, quantum computing-based data processing logic devices, etc., and are not limited to these.

[0111] In this application, all actions to acquire signals, information, or data are carried out in compliance with the relevant data protection laws and policies of the country where the location is situated, and with the authorization granted by the owner of the relevant device.

[0112] The technical features of the above embodiments can be combined in any way. For the sake of brevity, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0113] This document uses specific examples to illustrate the principles and implementation methods of this application. The descriptions of the above embodiments are only for the purpose of helping to understand the methods and core ideas of this application. Furthermore, those skilled in the art will recognize that, based on the ideas of this application, there will be changes in the specific implementation methods and application scope. Therefore, the content of this specification should not be construed as a limitation of this application.

Claims

1. A method for identifying key areas in an ecological security pattern, characterized in that, The method for identifying key areas of ecological security patterns includes: Obtain geospatial data of the area to be identified; the geospatial data includes: land use data, slope, elevation, water network, and road network; Based on geospatial data, morphological spatial pattern identification technology, area thresholds, and connectivity analysis were used to determine the ecological source areas; Based on geospatial data, resistance factor indicators are determined; and based on these indicators, the weights of the resistance factor indicators are determined using the analytic hierarchy process (AHP); then, a weighted superposition is performed to generate a comprehensive ecological resistance surface for the ecological source area; the resistance factor indicators include: land use type, slope, distance from railways and highways, elevation, and distance from water bodies; Based on the ecological source area and the comprehensive ecological resistance surface, the minimum cost path model is used to determine the ecological corridor; An ecological security pattern network diagram is constructed based on ecological source areas and ecological corridors; the ecological security pattern network diagram uses ecological source areas as network nodes and ecological corridors as edges connecting network nodes. The ecological security pattern network topology characteristics were analyzed based on the ecological security pattern network diagram, and the analysis results were obtained; the analysis results include the community structure and the core nodes and edge nodes within the community structure; Ecological pressure index factors of ecological source areas are obtained, and a comprehensive ecological pressure index of ecological source areas is determined based on ecological pressure index factors; the ecological pressure index factors include: soil erosion, population concentration pressure factor and land use intensity factor. Based on the network topology of the ecological security pattern network diagram and the comprehensive ecological pressure index, the ecological security pattern zoning of the area to be identified is determined; The process of determining the ecological source area based on geospatial data, using morphological spatial pattern recognition technology, area thresholding, and connectivity analysis, specifically includes: The types of ecological source areas are determined based on land use data; these types include: forest land, grassland, shrubland, water bodies, and wetlands. Based on the type of ecological source area, morphological spatial pattern identification technology is used to determine the core area patches in the ecological space; Core patches are obtained by screening core areas in the ecological space based on area thresholds; Connectivity analysis was performed on the core patches to determine the patch importance index; Ecological source areas are determined based on core patches with a patch importance index greater than 1.

2. The method for identifying key areas of ecological security pattern according to claim 1, characterized in that, The connectivity analysis of the core patches to determine the patch importance index specifically includes: Using formula Determine the plaque importance index dPC; Where PC is the potential connectivity index. , where n is the total number of core plaques. and A represents the area of ​​the i-th and j-th core patches, respectively. L P represents the total area of ​​the region to be identified. ij This represents the maximum product probability of all paths between the i-th and j-th core patches. PCremove The potential connectivity index of the core patches remaining after removing a single core patch.

3. The method for identifying key areas of ecological security pattern according to claim 1, characterized in that, The process involves determining ecological corridors based on ecological source areas and comprehensive ecological resistance surfaces using a minimum cost path model, specifically including: Using formula Determine the minimum cost path model; Where LCP is the final minimum cost path generated from the j-th ecological source to the i-th ecological source, and D ij R represents the spatial distance traversed from the j-th ecological source area to the i-th ecological source area. i Let represent the ecological resistance value of the i-th ecological source site in a certain direction in space, where m and n refer to the number of rows and columns traversed by the i-th to j-th ecological source sites in the ecological resistance raster data, where m is the number of rows and n is the number of columns. This is a function that takes the minimum value.

4. The method for identifying key areas of ecological security pattern according to claim 1, characterized in that, The analysis of the ecological security pattern network topology characteristics based on the ecological security pattern network diagram yields the following results: Based on the ecological security pattern network diagram, the GN algorithm is used to determine the community structure; Core-periphery structure analysis is performed on each community structure in the ecological security pattern network diagram to obtain the core nodes and peripheral nodes within the community structure.

5. The method for identifying key areas of ecological security pattern according to claim 1, characterized in that, The process of obtaining ecological pressure index factors of the ecological source area and determining the comprehensive ecological pressure index of the ecological source area based on these factors specifically includes: Using formula Determine the amount of soil erosion; among which, For the soil erosion modulus of the grid, For vegetation coverage, The rainfall erosivity index, The soil erodibility index. Factors related to soil and water conservation measures; Using formula Determine population agglomeration pressure factors Where pd is the population density; Using formula Determine land use intensity factors ;in, Let A be the area of ​​the i-th land use type in the grid, and let A be the area of ​​a single grid cell. Assign a value to the land use intensity of the i-th land use type; Using formula Determine the comprehensive ecological pressure index (EP) of the ecological source area; among which, This is the normalized value of soil erosion. This is the normalized value of the population agglomeration pressure factor. This is the normalized value of the land use intensity factor.

6. The method for identifying key areas of ecological security pattern according to claim 1, characterized in that, The process of determining the ecological security pattern zoning of the area to be identified based on the network topology of the ecological security pattern network diagram and the comprehensive ecological pressure index specifically includes: Based on the comprehensive ecological pressure index, the natural discontinuity method is used to classify the ecological pressure levels of ecological source areas. Based on the network topology of the ecological security pattern network diagram and the comprehensive ecological pressure index, the ecological source areas are divided into multiple ecological security clusters according to the ecological security pattern zoning principle. Within each ecological security cluster, an ecological security core area and an ecological security control area are defined. Based on the ecological pressure levels of the core ecological security area, the ecological security control area, and the ecological source area, ecologically fragile areas and ecologically sensitive areas are identified.

7. A computer device, comprising: A memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that the processor executes the computer program to implement the method for identifying key areas of ecological security patterns according to any one of claims 1-6.

8. A computer-readable storage medium having a computer program stored thereon, characterized in that, When executed by a processor, the computer program implements the method for identifying key areas of ecological security pattern as described in any one of claims 1-6.

9. A computer program product, comprising a computer program, characterized in that, When executed by a processor, the computer program implements the method for identifying key areas of ecological security pattern as described in any one of claims 1-6.

Images