Method and system for dynamically evaluating ecological value of astragalus sinicus based on multi-dimensional stripping algorithm

Abstract

The invention discloses a milk vetch ecological value dynamic evaluation method and system based on a multi-dimensional stripping algorithm, relates to the technical field of agricultural ecosystem evaluation, and effectively weakens background interference of different crop rotation conditions on ecological evaluation results by constructing a crop rotation background reference sequence and distinguishing continuous and discontinuous crop rotation states. The comparability and stability of ecological value evaluation are improved; astragalus sinicus turning and pressing operation is used as a disturbance stripping anchor point, and an alignment time window is constructed before and after turning and pressing, so that fine decoupling of an artificial farming operation disturbance effect and a long-term ecological change trend is realized, and the accuracy of turning and pressing ecological effect identification is enhanced; by performing parallel dimension stripping and unified time alignment on a gas absorption behavior and a gas emission behavior, positive and negative ecological effects are synthesized under the same time scale to form a gas regulation net effect sequence, so that ecological value bias caused by single-index analysis is avoided.

Description

Technical Field

[0001] This invention relates to the field of agricultural ecosystem assessment technology, and in particular to a method and system for dynamic assessment of the ecological value of milkvetch based on a multidimensional stripping algorithm. Background Technology

[0002] With the continuous advancement of green agricultural development and the "dual-carbon" strategy, the carbon sink function and ecological value assessment of farmland ecosystems have gradually become important research directions in agricultural ecology, agricultural resource utilization, and environmental science. Milk clover, a traditional green manure crop widely used in rice-growing areas of southern my country, can significantly improve soil physicochemical properties and increase organic matter content when incorporated into rice rotation systems, and it also has a significant impact on the absorption and emission of greenhouse gases from farmland. In recent years, research methods for assessing the ecological effects of milk clover rotation patterns have gradually expanded from single yield or soil nutrient indicators to quantitative analysis of emission characteristics of gases such as methane and nitrous oxide, as well as comprehensive ecological benefits. At the technical level, the application of ecological observation networks, long-term fixed-location experiments, and multi-source time series analysis methods has made it possible to obtain high-temporal-resolution ecological data for continuous multi-rotation cycles, laying a data foundation for conducting dynamic ecological value assessments.

[0003] However, existing ecological effect assessment techniques for milkvetch mostly employ periodic averaging or static comparative analysis methods, typically using changes in indicators over a fixed period before and after incorporation as the evaluation criterion. This makes it difficult to accurately characterize the temporal disturbance effect of incorporation, a key human intervention. Furthermore, ecological observation data under milkvetch rotation conditions often simultaneously incorporate multiple factors such as changes in rotation background, climate fluctuations, and agricultural operations, exhibiting significant non-stationarity and multidimensional coupling characteristics in their underlying driving mechanisms. Existing methods, when processing such data, generally lack the ability to distinguish between differences in rotation continuity and fail to effectively decouple incorporation disturbance from the long-term rotation background, making the assessment results susceptible to background noise interference. Summary of the Invention

[0004] In view of the problems existing in the ecological value assessment technology of milkvetch, this invention is proposed.

[0005] Therefore, the problem to be solved by this invention is how to extract and reconstruct the net effect of positive and negative ecological effects, and solve the problems of ecological value assessment results being greatly affected by background interference and discontinuous temporal expression in the prior art.

[0006] To solve the above-mentioned technical problems, the present invention provides the following technical solution: In a first aspect, this invention provides a method for dynamic assessment of the ecological value of milkvetch based on a multidimensional stripping algorithm. The method includes: collecting ecological observation records within multiple consecutive rotation cycles under milkvetch rotation conditions in rice-growing areas of southern China; classifying these records according to whether the rotation is continuous to form a rotation background baseline sequence; using the rotation background baseline sequence as input, setting the time point corresponding to the milkvetch overturning operation as the disturbance stripping anchor point; constructing aligned time windows before and after the overturning operation; segmenting and stripping the ecological observation records to generate an overturning disturbance response sequence; under the unified time alignment condition of the overturning disturbance response sequence, performing parallel dimensional stripping on the ecological observation records reflecting gas absorption and emission behaviors to obtain a gas regulation net effect sequence reflecting the synchronous changes of positive and negative ecological effects; temporally mapping the gas regulation net effect sequence and the rotation background baseline sequence; constraining and reconstructing the changes in ecological effects at each time period according to the rotation continuity status; and outputting a dynamic assessment result of the evolution of milkvetch's ecological value over time.

[0007] As a preferred embodiment of the dynamic assessment method for the ecological value of milkvetch based on a multidimensional stripping algorithm described in this invention, the formation of the crop rotation background benchmark sequence includes: according to the actual crop rotation arrangement of milkvetch in the southern rice-growing areas, edge computing nodes deployed in the fields collect and organize the ecological observation sensor data of the same plot in real time and upload it to the cloud or central server, arranging it according to the crop rotation cycle to construct the original observation sequence; in the original observation sequence, based on whether there is a lack of milkvetch planting or an interruption in the planting sequence between adjacent crop rotation cycles, the original observation sequence is judged for crop rotation continuity to form a continuous crop rotation observation sequence and a non-continuous crop rotation observation sequence; the corresponding time identifiers of the continuous crop rotation observation sequence and the non-continuous crop rotation observation sequence are retained respectively, and aligned according to a unified time axis to obtain the crop rotation background benchmark sequence.

