Urban ecological security pattern construction and space-time analysis method fusing ecological sensitivity

Abstract

The invention provides an ecological sensitivity fused urban ecological safety pattern construction and space-time analysis method. The method comprises the following steps: identifying a key ecological source and constructing a resistance surface; and extracting at least one ecological corridor from the key ecological source and the resistance surface input value MRC model, and constructing an ecological safety pattern based on the at least one ecological corridor, the ecological source and the resistance surface. According to the invention, social and economic data can be integrated into ecological pattern construction, and the rationality of an ecological security pattern is improved.

Description

Technical Field

[0001] The present invention relates to the technical field of ecological management, and in particular to a method for constructing an urban ecological security pattern and performing spatiotemporal analysis integrating ecological sensitivity. Background Art

[0002] A healthy and stable ecosystem is the most important material basis for the survival and development of human society. Especially in the process of urbanization, the essence of maintaining a healthy ecosystem is to meet the needs of human survival and development for the natural ecological environment, thereby achieving sustainable development of human society. With economic development and changes in population spatial structure, the rapid expansion of urban space has greatly squeezed the natural ecological space.

[0003] Ecological security pattern is an important natural solution to the ecological problems brought about by urbanization and industrialization. Ecological security pattern refers to a potential spatial pattern in the landscape, which is composed of some key parts, their locations and the spatial connections between the key parts. The main function of the ecological security pattern is to ensure that the ecological processes in certain specific areas can proceed safely and stably. By analyzing ecological sensitivity, constructing urban ecological security patterns, and analyzing spatial and temporal evolution trends, we can quantify the status of urban ecological security, deepen our understanding of urban ecosystems, and provide theoretical support and practical guidance for urban sustainable development and ecological environmental protection. In the existing technology, ecological sensitivity is used as an important part of urban ecological research. However, the sole use of ecological sensitivity does not make full use of social and economic data, which leads to the underestimation of landscape elements close to green areas, cultivated land, construction land and other areas in extracting ecological sources, resulting in the inability to reasonably construct an ecological security pattern strategy. Summary of the invention

[0004] In view of this, the purpose of the present invention is to provide a method for constructing an urban ecological security pattern and spatiotemporal analysis that integrates ecological sensitivity, in order to solve at least some of the above problems.

[0005] In the first aspect, the present application provides a method for constructing an urban ecological security pattern integrating ecological sensitivity, including: Identify key ecological sources and construct resistance surfaces; Extract at least one ecological corridor from the key ecological source and the resistance surface input value MRC model; constructing an ecological security pattern based on the at least one ecological corridor, ecological source and resistance surface; The key ecological sources are identified in the following manner: Access to physical geographic data; Determine the sensitivity factors based on the physical geographic data, and calculate the sensitivity level corresponding to each sensitivity factor; Determine the comprehensive sensitivity based on the sensitivity level corresponding to each sensitive factor; Setting a normalized vegetation threshold and calculating landscape connectivity, and identifying key ecological sources based on the normalized vegetation threshold, the comprehensive sensitivity and the landscape connectivity; The resistance surface is constructed as follows: Acquire social statistical data, and determine at least one influencing factor based on the social statistical data; Determining the weight of each influencing factor in the at least one influencing factor by using a hierarchical analysis method; The resistance surface is constructed based on the weight of each influencing factor.

[0006] Through this application, social and economic data can be integrated into the construction of the ecological pattern, thus improving the rationality of the ecological security pattern.

[0007] One possible design is that the sensitive factors include at least one of elevation, slope, aspect, water buffer zone, vegetation coverage and land use type.

