Construction method of ecological landscape corridor, electronic device and storage medium

Abstract

The embodiment of the present application provides a method for constructing an ecological landscape corridor, an electronic device, and a storage medium, the method comprising: obtaining sample plot data of a target area; calculating a habitat quality index of the target area through the sample plot data; determining a high habitat quality area from the target area according to the habitat quality index; constructing a primary ecological landscape corridor based on the high habitat quality area; using an insect pollination model of the target area to analyze the insect pollination service capacity in the primary ecological landscape corridor, and obtaining an analysis result of the insect pollination service capacity; optimizing the primary ecological landscape corridor according to the analysis result to obtain an ecological landscape corridor. The embodiment of the present application constructs a reasonable ecological landscape corridor based on the impact of the ecological landscape corridor on insect pollination, which can meet the ecological needs of insect pollinators and is beneficial to the diversity and genetic diversity of organisms in the ecological landscape corridor.

Description

Technical Field

[0001] The embodiments of the present application relate to the field of ecological construction technology, and in particular to a method for constructing an ecological landscape corridor, an electronic device, and a storage medium. Background Art

[0002] Ecological landscape corridors are connecting belts established to connect two or more fragmented natural habitats. They play a positive role in solving habitat fragmentation problems and maintaining biodiversity levels. Therefore, it is extremely important to scientifically and reasonably identify ecological landscape corridors.

[0003] The construction of existing ecological landscape corridors only considers how to improve the organization and distribution of plant populations. Insect pollination, as one of the important links in plant propagation in ecosystem services, helps to increase the genetic diversity of plant populations and maintain and enhance the biodiversity within ecological landscape corridors. It can improve the composition and distribution of plant populations, promote the exchange and migration of species between different habitats, and enhance the self-repair and recovery capacity of ecosystems.

[0004] Therefore, we are faced with the problem of being unable to improve biodiversity and protect ecosystem health and sustainable development by enhancing the service functions of insect pollination while building ecological landscape corridors. Summary of the invention

[0005] In view of this, a method for constructing an ecological landscape corridor is provided to at least partially solve the above problems.

[0006] According to a first aspect of an embodiment of the present application, a method for constructing an ecological landscape corridor is provided, comprising: obtaining sample plot data of a target area; calculating a habitat quality index of the target area through the sample plot data; determining a high habitat quality area from the target area according to the habitat quality index; constructing a primary ecological landscape corridor based on the high habitat quality area; using an insect pollination model of the target area to analyze the insect pollination service capacity in the primary ecological landscape corridor and obtain analysis results of the insect pollination service capacity; optimizing the primary ecological landscape corridor according to the analysis results to obtain an ecological landscape corridor.

[0007] According to a second aspect of an embodiment of the present application, there is provided an electronic device, comprising: a processor, a memory, a communication interface and a communication bus, wherein the processor, the memory and the communication interface communicate with each other via the communication bus; the memory is used to store at least one executable instruction, and the executable instruction enables the processor to execute the method described in the first aspect.

[0008] According to a third aspect of an embodiment of the present application, a computer storage medium is provided, on which a computer program is stored. When the program is executed by a processor, the method for constructing an ecological landscape corridor as described in the first aspect is implemented.

[0009] To sum up, in the embodiments of the present application, high habitat quality areas are determined through sample plot data of the target area, and ecological landscape corridors are constructed in combination with the insect pollination model. On the basis of solving the problem of habitat fragmentation and maintaining the level of biodiversity, this application further improves the service performance of insect pollination in biological diversity and protecting the health of the ecosystem, so as to enable the sustainable development of the ecosystem. BRIEF DESCRIPTION OF THE DRAWINGS

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

[0011] Figure 1 Shown is a flow chart of a method for constructing an ecological landscape corridor according to the present application;

[0012] Figure 2 Shown is a schematic structural diagram of an electronic device according to an embodiment of the present application. DETAILED DESCRIPTION

[0013] In order to enable those skilled in the art to better understand the technical solutions in the embodiments of the present application, the technical solutions in the embodiments of the present application will be clearly and completely described below in conjunction with the drawings in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, not all of the embodiments. All other embodiments obtained by ordinary technicians in this field based on the embodiments in the embodiments of the present application should fall within the scope of protection of the embodiments of the present application.

