Ecological restoration effect evaluation method for land and space comprehensive improvement area

Abstract

The invention discloses an ecological restoration effect evaluation method for a land and space comprehensive improvement area, and relates to the technical field of ecological environment monitoring. A terrain structure data set and a moisture dynamic data set are obtained after data analysis and image feature extraction; the method comprises the following steps: carrying out terrain disturbance evaluation by constructing a terrain disturbance response index Rtz, carrying out hydrological back-up state evaluation by constructing a hydrological back-up trend index Gsh when the terrain has an unsteady disturbance region, carrying out ecological adaptation evaluation by constructing an ecological adaptation comprehensive index SSZ when the hydrological back-up meets the standard, and generating a landform base map of a renovation region by calling DEM (Digital Elevation Model) data. And each evaluation result is subjected to visual partition display through GIS space rendering, and a repair effect report is automatically generated. According to the method, multi-dimensional quantitative evaluation of terrain disturbance, hydrological refilling and ecological adaptability of the renovation area is realized.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of ecological environment monitoring, in particular to a method for evaluating ecological restoration effect in a land space comprehensive regulation region. BACKGROUND

[0002] With the continuous advancement of land space comprehensive regulation work, ecological restoration has become one of the important means to support regional sustainable development. In typical ecological fragile areas such as mountains, hills and water conservation areas, the regulation engineering often involves the disturbance and reconstruction of the topographic structure, such as cutting slope filling, cutting and draining water system, artificial terrace construction and other operations, which further affect the regional natural hydrological process and ecological succession path. Especially after the implementation of the regulation, the regional ecosystem needs to adapt to the coupling relationship between the topographic pattern and the water supply condition, therefore, scientific evaluation of the ecological restoration effect after the regulation needs to consider the double factors of topographic structure change and hydrological dynamic process, and through comprehensive data extraction, model analysis and visual feedback means, a complete ecological adaptability evaluation framework is established.

[0003] The existing ecological restoration effect evaluation methods are mostly focused on single index analysis, such as vegetation coverage, ground humidity or biodiversity index, etc., lacking the analysis of the correlation response between topographic disturbance and hydrological recovery. In the actual regulation process, topographic change often causes changes in surface slope and catchment structure, thereby affecting runoff path and infiltration process, and the existing evaluation model often ignores this dynamic response effect, which is not easy to reflect the adjustment effect of topographic disturbance on hydrological pattern in real time. In addition, the hydrological analysis in the traditional method often relies on manual zoning, lacking adaptive evaluation mechanism based on actual data, resulting in lagging results and limited precision, which is not conducive to the dynamic tracking and regulation of regional ecological restoration. SUMMARY

[0004] In view of the deficiencies of the prior art, the present application provides a method for evaluating the ecological restoration effect in a land space comprehensive regulation region, which solves the problems in the background art.

[0005] To achieve the above purpose, the present application is realized by the following technical scheme: a method for evaluating the ecological restoration effect in a land space comprehensive regulation region, comprising the following steps: S1, collecting the topography before and after the regulation of the regulation region by digital elevation model DEM and remote sensing equipment, obtaining elevation data set and remote sensing image, performing data analysis and image feature extraction, and then performing dimensionless processing to obtain topographic structure data set and water dynamic data set; S2, performing disturbance response analysis on the hydrological and ecological system disturbance effect based on the topographic structure data set, and generating topographic disturbance evaluation according to the analysis result; S3, when the terrain disturbance evaluation is that there is a non-steady state disturbance area in the terrain, fitting the water dynamic data set with the disturbance response analysis result to perform hydrological back analysis on the treatment area, and generating a hydrological back state evaluation according to the analysis result; S4, when the hydrological back state evaluation is that the hydrological back is in line with the standard, fitting the disturbance response analysis result and the hydrological back analysis result, performing ecological adaptability analysis on the terrain disturbance and the hydrological back state of the treatment area, and generating an ecological adaptation evaluation according to the analysis result; S5, generating a topographic base map of the treatment area through DEM data and water flow path analysis, superimposing all analysis results to perform spatial partition visualization, and performing dynamic updating and comparison display.

[0006] Preferably, the S1 includes S11 and S12; S11, collecting the terrain of the treatment area before and after treatment through digital elevation model DEM and remote sensing equipment, obtaining an elevation data set and a remote sensing image, and performing data analysis based on the digital elevation model DEM and the elevation data set to obtain a geographic data set; The elevation data set includes a geographic elevation change ΔH(x,y); According to the geographic information system GIS, the terrain disturbance rate Dt(x,y) is calculated in units of pixels based on the geographic elevation change ΔH(x,y) of the corresponding grid unit before and after the treatment, and specifically , wherein H pre (x,y) is the elevation value before treatment; Based on the flow direction data and the basin accumulation in the digital elevation model DEM, the terrain catchment index Itp is calculated, specifically: Itp=ln(As×tanβ), wherein As represents the upstream catchment area of the grid unit, β represents the slope angle of the grid unit, ln represents the logarithmic function, and tan represents the tangent function; S12, image feature extraction is performed on the remote sensing image to obtain an image data set; The remote sensing image is processed through remote sensing image surface energy parameter inversion technology to obtain surface net radiation Rn, soil heat flux G and sensible heat flux H, and the latent heat flux λE is calculated according to the surface energy balance relationship λE=Rn-G-H, and then the vegetation evapotranspiration rate Ev is calculated based on the relationship between the latent heat flux λE and the vaporization latent heat constant lv, specifically , wherein the vaporization latent heat constant lv represents the water consumption capacity of vegetation, which is obtained by experiment and derivation, and ρw represents the density of water; The near-infrared band NIR and the short-wave infrared band SWIR in the remote sensing image are extracted by the band separation technology of the remote sensing image, and the difference and the sum of the reflectivity of the near-infrared band NIR and the short-wave infrared band SWIR are compared to obtain the normalized difference moisture index NDMI, and then the surface humidity change rate AMh is calculated according to the time sequence of the difference moisture index NDMI, specifically Wherein, NDMI(t) represents the difference moisture index NDMI at time t, NDMI(t+Δt) represents the difference moisture index NDMI at the next sampling time t+Δt, and Δt represents the sampling time interval.

[0007] Preferably, the S1 further comprises S13. S13, the obtained geographical data set and image data set are dimensionless processed to obtain a terrain structure data set and a water dynamic data set, and are written into a spatial database; The dimensionless processing is performed by the Max-Min maximum and minimum method to remove the dimensional influence of the environmental data. The terrain structure data set comprises a surface slope S(x,y), a terrain disturbance rate Dt(x,y) and a terrain catchment index Itp. The water dynamic data set comprises a vegetation evapotranspiration rate Ev(t) and a surface humidity change rate AMh(t).