[0008] One possible design is that, in the step of determining the sensitivity factor based on the natural geographical data and setting the sensitivity level corresponding to each sensitive factor, the sensitivity level corresponding to the vegetation coverage rate is determined in the following manner: Preprocessing the physical geographic data and calculating a plurality of normalized vegetation values; Setting a vegetation value confidence interval, and determining a maximum vegetation value and a minimum vegetation value, wherein the maximum vegetation value is the maximum value of the plurality of normalized vegetation values ​​in the vegetation value confidence interval, and the minimum vegetation value is the minimum value of the plurality of normalized vegetation values ​​in the vegetation value confidence interval; Determine a plurality of vegetation coverage values ​​based on the maximum vegetation value and the minimum vegetation value, wherein the plurality of vegetation coverage values ​​include: a vegetation coverage value corresponding to each normalized vegetation value in a plurality of normalized vegetation values, and a normalized vegetation value in a plurality of normalized vegetation values ​​corresponds to a vegetation coverage value; Determine at least one first sensitivity level, and determine a vegetation coverage value range corresponding to each of the at least one first sensitivity level, and determine a sensitivity level corresponding to each of the multiple vegetation coverage values ​​based on the vegetation coverage value range corresponding to each of the first sensitivity levels; Based on the sensitivity level corresponding to each vegetation coverage rate, the sensitivity level corresponding to the vegetation coverage rate is determined.

[0009] One possible design is that, in the step of determining the sensitivity factor based on the natural geographical data and setting the sensitivity level corresponding to each sensitive factor, the sensitivity level corresponding to the water buffer zone is determined in the following manner: Filling depressions, calculating flow direction and flow rate on the natural geographical data, and determining at least one river link; determining at least one second sensitivity level, and determining a river network level of each river link in the at least one river link based on the second sensitivity level; The sensitivity level corresponding to the water area buffer zone is determined based on the river network level of each river link.

[0010] One possible design is, In the step of setting the normalized vegetation threshold and calculating the landscape connectivity, and identifying the key ecological source based on the normalized vegetation threshold, the comprehensive sensitivity and the landscape connectivity, the landscape connectivity is calculated in the following manner: ; ; Among them, IIC—Integral Connectivity Index; PC—connectivity probability; i and j—different ecological patches; M—the total number of ecological patches in the landscape; and — patches i and j, is the number of connections between patches i and j; It is the biggest landscape attribute; It refers to the maximum product probability of all paths between source points i and j.

[0011] In a possible design, the step of determining the weight of each influencing factor in the at least one influencing factor by using the hierarchical analysis method comprises: Based on the relationship and characteristics of each factor among the at least one influencing factor, a hierarchical model is generated, wherein the hierarchical model includes a plurality of layers, the at least one influencing factor is allocated to the plurality of layers, and the influencing factors at the upper layer dominate the influencing factors at the lower layer or are affected by the factors at the lower layer; determining a comparison matrix for each of the plurality of strata; constructing a judgment matrix based on the comparison matrix of each layer in the plurality of layers; Outputting a weight value corresponding to each influencing factor in the at least one influencing factor based on the judgment matrix; Wherein, in the step of determining the comparison matrix of each layer in the multiple layers, the comparison matrix of the target layer in the multiple layers is determined as follows: Identify the influencing factors within the target layer; Based on the relative importance of each influencing factor in the target layer, the paired comparison method is adopted and the comparison matrix is ​​constructed using the ratio scale.

[0012] One possible design is: In the step of extracting at least one ecological corridor from the key ecological source and the resistance surface input value MRC model, the MRC model is as follows: ; Among them, MCR is the minimum cumulative resistance; —The positive correlation between the drag at any point in space and the drag at any other point; i and j—different ecological source points; —The distance from ecological source point i to ecological source point j; —iLand resistance coefficient during movement.

[0013] A possible design is that the method further includes: performing spatial syntax analysis on the ecological corridors to determine the priority of each ecological corridor, and optimizing the ecological security pattern based on the priority of each ecological corridor.

[0014] One possible design is to pre-treat the ecological corridor. Determine the axis map of each of the ecological corridors; Based on the axis map of each ecological corridor, the priority of each ecological corridor is determined by taking integration and selectivity as indicators; Among them, the axis map of an ecological corridor is determined in the following way: Determine all inflection points of the one ecological corridor axis; An axis map of the ecological corridor will be drawn based on all the inflection points.