[0014] It should be understood that the various steps described in the method embodiments of the present disclosure may be performed in different orders and / or in parallel. In addition, the method embodiments may include additional steps and / or omit the steps shown. The scope of the present disclosure is not limited in this respect.

[0015] The term "including" and its variations used in this document are open inclusions, that is, "including but not limited to". The term "based on" means "based at least in part on". The term "one embodiment" means "at least one embodiment"; the term "another embodiment" means "at least one other embodiment"; the term "some embodiments" means "at least some embodiments". Relevant definitions of other terms will be given in the description below. It should be noted that the concepts of "first", "second", etc. mentioned in this disclosure are only used to distinguish different devices, modules or units, and are not used to limit the order or interdependence of the functions performed by these devices, modules or units.

[0016] It should be noted that the modifications of "one" and "multiple" mentioned in the present disclosure are illustrative rather than restrictive, and those skilled in the art should understand that unless otherwise clearly indicated in the context, it should be understood as "one or more". The names of the messages or information exchanged between the multiple devices in the embodiments of the present disclosure are only for illustrative purposes and are not used to limit the scope of these messages or information.

[0017] The construction of existing ecological landscape corridors only considers how to improve the organization and distribution of plant populations, while insect pollination, as one of the important links in plant propagation in ecosystem services, is often overlooked. Based on this, an embodiment of the present application provides a method for constructing an ecological landscape corridor. The specific implementation of each embodiment of the present disclosure will be described in detail below in conjunction with the accompanying drawings.

[0018] For example, Figure 1 The flowchart of the method for constructing an ecological landscape corridor according to the present application is shown, which mainly includes the following steps:

[0019] Step S102, obtaining sample data of the target area.

[0020] In some embodiments, the sample data includes but is not limited to land cover (Land-Use and Land-CoverChange, LUCC) data, Digital Elevation Model (Digital Elevation Mode, DEM) data, Normalized Difference Vegetation Index (Normalized Difference Vegetation Index, NDVI) data, land use data, monthly average temperature data and monthly average sunshine data, plant species and quantity data, plant coverage data, pollinator species data, etc.

[0021] Exemplarily, land cover data include, but are not limited to: land cover type, land cover change analysis, land cover ratio, etc.

[0022] It can be understood that land cover types include but are not limited to: farmland, grassland, broad-leaved tree forest, coniferous tree forest, shrub forest, wasteland, hardened surface, wetland, etc.

[0023] In some embodiments, land cover changes are further analyzed based on the landscape pattern index around three types of plant landscapes, namely, tree forests, shrub forests and grasslands. Exemplarily, the analysis content includes: spatial configuration characteristics such as the number, shape and structure of patches of each type of plant landscape, landscape heterogeneity and degree of fragmentation, etc.

[0024] For example, the landscape pattern index includes, but is not limited to, the following: patch type and patch percentage (PLAND), patch density (NP), edge density (ED), landscape shape index (LSI), cohesion (COHESION), spread index (CONTAG) and Shannon's Diversity Index (SHDI). Referring to Table 1, in the embodiments of the present application, the landscape pattern index and its meaning are shown in Table 1 below:

[0025]

[0026] In some embodiments, the digital elevation model is used to analyze the elevation distribution characteristics of the target area and to calculate terrain indicators.

[0027] Exemplarily, the terrain indicators include, but are not limited to: the slope of the target area, the slope aspect of the target area, the soil moisture index of the target area, etc.

[0028] In some embodiments, the normalized vegetation index data is used to calculate the main indicators that affect the insect pollination service function, such as nesting suitability and nectar and pollen source plant abundance.

[0029] In some embodiments, the monthly average temperature data and the monthly average sunshine data are used in some embodiments, and the digital elevation model is used to analyze the elevation distribution characteristics of the target area and to calculate the terrain index.

[0030] Exemplarily, the terrain indicators include, but are not limited to: the slope of the target area, the slope aspect of the target area, the soil moisture index of the target area, etc.

[0031] In some embodiments, the normalized vegetation index data is used to calculate the main indicators that affect the insect pollination service function, such as nesting suitability and nectar and pollen source plant abundance.

[0032] In some embodiments, the monthly average temperature data and the monthly average sunshine data are used to calculate the activity coefficient of pollinating insects in different seasons.

[0033] For example, referring to Table 2, in the embodiments of the present application, some sample data are obtained through the public channels shown in Table 2 below:

[0034]

[0035] Illustratively, in some embodiments, the types and quantities of plants are obtained using a typical sampling method.

[0036] For example, in some embodiments, the types and numbers of pollinating insects are obtained by combining observation and netting methods.

[0037] Step S104, calculating the habitat quality index of the target area through the sample data.

[0038] In some embodiments, the threat source is mapped onto the land data grid obtained in step S102 to obtain the distance between the habitat grid and the threat source grid, and the influence of the distance between the habitat grid and the threat source grid on the habitat formation is calculated based on the distance between the habitat grid and the threat source grid. The threat level of the threat source is calculated based on the influence level, and the habitat quality index of each grid in the target area is calculated based on the threat level.

[0039] Threat sources are understood to refer to land use types that have threatening impacts on habitats, and such land use types generally include but are not limited to the following factors: interactions with native species, urbanization, agriculture, outdoor recreation and tourism development, livestock and ranching activities, reservoirs and other tap water modifications, water, air or soil pollution, mineral, natural gas, oil and geothermal extraction, or exploration industrial, institutional, or military activities, road construction and maintenance, loss of genetic variation, inbreeding depression or hybridization, aquifer depletion, wetland drainage or landfill, plant succession, background species interactions, disease, vandalism, etc.

[0040] Specifically, the influence degree value of the distance between the habitat grid and the threat source grid on the habitat formation is calculated based on the following formula 1 or the following formula 2.

[0041] , (Formula 1)

[0042] , (Formula 2)

[0043] Among them, in formula 1:

[0044] It is suitable to describe how the threat decays in space through a linear distance decay function. Furthermore, linear decay refers to the maximum effective distance reached by the linear distance between grid cells in space.

[0045] In formula 2:

[0046] Suitable for describing how the threat decays in space through an exponential distance decay function, illustratively, the exponential decay refers to a constant of 2.99, and further, in some embodiments, the impact of the threat is reduced by 95% (to 5%) at a specified maximum effective distance.

[0047] Used to represent the impact value. Used to characterize habitat grids Threat Source Grid The distance between Threat factors used to characterize threat sources, Threat factors used to characterize threat sources Habitat Grid The maximum impact distance.

[0048] In some embodiments, Habitat grid represented Threat Source Grid The distance between them is the linear distance.

[0049] Furthermore, is the effect of the distance between the habitat grid and the threat grid on the habitat formation. It can be understood that the degree of threat decreases as the distance between the grid and the threat source increases, so those grid cells closest to the threat will be more affected. Equation 1 or Equation 2 can be selected to describe how the threat is attenuated in space.

[0050] Furthermore, the threat level of the threat source is obtained based on the following formula.

[0051] , (Formula 3)

[0052] In formula 3:

[0053] Used to characterize the impact value, Used to characterize the threat source Habitat grids in habitat types The threat level It is used to characterize the number of threat factors included in the threat source. To characterize the The weight of the threat factor, To characterize the The number of grids occupied by the threat factor layer corresponding to each threat factor on the land data grid. The corresponding threat factor layer is used to represent the occupation of the land data raster grid. The number of threat factors per grid, Used to characterize habitat grids The accessibility level coefficient, To characterize the Habitat types and threatening factors The sensitivity coefficient.

[0054] Specifically, in some embodiments, Refers to the habitat within the target area. Grid The overall threat level is understandable, the threat level is habitat Grid The higher the threat level, the more degradation the habitat Grid The greater the degradation caused by the impact, the greater the degradation and vice versa.

[0055] Furthermore, the habitat quality index was obtained based on the following formula.

[0056] , (Formula 4)

[0057] In formula 4:

[0058] To characterize the Habitat grids in habitat types The habitat quality index, To characterize the The habitat suitability coefficient of the habitat type is Used to characterize the scale constant, Used to characterize the half-saturation constant.

[0059] Exemplarily, in some embodiments, the scaling constant takes a value of 2.5.

[0060] Step S106, determining a high habitat quality area from the target area according to the habitat quality index.

[0061] Specifically, the high habitat quality area in the target area is determined according to the habitat quality index in the target area obtained in step S104.

[0062] In some embodiments, the top 20% of the habitat quality index in the target area obtained in step S104 is taken as the high habitat quality area in the target area.

[0063] In some embodiments, the high habitat quality area is selected based on an area greater than two square kilometers.

