CN114580972B

CN114580972B - Ecological corridor construction method and device based on multi-element ecological source area

Info

Publication number: CN114580972B
Application number: CN202210458402.0A
Authority: CN
Inventors: 王健; 田园; 刘佳旭
Original assignee: Beijing Beilin Eco Pastoral Landscape Planning And Design Co ltd
Current assignee: Beijing Beilin Eco Pastoral Landscape Planning And Design Co ltd
Priority date: 2022-04-28
Filing date: 2022-04-28
Publication date: 2022-07-29
Anticipated expiration: 2042-04-28
Also published as: CN114580972A

Abstract

The invention relates to an ecological corridor construction method and device based on a multi-element ecological source area, relating to the field of landscape ecology and urban ecological planning, wherein the method comprises the following steps: obtaining multiband remote sensing data, underlying surface data, precipitation data, statistical data and night light index of a target area, performing index calculation to obtain a normalized vegetation index, a topographic relief degree, a target water source conservation factor, a carbon fixing factor and an oxygen releasing factor, then performing ecological source identification to obtain a grain supply source area, a water source conservation source area, a habitat maintenance source area, a carbon fixing and oxygen releasing source area and a soil maintenance source area, performing union operation to the obtained areas to obtain a comprehensive ecological source area distribution area as a source calculation ecological resistance cost value, and calculating an ecological corridor among the source areas for a target based on the minimum cost value. Therefore, the problem of reasonable protection and utilization of the urban ecological space is solved, a scientific analysis and demarcation method is provided for reservation and protection of the urban ecological space, and harmonious development of urban expansion and ecological protection is realized.

Description

Ecological corridor construction method and device based on multi-element ecological source area

Technical Field

The disclosure relates to the field of landscape ecology and urban ecological planning in the technical field of environmental science, in particular to an ecological corridor construction method and device based on a multi-element ecological source area.

Background

At present, the accelerated urbanization process promotes social progress and economic development, causes damage to the living environment, aggravates landscape fragmentation of the living habitat, and causes reduction of biological diversity and overall function degradation of an ecological system due to excessive fragmentation, artificial transformation and island transformation of the natural habitat, so that contradiction between human and nature is increasingly acute. Therefore, how to enhance the reasonable protection and utilization of ecological space in urban areas with scarce land resources to meet the multi-ecological and humanistic requirements of citizens and realize the ecological, flexible and harmonious transformation of cities becomes a new challenge for the sustainable development of cities.

Therefore, the concept of the ecological source and the ecological corridor comes into play, and it is very important how to identify the spatial distribution positions of the ecological source and the ecological corridor, and to scientifically protect and ecologically utilize the spatial distribution positions to ensure the urban ecological safety pattern in view of the special significance of the ecological source and the ecological corridor.

Disclosure of Invention

In order to solve the technical problems or at least partially solve the technical problems, the present disclosure provides an ecological corridor construction method and apparatus based on a multi-element ecological source.

The invention provides an ecological corridor construction method based on a multi-element ecological source area, which comprises the following steps:

obtaining multiband remote sensing data, underlying surface data, precipitation data, statistical data and night light index of a target area;

performing index calculation based on the multiband remote sensing data, the underlying surface data, the rainfall data, the statistical data and the night light index to obtain a normalized vegetation index, a topographic relief degree, a target water source conservation factor, a carbon fixation factor and an oxygen release factor;

performing ecological source land identification based on the normalized vegetation index, the topographic relief degree, the target water source conservation factor, the carbon-fixing factor and the oxygen-releasing factor to obtain a grain supply source land, a water source conservation source land, a habitat maintenance source land, a carbon-fixing oxygen-releasing source land and a soil conservation source land;

performing union operation on the grain supply source land, the water source conservation source land, the habitat maintenance source land, the carbon-fixing oxygen-releasing source land and the soil conservation source land to obtain a comprehensive ecological source land distribution area;

acquiring an ecological resistance value surface area layer corresponding to the comprehensive ecological source area distribution area, and calculating a cost consumption value based on the ecological resistance value surface area layer;

and acquiring the ecological corridor corresponding to the minimum cost value in the cost values as a target ecological corridor.

Optionally, the performing index calculation based on the multiband remote sensing data, the underlying surface data, the precipitation data, the statistical data, and the night light index to obtain a normalized vegetation index, a topographic relief degree, a target water source conservation factor, a carbon sequestration factor, and an oxygen release factor includes:

acquiring a near-infrared band gray value and an infrared band gray value of each grid based on the multiband remote sensing data, acquiring a difference value and a summation value between the near-infrared band gray value and the infrared band gray value, and calculating a ratio of the difference value and the summation value to obtain the normalized vegetation index of each grid;

acquiring a maximum elevation value and a minimum elevation value of each grid target distance in unit area based on the underlying surface data, and calculating a difference value between the maximum elevation value and the minimum elevation value to obtain the topographic relief degree of each grid;

acquiring a relief correction factor, acquiring the grid unit precipitation amount and the empirical runoff coefficient of each grid based on the precipitation data, and calculating based on the grid unit precipitation amount, the relief correction factor and the empirical runoff coefficient to obtain the target water source conservation factor;

calculating based on the normalized vegetation index of each grid and the average value of the carbon fixing rate of each grid to obtain the carbon fixing factor;

and calculating based on the normalized vegetation index of each grid and the average value of the oxygen release rate of each grid to obtain the oxygen release factor.

Optionally, the identifying an ecological source area based on the normalized vegetation index, the topographic relief degree, the target water source conserving factor, the carbon fixation factor and the oxygen release factor to obtain a grain supply source area, a water source conserving source area, a habitat maintenance source area, a carbon fixation and oxygen release source area and a soil conservation source area comprises:

carrying out average calculation on all normalized vegetation indexes in a farmland area to obtain a normalized vegetation index average value, carrying out calculation on the normalized vegetation index of each grid, the normalized vegetation index average value and the grain yield average value to obtain the grain yield of each grid, sequencing the grain yields of each grid from large to small on the basis of the grain yield of each grid, and obtaining grids with preset first values before sequencing as the grain supply source areas; wherein the total area of the grids with the preset first numerical value is larger than a preset first area threshold value;

sequencing the target water source conservation factors from large to small based on each grid, and acquiring grids with preset second values before sequencing as the water source conservation source; wherein the total area of the grid with the preset second numerical value is larger than a preset second area threshold;

