CN116579518A - Method for identifying protection and restoration targeting area of ecological space of water area - Google Patents

Abstract

The embodiment of the disclosure provides a method for identifying a target area for protecting and repairing a water ecological space. The method is applied to the technical field of water area ecological space protection. The method comprises the following steps: acquiring regional source data to be identified; according to preset indexes, identifying a high-value area of the aquatic environment and an ecological fragile area of the water, and screening the areas as ecological source areas; constructing a comprehensive aquatic ecological resistance surface, and calculating water rows Hong Langdao and aquatic life galleries among ecological source lands based on the comprehensive aquatic ecological resistance surface; identifying ecological protection pinch points and ecological restoration barrier points in the water row Hong Langdao and the water life corridor; the high-value area of the comprehensive aquatic environment, the water ecological fragile area, the ecological protection pinch point and the ecological restoration obstacle point are used as the target area for ecological management. The target area identified by the method can provide space guidance for ecological restoration engineering layout, and can realize the overall protection of the water ecological system, comprehensive treatment and restoration effect of water risk system defense.

Description

Method for identifying protection and restoration targeting area of ecological space of water area

Technical Field

The invention relates to the technical field of protection and restoration of ecological space of a water area, in particular to a method for identifying a target area for protection and restoration of ecological space of a water area.

Background

The water ecological space is various ecological spaces which provide places for ecological hydrologic process, maintain the health and stability of the water ecological system and ensure the safety of water resources. Protection and repair of water ecospace is a common challenge faced by all humans.

In the past, the ecological protection and restoration targets of the water body space are single in all countries. Most of their attention is directed to restoration and treatment of river water quality. The river ecological restoration method mainly comprising water quality restoration and treatment is generally adopted. These single-objective methods are generally directed to one or more types of ecological elements, and ignore the relevance and systematic nature of the ecological elements. These single target treatments are therefore difficult to meet for economic development. The method for comprehensively repairing the water body space is relatively lacking in multi-target, multi-element, multi-system and multi-system methods.

Disclosure of Invention

The disclosure provides a method, a device, equipment and a storage medium for identifying a water area ecological space protection restoration target area.

According to a first aspect of the present disclosure, a method for identifying a target area for protection and repair of a water ecological space is provided. The method comprises the following steps: acquiring regional source data to be identified; carrying out ecological space systematic evaluation on the region to be identified according to the region source data to be identified, and identifying a water ecological high-value region and/or a water ecological fragile region in the region to be identified as an ecological source land; constructing a comprehensive aquatic ecological resistance surface according to the regional source data to be identified, and calculating water rows Hong Langdao among ecological sources based on the comprehensive aquatic ecological resistance surface; identifying a regional river basin unit and constructing an aquatic life corridor; identifying ecological protection pinch points and ecological restoration barrier points in the water row Hong Langdao and the water life corridor; and taking the water environment high-value area, the water ecological fragile area and the ecological protection pinch point and the ecological restoration barrier point on the corridor in the source land as the water ecological space protection restoration targeting area of the area to be identified.

In the aspect and any possible implementation manner described above, there is further provided an implementation manner, where obtaining source data of an area to be identified includes: land use present data, digital elevation model, normalized vegetation index data, precipitation data, potential evapotranspiration data, soil texture type data, river water system data, population data, natural protected land distribution data, ecological red line and town development boundary data.

Aspects and any of the possible implementations as described above, further providing an implementation of calculating a water habitat high value area and/or a water ecology vulnerable area in an area to be identified: identifying a hydrologically regulated key zone and a selected aquatic indicative species distribution site; evaluating water conservation and water system connectivity; spatially superposing the hydrologic regulation key region, the water indicative species distribution region, the water source conservation and the water system connectivity high-value region to extract the water environment high-value region; calculating the pollution risk of the river basin and the flood disaster risk area; and extracting the river basin pollution high risk area and the flood disaster high risk area as water ecology fragile areas.

In the aspect and any possible implementation manner as described above, further providing an implementation manner, constructing a comprehensive aquatic ecological resistance surface according to the to-be-identified area source data includes: extracting a grade and an elevation from the digital elevation model; calculating a distance to a railway and a distance to a road according to the land property; the land type in the data land property is already a discretized attribute, so that only continuous values in five data of gradient, elevation, distance to a railway, distance to a road and normalized vegetation index are discretized; respectively endowing the discretized attributes with resistance system coefficients; giving different weights to the six types of data, including the discrete gradient, the discrete elevation, the discrete distance to the railway, the discrete distance to the road, the discrete normalized vegetation index and the land type; converting the six types of discretized data into a grid chart; and carrying out grid calculation of weighted summation on the six types of data according to the resistance system coefficient and the weight to obtain the comprehensive aquatic ecological resistance surface.

In aspects and any one of the possible implementations as described above, there is further provided an implementation, the obtaining a water flood corridor between ecological origins includes: calculating the weighted cost distance of any two water ecological source lands; creating a weighted cost distance surface, and calculating a minimum cost path from the weighted cost distances source by source; the water flood corridor is extracted from the comprehensive aquatic ecological resistance surface based on the circuit model principle by utilizing the minimum cost path and setting the accumulated resistance cut-off distance to determine the position and the shape of the corridor.

