Irrigation disaster monitoring method and system based on remote sensing monitoring

Abstract

The invention belongs to the technical field of irrigation disaster monitoring, and provides an irrigation disaster monitoring method and system based on remote sensing monitoring, and the method comprises the steps: carrying out the quantitative analysis of the growth data difference between an irrigation region and a surrounding region through a neighborhood matrix method, determining a suspected drought disaster region, and constructing a drought evaluation model of the suspected drought disaster region; determining whether the suspected drought disaster area is a drought disaster area; the irrigation area disaster monitoring efficiency is improved, the irrigation area disaster monitoring resource cost is saved, the drought time sequence is processed and analyzed by combining the self-correlation function and the partial self-correlation function, whether the drought occurrence of the drought disaster area has periodic changes or not is judged, the monitoring system of the drought disaster area is updated according to the periodic changes, and the drought disaster monitoring efficiency is improved. The irrigation area drought disaster monitoring efficiency is improved, the prediction model is constructed according to the analysis results of the autocorrelation function and the partial autocorrelation function, irrigation area drought prediction is achieved, irrigation area disaster comprehensive monitoring is achieved, and monitoring comprehensiveness is improved.

Description

Technical Field

[0001] The invention belongs to the technical field of irrigation disaster monitoring, and in particular relates to an irrigation disaster monitoring method and system based on remote sensing monitoring. Background Art

[0002] Disaster monitoring in irrigation areas can monitor and warn of natural disasters such as drought, floods, low temperature and frost in real time, provide timely disaster information for agricultural production, and effectively reduce the impact of disasters on agricultural production and ensure the safety and stability of agricultural production. However, there are still shortcomings in disaster monitoring in irrigation areas, such as low monitoring efficiency, inaccurate monitoring, and incomplete monitoring. Therefore, it is of great significance to study an irrigation disaster monitoring method and system based on remote sensing monitoring.

[0003] In the prior art, due to the limitations of monitoring means and the inadequacy of data processing methods, the monitoring results of drought disasters in irrigation areas are often inaccurate. For example, there is a lack of monitoring through remote sensing technology to obtain growth data of irrigation areas and adjacent areas and to conduct growth differentiation analysis, resulting in inaccurate monitoring results of drought disasters in irrigation areas. In the prior art, there is also a lack of combined monitoring through satellite remote sensing and unmanned aerial vehicle remote sensing technology. For example, there is a lack of preliminary disaster identification of irrigation areas through satellite remote sensing monitoring and analysis, and then the use of unmanned aerial vehicle remote sensing monitoring and analysis to finally identify disasters in irrigation areas and adjust the monitoring system of irrigation areas according to the identification and analysis results, resulting in low monitoring efficiency and waste of monitoring resource costs. In the prior art, there is also a lack of a comprehensive monitoring system. For example, while monitoring drought disasters in irrigation areas, a disaster prediction model is constructed according to the monitoring and analysis results.

[0004] To this end, the present invention provides an irrigation disaster monitoring method and system based on remote sensing monitoring. Summary of the invention

[0005] In order to make up for the deficiencies of the prior art, at least one technical problem raised in the background technology is solved.

[0006] The technical solution adopted by the present invention to solve its technical problem is: An irrigation disaster monitoring method based on remote sensing monitoring, comprising: S1: Identify the irrigation analysis area through preliminary comparison of the growth data of the irrigation area, introduce the neighborhood matrix method to perform differential analysis on the growth data of the irrigation analysis area, and identify the suspected drought disaster area; S2: Based on the historical weak growth data, historical temperature data and historical soil moisture data of the suspected drought-stricken areas, and using the entropy weight method to add data weights, a drought assessment model for the suspected drought-stricken areas is constructed; and based on the output results of the drought assessment model, it is determined whether the suspected drought-stricken areas are drought-stricken areas. If so, enter S3; S3: Screen and determine the drought results from multiple drought assessment results in the drought-stricken areas, extract the drought time series, and perform drought time series processing and analysis in combination with the autocorrelation function and the partial autocorrelation function to determine whether the drought in the drought-stricken areas has periodic changes. If so, enter S4; S4: Update the monitoring system of drought-affected areas according to periodic changes and build a drought prediction model for drought-affected areas.

[0007] As a further technical solution of the present invention, if at least one growth data in the growth data is lower than the minimum value of the normal range, the irrigation area is marked as an irrigation analysis area; Determine the total number of irrigation areas and the weights between irrigation areas, construct a neighborhood matrix, find the weights between the irrigation analysis area and other neighboring irrigation areas based on the neighborhood matrix, and calculate the growth data difference in combination with the growth data, where the growth data difference includes the normalized vegetation index difference, the vegetation cover difference, and the leaf area index difference; If any growth data difference is positive, the irrigation analysis area will be marked as a suspected drought disaster area.

[0008] As a further technical solution of the present invention, the weights between the irrigation areas are calculated by combining the center point distance between the irrigation areas with the Gaussian function.