[0015] In the second aspect, an embodiment of the present application provides a method for analyzing the spatiotemporal evolution of an urban ecological security pattern that integrates ecological sensitivity to construct an ecological security pattern for a target year, wherein the method described in the first aspect is used to determine the ecological security pattern for any year in the target year.

[0016] Based on the changes in the ecological security pattern in the target year, corresponding ecological management suggestions are put forward.

[0017] Other features and advantages of the present invention will be described in the following description, and partly become apparent from the description, or understood by practicing the present invention. The purpose and other advantages of the present invention are realized and obtained by the structures particularly pointed out in the description, claims and drawings.

[0018] In order to make the above-mentioned objects, features and advantages of the present invention more obvious and easy to understand, preferred embodiments are given below and described in detail with reference to the accompanying drawings. BRIEF DESCRIPTION OF THE DRAWINGS

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

[0020] Figure 1 A flow chart of a method for constructing an urban ecological security pattern integrating ecological sensitivity provided in an embodiment of the present invention; Figure 2 A flow chart of determining the sensitivity level corresponding to the vegetation coverage rate provided by an embodiment of the present invention; Figure 3 A flow chart for determining the sensitivity level corresponding to the water buffer zone provided in an embodiment of the present invention; Figure 4 This is an example diagram of the axis map provided in the embodiment of the present application. DETAILED DESCRIPTION

[0021] In order to better understand the technical solution of the embodiment of the present invention, the embodiment of the present invention is described in detail below with reference to the accompanying drawings.

[0022] It should be clear that the described embodiments are only part of the embodiments of the present invention, not all of them. Based on the embodiments in 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 embodiments of the present invention.

[0023] The terms used in the embodiments of the present invention are only for the purpose of describing specific embodiments, and are not intended to limit the embodiments of the present invention. The singular forms "a", "said" and "the" used in the embodiments of the present invention and the appended claims are also intended to include plural forms, unless the context clearly indicates other meanings.

[0024] At present, a healthy and stable ecosystem is the most important material basis for the survival and development of human society. Especially in the process of urbanization, the essence of maintaining the health of the ecosystem is to meet the needs of human survival and development for the natural ecological environment, thereby achieving sustainable development of human society. With the development of the economy and the changes in the spatial structure of the population, the rapid expansion of urban space has greatly squeezed the natural ecological space.

[0025] Ecological security pattern is an important natural solution to the ecological problems brought about by urbanization and industrialization. Ecological security pattern refers to a potential spatial pattern in the landscape, which is composed of some key parts, their locations, and the spatial connections between key parts. The main function of ecological security pattern is to ensure that ecological processes in certain specific areas can proceed safely and stably. By analyzing ecological sensitivity, constructing urban ecological security pattern, and analyzing the trend of spatiotemporal evolution, we can quantify the status of urban ecological security, deepen our understanding of urban ecosystems, and provide theoretical support and practical guidance for urban sustainable development and ecological environmental protection.

[0026] In the existing technology, ecological sensitivity is used as an important part of urban ecological research. However, the sole use of ecological sensitivity does not make full use of social and economic data, which leads to the underestimation of landscape elements close to green areas, cultivated land, construction land and other areas in extracting ecological sources, resulting in the inability to reasonably construct an ecological security pattern strategy.

[0027] In view of this, this application provides a method for constructing an urban ecological security pattern that integrates ecological sensitivity, so as to integrate social and economic data into the construction of the ecological pattern and enhance the rationality of the ecological security pattern.

[0028] In view of this, refer to Figure 1 This application provides a method for constructing an urban ecological security pattern that integrates ecological sensitivity, specifically including: S10: Identify key ecological sources and construct resistance surfaces.

[0029] S20: Extract at least one ecological corridor from the key ecological source and the resistance surface input value MRC model.

[0030] S30: Constructing an ecological security pattern based on the at least one ecological corridor, ecological source and resistance surface.