In aspects and any one of the possible implementations described above, there is further provided an implementation, a method of calculating an aquatic life corridor between ecological origins, including: acquiring a digital elevation model, and performing filling processing on the digital elevation model; carrying out flow direction analysis on the well-filled digital elevation model, and determining an 8 neighborhood direction flowing out of each pixel in the grid by using a D8 algorithm; calculating the total number of each flow direction grid in an accumulated mode by using a surface runoff model, calculating the accumulated confluence quantity of each water collecting area, and determining a confluence threshold value for generating surface runoffs; all grids with the flow rate larger than the threshold value are potential river networks, and then river, flow direction and water outlet are extracted by using a river network grading tool; vectorizing the first-level grid river network, and then utilizing a watershed tool to jointly determine the space range of a water collecting area, namely a river basin, according to the river, the flow direction and the water outlet; and (3) spatially connecting the regional main river with the river basin, and identifying a river basin unit of the main river as an aquatic life corridor.

In aspects and any one of the possible implementations described above, there is further provided an implementation, a method of identifying ecological protection pinch points in a water row Hong Langdao and a water life corridor, comprising: acquiring a water line Hong Langdao diagram, and calculating an accumulated current value of the landscape surface by iterative operation to identify important ecological elements by utilizing a circuit model based on a set corridor resistance threshold; the high value region of current density is a physiological protection pinch point within a selected threshold range.

In aspects and any one of the possible implementations as set forth above, there is further provided an implementation, a method for identifying an ecological restoration obstacle point, including: covering the area with a grid, and identifying obstacle points on the water row Hong Langdao and the water mission corridor graph by using a moving window searching method within the range of the selected searching radius and the step length; calculating an original value of the minimum cost distance between each search window and two ecological source grounds connected by a water row Hong Langdao and/or a water life corridor, namely, a resistance of the original value is a normal resistance value, and the original value indicates that barrier points exist; calculating a repair improvement coefficient of each search window reaching the minimum cost distance between two ecological source grounds connected by a water row Hong Langdao and/or a water life corridor, namely the minimum cost distance after the resistance becomes 1, and indicating that the barrier point is cleared; and calculating the difference value between the minimum cost distance and the original value when the resistance is 1, and taking the interval with the highest value of the difference value as a quantitative index for identifying the ecological restoration obstacle point.

According to a second aspect of the disclosure, an electronic device includes a memory and a processor, where the memory stores a computer program, and the processor implements a method for identifying a target area for protection and repair of a water ecological space when the processor executes the program.

According to a third aspect of the present disclosure, a computer readable storage medium, on which a computer program is stored, characterized in that the program, when executed by a processor, implements the method according to any one of claims 1-8. The device comprises a memory and a processor, wherein the memory stores instructions which when executed by the processor cause the processor to execute the method for identifying the water ecological space protection repair target area.

The method for accurately identifying the ecological protection restoration target area of the water body space by using the multi-target cooperation as a guide and through an ecological security pattern modeling technology is provided. Thus, the problem that the past single-target treatment and ecological protection restoration targeting region diagnosis method is difficult to meet the requirement of economic development is solved.

The ecological safety pattern-based water ecological space protection restoration target area is identified, space guidance can be provided for ecological restoration engineering layout, and the restoration effects of overall protection of an aquatic ecological system, comprehensive treatment of water resources and defense of a water risk system can be realized. Meanwhile, the land contradiction between regional development and ecological protection can be effectively relieved.

It should be understood that what is described in this summary is not intended to limit the critical or essential features of the embodiments of the disclosure nor to limit the scope of the disclosure. Other features of the present disclosure will become apparent from the following description.

Drawings

The above and other features, advantages and aspects of embodiments of the present disclosure will become more apparent by reference to the following detailed description when taken in conjunction with the accompanying drawings. For a better understanding of the present disclosure, and without limiting the disclosure thereto, the same or similar reference numerals denote the same or similar elements, wherein:

FIG. 1 is a block diagram of an apparatus for identifying a water volume ecospace protection repair target area in accordance with an embodiment of the present disclosure;

FIG. 2 is a flow chart of a method of identifying a water ecospace protection repair target zone in accordance with an embodiment of the present disclosure;

FIG. 3 is a computational flow diagram of an aquatic environment high value zone according to an embodiment of the present disclosure;

FIG. 4 is a computational flow diagram of a water ecology fragile zone according to an embodiment of the disclosure;

FIG. 5 is a computational flow diagram of a water row Hong Langdao and a water life corridor in accordance with an embodiment of the present disclosure;

FIG. 6 is a computational flow diagram of an ecological protection pinch point in accordance with an embodiment of the present disclosure;

FIG. 7 is a computational flow diagram of an ecological restoration obstacle point according to an embodiment of the present disclosure;

fig. 8 illustrates a block diagram of an exemplary electronic device capable of implementing embodiments of the present disclosure.

Detailed Description

For the purposes of making the objects, technical solutions and advantages of the embodiments of the present disclosure more apparent, the technical solutions of the embodiments of the present disclosure will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present disclosure, and it is apparent that the described embodiments are some embodiments of the present disclosure, but not all embodiments. All other embodiments, which can be made by one of ordinary skill in the art based on the embodiments in this disclosure without inventive faculty, are intended to be within the scope of this disclosure.