[0008] Preferably, the S2 comprises S21. S21, fitting the terrain structure data set, performing disturbance response analysis on the hydrological and ecological system disturbance effect caused by the terrain structure change according to the fitting result, and constructing a terrain disturbance response index Rtz to represent the influence intensity of the terrain disturbance on the catchment capacity and the slope structure, and reflect the change degree of the surface catchment path after the terrain morphology is changed, specifically Wherein, Rtz(x,y) represents the terrain disturbance response index at the digital elevation model DEM spatial grid (x,y) coordinate, is the first-order partial derivative of the terrain disturbance rate Dt(x,y) to the surface slope S(x,y), represents the response speed of the slope change to the terrain disturbance rate, and ln represents the logarithmic function.

[0009] Preferably, the S2 further comprises S22. S22, calculating the mean and the standard deviation of the terrain disturbance response index Rtz of the historical reclamation sample area by the statistical method, setting the sum of the mean and the standard deviation as a terrain disturbance critical threshold Ca, and performing terrain disturbance evaluation with the real-time obtained terrain disturbance response index Rtz, and the specific evaluation scheme is as follows: When the terrain disturbance response index Rtz≤the terrain disturbance critical threshold Ca, it indicates that the terrain tends to be stable, and the hydrological catchment structure is not affected by abnormal disturbance, at this time, it is marked as a terrain stable period, and normal monitoring is maintained. When the terrain disturbance response index Rtz is greater than the terrain disturbance critical threshold Ca, it indicates that the terrain has a non-steady-state disturbance area, at which time the hydrological backfilling backstepping analysis is started.

[0010] Preferably, the S3 comprises S31; S31, when the terrain disturbance evaluation is that the terrain has a non-steady-state disturbance area, performing hydrological backfilling backstepping analysis on the regulated area in the non-steady-state disturbance area according to the water dynamic data set, and constructing a hydrological backfilling trend index Gsh in combination with the terrain disturbance response index Rtz, for quantifying the difference relationship between the surface water infiltration and the ecological water consumption of the regulated area, reflecting the adaptability of the ecological water consumption and the groundwater dynamic under the terrain constraint, specifically , wherein ε represents a small constant to prevent the denominator from being zero, and the value is 10 -6 .

[0011] Preferably, the S3 further comprises S32; S32, calculating the average evapotranspiration rate before regulation and the average surface humidity change rate before regulation respectively according to the statistical method, performing difference processing and presetting as a hydrological backfilling critical threshold Tg, and performing hydrological backfilling state evaluation with the real-time obtained hydrological backfilling trend index Gsh, and the specific evaluation scheme is as follows; When the hydrological backfilling trend index Gsh is less than or equal to the hydrological backfilling critical threshold Tg, it indicates that the hydrological backfilling meets the standard and meets the ecological water demand, at which time the ecological adaptability analysis is triggered; When the hydrological backfilling trend index Gsh is greater than the hydrological backfilling critical threshold Tg, it indicates that the ecological hydrological relationship is unbalanced and the hydrological backfilling is insufficient, and there is evapotranspiration, at which time the sampling frequency is increased by 50% and the evaluation period is extended by 50%, and iterative evaluation is performed through S2, and if the ecological hydrological relationship is unbalanced for three times in a row, the ecological restoration effect is unqualified, and the restoration suggestion is regenerated.

[0012] Preferably, the S4 comprises S41; S41, when the hydrological backfilling state evaluation is that the hydrological backfilling meets the standard, fitting the terrain disturbance response index Rtz and the hydrological backfilling trend index Gsh, performing ecological adaptability analysis on the terrain disturbance and the hydrological backfilling state of the regulated area, and constructing an ecological adaptability comprehensive index SSZ, quantifying the collaborative response process of the ecological system to the terrain and hydrology double changes after regulation, and reflecting the self-regulation and self-balancing ability of the ecological system in the time dimension, specifically: , wherein t0 and t1 respectively represent the start and end times of the evaluation period after regulation, and dt represents the time change amount in the integral.

[0013] Preferably, the S4 further comprises S42; S42. Calculate the average value of the comprehensive ecological adaptation index SSZ in the past year when the ecosystem adapts to the current hydrological pattern using statistical methods, and set the average value as the ecological adaptation critical threshold Tc. Then, conduct an ecological adaptation assessment with the real-time acquired comprehensive ecological adaptation index SSZ. The specific assessment scheme is as follows. When the comprehensive ecological adaptation index SSZ is less than the critical threshold Tc for ecological adaptation, it indicates that the ecosystem is not yet adapted to the current hydrological pattern. At this time, the monitoring cycle is extended by 50%, the current remediation area is adjusted to the first assessment priority, and supplementary sampling is initiated to conduct iterative assessment through S2. When the comprehensive ecological adaptation index SSZ is greater than or equal to the ecological adaptation critical threshold Tc, it indicates that the ecosystem has adapted to the current hydrological pattern and entered a stable recovery period. At this time, the current ecological restoration area is marked as ecologically balanced, and an ecological restoration effect report is generated and sent to the land consolidation and management platform.

[0014] Preferably, S5 includes S51 and S52; S51. The DEM data of the remediation area is called to generate a basic geomorphological base map. Based on the water flow path analysis results, the water catchment network structure of the remediation area is depicted. The topographic disturbance response index Rtz, the hydrological replenishment trend index Gsh, and the ecological adaptation comprehensive index SSZ are mapped and superimposed onto the basic geomorphological base map according to spatial pixels. The disturbance intensity range of each pixel is obtained through GIS spatial rendering algorithm, and the results are automatically divided into stable area, unbalanced area, maladaptive area and adapted area. During the rendering process, areas with stable topography are visualized in green, areas with unbalanced ecological and hydrological relationships are visualized in red, areas where the ecosystem is not yet adapted to the current hydrological pattern are visualized in yellow, and areas where the ecosystem is adapted to the current hydrological pattern are visualized in blue. Then, the hierarchical layer and the geomorphological base map are integrated and displayed through the three-dimensional visualization interface. S52. The template-based automatic arrangement technology is used to integrate and format the topographic disturbance response index Rtz, the hydrological replenishment trend index Gsh, and the comprehensive ecological adaptation index SSZ, and store them in the land space management database. Then, the visualization refresh and interface push technology is used to update them in real time, and the restoration status of different management cycles can be dynamically queried and compared horizontally.