[0008] As a preferred embodiment of the dynamic assessment method for the ecological value of milkvetch based on a multidimensional stripping algorithm described in this invention, the setting of the disturbance stripping anchor point includes: for the rotation cycles contained in the continuous rotation observation sequence and the non-continuous rotation observation sequence respectively, when the milkvetch over-cultivation operation occurs, the field edge computing node captures or receives the event signal, immediately takes the time of occurrence as the disturbance stripping anchor point, and generates indices for the time window before over-cultivation and the time window after over-cultivation according to a preset fixed length.

[0009] As a preferred embodiment of the dynamic assessment method for the ecological value of milkvetch based on a multidimensional stripping algorithm described in this invention, the generation of the overburden disturbance response sequence includes: according to the indices of the time window before and after overburden, the cloud or central server segments the corresponding ecological observation records in the crop rotation background baseline sequence, and collects the ecological observation records before and after the overburden operation to form an overburden disturbance segmented sequence; the overburden disturbance segmented sequence is arranged according to the crop rotation cycle, and the time identifier corresponding to the crop rotation continuity state is retained to generate an overburden disturbance response sequence used to characterize the impact of the overburden operation disturbance.

[0010] As a preferred embodiment of the dynamic assessment method for the ecological value of milkvetch based on a multidimensional stripping algorithm described in this invention, the parallel dimensional stripping of ecological observation records reflecting gas absorption and emission behaviors includes: extracting ecological observation records of gas absorption behavior under the unified time alignment condition of the overturning disturbance response sequence to form a time-continuous gas absorption effect sequence; simultaneously, extracting ecological observation records of gas emission behavior to form a time-continuous gas emission effect sequence; standardizing the gas absorption effect sequence and the gas emission effect sequence based on their respective corresponding global warming potential equivalent values ​​of greenhouse gases; marking the standardized gas absorption effect sequence as a positive value sequence and the gas emission effect sequence as a negative value sequence, while maintaining strict synchronization of the two sequences on the time axis.

[0011] As a preferred embodiment of the dynamic assessment method for the ecological value of milkvetch based on a multidimensional stripping algorithm described in this invention, the parallel dimensional stripping of ecological observation records reflecting gas absorption and emission behaviors further includes: under the unified time alignment condition, using the gas absorption effect sequence marked with positive values ​​and the gas emission effect sequence marked with negative values ​​as parallel inputs, performing algebraic sum calculation of the positive and negative sequences at each identical time point to generate an initial gas regulation effect composite sequence; using the rotation continuity state corresponding to the overturning disturbance response sequence as a constraint, smoothing the initial gas regulation effect composite sequence to generate a net gas regulation effect sequence; and re-associating the time identifier of the net gas regulation effect sequence with the disturbance stripping anchor point and the corresponding rotation cycle to form a time evolution sequence.

[0012] As a preferred embodiment of the dynamic assessment method for the ecological value of milkvetch based on a multidimensional stripping algorithm described in this invention, the dynamic assessment result of the evolution of the ecological value of milkvetch over time includes: using the time axis of the crop rotation background baseline sequence as a unified benchmark, mapping the time identifier of each data point in the net gas regulation effect sequence to the corresponding disturbance time window position within the source crop rotation cycle; wherein, within each crop rotation cycle, the mapped time range is jointly defined by the pre-overturning time window and the post-overturning time window, and the mapped net gas regulation effect value is bound to the continuity state of the crop rotation cycle; based on the continuity state bound to each crop rotation cycle after mapping... Reconstruction constraints are applied to the net gas regulation effect values ​​within the corresponding time periods: for the rotation cycles corresponding to continuous cropping observation sequences, the net gas regulation effect value is directly used as the effective estimate of the ecological value change within the time period and marked as a high-confidence data segment; for the rotation cycles corresponding to non-continuous cropping observation sequences, the net gas regulation effect value needs to be multiplied by the reconstruction weight coefficient to generate a corrected estimate of the ecological value change and marked as a low-confidence data segment; the ecological value change estimates within all rotation cycles after constraint reconstruction are arranged and connected according to the mapped time order to form a complete time series, which is output as the dynamic assessment sequence of the ecological value of milkvetch.

[0013] Secondly, this invention provides a dynamic assessment system for the ecological value of *Astragalus membranaceus* based on a multidimensional stripping algorithm, comprising: The crop rotation baseline construction module is used to collect ecological observation records within multiple consecutive crop rotation cycles under the condition of purple clover rotation in the rice-growing areas of southern China, and classify the ecological observation records according to whether the crop rotation is continuous, forming a crop rotation background baseline sequence. The disturbance response stripping module is used to take the crop rotation background baseline sequence as input, set the time point corresponding to the milkvetch over-burden operation as the disturbance stripping anchor point, construct an aligned time window before and after the over-burden operation, segment and strip the ecological observation record, and generate an over-burden disturbance response sequence. The net effect analysis module is used to perform parallel dimensionality stripping on ecological observation records reflecting gas absorption and gas emission behavior under the unified time alignment condition of the overturning disturbance response sequence, so as to obtain a gas regulation net effect sequence reflecting the synchronous changes of positive and negative ecological effects. The ecological value assessment module is used to perform time mapping between the net gas regulation effect sequence and the crop rotation background baseline sequence, and to reconstruct the changes in ecological effects in each time period based on the crop rotation continuity status, and output the dynamic assessment results of the evolution of the ecological value of milkvetch over time.