acquiring a water source distance, a land condition and a weighting factor corresponding to the topographic relief degree, calculating based on the water source distance, the land condition, the topographic relief degree and the corresponding weighting factor to obtain the habitat maintenance values, sorting the habitat maintenance values from large to small based on the habitat maintenance values of each grid, and acquiring a grid with a preset third value before sorting as the habitat maintenance source; wherein the total grid area with the preset third numerical value is larger than a preset third area threshold;

respectively summing the carbon fixation factor and the oxygen release factor of each grid to obtain a comprehensive carbon fixation and oxygen release rate of each grid, and sequencing the comprehensive carbon fixation and oxygen release rates of each grid from large to small to obtain a grid with a preset fourth value before sequencing as the carbon fixation and oxygen release source; wherein the total area of the grid with the preset fourth numerical value is larger than a preset fourth area threshold;

acquiring rainfall erosion force and underlying surface data based on the rainfall data to acquire vegetation coverage factors, acquiring soil texture based on a preset soil database, analyzing based on the rainfall erosion force, the soil texture, the terrain relief and the vegetation coverage factors to obtain the soil sensitivity level of each grid, sorting the soil sensitivity levels of each grid from low to high based on the soil sensitivity level of each grid, and acquiring grids with a preset fifth value before sorting as the soil conservation source; and the total area of the grid with the preset fifth numerical value is larger than a preset fifth area threshold value.

Optionally, the calculating a cost value based on the ecological resistance value surface map layer includes:

acquiring a field distance and a target ecological resistance value based on the ecological resistance value surface area layer;

and calculating the field distance and the target ecological resistance value based on a preset formula to obtain the cost consumption value.

Optionally, the obtaining a target ecological resistance value based on the ecological resistance value surface area map layer includes:

acquiring a normalized vegetation index, an average value of ecological resistance equivalent and an average value of non-ecological space resistance equivalent of each ground cover type in the ecological resistance value surface area map layer;

carrying out average calculation based on the normalized vegetation indexes of the surface coverage types of the ecological space to obtain the average value of the corresponding normalized vegetation indexes;

calculating based on the normalized vegetation index, the ecological resistance equivalent average value and the normalized vegetation index average value of each surface coverage type of the ecological space to obtain a first ecological resistance value;

acquiring night light index values of surface coating types of each non-ecological space based on the night light index;

carrying out average calculation based on the night light indexes of the surface covering types of the non-ecological space to obtain corresponding average values of the night light indexes;

calculating based on the night light index numerical values of the surface covering types of the non-ecological space, the average value of the resistance equivalent of the non-ecological space and the average value of the night light index to obtain a second ecological resistance value;

and calculating based on the first ecological resistance value and the second ecological resistance value to obtain the target ecological resistance value.

The utility model provides an ecological corridor founds device based on many first ecological source ground, includes:

the data acquisition module is used for acquiring multiband remote sensing data, underlying surface data, precipitation data, statistical data and night light indexes of a target area;

the index calculation module is used for performing index calculation on the basis of the multiband remote sensing data, the underlying surface data, the precipitation data, the statistical data and the night light index to obtain a normalized vegetation index, a topographic relief degree, a target water source conservation factor, a carbon sequestration factor and an oxygen release factor;

a source area identification module for performing ecological source area identification based on the normalized vegetation index, the topographic relief degree, the target water source conservation factor, the carbon fixation factor and the oxygen release factor to obtain a grain supply source area, a water source conservation source area, a habitat maintenance source area, a carbon fixation and oxygen release source area and a soil maintenance source area;

the union calculation module is used for carrying out union operation on the grain supply source area, the water source conservation source area, the habitat maintenance source area, the carbon-fixing oxygen-releasing source area and the soil maintenance source area to obtain a comprehensive ecological source area distribution area;

the acquisition module is used for acquiring an ecological resistance value surface area layer corresponding to the comprehensive ecological source area distribution area;

the cost calculation module is used for calculating a cost consumption value based on the ecological resistance value surface area layer;

and the ecological corridor determining module is used for acquiring the ecological corridor corresponding to the minimum cost value in the cost values as a target ecological corridor.

Optionally, the index calculating module is specifically configured to:

Optionally, the source identification module is specifically configured to:

The present disclosure provides an electronic device, including: a processor; a memory for storing the processor-executable instructions; the processor is used for reading the executable instructions from the memory and executing the instructions to realize the ecological corridor construction method based on the multi-element ecological source area.

The present disclosure provides a computer-readable storage medium, wherein the storage medium stores a computer program for executing the method for constructing an ecological corridor based on multiple ecological sources according to the above embodiments.

Compared with the prior art, the technical scheme provided by the embodiment of the disclosure has the following advantages:

acquiring multiband remote sensing data, underlying surface data, precipitation data, statistical data and night light index of a target area, performing index calculation based on the multiband remote sensing data, the underlying surface data, the precipitation data, the statistical data and the night light index to obtain a normalized vegetation index, a topographic relief degree, a target water source conservation factor, a carbon fixation factor and an oxygen release factor, performing ecological source area identification based on the normalized vegetation index, the topographic relief degree, the target water source conservation factor, the carbon fixation factor and the oxygen release factor to obtain a grain supply source area, a water source conservation source area, a habitat maintenance source area, a carbon fixation and oxygen release source area and a soil conservation source area, performing union operation on the grain supply source area, the water source conservation source area, the habitat maintenance source area, the carbon fixation and oxygen release source area and the soil conservation source area to obtain a comprehensive ecological source area distribution area, and acquiring an ecological resistance value area map layer corresponding to the comprehensive ecological source area distribution area, and calculating a cost value based on the ecological resistance value surface area layer, and acquiring an ecological corridor corresponding to the minimum cost value in the cost values as a target ecological corridor. Therefore, a scientific analysis and demarcation method is provided for urban ecological space reservation and protection, and harmonious development of urban expansion and ecological protection is realized.

Drawings

The above and other features, advantages, and aspects of embodiments of the present disclosure will become more apparent by referring to the following detailed description when taken in conjunction with the accompanying drawings. Throughout the drawings, the same or similar reference numbers refer to the same or similar elements. It should be understood that the drawings are schematic and that elements and features are not necessarily drawn to scale.

Fig. 1 is a schematic flow chart of an ecological corridor construction method based on a multi-element ecological source area according to an embodiment of the present disclosure;

fig. 2 is a schematic flow chart of another method for constructing an ecological corridor based on a multi-element ecological source provided by the embodiment of the disclosure;

fig. 3 is a schematic view of an ecological source identification process based on quantitative data analysis and calculation according to an embodiment of the present disclosure;