In addition, the term "and/or" herein is merely an association relationship describing an association object, and means that three relationships may exist, for example, a and/or B may mean: a exists alone, A and B exist together, and B exists alone. In addition, the character "/" herein generally indicates that the front and rear associated objects are an "or" relationship.

The method for accurately identifying the ecological protection restoration target area of the water body space by using the multi-target cooperation as a guide and through an ecological security pattern modeling technology is provided. Thus, the problem that the traditional ecological protection restoration targeting area identification method of a single treatment target is difficult to meet systematic treatment and integral protection is solved.

The present disclosure may be used in the following scenarios: (1) providing spatial guidance for the ecological restoration project layout. (2) The whole protection of the water supply ecological system, the comprehensive treatment of water resources and the defense of the water risk system provide target areas, thereby achieving the purpose of improving the restoration effect. (3) When the regional development planning is formulated, the ecological protection target and the land are identified, so that the land contradiction between the regional development and the ecological protection is relieved.

Fig. 1 is a block diagram of an apparatus for identifying a water ecospace protection repair target area in accordance with an embodiment of the present disclosure.

As shown in fig. 1, the apparatus 100 includes an area source data storage and management apparatus 110; a source data processing and analyzing device 120, the device 120 including an ecological source land generating program 121, a water row Hong Langdao and an aquatic life corridor generating program 122, an ecological safety pattern generating and managing program 123 and an ecological protection pinch point and ecological restoration obstacle point generating program 124; a human-machine interaction interface 130.

The area source data storage and management device 110 includes a spatial database. The specific spatial data includes: land utilization status data, digital Elevation Model (DEM), normalized vegetation index (NDVI) data, precipitation data, potential evapotranspiration data, soil texture type data, river water system data, demographic data, natural protected land distribution data, ecological red line and town development boundary data. The database supports vector data and raster data.

The purpose of the ecological source generating program 121 is to identify a region of high water environmental values and a region of water ecological weakness using a selected variety of spatial data. The combination of the high-value area of the aquatic environment and the fragile area of the aquatic ecology is the water ecology source. The water ecological source is an area with important water ecological value and ecological vulnerability, and comprises important areas with ecological service functions such as water source conservation, water quality purification, biodiversity maintenance and the like, and vulnerable and sensitive areas such as river basin pollution risks, urban waterlogging points and the like.

The purpose of the water row Hong Langdao and water mission corridor generation routine 122 is to identify water rows Hong Langdao and water mission corridor for an area using selected various spatial data. The water flood passage is a passage which can ensure the migration of water birds, purify water quality, retain rain and flood and promote flood regulation; the water life corridor is a river basin unit through which a main river system and a lake wetland pass, is a key channel for guaranteeing migration of aquatic organisms such as fishes, amphibians and the like, and is a core area for maintaining sustainable development of water resources of a river pool and a city.

The purpose of the ecological safety pattern generation and management program 123 is to combine the ecological safety pattern with the water ecological source land and blue net corridor described above. Ecological safety patterns refer to landscape elements that have a critical meaning for maintaining the health and safety of ecological processes. Specifically, the ecological safety pattern is based on 'plaque-gallery-matrix', ecological elements are arranged in series on the whole, and an ecological network is built to improve the boundary integrity and structural connectivity of the ecological system.

The purpose of the ecological protection pinch and ecological restoration obstacle point generation program 124 is to identify ecological pinch and obstacle points for an area using an ecological security pattern. The ecological protection pinch point is a region with higher current density in the corridor, the ecological resistance value of the region is small, the position is critical, the ecological protection pinch point is a 'throat' region affecting the whole landscape connectivity, the smoothness of the occurrence of the ecological process of water is the largest, the ecological protection pinch point is irreplaceable, and the ecological protection pinch point is an important region for protecting in the future. Ecological restoration obstacle points refer to areas where movement of species between habitat patches is blocked, and removal of ecological obstacle points can increase connectivity between water ecological sources and places, and ecological restoration should be performed.

The man-machine interface 130 is responsible for controlling the program process, defining parameters, and outputting the results of the respective generated programs. The results are presented in the form of images. Icons, scales and legends are set for the images through a drawing tool. The mapping tool may extract attributes such as the area of the relevant map element. Likewise, the mapping tool may filter map elements based on attributes.

Fig. 2 is a flow chart of a method of identifying a water ecospace protection repair target zone in accordance with an embodiment of the present disclosure.

With reference to figure 2 of the drawings,

at 210, regional source data is collected;

in some embodiments, the specific data includes land use present data, digital Elevation Model (DEM) data, normalized vegetation index (NDVI) data, precipitation data, potential evapotranspiration data, soil texture type data, river water system data, demographic data, natural protectiveness distribution data, ecological red lines, and town development boundary data. The above data may be a vector diagram or a raster diagram. The resolution of the above data should be adjusted according to the area coverage. For region-level targeted region identification, the resolution of the data is typically 30 meters. If the data resolution is less than 30 meters, resampling the data by a resampling method, such as nearest neighbor method or spatial interpolation method, to unify the data resolution. If the data resolution is less than 30 meters, resampling the data by adopting a mean value aggregation method to unify the data resolution.