[0014] Thirdly, the present invention provides a computer device, including a memory and a processor, wherein the memory stores a computer program, wherein: when the computer program instructions are executed by the processor, they implement the steps of a dynamic assessment method for the ecological value of *Astragalus membranaceus* based on a multidimensional stripping algorithm as described in the first aspect of the present invention.

[0015] Fourthly, the present invention provides a computer-readable storage medium having a computer program stored thereon, wherein: when the computer program instructions are executed by a processor, they implement the steps of a dynamic assessment method for the ecological value of *Astragalus membranaceus* based on a multidimensional stripping algorithm as described in the first aspect of the present invention.

[0016] The beneficial effects of this invention are as follows: By constructing a crop rotation background benchmark sequence and distinguishing between continuous and discontinuous crop rotation states, this invention effectively weakens the background interference of different crop rotation conditions on the ecological assessment results, thereby improving the comparability and stability of ecological value assessment; by using the incorporation of milkvetch as a disturbance stripping anchor point and constructing an aligned time window before and after incorporation, it achieves a fine decoupling of the disturbance effect of human agricultural operations and the long-term ecological change trend, enhancing the accuracy of identifying the ecological effect of incorporation; by performing parallel dimensional stripping and unified time alignment of gas absorption behavior and gas emission behavior, it synthesizes positive and negative ecological effects on the same time scale to form a gas regulation net effect sequence, avoiding the ecological value bias caused by single-index analysis.

[0017] By further combining the continuous state of crop rotation to constrain and reconstruct the ecological effects, the assessment results are made to have confidence indicators and dynamic evolution characteristics, thereby realizing a continuous, objective and highly credible assessment of the ecological value of milkvetch over time, which significantly improves the ability of ecological management and decision support. Attached Figure Description

[0018] To more clearly illustrate the technical solutions of the embodiments of the present invention, the accompanying drawings used in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0019] Figure 1 This is a flowchart of a dynamic assessment method for the ecological value of milkvetch based on a multidimensional stripping algorithm.

[0020] Figure 2 This is a structural diagram of the dynamic assessment system for the ecological value of milkvetch based on a multidimensional stripping algorithm. Detailed Implementation

[0021] To make the above-mentioned objects, features and advantages of the present invention more apparent and understandable, the specific embodiments of the present invention will be described in detail below with reference to the accompanying drawings.

[0022] Many specific details are set forth in the following description in order to provide a full understanding of the invention. However, the invention may also be practiced in other ways different from those described herein, and those skilled in the art can make similar extensions without departing from the spirit of the invention. Therefore, the invention is not limited to the specific embodiments disclosed below.

[0023] Secondly, the term "one embodiment" or "embodiment" as used herein refers to a specific feature, structure, or characteristic that may be included in at least one implementation of the present invention. The phrase "in one embodiment" appearing in different places in this specification does not necessarily refer to the same embodiment, nor is it a single or selective embodiment that is mutually exclusive with other embodiments.

[0024] Figure 1 This is a flowchart illustrating a dynamic assessment method for the ecological value of milkvetch based on a multidimensional stripping algorithm, according to an embodiment of the present invention. Figure 1 As shown, a dynamic assessment method for the ecological value of milkvetch based on a multidimensional stripping algorithm includes: S1: Under the condition of purple clover rotation in the rice-growing areas of southern China, ecological observation records were collected for multiple consecutive rotation cycles, and the ecological observation records were classified according to whether the rotation was continuous to form a rotation background baseline sequence.

[0025] It should be noted that, due to the significant seasonal, plot-specific, and management differences in crop rotation patterns of milkvetch in southern rice-growing areas, it is difficult to distinguish between natural succession, the crop rotation system itself, and the ecological effects introduced by tilling operations based solely on ecological observation data within a short period or a single rotation cycle. In existing technologies, most ecological assessment methods typically use annual or quarterly statistical data directly as the background baseline, or simply smooth the ecological observation data using a simple time averaging method, without explicitly distinguishing the key agricultural system variable of whether crop rotation is continuous. This approach easily introduces background aliasing in the milkvetch rotation scenario, causing subsequent identification results of ecological effect changes to be interfered with by historical missing rotations or abnormal rotation cycles. Therefore, this invention proposes a background construction mechanism with crop rotation continuity as the core constraint. Through the time processing capabilities of edge computing nodes and crop rotation continuity discrimination rules, a crop rotation background baseline sequence can be directly used for multidimensional stripping algorithms. Specifically: S1.1: Based on the actual crop rotation arrangement of milkvetch in the southern rice-growing areas, edge computing nodes deployed in the fields collect ecological observation sensor data of the same plot in real time and organize it with local timestamps, and upload it to the cloud or central server, arrange it according to the crop rotation cycle, and construct the original observation sequence.

[0026] In existing agricultural ecological monitoring systems, ecological observation data is typically collected by multiple types of sensors and periodically uploaded to a central server, where it is then uniformly processed and sorted by time. This centralized processing method has two prominent problems in practical applications: First, the large number of field sensors and high sampling frequency result in a large volume of uploaded data and inconsistent timestamps, which can easily lead to time drift or data loss, especially in the unstable network environment of southern rice-growing areas. Second, time processing relies entirely on central processing, making it difficult to refine the organization based on the crop rotation cycle boundaries of specific plots, thus weakening the accuracy of subsequent analysis based on crop rotation cycles. To address these issues, this invention deploys edge computing nodes in the field as pre-processing units for ecological observation data.