[0064] Step S108, constructing a primary ecological landscape corridor based on the high habitat quality area.

[0065] A resistance surface is constructed using the high habitat quality areas in the target area obtained in step S106 using the sample plot data obtained in step S102, wherein the resistance surface is used to simulate and quantify the hindering or promoting effects of different landscape elements on species migration and ecological flow.

[0066] In the embodiment of the present application, The resistance surface of a high habitat quality area is calculated by the following formula:

[0067] , (Formula 5)

[0068] In formula 5:

[0069] No. The resistance surface of the high habitat quality area, To characterize the The sample data includes the terrain slope, Used to characterize terrain slope The weight coefficient of To characterize the The vegetation coverage rate included in the sample data is Used to characterize vegetation coverage The weight coefficient of To characterize the The land cover percentage included in the sample data is To characterize land cover The weight coefficient of To characterize the The ecological importance of the sample data included in each plot is To characterize ecological importance The weight coefficient of

[0070] Furthermore, according to The minimum cumulative resistance value is calculated from the resistance surface of the high habitat quality area, where the minimum cumulative resistance value is the minimum possible value of the cumulative resistance from the ecological source to the target area or another source. These paths are the routes with the least resistance for species to carry out seasonal migration or daily activities.

[0071] Furthermore, the minimum cumulative resistance is obtained based on the following formula.

[0072] , (Formula 6)

[0073] In formula 6:

[0074] Used to characterize the minimum cumulative resistance value, Used to characterize ecological origin to Distance to high habitat quality areas Used to characterize the minimum cumulative resistance value From the ecological source to the Distance to high habitat quality areas The positive correlation function between .

[0075] in, Acquired through step S106.

[0076] It can be understood that by calculating the minimum cumulative resistance value, the route with the minimum resistance value between the source points is obtained according to the minimum cumulative resistance value between the source points, and then the first generation ecological landscape corridor is constructed.

[0077] For example, the construction area of ​​the ecological landscape corridor may include land ecological routes, river ecological routes, road ecological routes, etc.

[0078] Step S110, using the insect pollination model of the target area, analyzing the insect pollination service capacity in the primary ecological landscape corridor, and obtaining the analysis result of the insect pollination service capacity.

[0079] Exemplarily, the insect pollination service capacity is affected by the intensity of insect activity. It is understandable that the intensity of insect activity is affected by external factors.

[0080] Exemplarily, in some embodiments, external influencing factors include but are not limited to one or more of the following: temperature, season, climate, etc.

[0081] Understandably, research in this field often uses the "three gullies" data as important reference data. Specifically, the ecological background circle is composed of three gullies: Changchengling Valley, Wanlong Valley and Taizicheng Valley. Due to the large difference in elevation among the three valleys, the difference in elevation between the highest and lowest points is 949m.

[0082] Specifically, elevation refers to the distance from a point along the plumb line to the absolute base surface, which is called absolute elevation, or elevation for short.

[0083] It is understandable that under environmental conditions at different altitudes, temperature is one of the most important environmental factors affecting the biodiversity pattern. It is understood in the art that the temperature will decrease by about 0.6°C for every 100m increase in altitude.

[0084] For example, in the same season, the maximum temperature difference between different spaces in the ecological landscape corridor is 5.70°C.

[0085] Furthermore, since the activities of pollinating insects are highly sensitive to temperature, and further, the temperature differences in the spatial level will affect the activities of pollinating insects, an elevation correction activity intensity parameter is introduced to correct the impact of elevation on temperature and thus affect the pollination activities of pollinating insects.

[0086] In some embodiments, the insect pollination model is constructed by the following steps:

[0087] Step 1: Calculate the temperature corresponding to the elevation of the habitat grid using the following formula:

[0088] , (Formula 7)

[0089] In formula 7:

[0090] Used to characterize habitat grids The elevation of your location Corresponding temperature , Used to represent the average elevation, Used to characterize habitat grids The elevation of the location, Used to represent the average temperature;

[0091] Step 2: Calculate the average pollinator body temperature of the habitat grid using the following formula:

[0092] , (Equation 8)

[0093] In formula 8:

[0094] Used to characterize habitat grids Average pollinator body temperature, R Used to characterize solar irradiance;

[0095] Step 3: Calculate the activity intensity of pollinators in the habitat grid for the month using the following formula:

[0096] , (Formula 9)

[0097] In formula 9:

[0098] Used to characterize the activity intensity of pollinating insects in that month.

[0099] Step 4: Calculate the relative activity intensity of pollinating insects on a monthly basis using the following formula:

[0100] , (Formula 10)

[0101] In formula 10:

[0102] Used to characterize habitat grids The relative activity intensity of pollinating insects in the region from month to month, Used to characterize habitat grids Pollinators Monthly activity intensity The maximum value of

[0103] Step 5: Calculate the elevation-corrected activity intensity of pollinators using the following formula:

[0104] , (Formula 11)

[0105] In formula 11:

[0106] Pollinators in habitat grids The intensity of the elevation correction activity within A raster representing the habitat in the study area Activity intensity corresponding to the lowest point in the interior .

[0107] The insect pollination model is constructed based on the activity intensity of the month and the elevation-corrected activity intensity. The insect pollination model is as follows:

[0108] )

[0109] in, Used to characterize the insect pollination service capacity index, Used to characterize habitats The intensity of the activity within the elevation correction; Used to characterize habitat grids Species within exist Monthly activity intensity, The number of reference months used to characterize the capacity of insect pollination services.

[0110] The relationship between pollinating insects and ecological landscape corridors is obtained according to the insect pollination model obtained in step S110, and the pollination ability analysis results of pollinating insects in the ecological landscape corridors are obtained according to the insect pollination model, and the pollination ability analysis results are used as important reference conditions for constructing ecological landscape corridors.

[0111] Step S112, optimizing the initial ecological landscape corridor according to the analysis result to obtain the ecological landscape corridor.

[0112] The insect pollination service capacity of the primary ecological landscape corridor is analyzed through the insect pollination model, and the primary ecological landscape corridor is further optimized to finally obtain the ecological landscape corridor.

[0113] In the embodiments of the present application, high habitat quality areas are determined through sample plot data of the target area, and an ecological landscape corridor is constructed in combination with an insect pollination model. On the basis of solving the problem of habitat fragmentation and maintaining the level of biodiversity, the present application further improves the service performance of insect pollination in biological diversity and protecting the health of the ecosystem, so as to enable the sustainable development of the ecosystem.

[0114] See also Figure 2 , shown is a schematic diagram of the structure of an electronic device according to an embodiment of the present application. The specific embodiment of the present application does not limit the specific implementation of the electronic device.

[0115] like Figure 2 As shown, the electronic device may include: a processor (processor) 202 , a communication interface (Communications Interface) 204 , a memory (memory) 206 , and a communication bus 208 .

[0116] in:

[0117] The processor 202 , the communication interface 204 , and the memory 206 communicate with each other via the communication bus 208 .

[0118] The communication interface 204 is used to communicate with other electronic devices or servers.

[0119] The processor 202 is used to execute the program 210, and specifically can execute the relevant steps in the above-mentioned ecological landscape corridor construction method embodiment.

[0120] Specifically, the program 210 may include program codes, which include computer operation instructions.

[0121] The processor 202 may be a CPU, or an application specific integrated circuit ASIC (Application Specific Integrated Circuit), or one or more integrated circuits configured to implement the embodiments of the present application. The one or more processors included in the smart device may be processors of the same type, such as one or more CPUs; or processors of different types, such as one or more CPUs and one or more ASICs.

[0122] The memory 206 is used to store the program 210. The memory 206 may include a high-speed RAM memory, and may also include a non-volatile memory (non-volatile memory), such as at least one disk memory.

[0123] The program 210 may include multiple computer instructions. Specifically, the program 210 may enable the processor 202 to execute operations corresponding to the method for constructing an ecological landscape corridor described in any of the aforementioned method embodiments through the multiple computer instructions.

[0124] The specific implementation of each step in the program 210 can refer to the corresponding description of the corresponding steps and units in the above method embodiment, and has corresponding beneficial effects, which will not be repeated here. Those skilled in the art can clearly understand that for the convenience and simplicity of description, the specific working process of the above-described devices and modules can refer to the corresponding process description in the above method embodiment, which will not be repeated here.

[0125] Through the electronic device of the embodiment of the present application, the high habitat quality area is determined through the sample data of the target area, and the ecological landscape corridor is constructed in combination with the insect pollination model. On the basis of solving the habitat fragmentation problem and maintaining the level of biodiversity, the present application further improves the service performance of insect pollination in biological diversity and protecting the health of the ecosystem, so as to enable the sustainable development of the ecosystem.