[0015] This invention provides a method for evaluating the ecological restoration effect of comprehensive land space management areas. It has the following beneficial effects: (1) This method constructs a topographic structure analysis system with the Digital Elevation Model (DEM) as its core in S1, and combines it with multispectral image data from remote sensing equipment to form a geographic dataset and an image dataset. Through band separation and energy inversion processing of remote sensing images, surface energy elements including net surface radiation Rn, soil heat flux G, and sensible heat flux H are extracted, and latent heat flux λE and vegetation evapotranspiration rate Ev are calculated. Then, normalized differential humidity (NDMI) is used to reflect the surface humidity change rate ΔMh. After Max-Min dimensionless processing, the data are written into a spatial database to form two types of datasets: topographic structure data set and water dynamics data set. This processing method realizes unified dimensional mapping and data standardization of topographic and hydrological features, making surface morphological parameters and ecological water consumption characteristics comparable and continuous, providing quantifiable and traceable basic data support for subsequent disturbance response and ecological adaptation assessment. By integrating multi-source data such as elevation, slope, energy, and humidity, synchronous perception of natural elements in the remediation area from both spatial and temporal dimensions was achieved, providing refined input conditions for dynamic analysis of ecological restoration.

[0016] (2) The method S2 calculates the topographic disturbance response index Rtz based on the topographic structure data set, quantifies the response intensity of the topography to the remediation activities in the form of the partial derivative of the slope change with the disturbance rate, and determines the critical threshold Ca of topographic disturbance through the historical sample area mean and standard deviation, so as to realize the identification of the topographic stable and non-steady-state areas. When there is non-steady-state disturbance in the topography, S3 calls the water dynamic data set to calculate the hydrological replenishment trend index Gsh, which is used to characterize the difference between surface infiltration and ecological water consumption; combined with the statistically obtained critical threshold Tg of hydrological replenishment, it can be determined whether the hydrology of the remediation area has reached the ecological balance state. If the ecological-hydrological relationship is unbalanced, the sampling frequency is automatically adjusted and the evaluation cycle is extended to form a self-correcting cyclic judgment. If the hydrological replenishment meets the standard, S4 constructs the ecological adaptation comprehensive index SSZ by fitting the integral of the topographic disturbance response index Rtz and the hydrological replenishment trend index Gsh, quantifies the comprehensive adaptability of the ecosystem to topographic and hydrological changes, and distinguishes the ecological balance and imbalance states based on the ecological adaptation critical threshold Tc. This hierarchical threshold system forms a progressive feedback loop from three dimensions: "topographic disturbance, hydrological replenishment, and ecological adaptation," ensuring the dynamism and continuity of the restoration effect assessment. It enables real-time diagnosis and quantitative determination of ecological stability under complex topography and multi-temporal conditions.

[0017] (3) In S5, this method generates a basic geomorphological base map of the remediation area by calling DEM data, and draws the catchment network structure based on the water flow path analysis results. Then, the topographic disturbance response index Rtz, the hydrological replenishment trend index Gsh, and the ecological adaptation comprehensive index SSZ are mapped to the geomorphological layer to achieve spatial pixel-level overlay rendering. Using GIS spatial rendering algorithms and three-dimensional visualization interfaces, different ecological states are coded in different zones: stable areas are green, unbalanced areas are red, maladaptive areas are yellow, and adapted areas are blue, so as to realize the intuitive spatial expression of the remediation effect. At the same time, using template-based automatic arrangement and data interface push technology, the calculation results of each indicator are automatically integrated, formatted, and written into the national land space remediation database, supporting dynamic refresh and horizontal comparative analysis within the remediation cycle. This step realizes the visualization and automated archiving of the evaluation results, enabling the remediation management department to grasp the progress, trend, and spatial distribution of ecological restoration in real time on the same platform, realizing the transformation from qualitative observation to quantitative decision-making, and providing traceable decision support for subsequent ecological management, policy evaluation, and optimization of restoration plans. Attached Figure Description

[0018] Figure 1 This is a schematic diagram illustrating the steps of a method for evaluating the ecological restoration effect of a comprehensive land space management area according to the present invention. Figure 2 This is a flowchart illustrating the steps of an evaluation method for the ecological restoration effect of a comprehensive land space management area according to the present invention. Detailed Implementation

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

[0020] Example 1: Please see Figure 1 This invention provides a method for evaluating the ecological restoration effect of comprehensive land space management areas. To achieve the above objectives, this invention is implemented through the following technical solution, including the following steps: S1. The topography of the treatment area before and after treatment is collected using digital elevation model (DEM) and remote sensing equipment to obtain elevation datasets and remote sensing images. After data analysis and image feature extraction, data preprocessing is performed to obtain topographic structure data sets and water dynamic data sets. S2. Based on the topographic structure data set, conduct disturbance response analysis on the hydrological and ecosystem disturbance effects, and generate a topographic disturbance assessment based on the analysis results; S3. When the terrain disturbance assessment indicates that there is an unsteady disturbance area in the terrain, the water dynamic data set is fitted with the disturbance response analysis results to conduct hydrological recharge back-inference analysis of the treatment area, and a hydrological recharge status assessment is generated based on the analysis results. S4. When the hydrological replenishment status assessment indicates that the hydrological replenishment meets the standards, the disturbance response analysis results and the hydrological replenishment back-inference analysis results are fitted together. An ecological adaptability analysis is conducted on the topographic disturbance and hydrological replenishment status of the treatment area, and an ecological adaptation assessment is generated based on the analysis results. S5. Generate a geomorphic base map of the remediation area through DEM data and water flow path analysis, and overlay all analysis results to visualize spatial partitioning, and perform dynamic updates and comparative displays.