At 220, screening the water ecological source land;

in some embodiments, a hydrologically regulated key zone and a selected aquatic indicative species distribution site are identified; calculating water conservation and water system connectivity; carrying out space superposition on a hydrologic regulation key region, a main water indicative species distribution region, a water source conservation region and a water system connectivity high-value region, and extracting the water environment high-value region; calculating the pollution risk of the river basin and the flood disaster risk area; extracting a river basin pollution high risk area and a flood disaster high risk area as water ecology fragile areas; the water habitat high value area and the water ecology fragile area are combined to generate the water ecology source land.

At 230, constructing an aquatic ecological resistance surface;

in some embodiments, multiple factors affecting the flow of the ecological process are selected, the factors are screened from the aspects of natural endowment and ecological disturbance, weights are given, and a resistance surface is comprehensively constructed.

At 240, identifying a water row Hong Langdao and an aquatic life corridor;

in some embodiments, the water row Hong Langdao and the water mission corridor are identified, respectively, using a circuit model and a surface runoff model based on the selected source data. The water row Hong Langdao and the water mission corridor are combined to form a regional blue net corridor.

At 250, identifying ecological protection pinch points and ecological restoration obstacle points;

in some embodiments, the ecological pinch points and the obstacle points are first identified based on the aquatic ecological resistance surface described above, and then using circuit theory.

At 260, identifying a water ecospace protection repair targeting region;

in some embodiments, the aquatic environment high value zone and the aquatic ecology fragile zone are superimposed as a aquatic ecology source land; and overlapping the water flood passage and the water mission passage to be used as a blue net passage together, so as to construct an ecological safety pattern. Based on the construction result of the ecological safety pattern, the high-value area of the water environment is used as an ecological protection source area, the water ecological fragile area is used as an ecological restoration source area, meanwhile, the space layout of elements such as ecological protection pinch points, ecological restoration barrier points and the like in the corridor is identified, and a target area to be protected and restored in the ecological safety pattern of the research area is comprehensively identified.

FIG. 3 is a computational flow diagram of an aquatic environment high value zone according to an embodiment of the present disclosure.

With reference to figure 3 of the drawings,

at 310, data is acquired;

in some embodiments, river water system data and land use status data are obtained.

At 320, identifying a hydrologic regulation key zone;

the hydrologic regulation key area is identified according to the water system type and the land type. The identification objects comprise a first-level river channel, a selected lake, a wetland, a reservoir and a selected water source protection area, wherein the selected lake, the wetland and the reservoir can supply groundwater. And extracting the area where the identification object is located, and outputting a hydrologic regulation key area.

At 330, identifying a selected aquatic specie distribution site;

in some embodiments, natural protectively distributed data and river water system data are obtained. The aquatic specie distribution sites are identified from a natural protected area vector diagram or grid diagram. Identification objects include those providing important habitat natural protection for fish, benthic organisms, amphibians and water birds. And selecting the identification object and establishing a buffer. The aquatic specie is distributed to include an area within 25 meters of the natural protected perimeter buffer.

At 340, water conservation is calculated;

in some embodiments, precipitation data, potential evapotranspiration data, soil texture type data, vegetation availability moisture content, land use status data, and river water system data are obtained. It may be a vector diagram or a grid diagram. The entire area is covered with a grid of selected size. Regional precipitation, potential evapotranspiration, soil depth, vegetation availability moisture content, land availability and basin boundary values within the grid are extracted for each grid. The grid has a plurality of values that can be averaged. The Water conservation evaluation value of each grid was calculated using the Water yield (Water yield) module in the InVEST model (Integrated Valuation of Ecosystem Services and Tradeoffs). Model inputs include precipitation, potential evapotranspiration, soil depth, vegetation availability moisture content, land utilization data, and basin boundary data. And outputting the water conservation evaluation value. The calculated water conservation evaluation values are divided into different levels by using a natural breakpoint method, and a high-value area with high extraction level is used as a water conservation functional source and ground.

At 350, water system connectivity is calculated;

in some embodiments, river water system data is acquired. In the case of a raster pattern, the vector pattern is converted. First, a river and a lake are used as foreground elements, and are assigned to 0, and other types of land are used as background, and are assigned to 1. The boundary width is set to 1. In the foreground element area, obtaining a water area landscape core area by a Morphological Space Pattern Analysis (MSPA) method; treating each core region as a patch; extracting the area of each plaque; meanwhile, calculating path connectivity values among the patches; a landscape potential connectivity index (probability of connectivity, PC) and a connectivity importance value (percentage variation in PC, dPC) for the core region are calculated based on the plaque area and connectivity values. Wherein the formula of the PC is as follows,

wherein: AL is the total area of the regional landscape; a, a _i Is the area of plaque i; a, a _j Area of plaque j; p (P) _*ij Is the maximum value in the connectivity of all paths between patches i and j; n is the total number of ecological plaques. Wherein the formula dPC is given by the formula,

wherein: PCi-remove is a possible connectivity index for the landscape after plaque i removal. A higher dPC value indicates a greater effect that the patch exhibits on landscape communication.