[0009] As a further technical solution of the present invention, the historical weak growth data, historical temperature data and historical soil moisture data of the suspected drought-affected area during multiple historical monitoring periods are obtained and normalized. The entropy weights of historical vulnerable growth data, historical temperature data, and historical soil moisture data are calculated as weight coefficients to construct a drought assessment model, and the drought index is output through the drought assessment model; If the drought index is greater than the drought index threshold, the suspected drought-hit area will be marked as a drought-hit area.

[0010] As a further technical solution of the present invention, drought assessment results with drought index greater than drought index threshold are marked as confirmed drought results, and all time points corresponding to confirmed drought results are extracted and integrated into drought time series.

[0011] As a further technical solution of the present invention, the ACF value and the PACF value of the drought time series are calculated, and the ACF graph and the PACF graph are drawn; In the ACF diagram, the point where the autocorrelation coefficient of the lag period exceeds the set threshold is marked as a significant autocorrelation peak; In the PACF diagram, the point where the partial autocorrelation coefficient of the lag period exceeds the set threshold is marked as a significant peak of partial autocorrelation; The significant peaks of autocorrelation and partial autocorrelation are processed and analyzed to obtain the periodicity value; If the periodicity value is greater than the periodicity threshold, it means that the drought in the drought-affected area has periodic changes.

[0012] As a further technical solution of the present invention, the following is provided: obtaining the number ratio of the significant autocorrelation peaks and the significant partial autocorrelation peaks; Obtain the time interval between every two adjacent autocorrelation significant peaks, and integrate them to obtain an autocorrelation interval group; Get the time interval between every two adjacent significant partial autocorrelation peaks, and integrate them to obtain the partial autocorrelation interval group; The variance and mean of the autocorrelation interval group and the partial autocorrelation interval group are calculated respectively, and the periodicity value is obtained by calculating the number proportion, variance and mean deviation of the significant autocorrelation peaks and the significant partial autocorrelation peaks.

[0013] As a further technical solution of the present invention, proportional weights are added to the autocorrelation interval group and the partial autocorrelation interval group respectively through the variance values ​​of the autocorrelation interval group and the partial autocorrelation interval group, and the output cycle update duration is calculated according to the proportional weights and the mean values ​​of the autocorrelation interval group and the partial autocorrelation interval group; The monitoring system for drought-stricken areas will be updated according to the periodic update duration.

[0014] As a further technical solution of the present invention, the stationarity test is performed on the drought time series; If the drought time series is unstable, the drought time series is differentiated once or multiple times until a stable series is obtained. The order of the difference is the difference order d; According to the truncation characteristics or attenuation trends of the ACF and PACF graphs, the values ​​of the autoregressive order p and the moving average order q are determined; The determined parameters p, d, and q are used in combination with the least squares estimation algorithm to fit the ARIMA model and obtain the drought prediction model for drought-stricken areas.

[0015] An irrigation disaster monitoring system based on remote sensing monitoring, comprising: Disaster area initial identification module: Identify the irrigation analysis area through preliminary comparison of the growth data of the irrigation area, introduce the neighborhood matrix method to perform differential analysis on the growth data of the irrigation analysis area, and identify suspected drought disaster areas; Disaster area determination module: Based on the historical weak growth data, historical temperature data and historical soil moisture data of the suspected drought disaster area, and using the entropy weight method to add data weights, a drought assessment model for the suspected drought disaster area is constructed; and based on the output results of the drought assessment model, it is determined whether the suspected drought disaster area is a drought disaster area; Disaster periodicity analysis module: Screen and determine drought results from multiple drought assessment results in drought-stricken areas, extract drought time series, and process and analyze drought time series in combination with autocorrelation function and partial autocorrelation function to determine whether drought in drought-stricken areas has periodic changes; Disaster monitoring adjustment and prediction module: Update the monitoring system of drought-stricken areas according to periodic changes, and build a drought prediction model for drought-stricken areas.

[0016] The beneficial effects of the present invention are as follows: 1. Through preliminary comparison of growth data, the irrigation analysis area is identified, and the neighborhood matrix method is introduced to perform differential analysis on the growth data of the irrigation analysis area and the adjacent irrigation areas around the irrigation analysis area, and the suspected drought disaster area is identified according to the analysis results; the neighborhood matrix method is used to quantitatively analyze the growth data differences between the irrigation area and the surrounding area, which can more keenly capture the characteristics of drought conditions and improve the accuracy of drought disaster monitoring in irrigation areas. The historical weak growth data, historical temperature data and historical soil moisture data of the suspected drought disaster area when drought occurred in multiple historical monitoring cycles are obtained, and the entropy weight method is used to add data weights, so as to construct a drought assessment model for the suspected drought disaster area; through the output results of the drought assessment model, it is determined whether the suspected drought disaster area is a drought disaster area; using The entropy weight method dynamically calculates the weights of data indicators according to the information content of historical data itself, so that the importance of parameters such as temperature and soil moisture is adjusted with the actual data distribution, which enhances the model's adaptive ability to drought characteristics in different irrigation areas, thereby further improving the accuracy of drought disaster monitoring in irrigation areas. Remote sensing satellite monitoring was first used to analyze the growth data of the irrigation analysis area and the surrounding adjacent irrigation areas to preliminarily analyze and judge whether there is a suspected drought in the irrigation area. In the event of suspected drought, unmanned aerial vehicle remote sensing monitoring was used to construct a drought assessment model by integrating parameters such as temperature and soil moisture, and the drought situation in the irrigation area was finally determined, which improved the efficiency of disaster monitoring in irrigation areas and saved the monitoring resource cost of disasters in irrigation areas.