[0021] In this embodiment, S1 achieves multi-dimensional perception of the topography and ecological status of the remediation area before and after remediation, and constructs a data acquisition system centered on the Digital Elevation Model (DEM) and remote sensing imagery. Through joint processing of the elevation dataset and remote sensing imagery, topographic structure features and water dynamic features are simultaneously acquired, forming topographic structure data sets and water dynamic data sets, providing a unified analytical basis for subsequent assessment. Unlike existing technologies that rely solely on a single data source for static topographic analysis, this step introduces a multi-source fusion and preprocessing mechanism during the data acquisition stage, achieving simultaneous modeling and structured organization of landform and hydrological status, ensuring data continuity and timeliness. This realizes a shift from "single static topographic analysis" to "dynamic temporal topographic and hydrological coupled analysis," which can more accurately reflect the real changes in the remediation area before and after remediation, providing comprehensive and quantifiable basic data support for the assessment of ecological restoration effects. S2 performs disturbance response analysis based on topographic structure data sets and generates a topographic disturbance assessment based on the analysis results, identifying changes in topographic structure stability. When an unsteady disturbance zone is detected, S3 automatically calls the water dynamic data set to conduct hydrological recharge back-inference analysis, assessing the matching relationship between ecological water consumption and topographic constraints from the perspective of infiltration and evapotranspiration. When hydrological recharge meets the standards, S4 performs a fitting analysis on the disturbance response analysis results and hydrological recharge back-inference analysis results to conduct an ecological adaptability analysis of the topographic disturbance and hydrological recharge status of the remediation area, and generates an ecological adaptability assessment based on the analysis results. Unlike traditional land remediation assessment methods that rely solely on qualitative indicators or single-period observations, this scheme introduces multiple parameters such as disturbance index, hydrological recharge trend, and comprehensive ecological adaptability index, improving the sensitivity and scientific rigor of ecological restoration status identification. It can detect potential ecological imbalance trends in the early stages of the remediation process and provide targeted remediation strategies for management. S5 generates a geomorphic base map and zonal rendering map of the remediation area by overlaying DEM data, water flow path analysis, and various assessment results, realizing spatial visualization and dynamic comparative display of the remediation area. By utilizing GIS rendering and 3D display interfaces, different ecological states are presented using color-coded zones, clearly showing the distribution characteristics of stable, unbalanced, maladaptive, and adapted zones. Combined with automatic data arrangement and periodic refresh mechanisms, the land consolidation platform enables dynamic updates and horizontal comparisons of results across multiple time periods. Compared to existing technologies that rely on manual mapping or static report output, this solution possesses data-driven automatic visualization and continuous display capabilities, enabling real-time monitoring, trend analysis, and report generation of consolidation effects. This forms an integrated and intelligent method for evaluating and displaying the effects of land spatial consolidation, allowing management units to make decisions based on dynamic data and improving the accuracy, timeliness, and operability of consolidation effect assessments.

[0022] Example 2: Please refer to Figure 2 Specifically: S1 includes S11 and S12; S11. Collect the topography of the treatment area before and after treatment using digital elevation model (DEM) and remote sensing equipment to obtain elevation datasets and remote sensing images, and perform data analysis based on the digital elevation model (DEM) and elevation datasets to obtain geographic datasets. The elevation dataset includes geographic elevation changes ΔH(x,y); Based on the geographic information system (GIS), the topographic disturbance rate Dt(x,y) is calculated in pixels from the geographic elevation change ΔH(x,y) of the corresponding raster units before and after the remediation. This represents the amount of landform change caused by the remediation activities. , where H pre (x,y) represents the elevation value before the remediation; The topographic catchment index Itp is calculated based on the flow direction data and watershed accumulation in the digital elevation model (DEM). It describes the ability of a surface unit to collect runoff. Specifically, Itp = ln(As × tanβ), where As represents the upstream catchment area of ​​the grid unit, β represents the slope angle of the grid unit, ln represents the logarithmic function, and tan represents the tangent function. S12. Extract image features from remote sensing images to obtain image datasets; Remote sensing images are processed using surface energy parameter inversion techniques to obtain net surface radiation Rn, soil heat flux G, and sensible heat flux H. Latent heat flux λE is calculated based on the surface energy balance equation λE=Rn-GH. Then, vegetation evapotranspiration rate Ev is calculated based on the relationship between latent heat flux λE and the latent heat of vaporization constant lv, representing the total amount of water evaporation and transpiration in the vegetation system per unit time. Among them, the latent heat of vaporization constant lv represents the water consumption capacity of vegetation, which is obtained through experimental derivation, and ρw represents the density of water; Near-infrared (NIR) and short-wave infrared (SWIR) bands in remote sensing images are extracted using band separation techniques. The difference and sum of reflectance between NIR and SWIR bands are compared to obtain the normalized differential humidity index (NDMI). The surface humidity change rate (ΔMh) is then calculated based on the NDMI time series to reflect the dynamic change in surface water content per unit time. , where NDMI(t) represents the differential humidity NDMI at time t, NDMI(t+Δt) represents the differential humidity NDMI at the next sampling time t+Δt, and Δt represents the sampling time interval.

[0023] S1 also includes S13; S13. Dimensionless processing is performed on the acquired geographic dataset and image dataset to obtain the topographic structure data set and the water dynamics data set, and then written into the spatial database. The dimensionless processing is performed by removing the dimensional influence of environmental data using the Max-Min method. The terrain structure data set includes surface slope S(x,y), terrain disturbance rate Dt(x,y), and terrain water catchment index Itp; The water dynamics data set includes vegetation evapotranspiration rate Ev(t) and surface humidity change rate ΔMh(t).

[0024] In this embodiment, a multi-source data acquisition mechanism for collaborative analysis of digital elevation model (DEM) and remote sensing images was established, enabling simultaneous acquisition and standardized processing of topographic structure and ecological moisture status before and after the remediation of the area. S11 uses the DEM and elevation dataset to extract and calculate geographic elevation change ΔH(x,y), surface slope S(x,y), topographic disturbance rate Dt(x,y), and topographic drainage index Itp, quantifying surface morphological changes and drainage characteristics. S12 utilizes energy inversion technology from remote sensing images to obtain net surface radiation Rn, soil heat flux G, and sensible heat flux H, and calculates latent heat flux λE and vegetation evapotranspiration rate Ev. Combined with the normalized differential humidity (NDMI) time series extracted by band separation, the surface humidity change rate ΔMh is calculated, thus constructing a complete dynamic moisture parameter system. After Max-Min dimensionless processing, all topographic and moisture characteristic data are uniformly written into the spatial database, forming a topographic structure data group and a dynamic moisture data group. This step achieves a unified description of topographic morphology parameters and eco-hydrological parameters in both time and space, solving the technical problems of inconsistent data dimensions and separation of topographic and water information in traditional topographic analysis. It enables the quantitative expression of the coupling relationship between geomorphic disturbance and water change, improving the accuracy, real-time performance, and comparability of subsequent disturbance response and ecological restoration assessments, and providing high-resolution, standardized basic data support for dynamic ecological assessment.

