Land space planning data analysis system for ecological restoration

Abstract

The invention relates to the technical field of data analysis, in particular to a territorial space planning data analysis system for ecological restoration, which comprises an interference parameter acquisition module, an ecological disturbance calculation module, a restoration area identification module, a priority level judgment module and a site selection stability screening module. According to the method, the accuracy and efficiency of ecological restoration are improved by capturing various disturbance factors of an ecological sensitive area and dynamically evaluating ecological degradation, and the signal of ecological degradation can be recognized earlier by comprehensively monitoring the change trend of earth surface coverage and the dynamic state of a water body and quickly responding to the subtle change of an ecological environment; ecological service function indexes are adopted for inter-year comparison, a long-term view angle about the health state of an ecological system is provided, prediction accuracy is improved, ecological restoration work can be implemented in a more targeted mode, and the lasting effect of the ecological restoration work is ensured through continuous monitoring of land utilization changes.

Description

Technical Field

[0001] The present invention relates to the field of data analysis technology, and in particular to a land space planning data analysis system for ecological restoration. Background Art

[0002] The field of data analysis technology involves the collection, processing, analysis, and mining of large amounts of structured or unstructured data through mathematical modeling, statistical methods, computer algorithms, and visualization techniques to reveal regularities, trends, and correlations. This technical field includes, but is not limited to, data collection and access, data warehousing and management, data modeling, statistical analysis, machine learning modeling, visualization, and intelligent decision support. The goal of data analysis is to quantitatively characterize the inherent structure of data, conduct associative reasoning, and predict trends through algorithms and models, thereby providing a verifiable quantitative basis for system operation optimization, strategy formulation, and resource allocation. In scenarios with dense multi-source spatial data, such as geographic space, land resources, and urban management, data analysis technology often integrates geographic information systems, remote sensing analysis, multidimensional databases, and spatial modeling tools to achieve unified processing and intelligent evaluation of multi-source heterogeneous spatial data.

[0003] The National Land Space Planning Data Analysis System for Ecological Restoration is an intelligent data processing and decision-making support system for ecological governance. It analyzes and processes planning data related to ecological and environmental elements within national land space. The system integrates and processes multiple data resources, including geographic information system data, current land use data, and data on the distribution of ecologically sensitive areas. Through spatial analysis and planning suitability assessments, it supports the site selection, prioritization, and configuration of restoration plans for ecological restoration areas, thereby enabling precise intervention in areas of ecological degradation and scientific optimization of restoration pathways.

[0004] Traditional analysis systems have obvious deficiencies in the dynamic monitoring of ecologically sensitive areas and the capture of long-term ecological change trends, especially in the early identification and continuous monitoring of ecologically degraded areas. Traditional systems lack sufficient sensitivity and predictive capabilities. For example, they do not have the ability to track and analyze surface cover change trends and water body dynamics in real time, making it impossible to identify and assess ecological interference factors in the early stages, delaying the optimal time for intervention. Traditional systems have failed to make full use of years of data to conduct trend analysis on ecological service functions, resulting in static assessments of ecosystem health status that lack the necessary predictability and operability, affecting the effective implementation of ecological restoration strategies and causing unstable or unsustainable results in ecological restoration work. Summary of the Invention

[0005] The purpose of the present invention is to solve the shortcomings of the existing technology and propose a land space planning data analysis system for ecological restoration.

[0006] To achieve the above objectives, the present invention adopts the following technical solution: a national land space planning data analysis system for ecological restoration, the system comprising:

[0007] The interference parameter acquisition module obtains remote sensing sequences of land cover in ecologically sensitive areas, extracts the trend of normalized vegetation index changes, calculates the frequency of land use changes and the land use type diversity index, and generates an information set of ecological disturbance sources;

[0008] The ecological disturbance calculation module sets the normalized vegetation index fluctuation threshold, the runoff variation warning standard, and the upper limit of the land use change frequency based on the ecological disturbance source information set, and determines whether the three disturbance factors in each spatial unit simultaneously meet the corresponding threshold standards. If so, it is marked as an ecological disturbance superposition unit and an ecological degradation disturbance judgment mark is generated;

[0009] The restoration area identification module extracts the number of buildings per unit area and the building function classification level within the spatial unit based on the ecological degradation interference judgment mark, determines whether there is a high-level building density area that overlaps with the ecological interference superposition unit, and generates an ecological restoration space screening result;

[0010] The priority determination module calls the ecological restoration space screening results, extracts the inter-annual series of the ecological service function index in the target unit, calculates the ratio compared to the regional average function degradation cycle as the urgency weight, and generates restoration priority classification information.

[0011] As a further solution of the present invention, the ecological disturbance source information set includes the normalized vegetation index change frequency index, the runoff change intensity level, the land use diversity ratio and the functional change behavior density; the ecological degradation interference judgment identifier is specifically a spatial disturbance level label, a judgment rule set and an interference factor cross-record table; the ecological restoration space screening result includes a restoration feasibility space mask, a development interference exclusion layer and a restoration adaptation unit number list; the restoration priority classification information specifically refers to the urgency level code, the ecological compression cycle ratio and the number of recovery buffer zone intervals.

[0012] As a further embodiment of the present invention, the system includes:

[0013] The interference parameter acquisition module obtains remote sensing sequences of land cover in ecologically sensitive areas, extracts the trend of normalized vegetation index changes, calculates the frequency of land use changes and the land use type diversity index, and generates an information set of ecological disturbance sources;

[0014] The ecological disturbance calculation module sets the normalized vegetation index fluctuation threshold, the runoff variation warning standard, and the upper limit of the land use change frequency based on the ecological disturbance source information set, and determines whether the three disturbance factors in each spatial unit simultaneously meet the corresponding threshold standards. If so, it is marked as an ecological disturbance superposition unit and an ecological degradation disturbance judgment mark is generated;

[0015] The restoration area identification module extracts the number of buildings per unit area and the building function classification level within the spatial unit based on the ecological degradation interference judgment mark, determines whether there is a high-level building density area that overlaps with the ecological interference superposition unit, and generates an ecological restoration space screening result;

[0016] The priority determination module calls the ecological restoration space screening results, extracts the inter-annual series of the ecological service function index in the target unit, calculates the ratio compared to the regional average function degradation cycle as the urgency weight, and generates restoration priority classification information.

[0017] As a further solution of the present invention, the ecological disturbance source information set includes the normalized vegetation index change frequency index, the runoff change intensity level, the land use diversity ratio and the functional change behavior density; the ecological degradation interference judgment identifier is specifically a spatial disturbance level label, a judgment rule set and an interference factor cross-record table; the ecological restoration space screening result includes a restoration feasibility space mask, a development interference exclusion layer and a restoration adaptation unit number list; the restoration priority classification information specifically refers to the urgency level code, the ecological compression cycle ratio and the number of recovery buffer zone intervals.