fig. 4 is a schematic structural diagram of an ecological corridor construction device based on a multi-element ecological source area provided by an embodiment of the disclosure.

Detailed Description

In order that the above objects, features and advantages of the present disclosure may be more clearly understood, aspects of the present disclosure will be further described below. It should be noted that the embodiments and features of the embodiments of the present disclosure may be combined with each other without conflict.

In the following description, numerous specific details are set forth in order to provide a thorough understanding of the present disclosure, but the present disclosure may be practiced in other ways than those described herein; it is to be understood that the embodiments disclosed in the specification are only a few embodiments of the present disclosure, and not all embodiments.

In general, ecological sources and ecological galleries are composed of green spaces with significant ecological functions, such as grasslands, wetlands, greenlands, woodlands, cultivated lands, gardens, water areas, etc. The ecological source land can effectively provide good foraging, inhabitation and breeding places for animals and plants, and can also play ecological functions of conserving water sources, regulating microclimate, eliminating pollutants and the like. The ecological corridor is used as a linear or strip-shaped ecological space with certain physical width and spatial continuity, and can provide high-quality foraging and migration channels for wild animals, promote the transfer and interaction of material flow, information flow and energy flow among multiple ecological systems, and realize the increase of biological diversity; in addition, the ecological corridor can also provide excellent leisure and rest space for urban residents.

Under the high-pressure driving of harmonious and sustainable development of urban expansion and ecological environment, the embodiment of the disclosure provides an ecological corridor construction method based on a multi-element ecological source area, provides a scientific analysis and demarcation method for urban ecological space reservation and protection, and realizes harmonious development of urban expansion and ecological protection.

Specifically, fig. 1 is a schematic flow chart of an ecological corridor construction method based on a multivariate ecological source area according to an embodiment of the present disclosure, which includes:

step 101, obtaining multiband remote sensing data, underlying surface data, precipitation data, statistical data and night light index of a target area.

The target area can be selected according to application scene requirements, such as XX province YY city or AA city BB area.

The multiband remote sensing data comprise near infrared band gray values, infrared band gray values and the like; underlying surface data such as land use, surface covering, terrain elevation, and the like; precipitation data such as precipitation volume, precipitation frequency, etc. for sites within the target area; statistical data such as average grain yield, average carbon fixation and oxygen release and the like facilitate subsequent spatial discretization; the night light index is, for example, night light index data obtained by using a Visible Infrared Imaging Radiometer Suite (VIIRS), and the data is subsequently used to perform spatial discretization processing on the ecological resistance value of the non-ecological space land.

And 102, performing index calculation based on the multiband remote sensing data, the underlying surface data, the precipitation data, the statistical data and the night light index to obtain a normalized vegetation index, a topographic relief degree, a target water source conservation factor, a carbon fixation factor and an oxygen release factor.

In the embodiment of the present disclosure, the calculation of the relevant digital indexes of the target area by using the mathematical model and using the collected multiband remote sensing data, underlying surface data, precipitation data, statistical data, night light index and other data of the target area includes: normalized Difference Vegetation Index (NDVI), topographic relief, target water source conservation factor, carbon sequestration factor, oxygen release factor, etc.

In the embodiment of the disclosure, the near-infrared band gray value and the infrared band gray value of each grid are obtained based on multi-band remote sensing data, the difference value and the sum value between the near-infrared band gray value and the infrared band gray value are obtained, the ratio of the difference value and the sum value is calculated, and the normalized vegetation index of each grid is obtained.

Specifically, the normalized vegetation index NDVI = (NIR-R) ÷ (NIR + R), where NIR is a near-infrared band gray scale value and R is an infrared band gray scale value.

In the embodiment of the disclosure, the maximum elevation value and the minimum elevation value of each grid target distance in unit area are obtained based on the underlying surface data, and the difference between the maximum elevation value and the minimum elevation value is calculated to obtain the topographic relief of each grid. The target distance can be set according to the application scene requirements, for example, the target distance is 1km in a target grid (the search window area is 1 km) ² ) Assigning the maximum elevation in the center area to a current grid in the center area; according to 1km ² Window (2)And range, traversing the universe, and assigning the most value (size and magnitude) to the corresponding center grid.

Specifically, relief degree D = H _max -H _min In the formula, H _max Is the maximum elevation per unit area, H _min Is the minimum elevation value in unit area.

In the embodiment of the disclosure, the waviness correction factor is obtained, the precipitation amount of the grid unit, the waviness correction factor and the empirical runoff coefficient of each grid are obtained based on precipitation data, and calculation is performed based on the precipitation amount of the grid unit and the empirical runoff coefficient to obtain the target water source conservation factor.

Specifically, the target water source conservation factor R = Px (1-alpha) xq, wherein P is the precipitation amount of the grid unit, alpha is the empirical runoff coefficient of different surface coatings, and q is the fluctuation correction factor.

In the embodiment of the disclosure, the carbon fixation factor is obtained by calculating based on the normalized vegetation index of each grid and the average value of the carbon fixation rate of each grid.

In particular, carbon sequestration factor C _i = NDVI _i ÷NDVI _{i are all} ×C _{i are all} In the formula, NDVI _i NDVI value, NDVI, for the ith formation overburden type _{i are all} Average value of NDVI as the ith formation overburden type, C _{i are all} The average value of the carbon fixing rate of the ith formation coating type is shown.

In the disclosed embodiment, the oxygen release factor is obtained by performing a calculation based on the normalized vegetation index of each grid and the average value of the oxygen release rate of each grid.

In particular, oxygen release factor O _i = NDVI _i ÷NDVI _{i are all} ×O _{i are all} In the formula, O _{i are all} The average value of the oxygen release rate of the ith formation coating type is shown.

And 103, identifying an ecological source area based on the normalized vegetation index, the topographic relief degree, the target water source conservation factor, the carbon fixation factor and the oxygen release factor to obtain a grain supply source area, a water source conservation source area, a habitat maintenance source area, a carbon fixation and oxygen release source area and a soil maintenance source area.

In the embodiment of the present disclosure, identifying an ecological source by discretization, weighting, and the like based on index data includes: a grain supply source place, a water source conservation source place, a habitat maintenance source place, a carbon-fixing oxygen-releasing source place and a soil maintenance source place.