[0025] Example 3: Please refer to Figure 2 Specifically: S2 includes S21; S21. Fit the topographic structure data set, and based on the fitting results, analyze the disturbance response of the hydrological and ecosystem disturbance effects caused by the change in topographic structure. Construct the topographic disturbance response index Rtz to represent the intensity of the impact of topographic disturbance on water catchment capacity and slope structure, reflecting the degree of change in surface water catchment paths after the remediation of the topographic morphology. Where Rtz(x,y) represents the terrain disturbance response index at the (x,y) coordinates of the DEM spatial grid. Let be the first partial derivative of the terrain disturbance rate Dt(x,y) with respect to the surface slope S(x,y), representing the response speed of slope change to the terrain disturbance rate, and ln denotes the logarithmic function.

[0026] S2 further includes S22; S22. Calculate the mean and standard deviation of the terrain disturbance response index Rtz of the historical remediation sample area using statistical methods. Set the sum of the mean and standard deviation as the terrain disturbance critical threshold Ca, and evaluate the terrain disturbance with the real-time acquired terrain disturbance response index Rtz. The specific evaluation scheme is as follows. When the topographic disturbance response index Rtz ≤ the topographic disturbance critical threshold Ca, it indicates that the topography tends to be stable and the hydrological runoff structure is not affected by abnormal disturbances. This period is marked as the topographic stability period and should be monitored normally. When the terrain disturbance response index Rtz > the terrain disturbance critical threshold Ca, it indicates that there is an unsteady disturbance zone in the terrain, and hydrological back-inference analysis is initiated at this time.

[0027] In this embodiment, S21 utilizes topographic structure data sets to analyze the disturbance response of hydrological and ecosystem disturbances caused by topographic structure changes, and constructs a topographic disturbance response index Rtz to quantify the impact of topographic disturbance on water catchment capacity and slope structure. The topographic disturbance response index Rtz reflects the reshaping effect of the landform on the water catchment path after the remediation activities. In partial differential dynamics and geomorphological evolution simulation, the response to disturbance is often represented by the first-order partial derivative of a disturbance field with respect to control variables. This represents the instantaneous response strength of the system's output variable F to the input variable G. In geomorphological analysis, this concept is used to describe the sensitive response of terrain disturbance (output variable terrain disturbance rate Dt(x,y)) to slope change (input variable surface slope S(x,y)), denoted as... In hydrogeomorphology, water catchment capacity is determined by the topographic water catchment index model. Defined as follows, where As is the upstream catchment area and β is the slope angle, used to measure the potential impact of topography on runoff convergence capacity. (Simply) This model only reflects the response of topography to slope but ignores the influence of hydrological processes. To address this, hydrodynamic influences are introduced, and Itp0 from the topographic catchment index model is incorporated into the disturbance response function. The nonlinear relationship between catchment and slope in the post-treatment area is also considered. Since As or tanβ may approach zero in some areas, a logarithmically smoothed form, ln(1+Itp0), is used to prevent logarithmic divergence. The sensitivity to topographic disturbance is then assessed. Multiplying the topographic disturbance response index Rtz by the water flow modulation term ln(1+Itp0) yields the topographic disturbance response index Rtz. In the formula, the surface slope S(x,y) and the topographic disturbance rate Dt(x,y) are dimensionless values, and the topographic water flow index Itp is a logarithmic function result, also dimensionless. Therefore, this formula is a dimensionless operation and has no dimensional influence. Based on the statistical characteristics of historical sample areas, S22 calculates the mean and standard deviation of the topographic disturbance response index Rtz, establishes a critical threshold Ca for topographic disturbance, and uses this as a real-time evaluation standard to achieve automatic identification and classification of the topographic stability state of the remediation area. When the topographic disturbance response index Rtz does not exceed the critical threshold Ca, the regional topographic structure is determined to be in a stable cycle, and the hydrological runoff is not disturbed; when the topographic disturbance response index Rtz exceeds the critical threshold Ca, it is automatically marked as an unsteady disturbance area and triggers hydrological back-calculation analysis. This technical solution introduces partial derivative functions to construct a response model and combines it with statistical threshold methods for dynamic judgment. It not only realizes the quantitative and temporal identification of terrain disturbance intensity, but also improves the accuracy and automation of geomorphological change monitoring. It overcomes the problems of traditional remediation assessment relying on manual interpretation and insufficient spatial resolution, thus enabling earlier detection of abnormal geomorphological change trends and providing a reliable dynamic early warning basis for subsequent hydrological replenishment and ecological restoration.

[0028] Example 4: Please refer to Figure 2 Specifically: S3 includes S31; S31. When the topographic disturbance assessment indicates the existence of an unsteady disturbance zone, a hydrological recharge back-analysis is conducted on the remediation area with the unsteady disturbance zone based on the water dynamics data set. A hydrological recharge trend index Gsh is constructed by combining the topographic disturbance response index Rtz. This index is used to quantify the difference between surface water infiltration and ecological water consumption in the remediation area, reflecting the dynamic adaptability of groundwater under topographic constraints and ecological water consumption. Specifically... Where ε represents a small constant to prevent the denominator from being zero, and its value is 10. -6 .

[0029] S3 further includes S32; S32. Calculate the average evapotranspiration rate and the average surface humidity change rate before treatment according to the statistical method. After the difference is processed, the value is preset as the critical threshold Tg for hydrological replenishment. The hydrological replenishment status is evaluated by comparing it with the real-time hydrological replenishment trend index Gsh. The specific evaluation scheme is as follows. When the hydrological recharge trend index Gsh ≤ the hydrological recharge critical threshold Tg, it indicates that the hydrological recharge meets the standard and satisfies the ecological water demand, and at this time, the ecological adaptability analysis is triggered. When the hydrological replenishment trend index Gsh > the hydrological replenishment critical threshold Tg, it indicates an imbalance in the eco-hydrological relationship, insufficient hydrological replenishment, and evapotranspiration. In this case, the sampling frequency is increased by 50%, the evaluation cycle is extended by 50%, and iterative evaluation is performed through S2. If the eco-hydrological relationship is imbalanced for three consecutive times, the ecological restoration effect is unqualified, and a new restoration recommendation is generated.

[0030] In this embodiment, when the topographic disturbance assessment results indicate that there is an unsteady disturbance zone in the remediation area, S31 conducts a hydrological recharge back-inference analysis on the area based on the water dynamic data group, and calculates the hydrological recharge trend index Gsh by combining the topographic disturbance response index Rtz to quantitatively characterize the difference between surface water infiltration and ecological water consumption, reflecting the dynamic adaptability of groundwater to topographic constraints.