[0027] The edge computing nodes maintain communication connections with multi-source ecological observation sensors deployed within the same plot, receiving ecological observation sensor data including but not limited to soil physicochemical indicators, vegetation cover status, and gas flux changes. The set of multi-source ecological observation records, indexed primarily by time and with time-based time processing and cycle attribution completed by the edge computing nodes within the same plot and crop rotation cycle, is recorded as the original observation sequence.

[0028] Edge computing nodes consist of a microcontroller, sensor interface, local storage unit, and power management module, and are deployed in field weather stations or data acquisition boxes.

[0029] S1.2: In the original observation sequence, based on whether there is a lack of milkvetch planting or an interruption in the planting sequence between adjacent crop rotation cycles, the original observation sequence is judged to determine the continuity of crop rotation, forming a continuous crop rotation observation sequence and a non-continuous crop rotation observation sequence.

[0030] Specifically, in the original observation sequence, the crop rotation cycle is used as the basic discrimination unit. The ecological observation records corresponding to complete milkvetch planting are checked sequentially within adjacent crop rotation cycles according to time sequence. For example, based on a typical farming calendar in southern rice-growing areas, a theoretical time range for milkvetch growth crop rotation cycles is defined (e.g., October to April of the following year). Using this time range as a template, multiple crop rotation cycle windows to be discriminated are divided by sliding along the time axis. Key ecological observation indicators that reflect vegetation cover, such as the Normalized Difference Vegetation Index (NDVI), are selected. An NDVI threshold representing effective milkvetch growth is set, for example, 0.3, and a minimum duration threshold representing growth continuity is set, for example, 30 days. For each crop rotation cycle window to be discriminated, the NDVI time-series data within that window is checked: if there is a continuous time period (length ≥ 30 days), and most NDVI observation values ​​within that time period are ≥ 0.3, then milkvetch planting is determined to exist in that crop rotation cycle; otherwise, milkvetch planting is determined to be missing in that crop rotation cycle.

[0031] If, between adjacent crop rotation cycles, there are consecutive ecological observation records of milkvetch covering the start and end times of that crop rotation cycle in the original observation sequence, and there are no missing milkvetch plantings or interruptions in the planting sequence, then the corresponding crop rotation cycles are determined to constitute a continuous crop rotation. If, within any crop rotation cycle, the original observation sequence lacks ecological observation records corresponding to milkvetch planting, or there is a significant interruption in the planting sequence, then the corresponding crop rotation cycle is determined to be a non-continuous crop rotation, and the ecological observation records within that crop rotation cycle are classified into the non-continuous crop rotation observation sequence.

[0032] S1.3: Retain the corresponding time identifiers for continuous crop rotation observation sequences and non-continuous crop rotation observation sequences respectively, and align them according to a unified time axis to obtain the crop rotation background benchmark sequence.

[0033] The alignment process according to a unified timeline includes: using the start date of the entire study period (covering all crop rotation cycles) as point 0, generating a continuous, equally spaced integer timeline index (e.g., Day 0, Day 1, ..., Day N) based on time units (e.g., day). Each ecological observation record is mapped to the corresponding point in the aforementioned timeline index according to its observation date; if there is no observation data at a certain time point, but data exists at the time points before and after it, no interpolation is performed. The background baseline sequence allows for missing data points, and continuity is ensured by the "crop rotation cycle" unit.

[0034] S2: Using the crop rotation background reference sequence as input, the time point corresponding to the milkvetch over-pressing operation is set as the disturbance stripping anchor point. An aligned time window is constructed before and after the over-pressing operation, and the ecological observation record is segmented and stripped to generate an over-pressing disturbance response sequence.

[0035] S2.1: For the rotation cycles contained in the continuous and non-continuous crop rotation observation sequences, when the milkvetch over-plowing operation occurs, the field edge computing node captures or receives the event signal, immediately uses the time of occurrence as the disturbance stripping anchor point, and generates indices for the time window before and after over-plowing according to the preset fixed length.

[0036] The fixed duration is 31 days, extending 15 days before and 15 days after the day of the tillage operation. It should be noted that the 2-3 weeks after tillage are a critical period for green manure decomposition, nutrient release, changes in soil microbial activity, and changes in greenhouse gas emission fluxes; the 2 weeks before tillage can be used to establish a stable local background baseline. This duration is sufficient to capture core disturbance signals without being so long as to obscure seasonal trends.

[0037] By generating time window indexes in real time at the edge computing nodes, not only is the repeated scanning of raw time data by the central server reduced, but the time windows before and after the disturbance also have a consistent time scale across different plots and different crop rotation cycles.

[0038] S2.2: The cloud or central server receives the disturbance stripping anchor points and corresponding time window index information uploaded by the edge computing nodes. According to the index of the time window before and after the overburdening, the cloud or central server segments the corresponding ecological observation records in the crop rotation background benchmark sequence by time. The ecological observation records before the overburdening operation and the ecological observation records after the overburdening operation are collected separately to form the overburdening disturbance segment sequence.

[0039] S2.3: Arrange the overturning disturbance segment sequence according to the rotation cycle, and retain the time identifier corresponding to the rotation continuity state to generate an overturning disturbance response sequence for characterizing the impact of overturning operation disturbance.