[0031] Here, the ecological sources can be identified in the following ways: first, obtain the natural geographical data. Specifically, the natural geographical data include: land utilization data, meteorological data, vegetation cover data, digital elevation data, soil data and vegetation net primary productivity (NPP) data; then determine the sensitive factors based on the natural geographical data, and calculate the sensitivity level corresponding to each sensitive factor.

[0032] After the sensitivity level corresponding to each sensitive factor is obtained, the comprehensive sensitivity can be determined based on the sensitivity level corresponding to each sensitive factor.

[0033] After the comprehensive sensitivity is obtained, a normalized vegetation threshold is set and the landscape connectivity is calculated, and key ecological sources are identified based on the normalized vegetation threshold, the comprehensive sensitivity and the landscape connectivity.

[0034] The key ecological sources are identified in the above manner.

[0035] Regarding the resistance surface, the resistance surface is constructed specifically using sociological statistical data. First, the social statistical data is obtained, and at least one influencing factor is determined based on the social statistical data.

[0036] Then, the weight of each influencing factor in at least one of the aforementioned influencing factors is determined by using the hierarchical analysis method. Finally, the resistance surface is constructed based on the weight of each influencing factor.

[0037] Specifically, at least one influencing factor corresponds to all influencing factors.

[0038] Through the above methods, social and economic data can be integrated into the construction of the ecological pattern to improve the rationality of the ecological security pattern.

[0039] As a preferred embodiment, the aforementioned sensitive factors include: at least one of elevation, slope, aspect, water buffer zone, vegetation coverage and land use type.

[0040] Reference Figure 2 As a preferred embodiment, the sensitivity level corresponding to the vegetation coverage rate can be determined in the following manner: S101: Preprocess the natural data and calculate a plurality of normalized vegetation values.

[0041] In this step, data preprocessing is achieved through operations such as unified projection coordinate system and radiation calibration. Specifically, the projection coordinate system and radiation calibration are first unified, and then the results of radiation calibration are mosaicked. The preprocessing is completed in the above manner.

[0042] For the normalized vegetation value (NDVI value), it can be calculated as follows: ; in, —Normalized vegetation value, is the reflectance value in the near-infrared band, is the reflectance value of the red light band.

[0043] After the normalized vegetation value is calculated, S102 is executed to set the vegetation value confidence interval and determine the maximum vegetation value and the minimum vegetation value. The maximum vegetation value is the maximum value of the multiple normalized vegetation values ​​in the vegetation value confidence interval.

[0044] Specifically, the vegetation confidence interval can be set to [5%, 95%]. After the vegetation confidence interval is set, the maximum vegetation value and the minimum vegetation value are calculated.

[0045] It should be noted that the maximum vegetation value is the maximum value of the multiple normalized vegetation values ​​in the vegetation value confidence interval, and the minimum vegetation value is the minimum value of the multiple normalized vegetation values ​​in the vegetation value confidence interval.

[0046] S103: Determine a plurality of vegetation coverage values ​​based on the maximum vegetation value and the minimum vegetation value.

[0047] In this step, the above-mentioned multiple vegetation coverage values ​​can be calculated in the following way: ; —Vegetation cover value; —Normalized vegetation value; —maximum vegetation value; —Minimum vegetation value.

[0048] S104: Determine at least one first sensitivity level, and determine a vegetation coverage value range corresponding to each of the at least one first sensitivity level, and determine a sensitivity level corresponding to each of the multiple vegetation coverage values ​​based on the vegetation coverage value range corresponding to each first sensitivity level; It should be noted that the first sensitivity level, i.e., the sensitivity level corresponding to the vegetation coverage rate, can be specifically divided into 5 sets, and each level increases by 0.2.

[0049] S105: Based on the sensitivity level corresponding to each vegetation coverage rate, determine the sensitivity level corresponding to the vegetation coverage rate.

[0050] In this step, the raster data is converted into surface data to facilitate the calculation of the area of ​​each sensitive area in the five different first sensitivity levels. An area column is added to the attribute table of the surface data, and the area of ​​each sensitive area is calculated using the geometry calculator.

[0051] Through the above method, the sensitivity level corresponding to the vegetation coverage rate is determined.