[0018] As a further solution of the present invention, the interference parameter acquisition module includes:

[0019] The vegetation index extraction submodule obtains remote sensing sequences of land cover in ecologically sensitive areas, extracts time series data of the normalized vegetation index, calculates the monthly normalized vegetation index change trend in the same area, summarizes the frequency of value changes between months within the year, establishes a spatial distribution layer based on the standard deviation of the annual normalized vegetation index change value, and obtains the fluctuation coefficient of vegetation change trend;

[0020] The hydrological fluctuation identification submodule uses the surface water dynamic simulation layer based on the vegetation change trend fluctuation coefficient to extract the seasonal unit time runoff change sequence of the main water systems in the basin, calculates the standard deviation of the four-season runoff values in each spatial unit and integrates them into a continuous change sequence table to obtain the seasonal amplitude coefficient of the hydrological change;

[0021] The land use statistics submodule collects annual remote sensing interpretation map data of land use types based on the seasonal amplitude coefficient of hydrological changes, constructs a land use change sequence of consecutive years in each spatial unit, counts the annual conversion times of functional types and the number of type changes, organizes them into a two-dimensional array to record the annual change trend indicators of each unit, and obtains the information set of ecological disturbance sources.

[0022] As a further solution of the present invention, the ecological disturbance calculation module includes:

[0023] The disturbance factor threshold setting module sets the normalized vegetation index fluctuation threshold, the runoff change ratio warning benchmark value, and the land use change frequency limit based on the ecological disturbance source information set, according to the vegetation change trend fluctuation coefficient, the hydrological change seasonal amplitude coefficient, and the land use change trend index, respectively, to establish the judgment baseline for the three types of disturbance factors and obtain the disturbance threshold parameter set;

[0024] The interference unit screening submodule calls the disturbance threshold parameter set to judge the data of each spatial unit in the ecological disturbance source information set, identifying whether the annual standard deviation of the normalized vegetation index is greater than the set fluctuation threshold, whether the runoff variation coefficient is greater than the set benchmark, and whether the frequency of land use change exceeds the limited value, screening the spatial units that meet all three conditions, and generating a list of composite disturbance spatial units;

[0025] The interference level generation submodule assigns an interference identification number to each unit according to the composite disturbance spatial unit list and marks it as an ecological disturbance superposition unit, records the corresponding three types of disturbance factor structures and maps them to the spatial data grid layer, and establishes an ecological degradation interference judgment mark.

[0026] As a further solution of the present invention, the repair area identification module includes:

[0027] The building density extraction submodule calls the regional building density layer based on the ecological degradation interference judgment mark, extracts the number of buildings per unit area in each spatial unit, calculates the ratio of the total number of buildings to the area in each unit, and summarizes the spatial building density information;

[0028] The function level comparison submodule matches the building function classification level layer according to the spatial building density information, extracts the corresponding building use level code in each type of spatial unit, determines whether the number of buildings per unit area is higher than the building density benchmark value, and whether its building use level is within the upper limit of the development level classification, and selects the corresponding units that meet the conditions to obtain the building function interference superposition unit set;

[0029] The suitable zoning submodule calls the building function interference superposition unit set and performs spatial superposition operation with the ecological degradation interference judgment mark to determine whether the two overlap in space, excludes the spatial units with overlapping relationships, assigns the remaining spatial units with ecological restoration suitability marks, and establishes the ecological restoration space screening results.

[0030] As a further solution of the present invention, the priority determination module includes:

[0031] The function sequence extraction submodule calls the ecological restoration space screening results, extracts the ecological service function index time series corresponding to the spatial unit, arranges the inter-annual value changes, establishes a function sequence distribution table under each unit, and obtains the inter-annual change trend set of ecological functions;

[0032] The decay cycle calculation submodule identifies the time periods of consecutive annual functional value decline based on the inter-annual variation trend set of ecological functions, calculates the functional succession urgency value of the spatial unit according to the duration of decline of each spatial unit and the ratio of decline in each period;

[0033] The urgency level assignment submodule sets the urgency baseline standard value according to the functional succession urgency value of the spatial unit, compares the urgency value of each spatial unit, identifies the spatial units that exceed the standard and marks them as high-level intervention objects, and obtains the restoration priority classification information.

[0034] As a further solution of the present invention, the formula for calculating the functional succession urgency value of the spatial unit is:

[0035]

[0036] Among them, U i represents the functional succession urgency value of the i-th spatial unit, T i represents the normalized value of the duration of the functional decline of the i-th spatial unit, represents the normalized value of the regional average functional degradation period, S i V represents the normalized value of the interannual variation slope of the function of the i-th spatial unit. i Represents the fluctuation variance of the ecological function value of the i-th spatial unit.

[0037] As a further embodiment of the present invention, the system further comprises:

[0038] The site stability screening module extracts the land use time series matrix of the corresponding spatial unit based on the restoration priority classification information, counts the annual number of land type transfers and the frequency of reverse changes in land function, and determines whether there are frequent direction reversals or the transfer frequency exceeds the set threshold. If so, it is determined that the long-term stability of the unit restoration is insufficient and a list of available spatial planning sites is generated;

[0039] The available list of spatial planning site selection includes a stability screening mask, a valid site selection area number, and an index table of corresponding indicators for eliminated units.

[0040] As a further solution of the present invention, the site stability screening module includes:

[0041] Using the trajectory extraction submodule to extract the land use time series data of the corresponding spatial unit according to the restoration priority classification information, construct the spatial land use type distribution information by year, obtain the land type change records of consecutive years in the unit, and establish the land use annual transfer sequence set;

[0042] The change frequency calculation submodule counts the number of land use type changes and the cumulative frequency of reverse land function changes in each spatial unit in adjacent years based on the land use annual transfer sequence set, and calculates the land function change intensity index of each unit;

[0043] The stability determination and screening submodule sets the land use stability benchmark standard according to the land function change intensity index of each unit, compares the intensity value of each unit, screens the spatial locations with frequent continuous fluctuations or concentrated reverse changes, eliminates the corresponding spatial units, and establishes an available list of spatial planning sites.

[0044] As a further solution of the present invention, the formula for calculating and obtaining the land function change intensity index of each unit is:

[0045]

[0046] Among them, X j represents the land function change intensity index of the jth spatial unit, n j represents the annual length included in the jth spatial unit, R j,t F represents the number of reverse changes that occurred in the j-th spatial unit in the t-th year, j,t represents the land function code of the jth spatial unit in year t, Z j,t It represents the total number of type changes that occurred in the j-th spatial unit in the t-th year.

[0047] Compared with the prior art, the advantages and positive effects of the present invention are:

[0048] In the present invention, by capturing various disturbance factors in ecologically sensitive areas and dynamically evaluating ecological degradation, the accuracy and efficiency of ecological restoration are improved. By comprehensively monitoring the changing trends of surface cover and water body dynamics, a rapid response is made to subtle changes in the ecological environment, and signals of ecological degradation can be identified earlier. Inter-annual comparison using the ecological service function index provides a long-term perspective on the health status of the ecosystem, increases the accuracy of prediction, and enables ecological restoration work to be implemented in a more targeted manner. Through continuous monitoring of land use changes, the lasting effect of ecological restoration work is ensured, the recurrence or regression that may occur in the ecological restoration process is reduced, the timeliness and accuracy of ecological restoration are improved, and the sustainability and stability of ecological restoration strategies are strengthened. BRIEF DESCRIPTION OF THE DRAWINGS

[0049] In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the following briefly introduces the drawings required for use in the description of the embodiments. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work.

[0050] Figure 1 is a system flow chart of the present invention;

[0051] Figure 2 Schematic diagram of the system framework of the present invention;

[0052] Figure 3 This is a flow chart of the interference parameter acquisition module of the present invention;

[0053] Figure 4 This is a flow chart of the ecological disturbance calculation module of the present invention;

[0054] Figure 5 This is a flowchart of the repair area identification module of the present invention;

[0055] Figure 6 is a flow chart of the priority determination module of the present invention;

[0056] Figure 7 This is a flow chart of the site stability screening module of the present invention. DETAILED DESCRIPTION

[0057] The technical solution of the present invention is described below in conjunction with the accompanying drawings.

[0058] In the embodiments of the present invention, words such as "exemplarily" and "for example" are used to indicate examples, illustrations, or explanations. Any embodiment or design described as an "exemplary" in the present invention should not be interpreted as being preferred or advantageous over other embodiments or designs. Rather, the use of the word "exemplary" is intended to present concepts in a concrete manner. Furthermore, in the embodiments of the present invention, "and / or" can mean both or either of the two.

[0059] In the embodiments of the present invention, the terms "image" and "picture" may be used interchangeably. It should be noted that, when the distinction between them is not emphasized, their intended meanings are the same. The terms "of," "corresponding," and "corresponding" may be used interchangeably. It should be noted that, when the distinction between them is not emphasized, their intended meanings are the same.

[0060] In the embodiments of the present invention, sometimes a subscript such as W1 may be written as a non-subscript such as W1. When the difference is not emphasized, the meanings to be expressed are the same.

[0061] In order to make the technical problems, technical solutions and advantages to be solved by the present invention clearer, a detailed description will be given below with reference to the accompanying drawings and specific embodiments.

[0062] See also Figure 1 , a national land space planning data analysis system for ecological restoration, the system includes:

[0063] The interference parameter acquisition module obtains remote sensing sequences of surface cover in ecologically sensitive areas, extracts the trend of the normalized vegetation index, calculates the degree of periodic fluctuation, obtains dynamic simulation layers of surface water bodies, extracts the seasonal variation intensity of the runoff variation coefficient, obtains annual interpretation data of land use types, calculates the frequency of land use changes and the land use type diversity index, and generates an information set of ecological disturbance sources;

[0064] The ecological disturbance calculation module is based on the ecological disturbance source information set. According to the annual standard deviation of the normalized vegetation index, the runoff variation coefficient, and the frequency of land use change, it sets the normalized vegetation index fluctuation threshold, the runoff variation warning standard, and the upper limit of the land use change frequency. It determines whether the three disturbance factors in each spatial unit simultaneously meet the corresponding threshold standards. If they do, the unit is marked as an ecological disturbance superposition unit and an ecological degradation disturbance judgment mark is generated.

[0065] The restoration area identification module uses the regional building density layer based on the ecological degradation interference judgment mark to extract the number of buildings per unit area and the building function classification level within the spatial unit. It then determines whether there is overlap between high-level building density areas and ecological interference superposition units. If there is no overlap, the corresponding unit is marked as suitable for ecological restoration, generating an ecological restoration space screening result.

[0066] The priority determination module calls the ecological restoration spatial screening results, extracts the interannual series of the ecological service function index within the target unit, analyzes the decline duration and decline slope, and calculates the ratio compared to the regional average function degradation cycle as the urgency weight. It determines whether it exceeds the urgency baseline. If so, it assigns a high priority level and generates restoration priority classification information.

[0067] The site stability screening module extracts the land use time series matrix of the corresponding spatial unit based on the restoration priority classification information, counts the annual number of land type transfers and the frequency of reverse land function changes, and determines whether there are frequent direction reversals or the transfer frequency exceeds the set threshold. If so, it is determined that the long-term stability of the unit restoration is insufficient and a list of available spatial planning sites is generated;

[0068] The frequency of reverse change of land use is the proportion of land development after restoration;

[0069] The ecological disturbance source information set includes the normalized vegetation index change frequency index, runoff change intensity level, land use diversity ratio and functional change behavior density. The ecological degradation interference judgment identification is specifically the spatial disturbance level label, the judgment rule set and the interference factor cross-record table. The ecological restoration spatial screening results include the restoration feasibility spatial mask, the development interference exclusion layer and the restoration adaptation unit number list. The restoration priority classification information specifically refers to the urgency level code, the ecological compression cycle ratio and the number of restoration buffer zone intervals. The available list of spatial planning site selection includes the stability screening mask, the site selection effective area number and the index table of the corresponding indicators of the exclusion unit.

[0070] See also Figure 2 and Figure 3 The interference parameter acquisition module includes the vegetation index extraction submodule, the hydrological fluctuation identification submodule, and the land use statistics submodule:

[0071] The vegetation index extraction submodule obtains remote sensing sequences of land cover in ecologically sensitive areas, extracts time series data of the normalized vegetation index, calculates the monthly normalized vegetation index change trend in the same area, summarizes the frequency of value changes between months within the year, establishes a spatial distribution layer based on the standard deviation of the annual normalized vegetation index change value, and obtains the fluctuation coefficient of vegetation change trend;

[0072] To obtain the remote sensing sequence of surface cover in ecologically sensitive areas, it is necessary to first specify the boundary range of the target ecological zone and call public remote sensing image data sources such as Sentinel-2 or Landsat 8. In terms of time span, monthly image data of 12 consecutive months are screened. The pixel values of the red band and near-infrared band of each month are extracted respectively using the NDVI calculation method. The normalized vegetation index value of each pixel is calculated using the formula NDVI = (NIR-RED) / (NIR+RED). Then, the spatially corresponding NDVI layer raster data matrix is generated. The NDVI raster images of all months are sorted by month to establish a time series structure. The monthly NDVI value sequence of each spatial unit pixel is differenced to determine the direction of monthly change and record the change frequency. Specifically, the number of positive and negative value changes in the whole year is counted with NDVI rising as positive value and falling as negative value. For example, the NDVI values of a pixel in 12 months are 0.71, 0.68, 0.65, 0.62, 0.60, 0.59, 0.61, 0.64, 0.67, 0.70, 0.72, and 0.73, indicating that the pixel experienced a downward trend 5 times and an upward trend 6 times in the year, with frequencies of 5 / 11 and 6 / 11, respectively. Subsequently, for each pixel's annual NDVI value series, its standard deviation was calculated to reflect its degree of fluctuation within the year. The larger the NDVI standard deviation, the more unstable the vegetation cover change in the area. NDVI fluctuation values > 0.1 were used as the medium fluctuation threshold, and > 0.2 as the severe fluctuation threshold. In specific areas, it can be seen that when the standard deviation of a wetland area reaches 0.28, it falls into the severe fluctuation range. The standard deviation values are then mapped to the corresponding spatial grid to construct a unified spatial distribution layer. This layer represents the distribution of NDVI changes in each region within the year, which is further used for spatial trend analysis of vegetation changes and the establishment of a vegetation change trend fluctuation coefficient.