[0040] S3: Under the condition of unified time alignment of the overturning disturbance response sequence, parallel dimensional stripping is performed on the ecological observation records reflecting gas absorption behavior and gas emission behavior to obtain the gas regulation net effect sequence reflecting the synchronous changes of positive and negative ecological effects.

[0041] S3.1: Under the condition of unified time alignment of the overturning disturbance response sequence, extract the ecological observation records of gas absorption behavior to form a time-continuous gas absorption effect sequence; at the same time, extract the ecological observation records of gas emission behavior to form a time-continuous gas emission effect sequence.

[0042] Among them, gas absorption behavior includes, but is not limited to, the net absorption of carbon dioxide during the growth and soil fixation of milkvetch; gas emission behavior includes, but is not limited to, the net emission of greenhouse gases such as methane and nitrous oxide produced by soil microbial activity after milkvetch is turned over.

[0043] The gas absorption effect sequence and the gas emission effect sequence are time-aligned to ensure that they have the exact same set of time points on the time axis. The specific alignment scheme is as follows: Based on the total time range covered by the aforementioned overpressure disturbance response sequence, a series of equally spaced common time points, termed the analysis time grid, is generated according to the highest time resolution required for the analysis (e.g., 1 hour or 1 day). The original data from the gas absorption effect sequence and the gas emission effect sequence are resampled onto the aforementioned analysis time grid using interpolation methods. If the original sequence has data both before and after a certain grid point, the effect value at that point is estimated using linear interpolation. Linear interpolation can effectively reflect trend changes within a short time window and has low computational complexity, making it suitable for rapid processing at the edge or in the cloud. If the grid point is outside the time range of the original data, the effect value at that point is marked as invalid or marginalized.

[0044] If a grid point has no valid data within a certain time range (such as two consecutive sampling periods), then that point will not be calculated or will be marked as missing data.

[0045] After the above processing, two sequences are obtained that are strictly aligned on the analysis time grid. It should be noted that this invention does not exclude the use of other interpolation methods (such as spline interpolation) or the implementation of an alignment method with equal periodic sampling through hardware synchronization.

[0046] S3.2: For the gas absorption effect sequence and gas emission effect sequence, based on their respective corresponding global warming potential equivalent values ​​of greenhouse gases, the effect intensity is standardized so that the ecological effects produced by different gas behaviors are comparable and additive in terms of intensity; this standardization process needs to be evaluated under the same time background.

[0047] Specifically, for each time record in the gas absorption effect sequence, the corresponding warming potential equivalent value is selected according to the corresponding gas type and used as the effect weight parameter to reflect the relative impact intensity of different gas emission or absorption behaviors at the ecological regulation level; the absorption effect quantity corresponding to the time record is weighted and converted to represent the positive ecological effect intensity under a unified effect scale.

[0048] In one optional implementation, the effect size at the corresponding time point is proportionally converted based on the effect weights, so that the effect size is adjusted proportionally as the weights change. It should be noted that the weighted conversion in this invention is not limited to a specific mathematical operation form; its core lies in uniformly mapping the intensity of the ecological effects of different gas behaviors through effect weights, and multiplication is only one possible implementation method.

[0049] Applying the same treatment to the gas emission effect sequence represents the intensity of negative ecological effects at a uniform effect scale.

[0050] S3.3: The standardized gas absorption effect sequence is labeled as a positive value sequence, and the gas emission effect sequence is labeled as a negative value sequence, while maintaining strict synchronization between the two sequences on the time axis.

[0051] S3.4: Under the unified time alignment condition, using the gas absorption effect sequence marked with positive values ​​and the gas emission effect sequence marked with negative values ​​as parallel inputs, at each same time point, perform the algebraic sum calculation of the positive and negative sequences to generate the initial gas regulation effect synthesis sequence.

[0052] The initial gas regulation effect synthesis sequence simultaneously reflects the combined regulatory results of two types of ecological behaviors—gas absorption and gas emission—at each time point, thus achieving a synchronous characterization of positive and negative ecological effects. Each sequence element corresponds to an instantaneous or stage-specific net ecological effect under a single overturning disturbance.

[0053] S3.5: Using the rotation continuity state corresponding to the overturning disturbance response sequence as a constraint, the initial gas regulation effect synthesis sequence is smoothed to eliminate short-term drastic fluctuations caused by discontinuous rotation or observational anomalies, forming a stable and continuous gas regulation net effect sequence.

[0054] The constraint of the rotation continuity state corresponding to the overturning disturbance response sequence means that when processing the initial gas regulation effect synthesis sequence, the continuous rotation time period and the non-continuous rotation time period are distinguished according to the rotation state corresponding to each time period, and different sequence processing strategies are adopted for the two types of time periods.

[0055] In this embodiment of the invention, the smoothing process employs a moving average method. Specifically, for data segments with a continuous state of 'continuous', a weaker smoothing constraint is applied. For example, a smaller moving time window (e.g., 3 or 5 days) is used. This effectively filters out minor noise while preserving as much detail as possible the true and rapid ecological response caused by overburden disturbance. For data segments with a discontinuous state of 'continuous', a stronger smoothing constraint is applied. For example, a larger moving time window (e.g., 7 or 10 days, or proportional to a 'fixed-length time window,' such as taking 1 / 4 of it) is used. This more effectively suppresses larger background disturbances and uncertainties caused by factors such as inconsistent planting history and large differences in soil background conditions, highlighting the general trend across cycles rather than individual fluctuations.