[0031] In land consolidation scenarios, precipitation (P) and runoff (R) are controlled by seasonal external conditions, while ecosystem restoration focuses on the difference between evapotranspiration (Ev) and surface moisture change (ΔMh), i.e., the deviation between ecological water consumption and actual water storage. A hydrological recharge difference function Fw(t) = Ev(t) − ΔMh(t) can be defined to represent the difference between ecological water consumption and surface water content change per unit time. If Fw > 0, it indicates that water consumption exceeds recharge; if Fw < 0, it indicates sufficient surface water recharge. The difference expression is a prototype function of hydrological difference derived from energy balance and water conservation. To eliminate numerical bias caused by geomorphic differences between different topographic units, this invention divides the hydrological difference function Fw(t) by the topographic disturbance response index Rtz(x,y) to construct a topographically normalized hydrological recharge trend index Gsh. In this formula, the vegetation evapotranspiration rate Ev(t) and the surface humidity change rate ΔMh(t) are dimensionless values ​​after dimensionless processing, and the topographic disturbance response index Rtz is a dimensionless value obtained from disturbance response analysis. Therefore, the entire formula is dimensionless, and there is no dimensional influence in the calculation process. S32 calculates the difference between the average evapotranspiration rate and the surface humidity change rate before remediation as the critical threshold Tg for hydrological recharge, and compares it with the real-time hydrological recharge trend index Gsh to determine whether the current eco-hydrological system is in equilibrium. If the hydrological recharge trend index Gsh exceeds the critical threshold Tg, a high-frequency sampling and period extension mechanism is automatically triggered, entering the iterative analysis and remediation suggestion generation stage. Real-time monitoring of the balance between ecological water consumption and hydrological infiltration, and dynamic capture of the hydrological regulation lag and groundwater response characteristics in the remediation area, compared with traditional static monitoring methods, accurately reveals the adaptation status and trend changes of the hydrological system in ecological restoration, improves the accuracy and timeliness of assessment, and provides a refined hydrological control basis for land space remediation.

[0032] Example 5: Please refer to Figure 2 Specifically: S4 includes S41; S41. When the hydrological recharge status assessment indicates that the hydrological recharge meets the standards, the topographic disturbance response index Rtz and the hydrological recharge trend index Gsh are fitted to conduct an ecological adaptability analysis of the topographic disturbance and hydrological recharge status of the remediation area, and an ecological adaptability comprehensive index SSZ is constructed to quantify the coordinated response process of the ecosystem to the dual changes in topography and hydrology after remediation, reflecting the ecosystem's self-regulation and self-balancing capacity in the time dimension, specifically: , where t0 and t1 represent the start and end times of the post-rectification assessment period, respectively, and dt represents the time change in the integral.

[0033] S4 further includes S42; S42. Calculate the average value of the comprehensive ecological adaptation index SSZ in the past year when the ecosystem adapts to the current hydrological pattern using statistical methods, and set the average value as the ecological adaptation critical threshold Tc. Then, conduct an ecological adaptation assessment with the real-time acquired comprehensive ecological adaptation index SSZ. The specific assessment scheme is as follows. When the comprehensive ecological adaptation index SSZ is less than the critical threshold Tc for ecological adaptation, it indicates that the ecosystem is not yet adapted to the current hydrological pattern. At this time, the monitoring cycle is extended by 50%, the current remediation area is adjusted to the first assessment priority, and supplementary sampling is initiated to conduct iterative assessment through S2. When the comprehensive ecological adaptation index SSZ is greater than or equal to the ecological adaptation critical threshold Tc, it indicates that the ecosystem has adapted to the current hydrological pattern and entered a stable recovery period. At this time, the current ecological restoration area is marked as ecologically balanced, and an ecological restoration effect report is generated and sent to the land consolidation and management platform.

[0034] In this embodiment, when the hydrological recharge status assessment results indicate that the hydrological conditions in the remediation area have reached the ecological recharge standard, S41 performs function fitting on the topographic disturbance response index Rtz and the hydrological recharge trend index Gsh to construct the comprehensive ecological adaptation index SSZ. This quantitatively analyzes the synergistic effect of topography and hydrology within the remediation area, reflecting the self-regulation and self-balancing capacity of the ecosystem over time. This formula originates from the time-varying energy conservation and disturbance response integral concepts of ecosystems, representing a combination of the "weighted integral model" in mathematics and the "time cumulative response function" in physics. The original mathematical prototype can be traced back to the time-varying system response formula. This form is used to describe the cumulative effect of the input function f(t) on the state variable g(t) of a system within the time interval t0, t1. In ecophysics and hydrology, it is often used to describe the "temporal coupling relationship between external driving forces and system response". In this scheme, the topographic disturbance response index Rtz(x,y) represents the spatial response intensity of the topographic structure to disturbance, and the hydrological replenishment trend index Gsh(t) represents the dynamic trend of the eco-hydrological replenishment process in the time dimension. The product of the two, Gsh(t)×Rtz(x,y), reflects the "coupling synergistic effect of topography and hydrology". The above integral model does not consider the inhibitory effect of humidity dynamics on the ecosystem response. The absolute value of the surface humidity change rate ΔMh(t) is introduced to reflect the "damping effect of water fluctuations" in the ecosystem, forming the topographic disturbance response index Rtz. In this formula, the topographic disturbance response index Rtz and the hydrological replenishment trend index Gsh are both dimensionless, and the surface humidity change rate ΔMh(t) is a dimensionless value after dimensionless processing. Both the numerator and denominator are dimensionless, and the integral is summed over time, resulting in a dimensionless quantity. S42 calculates the mean of the comprehensive ecological adaptation index SSZ based on historical annual data, sets it as the critical threshold Tc for ecological adaptation, and compares it with the real-time monitored comprehensive ecological adaptation index SSZ. When the comprehensive ecological adaptation index SSZ is lower than the critical threshold Tc, the system automatically extends the monitoring period and adjusts the sampling strategy, performing a second disturbance response analysis by returning to stage S2 to achieve self-feedback correction of the ecological imbalance state. When the comprehensive ecological adaptation index SSZ is greater than or equal to the critical threshold Tc, it is determined that the ecosystem has adapted to the current hydrological pattern, and an ecological restoration effect report is automatically generated and pushed to the land consolidation management platform. This implementation method enables dynamic quantification and thresholding of the adaptation process of the ecosystem in the remediation area, transforming ecological restoration monitoring from static results to continuous feedback assessment, ensuring the consistency and traceability of assessment results in both time and space dimensions, and improving the accuracy, responsiveness, and automation level of ecological restoration assessment.