[0052] Reference Figure 3 , for the sensitivity level corresponding to the water buffer zone: it can be calculated in the following way: S201: Fill depressions, calculate flow direction and flow rate on the natural geographical data, and determine at least one river link.

[0053] S202: Determine at least one second sensitivity level, and determine a river network level of each river link in the at least one river link based on the second sensitivity level.

[0054] S203: Determine the sensitivity level corresponding to the water area buffer zone based on the sensitivity level corresponding to each water area point.

[0055] Specifically, operations such as polygon to raster conversion and reclassification can be used to determine the sensitivity level corresponding to the water buffer zone, which will not be elaborated here.

[0056] As a preferred embodiment, in order to obtain and calculate the landscape connectivity, the following formula may be specifically used: ; ; Among them, IIC—Integral Connectivity Index; PC—connectivity probability; i and j—different ecological patches; M—the total number of ecological patches in the landscape; and — patches i and j, is the number of connections between patches i and j; It is the biggest landscape attribute; It refers to the maximum product probability of all paths between source points i and j.

[0057] The key ecological sources are identified in the above manner.

[0058] The following will specifically describe the steps of determining the weight of each influencing factor in the at least one influencing factor by using the analytic hierarchy process: First, a hierarchical structure model is generated based on the relationship and characteristics of each factor among the at least one influencing factor.

[0059] Specifically, the at least one influencing factor includes at least one of: elevation, distance to water source, distance to other buildings, distance to residence, slope and land use type.

[0060] Specifically, the hierarchical model includes multiple layers, the at least one influencing factor is distributed to the multiple layers, and in the multiple layers, the influencing factors at the upper layer dominate the influencing factors at the lower layer or are affected by the factors at the lower layer.

[0061] Then, a judgment matrix is ​​constructed based on the comparison matrix of each layer in the multiple layers.

[0062] Finally, a weight value corresponding to each influencing factor in the at least one influencing factor is output based on the judgment matrix.

[0063] It should be noted that the judgment matrix needs to be verified, specifically using hierarchical single sorting, total sorting and consistency check.

[0064] The so-called single sorting is to calculate the maximum eigenvalue and corresponding eigenvector of each pairwise comparison matrix. Total sorting refers to calculating the combined weight vector of the lowest layer for the target. The consistency test is performed using the consistency index, random consistency index and consistency ratio. If the test passes, the eigenvector is the weight vector; if not, the pairwise comparison matrix needs to be reconstructed.

[0065] For a target layer among multiple layers, the comparison matrix of the target layer among multiple layers can be determined in the following manner: Firstly, the influencing factors within the target layer are determined. Based on the relative importance of each influencing factor in the above target layer, the paired comparison method is adopted and the comparison matrix is ​​constructed using the comparison scale.

[0066] Specifically, the paired comparison method and the 1-9 comparison scale were used to construct a paired comparison matrix.

[0067] The above method can be used to analyze the weight of each influencing factor using the hierarchical analysis method.

[0068] In view of the above, at least one ecological corridor is extracted from the key ecological source and the resistance surface input value MRC model. In this step, the MRC model is as follows: ; Among them, MCR is the minimum cumulative resistance; —The positive correlation between the drag at any point in space and the drag at any other point; i and j—different ecological source points; —The distance from ecological source point i to ecological source point j; —iLand resistance coefficient during movement.

[0069] On the basis of the above-mentioned embodiment, the priority of each ecological corridor can be determined by performing space syntax analysis on the ecological corridor, and the ecological security pattern can be optimized based on the priority of each ecological corridor.

[0070] The basis of space syntax analysis is to reflect the relative accessibility of different locations in a city or region by measuring different characteristics of the network.

[0071] On the basis of constructing the ecological security pattern, the turning points of each ecological corridor are taken as discontinuity points to determine the preliminary corridors. Then a single corridor is divided into multiple axes to construct a corridor axis map. Corridors are prioritized according to the global integration and selectivity values.

[0072] Specifically, the ecological corridors were first preprocessed, and then the axis map of each ecological corridor was constructed.