[0017] 2. The drought results are screened and determined from multiple drought assessment results in the drought-stricken area, and the drought time series is extracted according to the time point of determining the drought results. The drought time series is processed and analyzed in combination with the autocorrelation function (ACF) and the partial autocorrelation function (PACF) to determine whether the drought in the drought-stricken area has periodic changes. If the drought in the drought-stricken area has periodic changes, the monitoring system of the drought-stricken area is updated according to the periodic changes, and the drought time series is processed by the autocorrelation function (ACF) and the partial autocorrelation function (PACF) to construct a drought prediction model for the drought-stricken area. The present invention adjusts the monitoring system according to the periodic changes in the drought-stricken area, thereby improving the monitoring efficiency of drought disasters in irrigation areas, and in the periodic change analysis, a prediction model is constructed according to the analysis results of the autocorrelation function (ACF) and the partial autocorrelation function (PACF), thereby realizing the prediction of drought in irrigation areas, realizing the comprehensive monitoring of disasters in irrigation areas, and improving the comprehensiveness of monitoring. BRIEF DESCRIPTION OF THE DRAWINGS

[0018] The present invention will be further described below in conjunction with the accompanying drawings.

[0019] Figure 1 is a flowchart of the steps of an irrigation disaster monitoring method based on remote sensing monitoring according to an embodiment of the present invention; Figure 2 It is a flowchart of an irrigation disaster monitoring system based on remote sensing monitoring according to an embodiment of the present invention. DETAILED DESCRIPTION

[0020] In order to make the technical means, creative features, objectives and effects achieved by the present invention easy to understand, the present invention is further explained below in conjunction with specific implementation methods. Example 1

[0021] See also Figure 1 As shown, an irrigation disaster monitoring method based on remote sensing monitoring according to an embodiment of the present invention comprises the following steps: S1: Use satellite sensors to monitor each irrigation area, obtain the growth data of the irrigation area, identify the irrigation analysis area through preliminary comparison of the growth data, introduce the neighborhood matrix method to perform differential analysis on the growth data of the irrigation analysis area and the adjacent irrigation areas around the irrigation analysis area, and identify the suspected drought disaster area based on the analysis results; The growth data of the irrigation area in S1 include normalized difference vegetation index (NDVI), vegetation coverage, leaf area index (LAI), etc.; Exemplarily, the process of acquiring the growth data of the irrigation area in S1 includes: During the monitoring period, remote sensing satellites are used to carry multi-spectral sensors to obtain image data in multiple bands, including near-infrared and red bands, which can meet the needs of calculating NDVI. Taking Landsat8 as an example, band 4 of its OLI (Operational LandImager) sensor corresponds to the red band, and band 5 corresponds to the near-infrared band. These band data can be used to calculate the NDVI value according to the formula: NDVI=(NIR−R) / (NIR+R), where NIR is the reflectance of the near-infrared band and R is the reflectance of the red band. During the monitoring period, medium-resolution satellite images such as Landsat and Sentinel-2 were used to estimate vegetation coverage based on the pixel binary model. First, the image was preprocessed, including radiometric calibration and atmospheric correction, to ensure data accuracy. Then, based on the differences in spectral characteristics between vegetation and non-vegetation in the image, auxiliary information such as the Normalized Difference Vegetation Index (NDVI) was used to divide the pixels into vegetation and non-vegetation parts, and then the vegetation coverage was calculated. During the monitoring period, the LAI value is inverted using data from medium-resolution satellites such as MODIS (Moderate-Resolution Imaging Spectroradiometer) combined with radiation transfer models (such as the PROSAIL model). The multi-band data obtained by the MODIS sensor can be used to calculate parameters such as vegetation index. These parameters are used as input to simulate the optical characteristics of the vegetation canopy through the radiation transfer model, and then matched with the actual observed satellite data to invert the LAI value; The process of identifying the irrigation analysis area through preliminary comparison of growth data in S1 is as follows: Based on the growth data of the irrigated area, the growth data are compared with the corresponding normal range; If in the growth data, the growth data are all within the normal range, no operation is performed; If at least one growth data in the growth data is lower than the minimum value of the normal range, the irrigation area is marked as an irrigation analysis area, and the growth data is marked as weak growth data; It should be noted that the NDVI value of normal vegetation in the irrigation area is usually between 0.2 and 0.8, and the LAI of normal vegetation is within a certain range. For example, the LAI of crops is generally between 2 and 6. During the vigorous growth period of crops, such as corn and wheat, the vegetation coverage is generally 60% to 90%. Shortly after sowing or about to be harvested, the coverage is lower, which may be only 10% to 30% in the early stage of sowing. Before harvest, as the leaves turn yellow and fall, the coverage may drop to 40% to 60%. The process of introducing the neighborhood matrix method in S1 to perform differential analysis on the growth data of the irrigation analysis area and the irrigation areas adjacent to the irrigation analysis area includes: Determine the total number of irrigation areas N and the weights between irrigation areas, and construct an N×N neighborhood matrix A:

[0022] Among them, Aij represents the weight between the i-th irrigation area and the j-th irrigation area; The calculation process of the weights between the irrigation areas includes: Determine the center point distance between the irrigation area and other adjacent irrigation areas, and integrate them to obtain a regional distance group (di1, di2, di3...dij), where dij represents the center point distance between the i-th irrigation area and the j-th irrigation area; In actual situations, areas closer to the center have a greater impact on vegetation growth. Therefore, a Gaussian function is used to calculate the weight. The specific calculation formula is as follows:

[0023] in, Represents the attenuation coefficient, the value is 1; For example, assuming that the center distance between irrigation area 1 and irrigation area 2 is d12, the center distance between irrigation area 1 and irrigation area 3 is d13, and d12 < d13, assuming d12 = 1, d13 = 2, then the weights for irrigation area 1 and irrigation area 2 are: , then the weights for irrigation area 1 and irrigation area 3 are ; The method for determining the distance between the center point of the irrigation area and other adjacent irrigation areas is: A1, Data acquisition and processing: Obtain geographical data of the irrigation area through satellite remote sensing images, aerial photogrammetry, etc., and process the data using GIS (Geographic Information System) software, including the unification of coordinate systems, data correction and registration, etc.; A2, determine the center point: In GIS software, for regular-shaped irrigation areas, the center point is determined by calculating the geometric center; for irregular-shaped areas, the centroid algorithm is used to calculate, that is, the center point is determined by calculating the weighted average of the coordinates of all points in the area; A3, Distance measurement: Use the distance measurement tool in the GIS software to directly measure the straight-line distance between the center points of the two irrigation areas; this tool will accurately calculate the actual distance between the two points based on the geographic coordinate system and projection method used; According to the N×N neighborhood matrix A, the weights between the irrigation analysis area and other adjacent irrigation areas are found, and the growth data difference is calculated in combination with the growth data. The growth data difference includes the normalized vegetation index difference, the vegetation coverage difference, and the leaf area index difference. The specific calculation includes: Normalized Difference Vegetation Index Difference: , where NDVIi represents the normalized vegetation index of the i-th irrigation area, and NDVIj represents the normalized vegetation index of the j-th irrigation area; Poor vegetation cover: , where FGi represents the vegetation coverage of the i-th irrigation area, and FGj represents the vegetation coverage of the j-th irrigation area; Leaf Area Index Difference: , where LAi represents the leaf area index of the i-th irrigation area, and LAj represents the leaf area index of the j-th irrigation area; For example, assuming that the irrigation analysis area is irrigation area 1, and the other adjacent irrigation areas are irrigation areas 2, 3, 4, and 5; then the normalized difference in vegetation index is: A12*(NDVI2-NDVI1)+A13*(NDVI3-NDVI1)+A14*(NDVI4-NDVI1)+A15*(NDVI5-NDVI1); It can be understood that the neighborhood matrix method can quantify the growth data of the irrigation analysis area and the surrounding areas, clearly present the differences between the areas in the form of a matrix, fully consider the spatial position relationship between the irrigation analysis area and the surrounding adjacent areas, and more accurately judge the growth status of the analysis area in the overall environment. Since the neighborhood matrix method considers multiple surrounding adjacent irrigation areas at the same time, when multiple surrounding areas have good growth conditions and the analysis area shows obvious growth lag, this comparative difference can more reliably indicate the possibility of drought disasters; Obtain the growth data difference corresponding to the weak growth data. If any growth data difference is positive, mark the irrigation analysis area as a suspected drought disaster area. S2: Obtain the historical weak growth data, historical temperature data and historical soil moisture data of the suspected drought-stricken areas in multiple historical monitoring periods, and use the entropy weight method to add data weights to build a drought assessment model for the suspected drought-stricken areas; use the weak growth data, temperature data and soil moisture data in the current monitoring period as input, and determine whether the suspected drought-stricken areas are drought-stricken areas through the output results of the drought assessment model; The historical weak growth data is the same as the weak growth data type in the current monitoring period; The historical vulnerable growth data, historical temperature data and historical soil moisture data when drought occurs in the multiple historical monitoring periods are obtained through historical monitoring reports of the irrigation area, the historical temperature data includes the soil temperature index (TI), and the historical soil moisture data includes the soil moisture index (SMI); The temperature data and soil moisture data in the current monitoring period are obtained through UAV remote sensing monitoring technology; The process of adding data weights using the entropy weight method includes: The historical weak growth data, historical temperature data, and historical soil moisture data of the suspected drought-affected areas during multiple historical monitoring periods were obtained, and normalized to map the data to the [0,1] interval. The specific normalization formula is:

[0024] Among them, x ab Represents the initial data, is the normalized data, a=1.2.3......n (n represents the number of samples, which is equal to the number of historical monitoring cycles), b=1.2.3......m (m represents the number of data types, i.e. the total number of types of historical weak growth data, historical temperature data and historical soil moisture data), max(x b ) and max(x b ) represents the maximum and minimum values ​​of item b’s data; Calculate the proportion p of the ath sample value under the bth data ab ; The specific formula is:

[0025] Calculate the information entropy e of item b data b , the specific formula is:

[0026] in, ; Calculate the entropy weight w of the b-th data b , the specific formula is:

[0027] The entropy weight of the data is used as the weight coefficient to construct a drought assessment model, which is as follows:

[0028] Among them, DI is the drought index, sz is the historical weak growth data, w1+w2+......+w b =1; This application provides an exemplary calculation description for the process of adding data weights using the entropy weight method and constructing a drought assessment model. For example, assuming that the normalized vegetation index NDVI, temperature index TI, and soil moisture index SMI within 5 historical monitoring cycles are obtained, as shown in Table 1 below; Table 1: Normalized difference vegetation index NDVI, temperature index TI and soil moisture index SMI data during the historical monitoring period; Historical monitoring period NDVI SMI TI 1 0.167 0.167 0.083 2 0.333 0.333 0 3 0 0 0.333 4 0.25 0.25 0.125 5 0.083 0.083 0.25 By formula The information entropy e of item b of data b ; It should be noted that there are 5 historical monitoring cycles. ; The calculation results are: (information entropy of normalized vegetation index) eNVDI = 0.918, (information entropy of soil moisture index) eSMI = 0.918, (information entropy of temperature index) eTI = 0.841; By formula Calculate the entropy weight w of the b-th data b ; The calculation result is: (Entropy weight of normalized difference vegetation index) wNVDI=0.24, (Entropy weight of soil moisture index) wSMI=0.24, (Entropy weight of temperature index) wTI=0.52; The drought assessment model constructed is:

[0029] The weak growth data, temperature data and soil moisture data in the current monitoring period are used as inputs to the drought assessment model, and the drought index DI is obtained as output; In some embodiments, the drought index DI is compared with a drought index threshold; If the drought index DI is greater than the drought index threshold, the suspected drought-hit area will be marked as a drought-hit area; If the drought index DI is less than or equal to the drought index threshold, no operation is performed; The technical solution of the embodiment of the present invention is as follows: by preliminary comparison of growth data, the irrigation analysis area is identified, and the neighborhood matrix method is introduced to perform differential analysis on the growth data of the irrigation analysis area and the irrigation areas adjacent to the irrigation analysis area, and the suspected drought disaster area is identified according to the analysis results; the neighborhood matrix method is used to quantitatively analyze the growth data differences between the irrigation area and the surrounding area, so as to more keenly capture the characteristics of the drought condition, improve the accuracy of drought disaster monitoring in the irrigation area, and obtain the historical weak growth data, historical temperature data and historical soil moisture data of the suspected drought disaster area when drought occurs in multiple historical monitoring cycles, and use the entropy weight method to add data weights, so as to construct a drought assessment model for the suspected drought disaster area; and through the output results of the drought assessment model, determine whether the suspected drought disaster area is drought-prone. disaster area; the entropy weight method is used to dynamically calculate the weights of data indicators according to the information content of the historical data itself, so that the importance of parameters such as temperature and soil moisture is adjusted according to the actual data distribution, which enhances the model's adaptive ability to the drought characteristics of different irrigation areas, thereby further improving the accuracy of drought disaster monitoring in irrigation areas. In addition, remote sensing satellite monitoring was first used to analyze the growth data of the irrigation analysis area and the surrounding adjacent irrigation areas to preliminarily analyze and judge whether there is a suspected drought in the irrigation area. In the event of suspected drought, unmanned aerial vehicle remote sensing monitoring was used to construct a drought assessment model by integrating parameters such as temperature and soil moisture, and the drought situation in the irrigation area was finally determined, which improved the efficiency of disaster monitoring in irrigation areas and saved the cost of monitoring resources for disasters in irrigation areas. Example 2