[0035] Example 6: Please refer to Figure 1 Specifically: S5 includes S51 and S52; S51. The DEM data of the remediation area is called to generate a basic geomorphological base map. Based on the water flow path analysis results, the water catchment network structure of the remediation area is depicted. The topographic disturbance response index Rtz, the hydrological replenishment trend index Gsh, and the ecological adaptation comprehensive index SSZ are mapped and superimposed onto the basic geomorphological base map according to spatial pixels. The disturbance intensity range of each pixel is obtained through GIS spatial rendering algorithm, and the results are automatically divided into stable area, unbalanced area, maladaptive area and adapted area. During the rendering process, areas with stable topography are visualized in green, areas with unbalanced ecological and hydrological relationships are visualized in red, areas where the ecosystem is not yet adapted to the current hydrological pattern are visualized in yellow, and areas where the ecosystem is adapted to the current hydrological pattern are visualized in blue. Then, the hierarchical layer and the geomorphological base map are integrated and displayed through the three-dimensional visualization interface. S52. The template-based automatic arrangement technology is used to integrate and format the topographic disturbance response index Rtz, the hydrological replenishment trend index Gsh, and the comprehensive ecological adaptation index SSZ, and store them in the land space management database. Then, the visualization refresh and interface push technology is used to update them in real time, and the restoration status of different management cycles can be dynamically queried and compared horizontally.

[0036] In this embodiment, S51 generates a basic geomorphic map based on the digital elevation model (DEM) data of the remediation area, and constructs a catchment network structure by combining the results of water flow path analysis. The topographic disturbance response index (Rtz), hydrological recharge trend index (Gsh), and ecological adaptation comprehensive index (SSZ) are then spatially mapped and overlaid at the pixel level. Using a GIS spatial rendering algorithm, different disturbance intensity ranges are automatically classified, dividing the geomorphic area into stable, unbalanced, maladaptive, and adapted zones. Different ecological states are visually displayed using green, red, yellow, and blue visual codes. A 3D visualization interface is then used to integrate and display the layers with the geomorphic map. S52 employs template-based automatic arrangement technology to integrate and format the results of the three types of indices, and writes the results into the national land space remediation database. Dynamic updates and horizontal comparative analysis for different remediation cycles are achieved through visualization refresh and interface push technologies. This implementation method enables visualization and automatic archiving of the entire process results from topographic disturbance and hydrological replenishment to ecological adaptation, allowing the spatial distribution and evolution trend of the remediation effect to be presented intuitively. It not only improves the readability of the remediation status assessment and the real-time nature of decision-making, but also overcomes the limitations of traditional manual mapping and static reports, forming a dynamic and continuously updated visualized assessment system for land space remediation.

[0037] Although embodiments of the invention have been shown and described, it will be understood by those skilled in the art that various changes, modifications, substitutions and alterations can be made to these embodiments without departing from the principles and spirit of the invention, the scope of which is defined by the appended claims and their equivalents.

Claims

1. A method for evaluating the ecological restoration effect of a land space comprehensive improvement area, characterized in that, The method comprises the following steps: S1, collecting the terrain before and after the regulation of the regulation area by a digital elevation model (DEM) and a remote sensing device, obtaining an elevation data set and a remote sensing image, performing data analysis and image feature extraction, and performing dimensionless processing to obtain a terrain structure data set and a water dynamic data set; S2, performing disturbance response analysis on the hydrological and ecological system disturbance effect based on the terrain structure data set, and generating a terrain disturbance evaluation according to the analysis result; S3, when the terrain disturbance evaluation is a non-steady state disturbance area, fitting the water dynamic data set and the disturbance response analysis result to perform hydrological backfilling reverse analysis on the regulation area, and generating a hydrological backfilling state evaluation according to the analysis result; S4, when the hydrological backfilling state evaluation is that the hydrological backfilling meets the standard, fitting the disturbance response analysis result and the hydrological backfilling reverse analysis result, performing ecological adaptability analysis on the terrain disturbance and the hydrological backfilling state of the regulation area, and generating an ecological adaptation evaluation according to the analysis result; S5, generating a regulation area topographic base map by DEM data and water flow path analysis, superimposing all analysis results for spatial partition visualization, and performing dynamic updating and comparison display.

2. The method according to claim 1, wherein, The S1 comprises S11 and S12; S11, collecting the terrain before and after the regulation of the regulation area by a digital elevation model (DEM) and a remote sensing device, obtaining an elevation data set and a remote sensing image, and performing data analysis based on the digital elevation model (DEM) and the elevation data set to obtain a geographic data set; The elevation data set comprises a geographic elevation change ΔH(x, y); According to the geographic information system GIS, the terrain disturbance rate Dt(x, y) is calculated in units of pixels based on the geographic elevation change AH(x, y) of the corresponding grid unit before and after the regulation, and specifically wherein, H pre (x, y) is the elevation value before the regulation; A terrain catchment index Itp is calculated based on the flow direction data in the digital elevation model (DEM) and the flow accumulation, specifically: Itp=ln(As×tanβ), wherein As represents the upstream catchment area of the grid cell, β represents the slope angle of the grid cell, ln represents the logarithmic function, and tan represents the tangent function; S12, image feature extraction is performed on the remote sensing image to obtain an image data set; The remote sensing image is processed by a remote sensing image surface energy parameter inversion technology to obtain a surface net radiation Rn, a soil heat flux G and a sensible heat flux H, and a latent heat flux λE is calculated according to a surface energy balance relationship formula λE = Rn - G - H, and then a vegetation evapotranspiration rate Ev is calculated based on a relationship between the latent heat flux λE and a vaporization latent heat constant lv, that is wherein the vaporization latent heat constant lv represents a vegetation water consumption capacity, which is obtained by experiment deduction, and ρw represents a water density. The near-infrared band NIR and the short-wave infrared band SWIR in the remote sensing image are extracted by a band separation technology of the remote sensing image, and the difference between the reflectivity of the near-infrared band NIR and the reflectivity of the short-wave infrared band SWIR is compared with the sum of the reflectivity, to obtain a normalized difference moisture index NDMI, and then the surface moisture change rate AMh is calculated according to a time sequence of the difference moisture index NDMI, specifically wherein NDMI(t) represents the difference moisture index NDMI at time t, NDMI(t+Δt) represents the difference moisture index NDMI at next sampling time t+Δt, and Δt represents a sampling time interval.

3. The method according to claim 2, wherein, The S1 further comprises S13; S13, dimensionless processing is performed on the obtained geographic data set and image data set to obtain a terrain structure data set and a water dynamic data set, and the data sets are written into a spatial database; The dimensionless processing is performed by a Max-Min maximum and minimum method to remove the dimensional influence of environmental data; The terrain structure data set comprises a surface slope S(x, y), a terrain disturbance rate Dt(x, y), and a terrain catchment index Itp; The water dynamic data set comprises a vegetation evapotranspiration rate Ev(t) and a surface humidity change rate ΔMh(t).