In the embodiment of the disclosure, average calculation is performed based on normalized vegetation indexes of all farmlands in a target area to obtain a normalized vegetation index average value, calculation is performed based on the normalized vegetation index and the normalized vegetation index average value of each grid and a grain yield average value to obtain a grain yield of each grid, the grain yields of each grid are sorted from large to small based on the grain yields of each grid, and a grid with a preset first value before sorting is obtained as a grain supply source; the total area of the grid with the preset first numerical value is larger than a preset first area threshold value. Therefore, effective identification of a high-yield region of a farmland can be accurately realized through a scientific NDVI-based grid discretization method, and a grain supply source place is further scientifically determined.

Specifically, the grain supply source area identification means that rasterized grain yield data is generated based on the average grain yield per unit area of a target area (for example, if the average grain yield of the target area is 500kg/ha, all arable lands are assigned with 500), and discretization is performed on the rasterized NDVI (one value of 500 is spatially discretized into different values of 50-900), and then, for example, the first 20% of high-yield areas and land blocks with the whole area larger than 50ha are selected as the grain supply source areas. The discretization method comprises the following steps: rasterized grain yield L _i = NDVI _Li ÷NDVI _{L is all} ×L _{Are all made of} In the formula, NDVI _Li For rasterised NDVI values of cultivated land, NDVI _{L is all} Is the average value of NDVI, L, of cultivated land _{Are all made of} Is the average grain yield.

In the embodiment of the disclosure, the target water source conservation factors based on each grid are sorted from large to small, and the grid with a preset second numerical value before sorting is obtained as a water source conservation source; and the total area of the grid with the preset second numerical value is larger than a preset second area threshold value. Therefore, the regional water source conservation quantity can be calculated simply and quickly, and a large amount of time can be saved compared with the existing model method; and the result is corrected through the topographic relief degree, so that the accuracy and precision of calculation can be improved.

Specifically, the identification of the water source conservation source area refers to that a land block with the whole area larger than 50ha, which is the first 20% of a high water source conservation area, is selected as a source area after the groundwater is supplied by calculation through an empirical runoff coefficient method and the topographic relief degree is adopted for correcting the water source conservation area based on a water balance principle.

More specifically, in order to simply and rapidly calculate the water source conservation quantity in the region and save a large amount of time, the accuracy and precision of calculation can be improved by correcting the result through the topographic relief degree. R = P x (1-alpha) x q, wherein P is the precipitation of the grid unit, alpha is the empirical runoff coefficient of different surface coatings, and q is the undulation correction factor.

In the embodiment of the disclosure, weight factors corresponding to a water source distance, a land condition and a topographic relief degree are obtained, calculation is performed based on the water source distance, the land condition and the topographic relief degree and the corresponding weight factors to obtain habitat maintenance values, sorting is performed from large to small based on the habitat maintenance values of each grid, and a grid with a preset third value before sorting is obtained as a habitat maintenance source; and the total area of the grid with the preset third numerical value is larger than a preset third area threshold value. Therefore, an index system for identifying the habitat maintenance source is scientifically, reasonably and comprehensively constructed, and the quality grade state of the animal habitat can be reflected.

Specifically, the habitat maintenance source identification means that four habitat maintenance factors such as a distance from a water source, a topographic condition, a distance from a construction site and a covertness (a geological condition) of a ground surface covering are integrated to perform weighted suitability analysis (for example, the weights are 0.3, 0.15, 0.25 and 0.3 in sequence) so as to reflect the habitat quality level of a space, for example, a land block with the high suitability area of the first 20% and the whole area larger than 50ha is selected as a habitat maintenance source.

In the embodiment of the disclosure, the carbon fixation factor and the oxygen release factor of each grid are summed respectively to obtain the comprehensive carbon fixation and oxygen release rate of each grid, and the comprehensive carbon fixation and oxygen release rates of each grid are sorted from large to small based on the comprehensive carbon fixation and oxygen release rate of each grid, so that the grid with a preset fourth value before sorting is obtained as a carbon fixation and oxygen release source; and the total area of the grid with the preset fourth numerical value is larger than a preset fourth area threshold value. Therefore, the ecological environmental benefits of carbon fixation and oxygen release are comprehensively considered, and the atmospheric environmental effect of the vegetation is comprehensively and quantitatively expressed.

Specifically, the carbon-fixing oxygen-releasing source identification refers to that on the basis of average carbon-fixing and oxygen-releasing rates of different earth surface vegetation types, discrete processing is respectively carried out on the average carbon-fixing and oxygen-releasing rates by utilizing NDVI of various land types, and after the sum of the comprehensive rates of carbon-fixing and oxygen-releasing, for example, a high-yield area of the first 20% and a land with the whole area larger than 50ha are selected as the carbon-fixing oxygen-releasing source.

Namely, the average value of the carbon fixation and oxygen release rates of the surface coating of the target area is obtained from the statistical data, and then discretization processing is carried out on the average value by combining the spatialized NDVI value (a spatial distribution map of the carbon fixation and oxygen release rates is obtained).

In the embodiment of the disclosure, rainfall erosion force and underlying surface data are acquired based on rainfall data to acquire vegetation coverage factors, soil texture is acquired based on a preset soil database, the rainfall erosion force, the soil texture, the topographic relief and the vegetation coverage factors are analyzed to obtain the soil sensitivity level of each grid, the soil sensitivity levels of each grid are sorted from low to high based on the soil sensitivity level of each grid, and the grid with a fifth value preset before sorting is acquired as a soil conservation source; and the total area of the grids with the preset fifth numerical value is larger than a preset fifth area threshold value. Therefore, the method can be simplified into five-level sensitivity indexes from four aspects of erosive power, soil texture, undulation degree, vegetation coverage and the like, quantitative evaluation can be conveniently and rapidly carried out on soil erosion sensitivity, and convenience and operability of evaluation are greatly improved.

Specifically, the identification of the soil conservation source land refers to a simplified soil erosion risk model based on a preset soil loss equation, and after the rainfall erosion force, the soil texture, the topographic relief degree, the vegetation coverage factor and other four factors are subjected to integrated evaluation, for example, a land block with the low sensitive area of the first 20% and the whole area larger than 50ha is selected as the soil conservation source land.

In particular, soil retention factors

Wherein R is rainfall erosion force; s is the soil texture; d is the relief degree of the terrain; c is vegetation coverage; the specific grading equivalent and score are detailed in the following table 1:

TABLE 1