[0073] The hydrological fluctuation identification submodule uses the surface water dynamic simulation layer based on the vegetation change trend fluctuation coefficient to extract the seasonal unit time runoff change sequence of the main water systems in the basin, calculates the standard deviation of the four-season runoff values in each spatial unit and integrates them into a continuous change sequence table to obtain the seasonal amplitude coefficient of hydrological change;

[0074] Based on the vegetation change trend fluctuation coefficient, it is necessary to extract the main water system range from the surface water dynamic simulation layer in the region, give priority to determining the boundaries of major hydrological units such as rivers, lakes, and wetlands, and then call the hydrological simulation data source such as SWAT or MODFLOW simulation results to extract the average runoff data of each water unit in each season to form a time series structure. The average flow rate corresponding to the spring, summer, autumn and winter seasons is calculated respectively to obtain the degree of runoff variation between the four seasons. The standard deviation calculation formula σ = √(∑(x_i-μ)2 / n) is used to calculate the fluctuation of the runoff values of each water unit in the four seasons. In this example, if the runoff values in spring, summer, autumn and winter are 48, 62, 45, and 51 respectively, the average value is 51.5 and the standard deviation is about 7.3, which constitutes the seasonal variation coefficient of the unit. The results of all water units are integrated to form a continuous change sequence, which is organized into a spatial distribution map of the standard deviation of runoff in the four seasons. All values are uniformly normalized to construct a spatial grid expression structure, which is used to merge with the aforementioned vegetation fluctuation spatial structure to obtain the seasonal amplitude coefficient of hydrological variation.

[0075] The land use statistics submodule collects annual remote sensing interpretation data of land use types based on the seasonal amplitude coefficient of hydrological changes, constructs a sequence of land use changes in each spatial unit, counts the annual conversion times of functional types and the number of type changes, and organizes them into a two-dimensional array to record the annual change trend indicators of each unit, thereby obtaining an information set of ecological disturbance sources;

[0076] According to the seasonal amplitude coefficient of hydrological variation, the interpretation results of annual land use remote sensing patches were collected in each spatial unit. A total of 11 data periods from 2013 to 2023 were used to construct an annual type identification time series matrix. The dominant land type in each spatial unit was recorded annually and coded in a consistent format, such as cultivated land as 1, forest land as 2, construction land as 3, wetland as 4, etc. The type change behavior was identified by year-on-year comparison. For example, if the annual land type sequence of a unit is 1→1→3→3→2→2→2→1→1→4→4, there are 7 type changes, and the change frequency is 7 / 10. Reverse change behavior refers to a type that recovers for a short period of time and then changes again. For example, if 1 appears in 1→3→1→3, recovers and then disappears, it is considered a reverse change. All reverse changes are counted and accumulated in the index sequence. At the same time, the diversity index is calculated based on the number of types of land type change codes in consecutive years. The Shannon diversity index formula H=-∑(p_i×ln p_i) is calculated, where p_i is the proportion of each type. After counting the change indicators of all unit land types, a two-dimensional array structure matrix is formed, and the field is linked with the seasonal amplitude coefficient of hydrological change to complete the spatial normalization integration and summary of all ecological disturbance related factors and obtain the ecological disturbance source information set.

[0077] See also Figure 2 and Figure 4,The ecological disturbance calculation module includes a disturbance factor threshold ,submodule, a disturbance unit screening submodule, and a disturbance level generation ,submodule;

[0078] The disturbance factor threshold setting module is based on the ecological disturbance source information set. According to the vegetation change trend fluctuation coefficient, the hydrological change seasonal amplitude coefficient and the land use change trend index, it sets the normalized vegetation index fluctuation threshold, the runoff change ratio warning benchmark value and the land use change frequency limit, establishes the judgment baseline of the three types of disturbance factors, and obtains the disturbance threshold parameter set.

[0079] Based on the ecological disturbance source information set, it is necessary to first extract three types of disturbance parameters: vegetation change trend fluctuation coefficient, hydrological change seasonal amplitude coefficient, and land use change trend index. In practical applications, taking a certain ecological unit as an example, its NDVI fluctuation standard deviation is 0.17, which belongs to the medium disturbance intensity range, its hydrological change amplitude standard deviation is 12.4, and the annual change frequency of land use type is 0.5 times / year. For these three data, the judgment thresholds need to be set separately. The setting method of the normalized vegetation index fluctuation threshold can be constructed by constructing a normalized expression through the maximum and average value of the NDVI standard deviation of the entire region, and the fluctuation threshold is set to 1.2 of the average value of the NDVI standard deviation of the entire region. times, setting the regional average to 0.12, the threshold is 0.144, the runoff change ratio warning benchmark value is set to the average of the four-season runoff standard deviation multiplied by 1.1 times. If the average is 11.0, it is set to 12.1. The land use change frequency limit is set to the 95th percentile value of the regional unit annual average change frequency, that is, all unit frequencies are sorted, and the 95th percentile position value is taken as the upper limit. For example, if the median of all frequencies is 0.36 and the 95th percentile is 0.68, 0.68 is used as the limit threshold. According to the above standards, a judgment baseline matrix is constructed, and the three values and their corresponding judgment boundaries are written into the parameter table. They are bound to their respective spatial indexes to generate a disturbance threshold parameter set.

[0080] The disturbance unit screening submodule calls the disturbance threshold parameter set to judge the data of each spatial unit in the ecological disturbance source information set, identifying whether the annual standard deviation of the normalized vegetation index is greater than the set fluctuation threshold, whether the runoff variation coefficient is greater than the set benchmark, and whether the frequency of land use change exceeds the limited value. It selects the spatial units that meet all three conditions and generates a list of composite disturbance spatial units.