[0056] In practice, after generating the initial gas regulation effect synthesis sequence, the data points are first marked according to the continuity of the crop rotation cycle to which the overburden event belongs. Then, the moving average algorithm with different window sizes is applied to smooth the data segments that belong to "continuous" and "discontinuous".

[0057] S3.6: Re-associate the time markers of the gas regulation net effect sequence with the disturbance stripping anchor point and the corresponding crop rotation cycle to form a time evolution sequence characterizing the net ecological effect that can be traced back to the specific overburden operation event.

[0058] Specifically, after forming a stable sequence of net gas regulation effects, the effects at each time point in the sequence are mapped and bound to the corresponding overburden disturbance stripping anchor point and the corresponding crop rotation cycle, so that each change in net ecological effect can be traced back to a specific overburden operation event. After re-correlation, it can be clearly seen which overburden operation caused the change in net effect.

[0059] S4: The net effect sequence of gas regulation and the baseline sequence of crop rotation are time-mapped, and the changes in ecological effects in each time period are reconstructed under constraints according to the continuity of crop rotation, so as to output the dynamic evaluation results of the ecological value of milkvetch evolving over time.

[0060] S4.1: Using the time axis of the crop rotation background reference sequence as a unified reference, map the time markers of each data point in the gas regulation net effect sequence to the corresponding disturbance time window position within the source crop rotation cycle.

[0061] Within each crop rotation cycle, the mapped time range is jointly defined by the time window before and after the overburdening, and the mapped net gas regulation effect value is bound to the continuity state of the crop rotation cycle.

[0062] S4.2: Based on the continuity state bound to each rotation cycle after mapping, apply reconstruction constraints to the net gas regulation effect value within the corresponding time period.

[0063] It should be noted that not all mapped net gas regulation effects are treated equally; instead, differentiated reconstruction rules are applied based on the continuity of crop rotation. The specific operation is as follows: For the time period corresponding to the rotation cycle with a continuous status of "continuous", the net effect value of gas regulation is directly used as the effective estimate of the change in ecological value during the time period and marked as a high-confidence data segment.

[0064] For the net effect value of gas regulation within the time period corresponding to the rotation cycle with a continuous state of "discontinuous", it is necessary to multiply it by the reconstruction weighting coefficient to characterize the estimation uncertainty introduced by the discontinuity of ecological processes caused by planting interruption, generate the corrected ecological value change estimate, and mark it as a low confidence data segment.

[0065] The reconstruction weighting coefficient is a numerical attenuation factor between 0 and 1, used to reduce the net gas regulation effect value calculated for discontinuous crop rotation cycles. For continuous crop rotation cycles, the reconstruction weighting coefficient is 1, meaning the net effect value is fully adopted with high confidence. For discontinuous crop rotation cycles, the reconstruction weighting coefficient is less than 1, meaning the weighting coefficient should reflect the degree of disruption to the continuity of ecological processes caused by the "interruption." For example, if other crops with significantly different ecological effects (such as vegetables) were planted during the interruption period, a lower reconstruction weighting coefficient (e.g., 0.6) is used; if it is natural fallow, a moderate value (e.g., 0.8) is used.

[0066] S4.3: The estimated ecological value changes within all crop rotation cycles, after constraint reconstruction, are arranged and connected according to the mapped time sequence to form a complete time series, which serves as the output of the dynamic assessment sequence for the ecological value of milkvetch. Specifically, using the crop rotation cycle sequence as the main sequence index, the estimated ecological value changes within the pre-plowing and post-plowing time windows of each crop rotation cycle are sequentially concatenated to form a continuous time series; during the concatenation process, the confidence level information corresponding to each time period is simultaneously retained.

[0067] The output of the dynamic evaluation results includes identifiers for high-confidence and low-confidence data segments to reflect the differences in the reliability of the evaluation results at different time periods.

[0068] Furthermore, such as Figure 2 As shown, this embodiment also provides a dynamic assessment system for the ecological value of milkvetch based on a multidimensional stripping algorithm, including: The crop rotation baseline construction module is used to collect ecological observation records within multiple consecutive crop rotation cycles under the condition of purple clover rotation in the rice-growing areas of southern China, and classify the ecological observation records according to whether the crop rotation is continuous, forming a crop rotation background baseline sequence. The disturbance response stripping module is used to take the crop rotation background baseline sequence as input, set the time point corresponding to the milkvetch over-burden operation as the disturbance stripping anchor point, construct an aligned time window before and after the over-burden operation, segment and strip the ecological observation record, and generate an over-burden disturbance response sequence. The net effect analysis module is used to perform parallel dimensionality stripping on ecological observation records reflecting gas absorption and gas emission behavior under the unified time alignment condition of the overturning disturbance response sequence, so as to obtain a gas regulation net effect sequence reflecting the synchronous changes of positive and negative ecological effects. The ecological value assessment module is used to perform time mapping between the net gas regulation effect sequence and the crop rotation background baseline sequence, and to reconstruct the changes in ecological effects in each time period based on the crop rotation continuity status, and output the dynamic assessment results of the evolution of the ecological value of milkvetch over time.