For example, the identification of the soil conservation source comprehensively considers four elements such as rainfall erosion force, soil texture, topographic relief degree and vegetation coverage, the grading value of each element is determined by a grading equivalent value assignment method, five grading value methods such as 1, 3, 5, 7 and 9 are adopted, and the grading threshold values of the rainfall erosion force are respectively 25, 100, 400 and 600; the soil texture threshold values are classified into gravel, sand, coarse sandy soil, fine sandy soil, clay, surface sandy soil, loam, sand loam, silt clay, soil clay, sand silt soil and silt soil in sequence; the topographic relief degree threshold values are respectively 20, 50, 100 and 300 m/km ² (ii) a The vegetation cover is divided mainly by the land utilization conditions, and the specific division threshold value land types are water bodies, swamps, rice fields, broad-leaved forests, coniferous forests, meadows, brush trees, sparse shrub grasslands, one-year-two-harvest, one-year-two-cropping-water-and-drought-and-desert, one-year-one-harvest and no vegetation.

It should be noted that the preset total area of the grid with the first value to the preset total area of the grid with the fifth value, and the preset threshold of the first area to the preset threshold of the fifth area may be selected and set according to an application scenario, and the disclosure is not limited specifically.

And 104, performing union operation on the grain supply source land, the water source conservation source land, the habitat maintenance source land, the carbon-fixing oxygen-releasing source land and the soil conservation source land to obtain a comprehensive ecological source land distribution area.

In the embodiments of the present disclosure, the integrated ecological source: based on five ecological functions of grain supply, water conservation, habitat maintenance, carbon fixation and oxygen release, soil maintenance and the like, a comprehensive ecological source distribution range is obtained through union operation, and diversified comprehensive ecological service functions are realized. That is, the generated 5 kinds of source region vector layers are subjected to union calculation to obtain a vector layer corresponding to the comprehensive ecological source region distribution area. The merging calculation is performed on the generated 5 source region vector layers to obtain the vector layer corresponding to the comprehensive ecological source region distribution area, which can be understood as performing spatial superposition on the 5 source region vector layers, merging different source region blocks with intersecting parts into a new large comprehensive unit, and directly bringing the source region unit blocks which are independent without intersection into the comprehensive ecological source region.

Therefore, the selected ecological source types are various and comprehensive, and the main ecological function characteristics of the natural environment can be reflected well, objectively and scientifically.

And 105, acquiring an ecological resistance value surface area layer corresponding to the comprehensive ecological source area distribution area, and calculating a cost value based on the ecological resistance value surface area layer.

And step 106, acquiring the ecological corridor corresponding to the minimum cost value in the cost values as a target ecological corridor.

The ecological resistance value surface area map layer refers to an ecological resistance map layer obtained by performing mean discretization processing by using a space discretization method based on a ground ecological resistance equivalent mean value raster file.

In the embodiment of the present disclosure, calculating the cost consumption value based on the surface area map layer of the ecological resistance value includes: and calculating the field distance and the target ecological resistance value based on a preset formula to obtain a cost consumption value. That is, the ecological source is equivalent to a circuit node, the non-ecological source is regarded as resistors with different resistance values, and the ecological corridor between the source and the ground is identified by adopting the minimum cumulative resistance model and the circuit model.

Therefore, the method for identifying the ecological corridor is scientific and reasonable, and the effective identification of the minimal resistance corridor for animal migration is realized by comprehensively superposing the field distance of the physical space and the resistance equivalent of the animal to the ground object perception.

In the embodiment of the disclosure, the normalized vegetation index, the average value of the ecological resistance equivalent and the average value of the non-ecological space resistance equivalent of each ground cover type in the ecological space in the ecological resistance value surface area image layer are obtained, the average calculation is performed based on the normalized vegetation index of each ground cover type in the ecological space to obtain the corresponding average value of the normalized vegetation index, the average value of the ecological resistance equivalent and the average value of the normalized vegetation index are calculated based on the normalized vegetation index, the average value of the ecological resistance equivalent and the average value of the normalized vegetation index of each ground cover type in the ecological space to obtain a first ecological resistance value, the night light index value of each ground cover type in the non-ecological space is obtained based on the night light index, the average calculation is performed based on the night light index of each ground cover type in the non-ecological space to obtain the corresponding average value of the night light index, the night light index value of each ground cover type in the non-ecological space, the average value of each ground cover type in the non-ecological space, And calculating the equivalent average value of the non-ecological space resistance and the average value of the light index at night to obtain a second ecological resistance value, and calculating based on the first ecological resistance value and the second ecological resistance value to obtain a target ecological resistance value.

Therefore, the NDVI and the night lamplight index are respectively selected for the ecological space and the non-ecological space to carry out space discretization processing on the ecological resistance value, the method is scientific and efficient, and the ecological resistance condition of the earth surface can be effectively reflected objectively.

Specifically, based on the landscape connectivity principle, the influence of any individual source area such as grain supply, water source conservation, habitat maintenance, carbon fixation and oxygen release, soil maintenance and the like on the connectivity of the whole landscape is evaluated, and quantitative expression is carried out by using a landscape connectivity importance index, so that the importance of the patches of each ecological source area is scientifically determined, and the importance of relevant ecological galleries is further facilitated to be judged.

Wherein, the importance of landscape connectivity refers to the importance of the source region patch to keep the connection with the ecological landscape, that is, the variation of the ecological landscape connectivity in the whole area broken (or removed) at this point, and the specific calculation formula is:

d _I (%)=100×（I-I _remove ）÷I；

wherein n represents the total number of plaques in the target region,

and

the areas of the patch i and the patch j are respectively represented,

is the total area of the target area and,

is the probability of the species spreading directly at plaque i and plaque j. I denotes the global index of all plaques in the landscape, I _remove Is the overall index value of the plaque remaining after removal of a single plaque. d _I Higher values indicate higher importance of the blob in landscape connectivity, also meaning d _I Higher value plaques are more prominent at the core of the target area.

Specifically, an ecological resistance equivalent method based on a land utilization type is adopted for assignment, NDVI and night light index are used for spatial discretization, and an ecological resistance accumulated cost grid map based on the ecological source space is constructed by integrating an ecological resistance value, an ecological source position, a spatial distance between the ecological source position and the source space and the like.

Specifically, rasterization assignment is performed based on the resistance equivalent mean value of the land utilization type (for example, different resistance value equivalents such as 1, 3, 50 and the like are respectively given to different land types such as grassland and forest land), wherein the ecological space (forest, grass, water, field, park, green land and the like) is subjected to space discretization processing based on the rasterized NDVI value; the non-ecological space (construction land, airport land, roads, home bases, rock gravel land, public utility land and the like) is subjected to spatial discretization processing based on the rasterized night light index value; and then merging the discretization grid image layers of different underlying surface ground conditions, and generating a space cost distance grid surface by combining the spatial position of the source ground, so as to obtain an ecological resistance accumulation cost grid image based on the position of the ecological source ground.