And sorting according to dPC from low to high, dividing the plaques into different levels to obtain a connectivity pattern of a core area, selecting a connectivity high-value threshold of the core area, and screening the plaques with the connectivity of more than 1 square kilometer as a water system connectivity high-value area. Here, the high value threshold may be selected based on the dPC level of plaque. For example, if divided into 5 stages, a 5 th stage section may be selected as the high value threshold.

At 360, calculating a aquatic environment high value zone;

in some embodiments, the above-mentioned critical area for hydrologic regulation, main water indicative species distribution, water conservation and high-value area for water system connectivity above the research area are uniformly converted into a grid chart, spatial superposition is carried out, the four values of each grid are calculated and added, and finally the chart is formed as a recognition result of the high-value area of the water habitat.

FIG. 4 is a computational flow diagram of a water ecology fragile zone according to an embodiment of the disclosure.

With reference to figure 4 of the drawings,

at 410, data is acquired;

in some embodiments, land use presence data and soil texture type data is obtained. In the case of vector diagrams, the vector diagrams are converted into raster diagrams.

At 420, calculating a risk of river basin contamination;

in some embodiments, a biophysical table of the property for selection is obtained from the correlation study. The biophysical table assigns three parameters for each land use category. The first is nutrient load, including nitrogen load (load_n). The nutrient load represents the nutrient input of the entire landscape. The second is the maximum trapping efficiency, which is a unitless floating point value between 0 and 1, denoted nitrogen (eff_n). Nutrient retention capacity for a particular vegetation type is expressed as the proportion of the amount of nutrient that is trapped upstream by the landscape. Thus, natural landscapes such as forests and wetlands typically have a higher interception capacity, while agricultural and urban land interception capacities are lower. The third parameter is the distance that a plaque of a particular locus type is assumed to hold nutrients at its maximum capacity, nitrogen being denoted crit-len n. After acquiring a biophysical table of the selected site type (e.g., table 1), the nitrogen (N) content within the grid cell is calculated using an NDR module in the invent model. And counting the output quantity of nitrogen (N) from the sub-basin scale, and extracting a high-value area as a basin pollution risk area according to the average nitrogen output quantity of the sub-basin.

Nature of the land use	load_n	eff_n	crit_len_n
				Cultivated land	29	0.25	25000
Arbor woodland	1.58	0.48	100000
				…	…	…	…
…	…	…	…
				Unused land	8.56	0.05	5000

Table 1: biophysical representation of geoproperty

At 430, a flood disaster area is calculated;

in some embodiments, perennial precipitation data and DEM are acquired. And filling the depression into the DEM, simulating the flow direction and the flow rate, and finding out the catchment basin area with runoff stagnation. And then the water collection point is determined as a strategic point for controlling the movement of water flow through the calculation of capturing the pouring point (Snap point) and the river connection (Stream link). And finally, according to urban annual flood data and flood control standards included in the precipitation data, respectively using flood flooding ranges of 1 meter, 2 meters and 2.5 meters to establish flood flooding ranges of 3 risk frequencies in 10 years, 20 years and 50 years, wherein the flood flooding ranges correspond to high, medium and low risk areas respectively, and the high risk areas are used as flood disaster high risk areas.

At 440, computing a water ecology vulnerable zone;

in some embodiments, the areas overlapping the river basin pollution high risk area and the flood disaster high risk area form the water ecology fragile area.

Fig. 5 is a computational flow diagram of a water row Hong Langdao and a water life corridor in accordance with an embodiment of the present disclosure.

With reference to figure 5 of the drawings,

at 510, a surface feature map is acquired;

in some embodiments, land utilization status is obtained, a digital elevation model, and NDVI vegetation coverage data. And acquiring a surface feature map, specifically comprising gradient and elevation, through a digital elevation model. Railway and road profiles are extracted from the current state of land use. The distance to the railway is calculated, the distance to the road. The following surface feature map is obtained: grade, elevation, NDVI, distance to rail, and distance to road. The map of the continuous values in the data includes grade, elevation, NDVI, distance to rail, and distance to road. Discretizing. And respectively endowing the discretized attribute and the selected land utilization type with resistance system coefficients according to the existing research. Each feature map is given a different weight. All the surface feature maps are converted into grid maps.

At 520, calculating a comprehensive aquatic ecological resistance surface;

and according to the resistance coefficient and the weight, carrying out weight addition and grid calculation on the gradient, the elevation, the NDVI, the distance from the railway and the distance from the road to obtain the comprehensive aquatic ecological resistance surface.

At 530, calculating a water flood gallery;

in some embodiments, the calculated aquatic ecology source land, and the integrated aquatic ecology resistance surface are obtained. And extracting the water flood passage from the water ecological source ground combined with the comprehensive water ecological resistance surface by using the circuit model. First, the cost weighted distances CWD (cost weighted distance) among the plurality of paired sources are calculated, a cost weighted distance surface is created to discriminate the least cost path LCP (drain-cost paths), then the weighted cost distances and the least cost path are calculated source by source, and finally, the cumulative resistance cutoff distances are set to determine the position and shape of the corridor, thereby identifying the water flood corridor.

At 540, calculating an aquatic life corridor;