[0081] After calling the disturbance threshold parameter set, each spatial unit in the ecological disturbance source information set is judged. The NDVI standard deviation, hydrological fluctuation coefficient, and annual land use change frequency of the unit need to be taken out respectively to determine whether the NDVI standard deviation is greater than the normalized fluctuation threshold of 0.144. If the value of a unit is 0.17, it is judged to be satisfied. If the hydrological fluctuation is 12.4, which is greater than the benchmark 12.1, it is also judged to be satisfied. If the annual land use change frequency is 0.52, which is less than the limit of 0.68, it is judged to be unsatisfied. Only when all three judgments are "satisfied" can the unit be marked as a composite disturbance unit. In actual operation, the spatial judgment logic nested expression method is used to set the rules, such as the logical condition expression NDVI_SD_i>0.144AND HY_VAR_i>12.1AND LU_FREQ_i>0.68, where NDVI_SD_i represents the standard deviation of NDVI for unit i, HY_VAR_i represents the hydrological fluctuation value, and LU_FREQ_i represents the annual land change frequency. Each unit in the spatial grid corresponds to three values. After performing cell-by-cell judgment on the grid, the result is marked as 1 or 0. All spatial units with a value of 1 are screened out to form the result grid. The corresponding spatial unit index number is further extracted to generate a list of composite disturbance spatial units.

[0082] The disturbance level generation submodule assigns a disturbance identification number to each unit based on the list of composite disturbance spatial units and marks it as an ecological disturbance superposition unit. It records the corresponding three types of disturbance factor structures and maps them to the spatial data grid layer to establish an ecological degradation disturbance judgment mark.

[0083] According to the list of composite disturbance spatial units, a disturbance identification number is assigned to each spatial unit that meets the conditions. The numbers are assigned in order of index. For example, if the list is G021, G045, and G103, they are numbered 001, 002, and 003, respectively. Each unit is bound to the three corresponding parameter values in the ecological disturbance source information set, including NDVI standard deviation, hydrological fluctuation amplitude coefficient, and annual land use frequency. The structural archiving of the three factor parameters is achieved through field connection to form a detailed table of unit disturbance data. At the same time, the numbering results are mapped to a raster layer consistent with the original spatial data structure. In the raster, each unit that meets the conditions is assigned a corresponding number, and the remaining units are assigned 0 to indicate that they are not identified. The number distribution is expressed in color-graded symbols in the layer, and a spatial visualization expression result is constructed, and finally an ecological degradation disturbance determination mark is established.

[0084] See also Figure 2 and Figure 5 ,The restoration area identification module includes the building density extraction submodule, the ,functionality level comparison submodule, and the suitable zoning determination submodule;

[0085] The building density extraction submodule uses the regional building density layer based on the ecological degradation interference identification, extracts the number of buildings per unit area in each spatial unit, calculates the ratio of the total number of buildings to the area in each unit, and summarizes the spatial building density information;

[0086] Based on the ecological degradation interference identification, it is necessary to call the regional building density layer. This layer is generated by superimposing the measured building vector data on the grid. The grid resolution is set to 100 meters × 100 meters. The space is divided according to this scale. The actual number of buildings in each grid cell is counted to calculate the building density per unit area. The formula D = N / A is used, where D is the building density, N is the number of buildings, and A is the grid area. For a standard area of 10,000 square meters, if there are 8 buildings, the density is 0.0008 buildings / square meter. Then, by summing the building outline areas, the total building area per unit S_total is calculated. S_total is then ratioed with the grid area to obtain the building coverage ratio B = S_total / A. Assuming the total building area is 4,200 square meters, the coverage ratio is 0.42. This ratio is summarized and written into the grid attribute table. A two-dimensional density matrix is constructed for all spatial units to form a spatial distribution building density expression structure, and spatial building density information is established based on this.

[0087] The functional level comparison submodule matches the building function classification level layer based on the spatial building density information, extracts the corresponding building use level code in each type of spatial unit, determines whether the number of buildings per unit area is higher than the building density benchmark value, and whether its building use level is within the upper limit of the development level classification, and selects the corresponding units that meet the conditions to obtain the building function interference superposition unit set;

[0088] According to the spatial building density information, it is necessary to match the building function classification level layer. Each building object in this layer contains a use field to indicate its function level. The field code is 1 for residential, 2 for commercial, 3 for industrial, 4 for office, and 5 for infrastructure. The levels are divided into low (12), medium (34), and high (5) according to the degree of development. If each spatial unit contains multiple-use buildings, its level proportion needs to be counted and the maximum use level needs to be extracted. In the judgment step, a building density benchmark value needs to be set. This value is based on the 85th percentile of the density value of all units in the area. The 85th percentile density value in the entire unit is set to 0.0011 buildings / square meter. This value is used as the density benchmark value. Then, it is judged whether the density of each unit is greater than this value, and whether its maximum use level is greater than or equal to 5. If the density of a unit is 0.0013 and the maximum use level is 5, then both conditions are met and recorded in the interference set. If the density of a unit is 0.0007, then the conditions are not met. The result set only contains spatial units that meet the two conditions, and finally the building function interference superposition unit set is obtained.

[0089] The suitable zoning submodule calls the building function interference superposition unit set and performs spatial superposition operation with the ecological degradation interference judgment mark to determine whether the two overlap in space. Spatial units with overlapping relationships are excluded, and the remaining spatial units are assigned ecological restoration suitability marks to establish the ecological restoration space screening results.

[0090] The building function interference superposition unit set is called and spatially superimposed with the ecological degradation interference judgment mark. Geometric coincidence judgment is performed in the spatial database. The operation method is to query whether each interference superposition unit exists within the boundary of the ecological degradation interference unit in each interference superposition unit. The spatial judgment expression intersects(geometry_A,geometry_B)=TRUE is used. If the spatial units of the two layers coincide, the unit is excluded; if not, it is retained. After eliminating all grid numbers that meet the spatial overlap, the remaining ones are the feasible restoration units in the non-interference overlapping area. The attribute tag "RESTORABLE=YES" is added to these units, and they are uniformly output as a single layer structure and assigned a unified code. Finally, the ecological restoration spatial screening results are established.

[0091] See also Figure 2 and Figure 6 ,The priority level determination module includes a function sequence extraction submodule, a decay period calculation submodule, and an urgency level assignment submodule;

[0092] The functional sequence extraction submodule calls the ecological restoration spatial screening results, extracts the ecological service function index time series corresponding to the spatial unit, arranges the inter-annual value changes, establishes the functional sequence distribution table under each unit, and obtains the inter-annual change trend set of ecological functions;