[0069] This embodiment also provides a computer device applicable to a dynamic assessment method for the ecological value of milkvetch based on a multidimensional stripping algorithm, including a memory and a processor; the memory is used to store computer-executable instructions, and the processor is used to execute the computer-executable instructions to implement the dynamic assessment method for the ecological value of milkvetch based on a multidimensional stripping algorithm as proposed in the above embodiment.

[0070] The computer device can be a terminal, comprising a processor, memory, communication interface, display screen, and input devices connected via a system bus. The processor provides computing and control capabilities. The memory includes non-volatile storage media and internal memory. The non-volatile storage media stores the operating system and computer programs. The internal memory provides an environment for the operation of the operating system and computer programs stored in the non-volatile storage media. The communication interface is used for wired or wireless communication with external terminals; wireless communication can be achieved through Wi-Fi, carrier networks, NFC (Near Field Communication), or other technologies. The display screen can be an LCD screen or an e-ink screen. The input devices can be a touch layer covering the display screen, buttons, a trackball, or a touchpad on the computer device's casing, or an external keyboard, touchpad, or mouse.

[0071] This embodiment also provides a storage medium storing a computer program that, when executed by a processor, implements a dynamic assessment method for the ecological value of milkvetch based on a multidimensional stripping algorithm, as proposed in the above embodiments.

[0072] In summary, the beneficial effects of this invention are as follows: By constructing a crop rotation background benchmark sequence and distinguishing between continuous and discontinuous crop rotation states, this invention effectively weakens the background interference of different crop rotation conditions on the ecological assessment results, improving the comparability and stability of ecological value assessment; by using the incorporation of milkvetch as a disturbance stripping anchor point and constructing an aligned time window before and after incorporation, it achieves a fine decoupling of the disturbance effect of human agricultural operations and the long-term ecological change trend, enhancing the accuracy of identifying the ecological effect of incorporation; by performing parallel dimensional stripping and unified time alignment of gas absorption behavior and gas emission behavior, it synthesizes positive and negative ecological effects on the same time scale to form a gas regulation net effect sequence, avoiding the ecological value bias caused by single-index analysis.

[0073] By further combining the continuous state of crop rotation to constrain and reconstruct the ecological effects, the assessment results are made to have confidence indicators and dynamic evolution characteristics, thereby realizing a continuous, objective and highly credible assessment of the ecological value of milkvetch over time, which significantly improves the ability of ecological management and decision support.

[0074] It should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and are not intended to limit it. Although the present invention has been described in detail with reference to preferred embodiments, those skilled in the art should understand that modifications or equivalent substitutions can be made to the technical solutions of the present invention without departing from the spirit and scope of the technical solutions of the present invention, and all such modifications or substitutions should be covered within the scope of the claims of the present invention.

Claims

1. A dynamic assessment method for the ecological value of milkvetch based on a multidimensional stripping algorithm, characterized in that: include: Under the condition of milkvetch rotation in the rice-growing areas of southern China, ecological observation records were collected over multiple consecutive rotation cycles. The ecological observation records were then classified according to whether the rotation was continuous, forming a baseline sequence for the rotation background. Using the aforementioned crop rotation background baseline sequence as input, the time point corresponding to the milkvetch over-plowing operation is set as the disturbance stripping anchor point. An aligned time window is constructed before and after the over-plowing operation, and the ecological observation records are segmented and stripped to generate an over-plowing disturbance response sequence. Under the condition of unified time alignment of the overturning disturbance response sequence, parallel dimensional stripping is performed on the ecological observation records reflecting gas absorption and gas emission behavior to obtain the gas regulation net effect sequence reflecting the synchronous changes of positive and negative ecological effects. The net effect sequence of gas regulation and the baseline sequence of crop rotation are time-mapped, and the changes in ecological effects in each time period are reconstructed under constraints based on the continuity of crop rotation, so as to output the dynamic evaluation results of the ecological value of milkvetch evolving over time.

2. The method for dynamic assessment of the ecological value of milkvetch based on a multidimensional stripping algorithm as described in claim 1, characterized in that: The formation of the crop rotation background benchmark sequence includes: According to the actual crop rotation arrangement of milkvetch in the southern rice-growing areas, edge computing nodes deployed in the fields collect ecological observation sensor data of the same plot in real time and organize it with local timestamps, and upload it to the cloud or central server, arrange it according to the crop rotation cycle, and construct the original observation sequence; In the original observation sequence, based on whether there is a lack of milkvetch planting or an interruption in the planting sequence between adjacent crop rotation cycles, the original observation sequence is judged to determine the continuity of crop rotation, forming a continuous crop rotation observation sequence and a non-continuous crop rotation observation sequence. The corresponding time markers of continuous and non-continuous crop rotation observation sequences are retained and aligned according to a unified time axis to obtain the crop rotation background benchmark sequence.

3. The method for dynamic assessment of the ecological value of milkvetch based on a multidimensional stripping algorithm as described in claim 2, characterized in that: The setting of the disturbance stripping anchor point includes: For the rotation cycles contained in the continuous and non-continuous crop rotation observation sequences, when the milkvetch over-plowing operation occurs, the field edge computing node captures or receives the event signal, immediately uses the time of occurrence as the disturbance stripping anchor point, and generates indices for the time window before and after over-plowing according to a preset fixed length.