[0030] See also Figure 1 As shown, based on Example 1, an irrigation disaster monitoring method based on remote sensing monitoring described in an embodiment of the present invention includes the following steps: S3: If the suspected drought-hit area is determined to be a drought-hit area, multiple drought assessment results of the drought-hit area within the preset monitoring period are obtained, and the drought results are screened and determined. The drought time series is extracted according to the time point of the drought results. The drought time series is processed and analyzed in combination with the autocorrelation function (ACF) and the partial autocorrelation function (PACF) to determine whether the drought in the drought-hit area has periodic changes. If so, enter S4; It is understandable that the time when drought occurs in irrigated areas is usually affected by many factors, with strong randomness and not stable enough. For example, there may be trend and seasonal changes. The ACF and PACF methods are applicable to various types of time series data, including time series of drought occurrence. Whether the data is stable or not, they can be analyzed. ACF and PACF can provide precise mathematical measures for quantifying the dependence in time series, which helps to more accurately determine the periodic changes in drought occurrence. The process of judging whether the drought in the drought-affected area has periodic changes includes: During the preset monitoring period, the drought index DI of the drought-affected area is obtained multiple times through the drought assessment model, and the drought assessment results with a drought index DI greater than the drought index threshold are marked as confirmed drought results. All time points corresponding to the confirmed drought results are extracted and integrated into a drought time series. Use statistical software or programming tools (such as Python's statsmodels library) to calculate the ACF and PACF values ​​of the drought time series and draw ACF and PACF graphs; In the ACF diagram, the point where the autocorrelation coefficient of the lag period exceeds the set threshold is marked as a significant autocorrelation peak; In the PACF diagram, the point where the partial autocorrelation coefficient of the lag period exceeds the set threshold is marked as a significant peak of partial autocorrelation; Obtain the number of significant autocorrelation peaks and partial autocorrelation peaks, which are A1 and A2; Obtain the time interval between every two adjacent autocorrelation significant peaks, and integrate them to obtain an autocorrelation interval group (x1, x2, x3...xn), where xn represents the nth time interval in the autocorrelation interval group; Get the time interval between every two adjacent significant partial autocorrelation peaks, and integrate them to obtain a partial autocorrelation interval group (y1, y2, y3...yz), where yz represents the zth time interval in the partial autocorrelation interval group; By formula: The periodicity value ZC is obtained, where s1, s2, s3, s4, and s5 are all weight coefficients, s1 is 0.19, s2 is 0.21, s3 is 0.17, s4 is 0.23, and s5 is 0.2; JC represents the preset monitoring period, FC1 represents the variance value of the autocorrelation interval group, and FC2 represents the variance value of the partial autocorrelation interval group; It can be understood that by calculating and analyzing the proportion of the number of significant autocorrelation peaks and partial autocorrelation peaks, the variance of the time interval and the difference in the average time interval, it is possible to judge whether the significant peaks appear in multiple lag periods, the stability of the time interval and the deviation of the time interval between the significant autocorrelation peaks and partial autocorrelation peaks in the ACF and PACF diagrams, so as to comprehensively judge whether the drought in the drought-stricken areas has periodic changes; In some embodiments, the periodicity value ZC is compared to a periodicity threshold; If the periodicity value ZC is greater than the periodicity threshold, it means that the drought in the drought-affected area has periodic changes; If the periodicity value ZC is less than or equal to the periodicity threshold, it means that the drought in the drought-affected area does not have periodic changes; then no operation is performed; S4: If the drought in the drought-stricken area has periodic changes, the monitoring system of the drought-stricken area is updated according to the periodic changes, and the drought time series is processed through the autocorrelation function (ACF) and partial autocorrelation function (PACF) to build a drought prediction model for the drought-stricken area; The step S4 of updating the monitoring system of drought-affected areas according to periodic changes specifically includes: The variance values ​​of the autocorrelation interval group and the partial autocorrelation interval group are used to add proportional weights to the autocorrelation interval group and the partial autocorrelation interval group, and the period update duration is output according to the proportional weights. The specific process includes: The proportional weights for the autocorrelation interval groups are:

[0031] The proportional weights for the partial autocorrelation interval groups are: ; Calculate the output cycle update duration GX: ; For drought-stricken areas, the monitoring interval of the drought-stricken areas is adjusted by the periodic update time; the drought-stricken areas are monitored according to the interval of the periodic update time; In S4, the process of processing the drought time series by using the autocorrelation function (ACF) and the partial autocorrelation function (PACF) to construct the drought prediction model for the drought-affected area includes: A1, use the ADF (Augmented Dickey-Fuller) test statistical method to test the stability of the drought time series. If the drought time series is stable, no treatment is performed (d is 0); if the drought time series is unstable, the drought time series is differentiated once or multiple times until a stable series is obtained. The order of the difference is the d (difference order) parameter in the ARIMA model; A2, determine the values ​​of the autoregressive order p and the moving average order q according to the truncation characteristics or attenuation trend of the ACF and PACF graphs; Exemplarily, the process of determining the values ​​of p and q is as follows: If the ACF graph decays rapidly to near zero after a certain lag order, or shows truncation characteristics, it indicates the order p of the autoregressive term (AR term); The PACF graph shows the partial correlation between the time series data and its lagged version after excluding the influence of other lagged terms. The PACF graph decays rapidly to near zero after a certain lag order, or shows truncation characteristics, indicating the order q of the moving average term (MA term); A3, using the determined parameters p, d, q and combining with the least squares estimation algorithm to fit the ARIMA model, the drought prediction model for the drought-stricken area is obtained.

[0032] The technical solution of the embodiment of the present invention is: screening and determining drought results from multiple drought assessment results in a drought-stricken area, extracting a drought time series according to the time point of determining the drought result, processing and analyzing the drought time series in combination with an autocorrelation function (ACF) and a partial autocorrelation function (PACF), judging whether the drought in the drought-stricken area has periodic changes, and if the drought in the drought-stricken area has periodic changes, updating the monitoring system of the drought-stricken area according to the periodic changes, and processing the drought time series through an autocorrelation function (ACF) and a partial autocorrelation function (PACF), and constructing a drought prediction model for the drought-stricken area. The present invention adjusts the monitoring system according to the periodic changes in the drought-stricken area, thereby improving the monitoring efficiency of drought disasters in irrigation areas, and in the periodic change analysis, constructs a prediction model according to the analysis results of the autocorrelation function (ACF) and the partial autocorrelation function (PACF), realizes the prediction of drought in irrigation areas, realizes the comprehensive monitoring of disasters in irrigation areas, and improves the comprehensiveness of monitoring. Example 3

[0033] See also Figure 2 As shown, an irrigation disaster monitoring system based on remote sensing monitoring according to an embodiment of the present invention includes: Disaster area initial identification module: Use satellite sensors to monitor each irrigation area, obtain the growth data of the irrigation area, identify the irrigation analysis area through preliminary comparison of the growth data, introduce the neighborhood matrix method to perform differential analysis on the growth data of the irrigation analysis area and the adjacent irrigation areas around the irrigation analysis area, and identify the suspected drought disaster area based on the analysis results; Disaster area determination module: obtain the historical weak growth data, historical temperature data and historical soil moisture data of the suspected drought disaster area in multiple historical monitoring periods, and use the entropy weight method to add data weights to build a drought assessment model for the suspected drought disaster area; use the weak growth data, temperature data and soil moisture data in the current monitoring period as input, and determine whether the suspected drought disaster area is a drought disaster area through the output results of the drought assessment model; Disaster periodicity analysis module: If the suspected drought-hit area is confirmed to be a drought-hit area, multiple drought assessment results of the drought-hit area within the preset monitoring period are obtained, and the drought results are screened and determined. The drought time series is extracted according to the time point of the drought results. The drought time series is processed and analyzed in combination with the autocorrelation function (ACF) and partial autocorrelation function (PACF) to determine whether the drought in the drought-hit area has periodic changes; Disaster monitoring adjustment and prediction module: If the drought in the drought-stricken area has periodic changes, the monitoring system of the drought-stricken area will be updated according to the periodic changes, and the drought time series will be processed through the autocorrelation function (ACF) and partial autocorrelation function (PACF) to build a drought prediction model for the drought-stricken area.

[0034] The above shows and describes the basic principles, main features and advantages of the present invention. It should be understood by those skilled in the art that the present invention is not limited to the above embodiments. The above embodiments and descriptions are only for explaining the principles of the present invention. Without departing from the spirit and scope of the present invention, the present invention may have various changes and improvements, which fall within the scope of the present invention. The scope of protection of the present invention is defined by the attached claims and their equivalents.

Claims

1. A method for monitoring irrigation disasters based on remote sensing monitoring, characterized in that: include: S1: Identify the irrigation analysis area through preliminary comparison of the growth data of the irrigation area, introduce the neighborhood matrix method to perform differential analysis on the growth data of the irrigation analysis area, and identify the suspected drought disaster area; S2: Based on the historical weak growth data, historical temperature data and historical soil moisture data of the suspected drought-stricken areas, and using the entropy weight method to add data weights, a drought assessment model for the suspected drought-stricken areas is constructed; and based on the output results of the drought assessment model, it is determined whether the suspected drought-stricken areas are drought-stricken areas. If so, enter S3; S3: Screen and determine the drought results from multiple drought assessment results in the drought-stricken areas, extract the drought time series, and perform drought time series processing and analysis in combination with the autocorrelation function and the partial autocorrelation function to determine whether the drought in the drought-stricken areas has periodic changes. If so, enter S4; S4: Update the monitoring system of drought-affected areas according to periodic changes and build a drought prediction model for drought-affected areas.

2. The irrigation disaster monitoring method based on remote sensing monitoring according to claim 1, characterized in that: If, among the growth data, at least one growth data is below the minimum value of the normal range, the irrigation area is marked as an irrigation analysis area; Determine the total number of irrigation areas and the weights between irrigation areas, construct a neighborhood matrix, find the weights between the irrigation analysis area and other neighboring irrigation areas based on the neighborhood matrix, and calculate the growth data difference in combination with the growth data, where the growth data difference includes the normalized vegetation index difference, the vegetation cover difference, and the leaf area index difference; If any growth data difference is positive, the irrigation analysis area will be marked as a suspected drought disaster area.

3. The irrigation disaster monitoring method based on remote sensing monitoring according to claim 2 is characterized in that: The weights between irrigation areas are calculated by combining the center point distance between irrigation areas with the Gaussian function.

4. The irrigation disaster monitoring method based on remote sensing monitoring according to claim 1, characterized in that: Obtain historical weak growth data, historical temperature data, and historical soil moisture data of suspected drought-affected areas during multiple historical monitoring periods, and perform normalization processing. The entropy weights of historical vulnerable growth data, historical temperature data, and historical soil moisture data are calculated as weight coefficients to construct a drought assessment model, and the drought index is output through the drought assessment model; If the drought index is greater than the drought index threshold, the suspected drought-hit area will be marked as a drought-hit area.

5. The irrigation disaster monitoring method based on remote sensing monitoring according to claim 4 is characterized in that: Drought assessment results with drought indices greater than the drought index threshold are marked as confirmed drought results, and all time points corresponding to confirmed drought results are extracted and integrated into drought time series.

6. The irrigation disaster monitoring method based on remote sensing monitoring according to claim 5, characterized in that: Calculate the ACF and PACF values ​​of the drought time series and draw the ACF and PACF graphs; In the ACF diagram, the point where the autocorrelation coefficient of the lag period exceeds the set threshold is marked as a significant autocorrelation peak; In the PACF diagram, the point where the partial autocorrelation coefficient of the lag period exceeds the set threshold is marked as a significant peak of partial autocorrelation; The significant peaks of autocorrelation and partial autocorrelation are processed and analyzed to obtain the periodicity value; If the periodicity value is greater than the periodicity threshold, it means that the drought in the drought-affected area has periodic changes.

7. The irrigation disaster monitoring method based on remote sensing monitoring according to claim 5, characterized in that: Get the number and proportion of significant autocorrelation peaks and partial autocorrelation peaks; Obtain the time interval between every two adjacent autocorrelation significant peaks, and integrate them to obtain an autocorrelation interval group; Get the time interval between every two adjacent significant partial autocorrelation peaks, and integrate them to obtain the partial autocorrelation interval group; The variance and mean of the autocorrelation interval group and the partial autocorrelation interval group are calculated respectively, and the periodicity value is obtained by calculating the number proportion, variance and mean deviation of the significant autocorrelation peaks and the significant partial autocorrelation peaks.

8. The irrigation disaster monitoring method based on remote sensing monitoring according to claim 1, characterized in that: The variance values ​​of the autocorrelation interval group and the partial autocorrelation interval group are used to add proportional weights to the autocorrelation interval group and the partial autocorrelation interval group, and the output cycle update duration is calculated according to the proportional weights and the means of the autocorrelation interval group and the partial autocorrelation interval group; The monitoring system for drought-stricken areas will be updated according to the periodic update duration.

9. The irrigation disaster monitoring method based on remote sensing monitoring according to claim 1, characterized in that: Conduct stationarity tests on drought time series; If the drought time series is unstable, the drought time series is differentiated once or multiple times until a stable series is obtained. The order of the difference is the difference order d; According to the truncation characteristics or attenuation trends of the ACF and PACF graphs, the values ​​of the autoregressive order p and the moving average order q are determined; The determined parameters p, d, and q are used in combination with the least squares estimation algorithm to fit the ARIMA model and obtain the drought prediction model for drought-stricken areas.

10. An irrigation disaster monitoring system based on remote sensing monitoring, characterized in that: include: Disaster area initial identification module: Identify the irrigation analysis area through preliminary comparison of the growth data of the irrigation area, introduce the neighborhood matrix method to perform differential analysis on the growth data of the irrigation analysis area, and identify suspected drought disaster areas; Disaster area determination module: Based on the historical weak growth data, historical temperature data and historical soil moisture data of the suspected drought disaster area, and using the entropy weight method to add data weights, a drought assessment model for the suspected drought disaster area is constructed; and based on the output results of the drought assessment model, it is determined whether the suspected drought disaster area is a drought disaster area; Disaster periodicity analysis module: Screen and determine drought results from multiple drought assessment results in drought-stricken areas, extract drought time series, and process and analyze drought time series in combination with autocorrelation function and partial autocorrelation function to determine whether drought in drought-stricken areas has periodic changes; Disaster monitoring adjustment and prediction module: Update the monitoring system of drought-stricken areas according to periodic changes, and build a drought prediction model for drought-stricken areas.

Images