[0093] To call the ecological restoration spatial screening results, it is necessary to first clarify the unique index number of each spatial unit in the screening results, and combine it with the Ecological Service Function Index (ESFI) database to obtain its corresponding inter-annual time series value. The index usually represents the functional integrity level as a decimal between 0 and 1. The year coverage is set to 2010 to 2023, a total of 14 years. In the GIS platform, the screening unit is linked to the ESFI data field through spatial linkage. The functional index value of each unit is extracted year by year, and the functional value table data with the structure of {unit ID, year 1 value, year 2 value, ... year 14 value} is constructed. It is arranged in chronological order and stored as a two-dimensional array structure. Then, the function index is executed for each unit. Identify the change trend, determine whether there is an increase or decrease in function value between adjacent years, record the change direction and magnitude, and calculate the interannual change rate at the same time to construct a function sequence distribution table for each unit. For example, if the ESFI of a unit from 2010 to 2023 is 0.83, 0.81, 0.79, 0.77, 0.75, 0.76, 0.74, 0.72, 0.70, 0.69, 0.70, 0.71, 0.72, and 0.74, then the unit has continued to decline in the first five years and then fluctuated and rebounded. The length of the decline period can be recorded as 5 years, and the maximum interannual decline is 0.02. The statistical results are bound to the unit ID to establish the interannual change trend set of ecological function.

[0094] The decay cycle calculation submodule identifies the time periods of consecutive annual functional value decline based on the interannual change trend set of ecological functions. It calculates the functional succession urgency value of the spatial unit according to the duration of the decline of each spatial unit and the ratio of the decline in each period.

[0095] The formula for calculating the functional succession urgency value of a spatial unit is:

[0096]

[0097] Among them, U i It represents the urgency of functional succession of the i-th spatial unit, which is used to reflect the urgency of ecological function decline of the spatial unit. i represents the normalized value of the duration of the functional decline of the i-th spatial unit, represents the normalized value of the regional average functional degradation period, S i V represents the normalized value of the interannual variation slope of the function of the i-th spatial unit. i represents the fluctuation variance of the ecological function value of the i-th spatial unit;

[0098] This indicator is used to assess whether the severity of degradation and degradation rate of the ecological function succession trend of each spatial unit constitutes an urgent need for restoration. It is a key supporting quantitative factor for constructing a restoration priority determination model. Indicates whether the degradation time of the unit is longer than the average degradation period. A larger value indicates a more persistent degradation. It represents the comprehensive weight of the interannual decline rate and stability. The faster the decline and the lower the fluctuation, the larger the item, which means the function is deteriorating rapidly and continuously. i The larger the value, the longer-term, faster, and more stable the ecological function of the spatial unit is declining, and it needs to be prioritized in the restoration target. It is a comprehensive indicator to measure the "urgency of ecological restoration."

[0099] Based on the interannual trend set of ecological function changes, it is necessary to identify the functional decline period of each spatial unit within consecutive years. By judging whether the difference between the functional index of each year and the previous year is negative, the number of consecutive negative years is counted and recorded as the decline duration period T. i , and then calculate the average annual functional decline S during the decline period i Variance V i , calculate the functional succession urgency value U of the spatial unit i , using the formula:

[0100]

[0101] Among them, assuming that the functional index of a unit dropped from 0.85 in 2011 to 0.73 in 2016, then T i =6, the average regional degradation cycle is 7 years, that is, Therefore The average annual decrease is (0.85-0.73) / 6≈0.02, and the standard deviation of the function value is calculated using the formula Calculate, set to 0.015, then V i =0.000225, substituting into the formula:

[0102]

[0103] In this way, the urgency scores of all spatial units are completed to obtain the urgency values of functional succession;

[0104] The urgency level assignment submodule sets the urgency baseline standard value based on the functional succession urgency value of the spatial unit, compares the urgency value of each spatial unit, identifies the spatial units that exceed the standard and marks them as high-level intervention targets, and obtains the restoration priority classification information;

[0105] According to the functional succession urgency value, the urgency baseline standard value needs to be set. This value is set as the 80th percentile of the urgency values of all spatial units in the entire area. That is, after the urgency values of all units are arranged in ascending order, the value at the 80th position is taken as the judgment standard. For example, there are 1,000 units in the area, and the 800th value after sorting is 1.42. This value is used as the urgency benchmark, and the urgency value of each spatial unit is judged in turn to see whether it is greater than 1.42. If the urgency of a unit is 1.66, it is determined to be a high urgency object, and its level is marked as "Priority Level I" in the attribute field. If it is less than or equal to the benchmark value, it is "Non-priority" or "Level II". Finally, the urgency level information of each spatial unit is written into the spatial layer and output into the classification field to establish the restoration priority classification information.

[0106] See also Figure 2 and Figure 7 ,The site selection stability screening module includes the utilization trajectory extraction submodule, the change frequency calculation submodule, and the ,stability determination and screening submodule;

[0107] The trajectory extraction submodule is used to extract the land use time series data of the corresponding spatial unit according to the restoration priority classification information, construct the spatial land use type distribution information by year, obtain the land type change records of consecutive years within the unit, and establish the land use annual transfer sequence set;

[0108] According to the restoration priority classification information, the spatial units marked as "Priority I" are first screened out, and their corresponding spatial codes are extracted. The dominant land class code information of the corresponding units in consecutive years is queried in the land use time series database. The data is usually constructed into a matrix with year as the horizontal axis and unit ID as the vertical axis. For example, from 2010 to 2023, there are 14 years. The spatial unit number is G_048. In each year, the land class codes are: 2, 2, 3, 3, 1, 1, 3, 2, 2, 2, 3, 3, 2, 1. Then, the annual land use sequence of the unit is established in sequence. Then, for each unit, the land class codes of adjacent years in the sequence are scanned year by year to see if they have changed. All changes are marked as 1 and no changes are marked as 0. These are recorded as change trajectory values. At the same time, the change type pairs between each year, such as (2→3) and (3→1), are retained to form a two-dimensional transfer record array. Subsequently, the annual change trajectories and transfer pairs of all units are constructed into a database table. Each row corresponds to a spatial unit, and each column records the land class change information, thereby establishing a land use annual transfer sequence set.

[0109] The change frequency calculation submodule is based on the land use annual transfer sequence set, and counts the number of land use type changes and the cumulative frequency of reverse land function changes in each spatial unit in adjacent years, and calculates the land function change intensity index of each unit;

[0110] The formula for calculating the land function change intensity index of each unit is:

[0111]

[0112] Among them, X j represents the land function change intensity index of the jth spatial unit, n j represents the annual length included in the jth spatial unit, R j,t F represents the number of reverse changes that occurred in the j-th spatial unit in the t-th year, j,t represents the land function code of the jth spatial unit in year t, Z j,t represents the total number of type changes that occurred in the jth spatial unit in the tth year;

[0113] This value is used to measure the degree of land use instability of a spatial unit over time, that is, whether there are frequent reverse changes or use substitution behaviors, reflecting its "long-term planning sustainability" as a restoration area. The numerator reflects two unstable factors: the frequency of reverse change behaviors and the degree of inter-annual jump in land use functions; the denominator takes into account the overall change frequency (to prevent the amplification of the impact of a single change); the formula is an annual average and has the ability to normalize over time. j The larger the value, the more frequent changes in the land use, the difficulty in implementing the policy, and the lack of stability in the use of the space after ecological restoration. Therefore, it is used to eliminate unstable sites.