Specifically, the discretization method of the resistance value of the ecological space comprises the following steps: rasterized resistance equivalent F _{Raw i} = NDVI _{Raw i} ÷NDVI _{All of them are} ×T _{All of them are} In the formula, NDVI _{Raw i} NDVI value for the ith ground cover type of the ecological space, NDVI _{All of them are} Average value of NDVI as the ith formation overburden type, T _{All of them are} The ecological resistance equivalent average value of the ith ground coating type is shown.

The discretization method of the resistance value of the non-ecological space comprises the following steps: rasterized resistance equivalent F _{Non i} = G _{Non i} ÷G _{Is not all i} ×Y _{Is not all i} In the formula, G _{Non i} Night light index value of the ith non-ecological space ground surface coating type G _{Is not all i} Mean night light index, Y, for the ith non-ecological space surface covering type _{Is not all i} The ecological resistance equivalent average value of the ith ground coating type is shown.

Therefore, the space discretization method of the ecological resistance value can effectively realize the spatialization, differentiation and differentiation of the equivalent ecological resistance value.

Thereby accumulating the cost value

(i =1, 2, … n; j =1, 2, … m), wherein E is _i The field distance from a certain unit i in the space to the ecological source place; f _j The ecological resistance value of a certain landscape unit j in the space range; m _i And accumulating the cost value from the landscape unit i to the source. n is the total number of basic ecological landscape units.

The minimum cumulative cost-of-consumption path (i.e., the ecological corridor) can then be expressed by a mathematical model as:

；

wherein MCR represents a minimum cumulative cost-of-consuming path; min () means finding the minimum value; m _i And accumulating the cost value from the landscape unit i to the source.

The ecological corridor belongs to a communication type corridor, is used as a channel for ensuring the flow of energy and materials between ecological source areas, and is a key ecological land for promoting the communication of material flow, ecological flow and ecological process and realizing the functional integrity of a regional ecological system. In the embodiment of the disclosure, the ecological corridor is identified by adopting a minimum cumulative resistance model and a circuit model, and the minimum cumulative resistance model considers that the ecological process horizontal flow needs to overcome corresponding landscape resistance, so that the ecological resistance reflects the degree of resistance of landscape internal migration. The ecological source area can be equivalent to a circuit node, the non-ecological source area is regarded as resistors with different resistance values, and the ecological corridor is identified based on the current value passing through the non-ecological source area. Therefore, the orderly expansion of urban development is realized, the scientific reservation of urban ecological source areas and corridor spaces can be achieved, and the urban and ecological symbiosis is promoted.

In summary, the ecological corridor construction method based on the multi-ecological source area of the embodiment of the disclosure obtains the normalized vegetation index, the topographic relief, the target water source conservation factor, the carbon fixation factor and the oxygen release factor by obtaining the multiband remote sensing data, the underlying surface data, the precipitation data, the statistical data and the night light index of the target area, performs index calculation based on the multiband remote sensing data, the underlying surface data, the precipitation data, the statistical data and the night light index, performs ecological source area identification based on the normalized vegetation index, the topographic relief, the target water source conservation factor, the carbon fixation factor and the oxygen release factor, obtains the grain supply source area, the water source conservation source area, the habitat maintenance source area, the carbon fixation oxygen release source area and the soil maintenance source area, performs union operation on the grain supply source area, the water source conservation source area, the habitat maintenance source area, the carbon fixation oxygen release source area and the soil maintenance source area to obtain the comprehensive ecological source area distribution area, and acquiring an ecological resistance value surface area layer corresponding to the comprehensive ecological source ground distribution area, calculating a consumption cost value based on the ecological resistance value surface area layer, and acquiring an ecological corridor corresponding to the minimum consumption cost value in the consumption cost values as a target ecological corridor. Therefore, a scientific analysis and demarcation method is provided for urban ecological space reservation and protection, and harmonious development of urban expansion and ecological protection is realized.

Based on the description of the above embodiment, as shown in fig. 2, the method includes data collection (remote sensing data, underlying surface data, precipitation data, statistical data and lighting data), index calculation (normalized vegetation index, terrain factor, land element, land distance factor, runoff factor, fixed carbon, oxygen release factor and erosion factor), source and ground identification (grain supply, water conservation, habitat maintenance, fixed carbon oxygen release and soil conservation), resistive surface construction (land resistance equivalent and grid discretization) and corridor construction according to a minimum resistance model method, so as to obtain a minimum resistance ecological corridor connecting ecological source and ground.

As shown in fig. 3, a data layer (remote sensing data, DEM (Digital Elevation Model) data, land utilization data, precipitation data, soil data, statistical data, and light data) connection element layer (field yield, NDVI, land type element, waviness, precipitation, runoff coefficient, fixed carbon rate, oxygen release rate, water source distance, construction land distance, light index, hidden condition, rainfall erosion force, soil texture, and vegetation cover) connection analysis layer (grain supply, water source conservation, habitat maintenance, fixed carbon oxygen release, and soil maintenance) is included.

Therefore, quantitative indexes such as normalized vegetation indexes and erosion factors are calculated scientifically through collection and analysis of data such as remote sensing, bedding surface and rainfall, and ecological functions such as grain supply, water source conservation, habitat maintenance, carbon and oxygen fixation and release, soil conservation and the like, such as a land block which is located in the top 20% of superior areas and has a continuous whole area larger than 50ha, are selected as ecological source areas. And analyzing the importance of the ecological source areas based on the landscape connectivity importance index so as to identify the protection level of each ecological source area. And extracting the optimal ecological corridor among all ecological source places by utilizing the minimum accumulated resistance model and the circuit model and based on the discretized grid ecological resistance surface, thereby providing effective guarantee for the urban ecological safety pattern.

As a scene example, taking T city of S province as an example, combining a land utilization database of T city, screening out farmland attribute land types including paddy fields, dry lands, irrigated lands and the like, inquiring statistical data to know that the annual average grain yield of T city is 6712.14kg/ha, the average NDVI value of a farmland area is 0.38, and carrying out spatial discretization on the grain yield of the farmland area by utilizing a rasterized NDVI value. And selecting a land with grain yield of the first 20% in the high-yield area and continuous whole area of more than 50ha as a grain supply source. For example, T city includes 11 grain supply source areas, the total area is 976ha, and the largest source area, 257ha, accounts for 26% of the total grain supply source area.