[0126] The present application also provides a computer-readable storage medium storing instructions for causing a machine to execute the method for constructing an ecological landscape corridor as described herein. Specifically, a system or device equipped with a storage medium may be provided, on which a software program code for implementing the functions of any of the above embodiments is stored, and a computer (or CPU or MPU) of the system or device reads and executes the program code stored in the storage medium.

[0127] In this case, the program code read from the storage medium itself can implement the function of any one of the above embodiments, so the program code and the storage medium storing the program code constitute part of the present application.

[0128] The storage medium embodiments for providing the program code include a floppy disk, a hard disk, a magneto-optical disk, an optical disk (such as CD-ROM, CD-R, CD-RW, DVD-ROM, DVD-RAM, DVD-RW, DVD+RW), a magnetic tape, a non-volatile memory card, and a ROM. Alternatively, the program code can be downloaded from a server computer via a communication network.

[0129] An embodiment of the present application also provides a computer program product, including computer instructions, which instruct a computing device to perform any corresponding operation in the above-mentioned multiple method embodiments.

[0130] It should be pointed out that, according to the needs of implementation, the various components / steps described in the embodiments of the present application can be split into more components / steps, or two or more components / steps or partial operations of components / steps can be combined into new components / steps to achieve the purpose of the embodiments of the present application.

[0131] The above-mentioned method according to the embodiment of the present application can be implemented in hardware, firmware, or implemented as software or computer code that can be stored in a recording medium (such as a CD ROM, RAM, floppy disk, hard disk or magneto-optical disk), or implemented as a computer code originally stored in a remote recording medium or a non-temporary machine-readable medium downloaded through a network and to be stored in a local recording medium, so that the method described herein can be stored in such software processing on a recording medium using a general-purpose computer, a dedicated processor or programmable or dedicated hardware (such as an ASIC or FPGA). It can be understood that a computer, a processor, a microprocessor controller or programmable hardware includes a storage component (e.g., RAM, ROM, flash memory, etc.) that can store or receive software or computer code, and when the software or computer code is accessed and executed by a computer, a processor or hardware, the navigation method described herein is implemented. In addition, when a general-purpose computer accesses the code for implementing the navigation method shown herein, the execution of the code converts the general-purpose computer into a dedicated computer for executing the navigation method shown herein.

[0132] Those of ordinary skill in the art will appreciate that the units and method steps of each example described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, or a combination of computer software and electronic hardware. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professional and technical personnel can use different methods to implement the described functions for each specific application, but such implementation should not be considered to be beyond the scope of the embodiments of the present application.

[0133] The above implementation methods are only used to illustrate the embodiments of the present application, and are not limitations on the embodiments of the present application. Ordinary technicians in the relevant technical field can make various changes and modifications without departing from the spirit and scope of the embodiments of the present application. Therefore, all equivalent technical solutions also belong to the scope of the embodiments of the present application. The scope of patent protection of the embodiments of the present application should be limited by the claims.

Claims

1. A method for constructing an ecological landscape corridor, characterized in that: include: Obtain sample data of the target area; Calculating the habitat quality index of the target area through the sample data; determining a high habitat quality area from the target area according to the habitat quality index; According to the sample data, the resistance surface of the ith high habitat quality area is calculated by the following formula: R i The resistance surface used to characterize the i-th high habitat quality area, S i Used to characterize the terrain slope included in the i-th sample plot data, α i Used to characterize the terrain slope S i The weight coefficient, V i Used to characterize the vegetation coverage rate included in the i-th sample data, β i Used to characterize vegetation coverage V i The weight coefficient, P i It is used to characterize the land cover ratio included in the i-th sample data, γ i Used to characterize land cover P i The weight coefficient, E i It is used to characterize the ecological importance of the i-th sample data, δ i Used to characterize ecological importance i The weight coefficient of According to the resistance surface of the i-th high habitat quality area, the minimum cumulative resistance value is calculated by the following formula: MCR is used to characterize the minimum cumulative resistance value, D ij It is used to characterize the distance from the ecological source to the i-th high habitat quality area, f min It is used to characterize the minimum cumulative resistance value MCR and the distance D from the ecological source to the i-th high habitat quality area. ij The positive correlation function between According to the minimum cumulative resistance value, constructing a first generation ecological landscape corridor; The insect pollination model of the target area is used to analyze the insect pollination service capacity in the primary ecological landscape corridor, and the analysis results of the insect pollination service capacity are obtained; wherein the insect pollination model is shown in the following formula: C is used to represent the insect pollination service capacity index; A hx It is used to characterize the elevation-corrected activity intensity of pollinating insects in the habitat grid x. The elevation-corrected activity intensity is calculated using the following formula: A cj Used to represent the relative activity intensity of pollinators in habitat grid x from month to month, A jm It is used to characterize the activity intensity A corresponding to the lowest point in the habitat grid m within the study area. j ; A xjλ It is used to characterize the activity intensity of species λ in habitat grid x in month j. The activity intensity of pollinating insects in the habitat grid in that month is calculated by the following formula: j =-39.3+4.01T bx , A j Used to characterize the activity intensity of pollinating insects in a month, T bx The average pollinator body temperature used to characterize habitat grid x; n is used to characterize the reference month for the analysis of insect pollination service capacity; The initial ecological landscape corridor is optimized according to the analysis results to obtain an ecological landscape corridor.