[0114] Based on the annual land use transfer sequence set, the number of land use type changes and the frequency of reverse land function changes in adjacent years in each spatial unit are counted respectively. The statistical logic of the number of changes is to determine whether the land class codes in adjacent years are unequal. If they are not equal, it is regarded as a type change behavior. For example, if the sequence is 2, 2, 3, 3, 1, 1, 3, 2, 2, 2, 3, 3, 2, 1, there are 9 changes in total, and the cumulative number is Z. j,t =9, and then determine whether there is a reverse behavior from one type to another and then back to the previous type. For example, 2→3→2 and 3→1→3 are typical reverse change structures. 2→3→2 and 3→1→3 are identified once each in the sequence, and the cumulative R j,t =2, then according to the absolute value of the difference in land function code for each year |F j,t -F j,t-1 To calculate the magnitude of the change, set the year of the reverse change as t, and take the difference between the land classification codes of that year and the previous year. If the codes are 3 and 1, the difference is 2. Substitute this into the formula:

[0115]

[0116] Set n of space unit number G_048 j =13, that is, 13 pairs of adjacent years, Rj,t =2, the sum of the differences in each year is 18, Z j,t =9, substituting into After the operation is completed, the corresponding spatial unit field is written to obtain the land function change intensity index;

[0117] The stability determination and screening submodule sets the land use stability benchmark standard based on the land function change intensity index of each unit, compares the intensity value of each unit, screens the spatial locations with frequent continuous fluctuations or concentrated reverse changes, eliminates the corresponding spatial units, and establishes a list of available spatial planning sites;

[0118] According to the land function change intensity index, the land use stability benchmark standard value is set, and the 85th percentile value of the change intensity value of the global spatial unit is used as the judgment baseline. For example, all units X j Sort from small to large, extract the value of the 85th position as the stability benchmark, assuming that the value is 0.137, and all units with a change intensity index greater than 0.137 are judged as areas with weak stability. Then, count whether these spatial units have multiple consecutive reverse change behaviors. The judgment rule is that if there are more than two reverse behaviors within three years, it is an unstable trend. Such units are removed from the spatial restoration candidates and marked as "UNSTABLE" in the result layer. The remaining spatial units are retained and assigned the value "USABLE", and all available spatial units are summarized to form an output list, and finally an available list for spatial planning site selection is established.

[0119] It should be understood that the term "and / or" as used herein simply describes a relationship between associated objects, indicating that three possible relationships exist. For example, "A and / or B" can represent: A alone, A and B together, or B alone. A and B can be singular or plural. Furthermore, the character " / " as used herein generally indicates an "or" relationship between the associated objects, but it may also indicate an "and / or" relationship. For specific understanding, please refer to the context.

[0120] In this disclosure, "at least one" means one or more, and "plurality" means two or more. "At least one of the following" or similar expressions refers to any combination of these items, including any combination of single or plural items. For example, "at least one of a, b, or c" can mean: a, b, c, ab, ac, bc, or abc, where a, b, and c can be single or plural.

[0121] It should be understood that in various embodiments of the present invention, the size of the serial numbers of the above-mentioned processes does not mean the order of execution. The execution order of each process should be determined by its function and internal logic, and should not constitute any limitation on the implementation process of the embodiments of the present invention.

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

[0123] Those skilled in the art will clearly understand that, for the convenience and brevity of description, the specific working processes of the above-described equipment, devices and units can refer to the corresponding processes in the aforementioned method embodiments and will not be repeated here.

[0124] In the several embodiments provided by the present invention, it should be understood that the disclosed devices, apparatuses and methods can be implemented in other ways. For example, the device embodiments described above are merely illustrative. For example, the division of the units is merely a logical function division. In actual implementation, there may be other division methods, such as multiple units or components can be combined or integrated into another device, or some features can be ignored or not executed. Another point is that the mutual coupling or direct coupling or communication connection shown or discussed can be through some interface, indirect coupling or communication connection of the device or unit, which can be electrical, mechanical or other forms.

[0125] The units described as separate components may or may not be physically separate, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed across multiple network units. Some or all of these units may be selected to achieve the purpose of this embodiment according to actual needs.

[0126] In addition, each functional unit in each embodiment of the present invention may be integrated into one processing unit, or each unit may exist physically separately, or two or more units may be integrated into one unit.

[0127] If the functions are implemented in the form of software functional units and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, or the part that contributes to the prior art, or the part of the technical solution, can be embodied in the form of a software product. The computer software product is stored in a storage medium and includes several instructions for enabling a computer device (which can be a personal computer, server, or network device, etc.) to execute all or part of the steps of the method described in each embodiment of the present invention. The aforementioned storage medium includes various media that can store program codes, such as a USB flash drive, a mobile hard disk, a read-only memory (ROM), a random access memory (RAM), a magnetic disk, or an optical disk.

[0128] The above description is merely a specific embodiment of the present invention, but the scope of protection of the present invention is not limited thereto. Any modifications or substitutions that can be easily conceived by a person skilled in the art within the technical scope disclosed in the present invention should be included in the scope of protection of the present invention. Therefore, the scope of protection of the present invention should be based on the scope of protection of the claims.

Claims

1. A national land space planning data analysis system for ecological restoration, characterized in that: The system comprises: The interference parameter acquisition module obtains remote sensing sequences of land cover in ecologically sensitive areas, extracts the trend of normalized vegetation index changes, calculates the frequency of land use changes and the land use type diversity index, and generates an information set of ecological disturbance sources; The ecological disturbance calculation module sets the normalized vegetation index fluctuation threshold, the runoff variation warning standard, and the upper limit of the land use change frequency based on the ecological disturbance source information set, and determines whether the three disturbance factors in each spatial unit simultaneously meet the corresponding threshold standards. If so, it is marked as an ecological disturbance superposition unit and an ecological degradation disturbance judgment mark is generated; The restoration area identification module extracts the number of buildings per unit area and the building function classification level within the spatial unit based on the ecological degradation interference judgment mark, determines whether there is a high-level building density area that overlaps with the ecological interference superposition unit, and generates an ecological restoration space screening result; The priority determination module calls the ecological restoration space screening results, extracts the inter-annual series of the ecological service function index in the target unit, calculates the ratio compared to the regional average function degradation cycle as the urgency weight, and generates restoration priority classification information.

2. The national land space planning data analysis system for ecological restoration according to claim 1 is characterized in that: The ecological disturbance source information set includes the normalized vegetation index change frequency index, runoff change intensity level, land use diversity ratio and functional change behavior density; the ecological degradation interference judgment identifier is specifically a spatial disturbance level label, a judgment rule set and an interference factor cross-record table; the ecological restoration spatial screening result includes a restoration feasibility space mask, a development interference exclusion layer and a restoration adaptation unit number list; the restoration priority classification information specifically refers to the urgency level code, the ecological compression cycle ratio and the number of recovery buffer zone intervals.