in some embodiments, DEM and river water system data is acquired. And (3) filling the depression of the DEM data, carrying out flow direction analysis on the filled DEM data based on a surface runoff model, determining the 8 neighborhood direction flowing out of each pixel in the grid by using a D8 algorithm, accumulating the total number of each flow-direction grid, and calculating the accumulated flow rate of each water collecting area. A confluence threshold is determined at which surface runoffs can be produced, and all grids with confluence volumes greater than the threshold are potential river networks. And (3) dividing the river network into three stages by using a river network grading tool, and extracting the river, the flow direction and the water outlet. Vectorizing the first-level grid river network, and then utilizing a watershed tool to jointly determine the space range of the water collecting area, namely the river basin, according to the river, the flow direction and the water outlet, so as to obtain a water collecting area diagram. Finally, the main river and the river basin in the area are spatially linked, and the river basin unit of the main river is identified and used as the water life corridor.

At 550, calculate a blue net corridor;

in some embodiments, the row Hong Langdao and the water life gallery are acquired, superimposed, and integrated to form the blue net gallery of the research area.

Fig. 6 is a computational flow diagram of an ecological protection pinch point in accordance with an embodiment of the present disclosure.

With reference to figure 6 of the drawings,

at 610, data is acquired;

in some embodiments, a water row Hong Langdao map, an aquatic life corridor map, integrated aquatic life resistance surface and natural protection land distribution data, and land utilization status data are obtained.

At 620, calculating an accumulated current value for the ecological surface;

in some embodiments, based on the water line Hong Langdao map and the comprehensive aquatic ecological resistance surface, in combination with setting a resistance threshold, an electrical circuit model is utilized to calculate the cumulative current value of the water line flood gallery surface by iterative operations to identify important landscape elements. The output is a grid plot with the accumulated current value of the water flood gallery surface.

At 630, ecological protection pinch points are identified;

in some embodiments, regions of grid point current density values above a threshold are identified within a selected threshold range based on a grid map of accumulated current values. These high value areas are ecological protection pinch points.

At 640, supplemental ecological pinch points are identified;

In some embodiments, important habitat patches for natural protection areas, national parks, etc. are extracted from land use present data. And (3) carrying out superposition analysis on the aquatic life corridor and important habitat plaques of natural protection areas, national parks and the like, and extracting the superposition area to serve as a complementary ecological pinch point.

Fig. 7 is a flowchart of the calculation of an ecological restoration obstacle point according to an embodiment of the present disclosure.

With reference to figure 7 of the drawings,

at 710, data is acquired;

in some embodiments, a blue net corridor map and a synthetic aquatic ecological resistance surface are first obtained.

At 720, an original value of the minimum cost distance is calculated;