2. The method according to claim 1, characterized in that The step of calculating the habitat quality index of the target area through the sample data includes: Map the threat sources onto the land data grid and obtain the distance between the habitat grid and the threat source grid; Calculate the influence of the distance between the habitat grid and the threat source grid on the habitat formation; According to the impact degree value, the threat level of the threat source is calculated by the following formula: i rxy Used to characterize the impact value, D xj is used to characterize the threat level of the threat source to the habitat grid x in the jth habitat type, R is used to characterize the number of threat factors included in the threat source, and w r The weight used to characterize the rth threat factor, Y r The number of grids occupied by the threat factor layer corresponding to the rth threat factor on the land data grid, r y It is used to characterize the number of threat factors of the corresponding threat factor layer occupying the yth grid on the land data grid, β x The coefficient used to characterize the accessibility level of the habitat grid x, S jr It is used to characterize the sensitivity coefficient of the jth habitat type to the threat factor r; According to the threat level of the threat source, the habitat quality index is calculated by the following formula: Q xj The habitat quality index used to characterize the habitat grid x in the jth habitat type, H j It is used to characterize the habitat suitability coefficient of the jth habitat type, z is used to characterize the scale constant, and k is used to characterize the half-saturation constant.

3. The method according to claim 2, characterized in that The calculation of the influence degree value of the distance between the habitat grid and the threat source grid on the habitat formation includes: The influence degree value of the distance between the habitat grid and the threat source grid on the habitat formation is calculated by the following formula 1 or formula 2; Formula 1 includes: Formula 2 includes: d xy It is used to represent the distance between the habitat grid x and the threat source grid y, d rmax It is used to characterize the maximum impact distance between the threat factor r included in the threat source and the habitat grid x.

4. The method according to claim 1, characterized in that: The temperature corresponding to the elevation of the habitat grid x is calculated using the following formula: T x =0.006(H0-H x )+T0 Among them, T x The height H used to represent the location of the habitat grid x x The corresponding temperature T x , H0 is used to characterize the average elevation, H x It is used to represent the elevation of the location of the habitat grid x, and T0 is used to represent the average temperature; The average pollinator body temperature of habitat grid x is calculated by the following formula: T bx =-0.62+1.027T x +0.006R Among them, T bx is used to characterize the average pollinator body temperature of the habitat grid x, and R is used to characterize the solar irradiance.

5. The method according to claim 4, characterized in that The relative activity intensity of pollinating insects on a monthly basis is calculated using the following formula: Among them, A cj It is used to characterize the relative activity intensity of pollinating insects in the habitat grid x from month to month, A j0 A is used to characterize the activity intensity A of pollinators in the habitat grid x in the jth month. j The maximum value of .

6. An electronic device comprising: A processor, a memory, a communication interface and a communication bus, wherein the processor, the memory and the communication interface communicate with each other via the communication bus; The memory is used to store at least one executable instruction, and the executable instruction enables the processor to perform an operation corresponding to the method according to any one of claims 1 to 5.

7. A computer storage medium having a computer program stored thereon, wherein when the program is executed by a processor, the method according to any one of claims 1 to 5 is implemented.

Images