[0073] Reference Figure 4 ,For the axis map, the following rules need to be followed: ensure that the connection relationship between convex spaces remains unchanged during the transformation process, and at the same time, the longest and smallest number of axes should be selected to run through all convex spaces.

[0074] Finally, based on the axis map of each ecological corridor, the priority of each ecological corridor is determined with integration and selectivity as indicators.

[0075] Integration describes the state of aggregation or dispersion between an element and other elements in a spatial system. The higher the integration, the more accessible and central it is, making it easier to attract and aggregate species flows. The formula is as follows. RRA describes the deviation from the standard value.

[0076] ; —Ecological Corridor Deviation from the standard value; —Ecological Corridor degree of integration.

[0077] Choice is used to measure the frequency with which an element in a spatial system is selected as the shortest topological distance between two nodes. The higher the choice, the more likely it is that the space will be traversed by biological flows.

[0078] Selectivity is closely related to global integration, which indicates the relative importance and accessibility of a spatial element in the overall spatial structure. Spatial elements with high selectivity usually have higher global integration values, meaning they are easier to access and connect.

[0079] As a preferred embodiment, the axis map of an ecological corridor is determined in the following manner: First, all the inflection points of the axis of the ecological corridor are determined, and then, an axis map of the ecological corridor is drawn based on all the inflection points.

[0080] On the basis of the above-mentioned embodiments, the present application also provides a method for analyzing the spatiotemporal evolution of urban ecological security patterns integrating ecological sensitivity, specifically using the above-mentioned method to construct the ecological security pattern of any year in the target year, and performing analysis based on the ecological security pattern of the target year; Through the above methods, by integrating the ecological security pattern over many years, corresponding ecological management suggestions can be put forward.

[0081] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present invention should be included in the scope of protection of the present invention.

Claims

1. A method for constructing an urban ecological security pattern integrating ecological sensitivity, characterized in that: include: Identify key ecological sources and construct resistance surfaces; Extract at least one ecological corridor from the key ecological source and the resistance surface input value MRC model; constructing an ecological security pattern based on the at least one ecological corridor, ecological source and resistance surface; The key ecological sources are identified in the following manner: Access to physical geographic data; Determine the sensitivity factors based on the physical geographic data, and calculate the sensitivity level corresponding to each sensitivity factor; Determine the comprehensive sensitivity based on the sensitivity level corresponding to each sensitive factor; Setting a normalized vegetation threshold and calculating landscape connectivity, and identifying key ecological sources based on the normalized vegetation threshold, the comprehensive sensitivity and the landscape connectivity; The resistance surface is constructed as follows: Acquire social statistical data, and determine at least one influencing factor based on the social statistical data; Determining the weight of each influencing factor in the at least one influencing factor by using a hierarchical analysis method; The resistance surface is constructed based on the weight of each influencing factor.

2. The method according to claim 1, characterized in that The sensitive factors include at least one of elevation, slope, aspect, water buffer zone, vegetation coverage and land use type.

3. The method according to claim 2, characterized in that: In the step of determining the sensitivity factor based on the physical geographic data and setting the sensitivity level corresponding to each sensitive factor, the sensitivity level corresponding to the vegetation coverage rate is determined in the following manner: Preprocessing the physical geographic data and calculating a plurality of normalized vegetation values; Setting a vegetation value confidence interval, and determining a maximum vegetation value and a minimum vegetation value, wherein the maximum vegetation value is the maximum value of the plurality of normalized vegetation values ​​in the vegetation value confidence interval, and the minimum vegetation value is the minimum value of the plurality of normalized vegetation values ​​in the vegetation value confidence interval; Determine a plurality of vegetation coverage values ​​based on the maximum vegetation value and the minimum vegetation value, wherein the plurality of vegetation coverage values ​​include: a vegetation coverage value corresponding to each normalized vegetation value in a plurality of normalized vegetation values, and a normalized vegetation value in a plurality of normalized vegetation values ​​corresponds to a vegetation coverage value; Determine at least one first sensitivity level, and determine a vegetation coverage value range corresponding to each of the at least one first sensitivity level, and determine a sensitivity level corresponding to each of the multiple vegetation coverage values ​​based on the vegetation coverage value range corresponding to each of the first sensitivity levels; Based on the sensitivity level corresponding to each vegetation coverage rate, the sensitivity level corresponding to the vegetation coverage rate is determined.