3. The national land space planning data analysis system for ecological restoration according to claim 2 is characterized in that: The interference parameter acquisition module includes: The vegetation index extraction submodule obtains remote sensing sequences of land cover in ecologically sensitive areas, extracts time series data of the normalized vegetation index, calculates the monthly normalized vegetation index change trend in the same area, summarizes the frequency of value changes between months within the year, establishes a spatial distribution layer based on the standard deviation of the annual normalized vegetation index change value, and obtains the fluctuation coefficient of vegetation change trend; The hydrological fluctuation identification submodule uses the surface water dynamic simulation layer based on the vegetation change trend fluctuation coefficient to extract the seasonal unit time runoff change sequence of the main water systems in the basin, calculates the standard deviation of the four-season runoff values in each spatial unit and integrates them into a continuous change sequence table to obtain the seasonal amplitude coefficient of the hydrological change; The land use statistics submodule collects annual remote sensing interpretation map data of land use types based on the seasonal amplitude coefficient of hydrological changes, constructs a land use change sequence of consecutive years in each spatial unit, counts the annual conversion times of functional types and the number of type changes, organizes them into a two-dimensional array to record the annual change trend indicators of each unit, and obtains the information set of ecological disturbance sources.

4. The national land space planning data analysis system for ecological restoration according to claim 3 is characterized in that: The ecological disturbance calculation module includes: The disturbance factor threshold setting module sets the normalized vegetation index fluctuation threshold, the runoff change ratio warning benchmark value, and the land use change frequency limit based on the ecological disturbance source information set, according to the vegetation change trend fluctuation coefficient, the hydrological change seasonal amplitude coefficient, and the land use change trend index, respectively, to establish the judgment baseline for the three types of disturbance factors and obtain the disturbance threshold parameter set; The interference unit screening submodule calls the disturbance threshold parameter set to judge the data of each spatial unit in the ecological disturbance source information set, identifying whether the annual standard deviation of the normalized vegetation index is greater than the set fluctuation threshold, whether the runoff variation coefficient is greater than the set benchmark, and whether the frequency of land use change exceeds the limited value, screening the spatial units that meet all three conditions, and generating a list of composite disturbance spatial units; The interference level generation submodule assigns an interference identification number to each unit according to the composite disturbance spatial unit list and marks it as an ecological disturbance superposition unit, records the corresponding three types of disturbance factor structures and maps them to the spatial data grid layer, and establishes an ecological degradation interference judgment mark.

5. The national land space planning data analysis system for ecological restoration according to claim 4 is characterized in that: The repair area identification module includes: The building density extraction submodule calls the regional building density layer based on the ecological degradation interference judgment mark, extracts the number of buildings per unit area in each spatial unit, calculates the ratio of the total number of buildings to the area in each unit, and summarizes the spatial building density information; The function level comparison submodule matches the building function classification level layer according to the spatial building density information, extracts the corresponding building use level code in each type of spatial unit, determines whether the number of buildings per unit area is higher than the building density benchmark value, and whether its building use level is within the upper limit of the development level classification, and selects the corresponding units that meet the conditions to obtain the building function interference superposition unit set; The suitable zoning submodule calls the building function interference superposition unit set and performs spatial superposition operation with the ecological degradation interference judgment mark to determine whether the two overlap in space, excludes the spatial units with overlapping relationships, assigns the remaining spatial units with ecological restoration suitability marks, and establishes the ecological restoration space screening results.

6. The national land space planning data analysis system for ecological restoration according to claim 5 is characterized in that: The priority determination module includes: The function sequence extraction submodule calls the ecological restoration space screening results, extracts the ecological service function index time series corresponding to the spatial unit, arranges the inter-annual value changes, establishes a function sequence distribution table under each unit, and obtains the inter-annual change trend set of ecological functions; The decay cycle calculation submodule identifies the time periods of consecutive annual functional value decline based on the inter-annual variation trend set of ecological functions, calculates the functional succession urgency value of the spatial unit according to the duration of decline of each spatial unit and the ratio of decline in each period; The urgency level assignment submodule sets the urgency baseline standard value according to the functional succession urgency value of the spatial unit, compares the urgency value of each spatial unit, identifies the spatial units that exceed the standard and marks them as high-level intervention objects, and obtains the restoration priority classification information.

7. The national land space planning data analysis system for ecological restoration according to claim 6 is characterized in that: The formula for calculating the functional succession urgency value of the spatial unit is: Among them, U i represents the functional succession urgency value of the i-th spatial unit, T i represents the normalized value of the duration of the functional decline of the i-th spatial unit, represents the normalized value of the regional average functional degradation period, S i V represents the normalized value of the interannual variation slope of the function of the i-th spatial unit. i Represents the fluctuation variance of the ecological function value of the i-th spatial unit.

8. The national land space planning data analysis system for ecological restoration according to claim 7 is characterized in that: The system further comprises: The site stability screening module extracts the land use time series matrix of the corresponding spatial unit based on the restoration priority classification information, counts the annual number of land type transfers and the frequency of reverse changes in land function, and determines whether there are frequent direction reversals or the transfer frequency exceeds the set threshold. If so, it is determined that the long-term stability of the unit restoration is insufficient and a list of available spatial planning sites is generated; The available list of spatial planning site selection includes a stability screening mask, a valid site selection area number, and an index table of corresponding indicators for eliminated units.

9. The national land space planning data analysis system for ecological restoration according to claim 8, characterized in that: The site stability screening module includes: Using the trajectory extraction submodule to extract the land use time series data of the corresponding spatial unit according to the restoration priority classification information, construct the spatial land use type distribution information by year, obtain the land type change records of consecutive years in the unit, and establish the land use annual transfer sequence set; The change frequency calculation submodule counts the number of land use type changes and the cumulative frequency of reverse land function changes in each spatial unit in adjacent years based on the land use annual transfer sequence set, and calculates the land function change intensity index of each unit; The stability determination and screening submodule sets the land use stability benchmark standard according to the land function change intensity index of each unit, compares the intensity value of each unit, screens the spatial locations with frequent continuous fluctuations or concentrated reverse changes, eliminates the corresponding spatial units, and establishes an available list of spatial planning sites.

10. The national land space planning data analysis system for ecological restoration according to claim 9, characterized in that: The formula for calculating the land function change intensity index of each unit is: Among them, X j represents the land function change intensity index of the jth spatial unit, n j represents the annual length included in the jth spatial unit, R j,t F represents the number of reverse changes that occurred in the j-th spatial unit in the t-th year, j,t represents the land function code of the jth spatial unit in year t, Z j,t It represents the total number of type changes that occurred in the j-th spatial unit in the t-th year.

Images