Further, screening different underlying surface types based on a T city land utilization database, utilizing a water balance principle, adopting empirical runoff coefficients and rasterized precipitation data after interpolation, correcting the data by utilizing topographic relief, and calculating the supply of precipitation to underground water as water source conservation quantity; and selecting the land with the whole area larger than 50ha in the first 20 percent of high water source conservation area as a water source conservation source land. For example, the T city comprises 13 water source conservation source areas, the total area is 7763ha, the maximum source area is 4406ha, and the ratio of the source areas to the total water source conservation source area is 57%; the area of the source area is 1580ha, and the proportion is 20%.

Further, the habitat maintains the identification of origin. Comprehensively considering four factors such as the distance from a water source, topographic conditions, the distance from a construction site, the concealment of a ground cover and the like, and carrying out weighted quantitative analysis (the weights are 0.3, 0.15, 0.25 and 0.3 in sequence) to reflect the habitat quality level of a land space; the threshold division of each habitat maintenance factor is scientifically set by combining the activity capability and range of animals, and the threshold and equivalent assignment conditions of each factor are detailed in the following table. The land with the high suitability area of the first 20 percent and the whole area more than 50ha is selected as the source land. For example, T city includes 18 habitat maintenance sources and the total area is 109 km ² Maximum source area 23 km ² The proportion of the area of the total habitat maintenance source is 21 percent; secondly, the source area is 17 km ² The proportion is up to 15%.

Specifically, the habitat maintenance source identification element index division table is shown in table 2:

TABLE 2

And further, identifying a carbon-fixing oxygen-releasing source area, performing discretization treatment on the carbon-fixing oxygen-releasing source area and the carbon-fixing oxygen-releasing source area by utilizing the gridded NDVI based on the average carbon-fixing and oxygen-releasing rates of different surface vegetation types to obtain a carbon-fixing and oxygen-releasing distribution map of T city, then performing spatial summation on the carbon-fixing and oxygen-releasing capacities, and selecting the land blocks with the first 20% of high-yield areas and the whole area larger than 50ha as the carbon-fixing and oxygen-releasing source areas. For example, T city includes 13 carbon-fixing oxygen-releasing source places, and the total area is 102 km ² Maximum source area of about 77 km ² The area proportion of the total carbon-fixing oxygen release source is 75 percent; the second is a source area of about 7 km ² The proportion is 7 percent.

Further, identifying a soil conservation source, comprehensively considering rainfall erosion force, soil texture, topographic relief degree, vegetation coverage factor and the like to perform integrated evaluation based on a soil erosion risk model simplified by a general soil loss equation, and carrying out integrated evaluation on the soil conservation factor B = (R multiplied by S multiplied by D multiplied by C) ^1/4 Wherein R is rainfall erosion force; s is the soil texture; d is the relief degree of the terrain; c is vegetation coverage; the specific classification criteria and score equivalents are detailed in table 1. And finally, selecting the land blocks with the first 20 percent of low-sensitivity areas with better soil retention function and the whole continuous area larger than 50ha as soil retention source lands. For example, T city comprises 10 soil conservation source areas with total area of 57 km ² Maximum source area of about 21 km ² The proportion of the total soil to the area of the source land is 36 percent; the second is a source area of about 16 km ² The ratio is 27%.

Furthermore, the comprehensive ecological source area, namely the comprehensive ecological source area distribution range in T city is obtained through union operation based on five ecological functions of grain supply, water source conservation, habitat maintenance, carbon fixation and oxygen release, soil maintenance and the like, and diversified comprehensive ecological service functions are realized. For example, the area of the T city integrated ecological source is 268.9 km ² The area of the ecological source which can provide 3 ecological functions is 6.4 km ² 2.4% of the total source area, canThe ecological source area capable of providing 2 ecological functions is 74.5 km ² 27.7%, and 188.0 km of ecological source area only capable of providing 1 ecological function ² And accounts for 69.9 percent.

Further, recognizing the importance of the ecological source, specifically, evaluating the influence on the overall landscape connectivity under the condition that any individual unit of the comprehensive ecological source is absent based on the landscape connectivity principle, and quantitatively expressing the influence by using a landscape connectivity importance index so as to scientifically determine the importance of the plaques of each ecological source. For example, based on the natural breakpoint method, the landscape has the importance of more than 20 percent and is used as an extremely important protection area of an ecological source ground with the area of 171.9 km ² The type of source should be strictly protected, and a target distance buffer zone is reasonably constructed; the landscape connectivity importance is 10% -20%, the landscape is used as an important protection area of an ecological source ground, and the area is 52.7 km ² The type of source should enhance protection and orderly ecological utilization; the landscape connectivity importance is less than 10 percent and is used as an ecological source area key restoration area with the area of 44.2 km ² The type of source should enhance the restoration and protection work and improve the ecological environment quality orderly.

Further, constructing a grid resistance surface, specifically, based on animal life habit characteristics and combined with existing research results, assigning resistance equivalent mean values to land types with different underlying surface conditions, and selecting a rasterized NDVI value for the resistance equivalent mean values of ecological spaces (forests, grasses, water, fields, parks, greenbelts and the like) to perform discretization treatment; selecting a rasterized night light index value for the mean value of the resistance equivalent of a non-ecological space (airport land, commercial service facility land, industrial land, roads, home bases, public facility land and the like) to perform spatial discretization; and then, carrying out spatial combination on the discretization grid image layers of different underlying surface land conditions to obtain a global ecological resistance surface grid image of T city. For example, after discretization, the resistance values of different places are spatially distributed, and the regions with large resistance equivalent values are mainly distributed in regions with relatively high urbanization level.

And finally, recognizing the ecological corridor, specifically, recognizing the ecological corridor by adopting a minimum cumulative resistance model and a circuit model, wherein the minimum cumulative resistance model considers that the horizontal flow of the ecological process needs to overcome the corresponding landscape resistance, so that the ecological resistance can reflect the resistance degree of the internal migration of the landscape. The ecological source area can be equivalent to a circuit node, the non-ecological source area is regarded as resistors with different resistance values, and the ecological corridor is identified based on the current value passing through the non-ecological source area. For example, 38 optimal ecological galleries exist among 25 ecological source places in T city, and the total length is 256 km; the longest corridor therein is about 32 km long. For galleries with longer distance (more than 10 km), emphasis needs to be placed on strengthening protection and restrictive development within the range of 50-100 m of target distance so as to ensure smooth and effective transmission of material flow and energy flow among ecological systems; for galleries with short distance (within 3 km), strengthening protection and gradual widening are needed, and effective connection among different sources and places is realized, so that the biological migration cost is reduced, and the range and the quality of animal habitats are improved.