4. The method for dynamic assessment of the ecological value of milkvetch based on a multidimensional stripping algorithm as described in claim 3, characterized in that: The generated overturning perturbation response sequence includes: According to the indexes of the time window before and after the overburdening operation, the cloud or central server segments the corresponding ecological observation records in the crop rotation background baseline sequence by time, and collects the ecological observation records before and after the overburdening operation to form the overburdening disturbance segment sequence. The overturning disturbance segment sequence is arranged in the order of the rotation cycle, and the time identifier corresponding to the rotation continuity state is retained to generate an overturning disturbance response sequence that characterizes the impact of overturning operation disturbance.

5. The method for dynamic assessment of the ecological value of milkvetch based on a multidimensional stripping algorithm as described in claim 4, characterized in that: The parallel dimensional stripping of ecological observation records reflecting gas absorption and emission behaviors includes: Under the condition of uniform time alignment of the overturning disturbance response sequence, ecological observation records of gas absorption behavior are extracted to form a time-continuous gas absorption effect sequence; at the same time, ecological observation records of gas emission behavior are extracted to form a time-continuous gas emission effect sequence. The gas absorption effect sequence and the gas emission effect sequence are standardized based on their respective global warming potential equivalent values ​​of greenhouse gases. The standardized gas absorption effect sequence was labeled as a positive value sequence, and the gas emission effect sequence was labeled as a negative value sequence, while maintaining strict synchronization between the two sequences on the time axis.

6. The method for dynamic assessment of the ecological value of milkvetch based on a multidimensional stripping algorithm as described in claim 5, characterized in that: The parallel dimensionality stripping of ecological observation records reflecting gas absorption and emission behaviors also includes: Under the unified time alignment condition, the gas absorption effect sequence marked with positive values ​​and the gas emission effect sequence marked with negative values ​​are used as parallel inputs. At each same time point, the algebraic sum calculation of the positive and negative sequences is performed to generate the initial gas regulation effect synthesis sequence. Using the rotation continuity state corresponding to the overturning disturbance response sequence as a constraint, the initial gas regulation effect synthesis sequence is smoothed to generate the net gas regulation effect sequence. The time markers of the gas regulation net effect sequence are re-associated with the disturbance stripping anchor point and the corresponding rotation cycle to form a time evolution sequence.

7. The method for dynamic assessment of the ecological value of milkvetch based on a multidimensional stripping algorithm as described in claim 6, characterized in that: The dynamic assessment results of the ecological value of milkvetch evolving over time include: Using the time axis of the crop rotation background reference sequence as a unified reference, the time markers of each data point in the gas regulation net effect sequence are mapped to the corresponding disturbance time window positions within the source crop rotation cycle. Within each crop rotation cycle, the mapped time range is jointly defined by the time window before and after the overburdening, and the mapped net gas regulation effect value is bound to the continuity state of the crop rotation cycle. Based on the continuous state bound to each rotation cycle after mapping, a reconstruction constraint is applied to the net gas regulation effect value within the corresponding time period: For the rotation cycle corresponding to the continuous crop rotation observation sequence, the net effect value of gas regulation is directly used as an effective estimate of the ecological value change in the corresponding time period and marked as a high-confidence data segment; For the rotation cycle corresponding to the non-continuous crop rotation observation sequence, the net effect value of gas regulation needs to be multiplied by the reconstruction weighting coefficient to generate the corrected estimate of the change in ecological value, and marked as a low-confidence data segment; The estimated values ​​of ecological value changes in all crop rotation cycles, after being reconstructed under constraints, are arranged and connected in the mapped time sequence to form a complete time series, which is then output as the dynamic assessment sequence of the ecological value of milkvetch.

8. A dynamic ecological value assessment system for milkvetch based on a multidimensional stripping algorithm, based on the dynamic ecological value assessment method for milkvetch based on a multidimensional stripping algorithm as described in any one of claims 1 to 7, characterized in that: Also includes: The crop rotation baseline construction module is used to collect ecological observation records within multiple consecutive crop rotation cycles under the condition of purple clover rotation in the rice-growing areas of southern China, and classify the ecological observation records according to whether the crop rotation is continuous, forming a crop rotation background baseline sequence. The disturbance response stripping module is used to take the crop rotation background baseline sequence as input, set the time point corresponding to the milkvetch over-burden operation as the disturbance stripping anchor point, construct an aligned time window before and after the over-burden operation, segment and strip the ecological observation record, and generate an over-burden disturbance response sequence. The net effect analysis module is used to perform parallel dimensionality stripping on ecological observation records reflecting gas absorption and gas emission behavior under the unified time alignment condition of the overturning disturbance response sequence, so as to obtain a gas regulation net effect sequence reflecting the synchronous changes of positive and negative ecological effects. The ecological value assessment module is used to perform time mapping between the net gas regulation effect sequence and the crop rotation background baseline sequence, and to reconstruct the changes in ecological effects in each time period based on the crop rotation continuity status, and output the dynamic assessment results of the evolution of the ecological value of milkvetch over time.

9. A computer device comprising a memory and a processor, wherein the memory stores a computer program, characterized in that: When the processor executes the computer program, it implements the steps of the dynamic assessment method for the ecological value of milkvetch based on a multidimensional stripping algorithm as described in any one of claims 1 to 7.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by the processor, it implements the steps of the dynamic assessment method for the ecological value of milkvetch based on a multidimensional stripping algorithm as described in any one of claims 1 to 7.

Images