4. The method according to claim 3, characterized in that In the step of determining the sensitivity factor based on the natural geographical data and setting the sensitivity level corresponding to each sensitive factor, the sensitivity level corresponding to the water area buffer zone is determined in the following manner: Filling depressions, calculating flow direction and flow rate on the natural geographical data, and determining at least one river link; determining at least one second sensitivity level, and determining a river network level of each river link in the at least one river link based on the second sensitivity level; The sensitivity level corresponding to the water area buffer zone is determined based on the river network level of each river link.

5. The method according to claim 3, characterized in that: In the step of setting the normalized vegetation threshold and calculating the landscape connectivity, and identifying the key ecological source based on the normalized vegetation threshold, the comprehensive sensitivity and the landscape connectivity, the landscape connectivity is calculated in the following manner: ; ; Among them, IIC—Integral Connectivity Index; PC—connectivity probability; i and j—different ecological patches; M—the total number of ecological patches in the landscape; and — patches i and j, is the number of connections between patches i and j; It is the biggest landscape attribute; It refers to the maximum product probability of all paths between source points i and j.

6. The method according to claim 1, characterized in that The step of determining the weight of each influencing factor in the at least one influencing factor by using the hierarchical analysis method comprises: Based on the relationship and characteristics of each factor among the at least one influencing factor, a hierarchical model is generated, wherein the hierarchical model includes a plurality of layers, the at least one influencing factor is allocated to the plurality of layers, and the influencing factors at the upper layer dominate the influencing factors at the lower layer or are affected by the factors at the lower layer; determining a comparison matrix for each of the plurality of strata; constructing a judgment matrix based on the comparison matrix of each layer in the plurality of layers; Outputting a weight value corresponding to each influencing factor in the at least one influencing factor based on the judgment matrix; Wherein, in the step of determining the comparison matrix of each layer in the multiple layers, the comparison matrix of the target layer in the multiple layers is determined as follows: Identify the influencing factors within the target layer; Based on the relative importance of each influencing factor in the target layer, the pairwise comparison method is adopted to construct a comparison matrix using the ratio scale; The at least one influencing factor includes at least one of elevation, distance to a water source, distance to other buildings, distance to a residence, slope and land use type.

7. The method according to claim 4, characterized in that The method further comprises: In the step of extracting at least one ecological corridor from the key ecological source and the resistance surface input value MRC model, the MRC model is as follows: ； Among them, MCR is the minimum cumulative resistance; —The positive correlation between the drag at any point in space and the drag at any other point; i and j—different ecological source points; —The distance from ecological source point i to ecological source point j; —iLand resistance coefficient during movement.

8. The method according to any one of claims 1 to 7, characterized in that: The method further includes: performing spatial syntax analysis on the ecological corridors to determine the priority of each ecological corridor, and optimizing the ecological security pattern based on the priority of each ecological corridor.

9. The method according to any one of claims 1 to 7, characterized in that: The ecological corridor is pre-processed, Determine the axis map of each of the ecological corridors; Based on the axis map of each ecological corridor, the priority of each ecological corridor is determined by taking integration and selectivity as indicators; Among them, the axis map of an ecological corridor is determined in the following way: Determine all inflection points of the one ecological corridor axis; An axis map of the ecological corridor will be drawn based on all the inflection points.

10. A method for analyzing the spatiotemporal evolution of urban ecological security pattern integrating ecological sensitivity, characterized in that: Constructing an ecological security pattern for a target year, wherein the ecological security pattern for any year in the target year is determined by using the method described in any one of claims 1 to 9; Based on the changes in the ecological security pattern in the target year, corresponding ecological management suggestions are put forward.

Images