4. The method according to claim 3, wherein, The S2 comprises S21; S21, fitting the terrain structure data set, and according to the fitting result, analyzing the disturbance response of the hydrological and ecological system disturbance effect caused by the terrain structure change, and constructing a terrain disturbance response index Rtz, indicating the influence intensity of the terrain disturbance on the catchment capacity and slope structure, reflecting the change degree of the terrain morphology after the regulation on the surface catchment path, specifically wherein Rtz(x, y) represents the terrain disturbance response index at the spatial grid (x, y) coordinates of the digital elevation model DEM, is the first-order partial derivative of the terrain disturbance rate Dt(x, y) to the ground slope S(x, y), indicating the response speed of the slope change to the terrain disturbance rate, and ln represents the logarithmic function.

5. The method according to claim 4, wherein, The S2 further comprises S22; S22, the mean and standard deviation of the terrain disturbance response index Rtz of the historical regulation sample area are calculated by a statistical method, the sum of the mean and standard deviation is set as a terrain disturbance critical threshold Ca, and a terrain disturbance evaluation is performed on the real-time obtained terrain disturbance response index Rtz, and the specific evaluation scheme is as follows; When the terrain disturbance response index Rtz is less than or equal to the terrain disturbance critical threshold Ca, it indicates that the terrain tends to be stable, and the hydrological catchment structure is not affected by abnormal disturbance, at this time, it is marked as a terrain stable period, and normal monitoring is maintained; When the terrain disturbance response index Rtz is greater than the terrain disturbance critical threshold Ca, it indicates that the terrain has a non-steady-state disturbance area, at which time the hydrological back analysis is started.

6. The method according to claim 5, wherein, The S3 includes S31; S31, when the topographic disturbance evaluation is that there is a non-steady disturbance area in the topography, performing hydrological backfilling and backstepping analysis on the treatment area with the non-steady disturbance area according to the water dynamic data set, and combining the topographic disturbance response index Rtz to construct a hydrological backfilling trend index Gsh, which is used to quantify the difference between the surface water infiltration and the ecological water consumption in the treatment area, and reflects the adaptability of the ecological water consumption and the groundwater dynamic under the constraint of the topography, specifically wherein ε represents a small constant to prevent the denominator from being zero, and the value is 10 -6 .

7. The method according to claim 6, wherein the method is characterized by, The S3 further includes S32; S32, the average evapotranspiration rate before the regulation and the average surface humidity change rate before the regulation are calculated according to the statistical method, difference processing is performed, and the hydrological back-up critical threshold Tg is preset, and the hydrological back-up trend index Gsh obtained in real time is evaluated, and the specific evaluation scheme is as follows: When the hydrological back-up trend index Gsh is less than or equal to the hydrological back-up critical threshold Tg, it indicates that the hydrological back-up meets the standard and meets the ecological water demand, at which time the ecological adaptability analysis is triggered; When the hydrological back-up trend index Gsh is greater than the hydrological back-up critical threshold Tg, it indicates that the ecological hydrological relationship is unbalanced and the evapotranspiration exists due to insufficient hydrological back-up, at which time the sampling frequency is increased by 50% and the evaluation period is extended by 50%, and the iterative evaluation is performed through S2, if the ecological hydrological relationship is unbalanced for three times in succession, the ecological restoration effect is unqualified, and the restoration suggestion is regenerated.

8. The method according to claim 7, wherein, The S4 includes S41; S41、in the hydrological backfill state evaluation is hydrological backfill accord with standard, will topographic disturbance response index Rtz and hydrological backfill trend index Gsh carry out fitting, the topographic disturbance and hydrological backfill state of the harness area is carried out ecological adaptability analysis, and constructs ecological adaptation comprehensive index SSZ, quantifies the collaborative response process of the ecosystem to the double changes of terrain and hydrology after the harness, reflects the self-regulation and self-balancing ability of the ecosystem in time dimension, specifically: Wherein, t0 and t1 respectively represent the start and end time of the evaluation period after the harness, and dt represents the time change in the integral.

9. The method according to claim 8, wherein, The S4 further includes S42; S42, the ecological adaptation comprehensive index SSZ mean value when the ecological system adapts to the current hydrological pattern in the past year is calculated by the statistical method, and the mean value is preset as the ecological adaptation critical threshold Tc, and then the ecological adaptation comprehensive index SSZ is evaluated, and the specific evaluation scheme is as follows: When the ecological adaptation comprehensive index SSZ is less than the ecological adaptation critical threshold Tc, it indicates that the ecological system is not adapted to the current hydrological pattern, at which time the monitoring period is extended by 50%, the current regulation area is adjusted to the first evaluation priority, and the supplementary sampling is started to perform the iterative evaluation through S2; When the ecological adaptation comprehensive index SSZ is greater than or equal to the ecological adaptation critical threshold Tc, it indicates that the ecological system is adapted to the current hydrological pattern, and the ecological system enters the stable recovery period, at which time the state of the current ecological regulation area is marked as ecological balance, and the ecological restoration effect report is generated and sent to the land consolidation management platform.

10. The method according to claim 1, wherein the method is characterized by, The S5 includes S51 and S52; S51, the DEM data of the regulation area is called to generate a geomorphological basic map, the water catchment network structure of the regulation area is depicted according to the water flow path analysis result, the terrain disturbance response index Rtz, the hydrological back-up trend index Gsh and the ecological adaptation comprehensive index SSZ are mapped and superimposed on the geomorphological basic map according to the space pixel, the disturbance intensity interval of each pixel is obtained through the GIS spatial rendering algorithm, and the result is automatically divided into stable area, unbalanced area, non-adaptation area and adaptation area, in the rendering process, the terrain tends to be stable area is visualized by green coding, the ecological hydrological relationship imbalance area is visualized by red coding, the ecological system is not adapted to the current hydrological pattern area is visualized by yellow coding, and the ecological system is adapted to the current hydrological pattern area is visualized by blue coding, and then the hierarchical layers and the geomorphological base map are fused and displayed through the three-dimensional visualization interface; S52, using templated automatic arrangement technology to integrate and format the terrain disturbance response index Rtz, hydrological backfill trend index Gsh and ecological adaptation comprehensive index SSZ, and store them in the land space regulation database. Then, using visual refresh and interface pushing technology for real-time updating, the dynamic query and horizontal comparison of the repair status of different regulation periods are carried out.

Images