Fig. 4 is a schematic structural diagram of an ecological corridor constructing apparatus based on a multi-ecological source area, which may be implemented by software and/or hardware, and may be generally integrated in an electronic device.

As shown in fig. 4, the apparatus includes:

the data acquisition module 201 is used for acquiring multiband remote sensing data, underlying surface data, precipitation data, statistical data and night light indexes of a target area;

the index calculation module 202 is configured to perform index calculation based on the multiband remote sensing data, the underlying surface data, the precipitation data, the statistical data, and the night light index to obtain a normalized vegetation index, a topographic relief degree, a target water source conservation factor, a carbon sequestration factor, and an oxygen release factor;

a source area identification module 203, configured to perform ecological source area identification based on the normalized vegetation index, the topographic relief degree, the target water source conservation factor, the carbon fixation factor, and the oxygen release factor, so as to obtain a grain supply source area, a water source conservation source area, a habitat maintenance source area, a carbon fixation and oxygen release source area, and a soil conservation source area;

a union calculation module 204, configured to perform union operation on the grain supply source, the water source conservation source, the habitat maintenance source, the carbon-fixing and oxygen-releasing source, and the soil conservation source to obtain a comprehensive ecological source distribution area;

an obtaining module 205, configured to obtain an ecological resistance value surface area map layer corresponding to the comprehensive ecological source area distribution area;

a cost calculation module 206, configured to calculate a cost consumption value based on the ecological resistance value surface area layer;

and the ecological corridor determining module 207 is used for acquiring the ecological corridor corresponding to the minimum cost value in the cost values as a target ecological corridor.

Optionally, the index calculating module 202 is specifically configured to:

Optionally, the source identification module 203 is specifically configured to:

Optionally, the calculation cost module 206 is configured to include:

the obtaining unit is used for obtaining a field distance and a target ecological resistance value based on the ecological resistance value surface area layer;

and the calculating unit is used for calculating the field distance and the target ecological resistance value based on a preset formula to obtain the cost consumption value.

The obtaining unit is specifically configured to:

acquiring a normalized vegetation index, an ecological resistance equivalent average value and a non-ecological space resistance equivalent average value of each ground cover type of ecological space in the ecological resistance value surface area map layer;

The ecological corridor construction device based on the multi-element ecological source area, provided by the embodiment of the disclosure, can execute the ecological corridor construction method based on the multi-element ecological source area, provided by any embodiment of the disclosure, and has corresponding functional modules and beneficial effects of the execution method.

In accordance with one or more embodiments of the present disclosure, there is provided an electronic device including: a processor; a memory for storing the processor-executable instructions; the processor is used for reading the executable instructions from the memory and executing the instructions to realize the ecological corridor construction method based on the multivariate ecological source region provided by the disclosure.

According to one or more embodiments of the present disclosure, there is provided a computer-readable storage medium storing a computer program for executing any one of the multivariate ecological source based ecological corridor construction methods provided by the present disclosure.

It is noted that, in this document, relational terms such as "first" and "second," and the like, may be used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions. Also, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising an … …" does not exclude the presence of other identical elements in a process, method, article, or apparatus that comprises the element.

The foregoing are merely exemplary embodiments of the present disclosure, which enable those skilled in the art to understand or practice the present disclosure. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the generic principles defined herein may be applied to other embodiments without departing from the spirit or scope of the disclosure. Thus, the present disclosure is not intended to be limited to the embodiments shown herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.

Claims

1. An ecological corridor construction method based on a multi-element ecological source area is characterized by comprising the following steps:

2. The ecological corridor construction method based on the multi-ecological source area according to claim 1, wherein index calculation is performed based on the multi-band remote sensing data, the underlying surface data, the precipitation data, the statistical data and the night light index to obtain a normalized vegetation index, a topographic relief degree, a target water source conservation factor, a carbon sequestration factor and an oxygen release factor, and the method comprises the following steps:

3. The ecological corridor construction method based on multi-ecological source land of claim 1, wherein the ecological source land identification based on the normalized vegetation index, the topographic relief degree, the target water source conservation factor, the carbon fixation factor, the oxygen release factor and other related indexes to obtain a grain supply land, a water source conservation land, a habitat maintenance land, a carbon fixation and oxygen release source land and a soil conservation land comprises:

acquiring a water source distance, a land condition and a weighting factor corresponding to the topographic relief degree, calculating based on the water source distance, the land condition, the topographic relief degree and the corresponding weighting factor to obtain habitat maintenance values, sequencing from large to small based on the habitat maintenance values of each grid, and acquiring a grid with a preset third value before sequencing as a habitat maintenance source; wherein the total grid area with the preset third numerical value is larger than a preset third area threshold;

4. The ecological corridor construction method based on multi-element ecological source area according to claim 1, wherein the calculating cost value based on the ecological resistance value surface area map layer comprises:

5. The ecological corridor construction method based on multi-element ecological source area according to claim 4, wherein the obtaining of the target ecological resistance value based on the ecological resistance value surface area map layer comprises:

6. The utility model provides an ecological corridor construction equipment based on many first ecological sources ground which characterized in that includes:

the data acquisition module is used for acquiring multiband remote sensing data, underlying surface data, precipitation data, statistical data and night light index of a target area;

the union calculation module is used for carrying out union operation on the grain supply source land, the water source conservation source land, the habitat maintenance source land, the carbon-fixing oxygen-releasing source land and the soil maintenance source land to obtain a comprehensive ecological source land distribution area;

7. The ecological corridor construction device based on multi-element ecological source area according to claim 6, wherein the index calculation module is specifically configured to:

8. The ecological corridor construction device based on multivariate ecological source areas according to claim 6, wherein the source area identification module is specifically configured to:

sequencing the target water source conservation factors from large to small based on each grid, and acquiring grids with preset second values before sequencing as the water source conservation source; the total area of the grids with the preset second numerical value is larger than a preset second area threshold value;

9. An electronic device, characterized in that the electronic device comprises:

a processor;

a memory for storing the processor-executable instructions;

the processor is used for reading the executable instructions from the memory and executing the instructions to realize the multivariate ecological source based ecological corridor construction method as claimed in any one of claims 1-5.

10. A computer-readable storage medium, characterized in that the storage medium stores a computer program for executing the ecological corridor construction method based on multi-ecological source area according to any one of claims 1-5.