in some embodiments, a minimum cost distance (Least cost distance, LCD') for each search window to reach between two ecological origins connected by a blue grid corridor is calculated on the blue grid corridor map using a moving window search method within a selected search radius and step size range. The formula is as follows,

LCD′＝CWD1 _MIN +CWD2 _MIN +(L×R′)

where CWD1MIN and CWD2MIN are the minimum cumulative resistance values of the search window to source 1 and source 2, respectively, L is the longest axis length of the search window, and R' is the characteristic resistance value of the surrogate barrier. The resistance is a normal resistance value obtained from the comprehensive aquatic ecological resistance surface when the original value is calculated, and represents that an obstacle point exists.

At 730, a repair improvement coefficient for the minimum cost distance is calculated;

in some embodiments, the minimum cost distance is calculated also using the minimum cost distance formula, assuming the obstacle point is cleared, after the resistance becomes 1. The result of dividing the minimum cost distance by the search diameter is a repair improvement factor.

At 740, calculating a quantitative indicator identifying the obstacle point;

in some embodiments, the minimum cost distance after a resistance value of 1 is calculated as a reduction from the original value, and the difference is used as a quantitative index for identifying the obstacle point.

At 750, identifying an ecological restoration obstacle point;

in some embodiments, the quantitative index is divided into three types by adopting a natural breakpoint method, and the highest value interval is used as an ecological restoration barrier point.

It should be noted that, for simplicity of description, the foregoing method embodiments are all described as a series of acts, but it should be understood by those skilled in the art that the present disclosure is not limited by the order of acts described, as some steps may be performed in other orders or concurrently in accordance with the present disclosure. Further, those skilled in the art will also appreciate that the embodiments described in the specification are all alternative embodiments, and that the acts and modules referred to are not necessarily required by the present disclosure.

The foregoing is a description of embodiments of the method, and the following further describes embodiments of the present disclosure through examples of apparatus.

Fig. 8 shows a schematic block diagram of an electronic device 800 that may be used to implement embodiments of the present disclosure. Electronic devices are intended to represent various forms of digital computers, such as laptops, desktops, workstations, personal digital assistants, servers, blade servers, mainframes, and other appropriate computers. The electronic device may also represent various forms of mobile devices, such as personal digital processing, cellular telephones, smartphones, wearable devices, and other similar computing devices. The components shown herein, their connections and relationships, and their functions, are meant to be exemplary only, and are not meant to limit implementations of the disclosure described and/or claimed herein.

The electronic device 800 includes a computing unit 801 that can perform various appropriate actions and processes according to a computer program stored in the ROM802 or a computer program loaded from a storage unit 808 into the RAM 803. In the RAM803, various programs and data required for the operation of the electronic device 800 can also be stored. The computing unit 801, the ROM802, and the RAM803 are connected to each other by a bus 804. An I/O interface 805 is also connected to the bus 804.

Various components in electronic device 800 are connected to I/O interface 805, including: an input unit 806 such as a keyboard, mouse, etc.; an output unit 807 such as various types of displays, speakers, and the like; a storage unit 808, such as a magnetic disk, optical disk, etc.; and a communication unit 809, such as a network card, modem, wireless communication transceiver, or the like. The communication unit 809 allows the electronic device 800 to exchange information/data with other devices through a computer network such as the internet and/or various telecommunication networks.

The computing unit 801 may be a variety of general and/or special purpose processing components having processing and computing capabilities. Some examples of computing unit 801 include, but are not limited to, a Central Processing Unit (CPU), a Graphics Processing Unit (GPU), various specialized Artificial Intelligence (AI) computing chips, various computing units running machine learning model algorithms, a Digital Signal Processor (DSP), and any suitable processor, controller, microcontroller, etc. The computing unit 801 performs the various methods and processes described above, such as the method of identifying a water ecospace protection repair target area. For example, in some embodiments, the method of identifying a water ecospace protection repair target area may be implemented as a computer software program tangibly embodied on a machine-readable medium, such as storage unit 808. In some embodiments, part or all of the computer program may be loaded and/or installed onto the electronic device 800 via the ROM802 and/or the communication unit 809. When the computer program is loaded into RAM803 and executed by computing unit 801, one or more steps of the method of identifying a water ecospace protection repair target area described above may be performed. Alternatively, in other embodiments, the computing unit 801 may be configured to perform the method of identifying a water ecospace protection repair target area by any other suitable means (e.g., by means of firmware).

Various implementations of the systems and techniques described here above may be implemented in digital electronic circuitry, integrated circuitry, field Programmable Gate Arrays (FPGAs), application Specific Integrated Circuits (ASICs), application Specific Standard Products (ASSPs), systems-on-chips (SOCs), load programmable logic devices (CPLDs), computer hardware, firmware, software, and/or combinations thereof. These various embodiments may include: implemented in one or more computer programs, the one or more computer programs may be executed and/or interpreted on a programmable system including at least one programmable processor, which may be a special purpose or general-purpose programmable processor, that may receive data and instructions from, and transmit data and instructions to, a storage system, at least one input device, and at least one output device.

Program code for carrying out methods of the present disclosure may be written in any combination of one or more programming languages. These program code may be provided to a processor or controller of a general purpose computer, special purpose computer, or other programmable data processing apparatus such that the program code, when executed by the processor or controller, causes the functions/operations specified in the flowchart and/or block diagram to be implemented. The program code may execute entirely on the machine, partly on the machine, as a stand-alone software package, partly on the machine and partly on a remote machine or entirely on the remote machine or server.

In the context of this disclosure, a machine-readable medium may be a tangible medium that can contain, or store a program for use by or in connection with an instruction execution system, apparatus, or device. The machine-readable medium may be a machine-readable signal medium or a machine-readable storage medium. The machine-readable medium may include, but is not limited to, an electronic, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any suitable combination of the foregoing. More specific examples of a machine-readable storage medium would include an electrical connection based on one or more wires, a portable computer diskette, a hard disk, a Random Access Memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disc read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing.

To provide for interaction with a user, the systems and techniques described here can be implemented on a computer having: display means for displaying information to a user; and a keyboard and pointing device (e.g., a mouse or trackball) by which a user can provide input to the computer. Other kinds of devices may also be used to provide for interaction with a user; for example, feedback provided to the user may be any form of sensory feedback (e.g., visual feedback, auditory feedback, or tactile feedback); and input from the user may be received in any form, including acoustic input, speech input, or tactile input.

The systems and techniques described here can be implemented in a computing system that includes a background component (e.g., as a data server), or that includes a middleware component (e.g., an application server), or that includes a front-end component (e.g., a user computer having a graphical user interface or a web browser through which a user can interact with an implementation of the systems and techniques described here), or any combination of such background, middleware, or front-end components. The components of the system can be interconnected by any form or medium of digital data communication (e.g., a communication network). Examples of communication networks include: local Area Networks (LANs), wide Area Networks (WANs), and the internet.

The computer system may include a client and a server. The client and server are typically remote from each other and typically interact through a communication network. The relationship of client and server arises by virtue of computer programs running on the respective computers and having a client-server relationship to each other. The server may be a cloud server, a server of a distributed system, or a server incorporating a blockchain.

It should be appreciated that various forms of the flows shown above may be used to reorder, add, or delete steps. For example, the steps recited in the present disclosure may be performed in parallel, sequentially, or in a different order, provided that the desired results of the disclosed aspects are achieved, and are not limited herein.

The above detailed description should not be taken as limiting the scope of the present disclosure. It will be apparent to those skilled in the art that various modifications, combinations, sub-combinations and alternatives are possible, depending on design requirements and other factors. Any modifications, equivalent substitutions and improvements made within the spirit and principles of the present disclosure are intended to be included within the scope of the present disclosure.

Claims (10)

1. A method for identifying a water ecological space protection repair targeting area comprises the following steps:

acquiring regional source data to be identified;

carrying out ecological space systematic evaluation on the region to be identified according to the region source data to be identified, and identifying a water ecological high-value region and/or a water ecological fragile region in the region to be identified as an ecological source land;

constructing a comprehensive aquatic ecological resistance surface according to the regional source data to be identified, and calculating water rows Hong Langdao among ecological sources based on the comprehensive aquatic ecological resistance surface;

Identifying a regional river basin unit and constructing an aquatic life corridor;

identifying ecological protection pinch points and ecological restoration barrier points in the water row Hong Langdao and the water life corridor;

and taking the water environment high-value area, the water ecological fragile area, the ecological protection pinch point and the ecological restoration barrier point as the water ecological space protection restoration targeting area of the area to be identified.

2. The method of claim 1, wherein the region source data to be identified comprises: land use present data, digital elevation model, normalized vegetation index data, precipitation data, potential evapotranspiration data, soil texture type data, river water system data, population data, natural protected land distribution data, ecological red line and town development boundary data.

3. The method of claim 2, wherein the water-habitat high value area and/or water ecology fragile area in the area to be identified comprises:

identifying a hydrologically regulated key zone and a selected aquatic indicative species distribution site;

evaluating water conservation function and water system connectivity;

spatially superposing the hydrologic regulation key region, the water indicative species distribution region, the water source conservation and the water system connectivity high-value region to extract the water environment high-value region;

Calculating the pollution risk of the river basin and the flood disaster risk area;

and extracting the river basin pollution high risk area and the flood disaster high risk area as water ecology fragile areas.

4. The method of claim 2, wherein constructing a comprehensive aquatic ecological resistance surface from the regional source data to be identified comprises:

extracting a grade and an elevation from the digital elevation model;

calculating a distance to a railway and a distance to a road according to the land property;

the land type in the data land property is already a discretized attribute, so that only continuous values in five data of gradient, elevation, distance to a railway, distance to a road and normalized vegetation index are discretized;

respectively endowing the discretized attributes with resistance system coefficients;

giving different weights to the six types of data, including the discrete gradient, the discrete elevation, the discrete distance to the railway, the discrete distance to the road, the discrete normalized vegetation index and the land type;

converting the six types of discretized data into a grid chart;

and carrying out grid calculation of weighted summation on the six types of data according to the resistance system coefficient and the weight to obtain the comprehensive aquatic ecological resistance surface.

5. The method of claim 1, wherein the water flood corridor between ecological sources comprises:

calculating the weighted cost distance of any two water ecological source lands;

creating a weighted cost distance surface, and calculating a minimum cost path from the weighted cost distances source by source;

and utilizing the minimum cost path, setting a cumulative resistance cut-off distance to determine the position and the shape of the corridor, and extracting the water flood corridor from the comprehensive aquatic ecological resistance surface based on a circuit model principle.

6. The method of claim 1, wherein the aquatic life corridor between ecological sources comprises:

acquiring a digital elevation model, and performing filling processing on the digital elevation model;

carrying out flow direction analysis on the well-filled digital elevation model, and determining an 8 neighborhood direction flowing out of each pixel in the grid by using a D8 algorithm;

calculating the total number of each flow direction grid in an accumulated way by using a surface runoff model, calculating the accumulated confluence quantity of each water collecting area, and determining a confluence threshold value capable of generating surface runoffs;

all grids with the flow rate larger than the threshold value are potential river networks, and then river, flow direction and water outlet are extracted by using a river network grading tool;

Vectorizing the first-level grid river network, and then utilizing a watershed tool to jointly determine the space range of a water collecting area, namely a river basin, according to the river, the flow direction and the water outlet;

and (3) spatially connecting the regional main river with the river basin, and identifying a river basin unit of the main river as an aquatic life corridor.

7. The method of claim 1, wherein the method of identifying ecological protection pinch points in the water row Hong Langdao and water life corridor comprises:

acquiring a water line Hong Langdao diagram, and calculating an accumulated current value of the landscape surface by iterative operation to identify important ecological elements by utilizing a circuit model based on a set corridor resistance threshold;

the high value region of current density is a physiological protection pinch point within a selected threshold range.

8. The method of claim 1, wherein the method of identifying an ecological restoration obstacle point comprises:

covering the area with a grid, and identifying obstacle points on the water row Hong Langdao and the water mission corridor graph by using a moving window searching method within the range of the selected searching radius and the step length;

calculating an original value of the minimum cost distance between each search window and two ecological source grounds connected by a water row Hong Langdao and/or a water life corridor, namely, a resistance of the original value is a normal resistance value, and the original value indicates that barrier points exist;

Calculating a repair improvement coefficient of each search window reaching the minimum cost distance between two ecological source grounds connected by a water row Hong Langdao and/or a water life corridor, namely the minimum cost distance after the resistance becomes 1, and indicating that the barrier point is cleared;

and calculating the difference value between the minimum cost distance and the original value when the resistance value is 1, and taking the interval with the highest value of the difference value as a quantitative index for identifying the ecological restoration obstacle point.

9. An electronic device comprising a memory and a processor, the memory having stored thereon a computer program, characterized in that the processor, when executing the program, implements the method according to any of claims 1-8.

10. A computer readable storage medium, on which a computer program is stored, characterized in that the program, when being executed by a processor, implements the method according to any one of claims 1-8.