WO2023282787A1 - Method of generating soil degradation maps - Google Patents

Abstract

The present method envisages downloading, for each point of the Earth's surface to be studied, Landsat satellite images from different satellite observation frames dating back over a period of many years. With the aid of a trained neural network, images are selected which do not contain defects that prevent the calculation of vegetation indices. An NDVI vegetation index is calculated for each usable image. The mean occurrence of low NDVI values is calculated by applying the principle of binary logic and measuring the frequency of occurrence of a binary attribute. The values of binary maps for each pixel are summed and divided by the number of Landsat scenes selected for an agricultural field. Values greater than 0.5 are taken to represent a zone of soil degradation, and binary maps of soil degradation are created using a 0.5 threshold. Empirical field measurements of the depth of the humus horizon are performed to verify the obtained results. The method makes it possible to create a reliable and maximally realistic cartographic representation of soil degradation.

Description

Method for generating soil degradation maps

The technical solution claimed as an invention relates to the field of agriculture, namely to precise precision coordinate farming, and is intended for the use of deep machine learning (neural networks) in the selection of satellite images over a long period in order to create soil degradation maps that reflect degradation areas within agricultural fields.

The use of satellite information is now widespread in agricultural work. Satellite images help to obtain dynamic information about agricultural land and respond in time to emerging deviations. At present, the development of space monitoring technologies has expanded the possibilities of operational tracking of the state of agricultural land. The use of satellite information data makes it possible to increase the accuracy of measuring biometric indicators of fields, which, using neural networks, makes it possible to interpolate the values of these indicators over a long period, and in the shortest possible time. The task of obtaining an integral characteristic of soil changes over decades is relevant. The integral characteristics are the values of the coefficients of the multitemporal line of soils. Reconciliation of computational data with empirical data obtained as a result of inspecting a particular field on the spot makes it possible to avoid obtaining unspecified information as much as possible.

Mapping areas of soil degradation is a complex and time-consuming procedure. Degradation can be predicted by modeling degradation processes and obtaining terrain characteristics from topographic maps and digital elevation models. Degradation sites predicted by modeling require further confirmation, as degradation may or may not occur under the same topography and climate. Degradation can be detected based on the analysis of a large number of satellite images in the manual interpretation mode. In this case, the accuracy of detection and determination of the boundaries of degraded territories is higher than in land-based surveys. The main disadvantages of manual interpretation of a large number of satellite images include the high cost of highly skilled human labor, which cannot be automated, and low productivity. Modern trends in automating the identification and mapping of areas of soil degradation are based on the analysis of vegetation indices of vegetation, the analysis of which causes problems if they are carried out in the automatic mode of selecting images suitable for calculations. For example, cloud masks are found in archives of satellite images, but usually these masks are not enough to diagnose the presence or absence of clouds on satellite images, and including even one unsuitable image in the processing of multi-time satellite imagery completely distorts the calculation results. The current trend in selecting satellite images is based on the use of deep machine learning, which allows you to select the necessary satellite images.

Pattern recognition is one of the most demanded tasks of artificial intelligence and can be used in a wide variety of application areas. Neural networks are currently used most widely for pattern recognition. In some cases, they show results comparable or superior to human capabilities in solving similar problems. But in a number of cases, the expansion of their use is hindered by the difficulty of interpreting the resulting solution. Indeed, for a number of the most critical applications, verification and human approval/rejection of the automatic solution is required or naturally expected.

Deep machine learning in the form of convolutional neural networks is becoming more widespread in various fields of scientific and technical activity. Neural networks are used to calculate and evaluate changes in land use over long periods of time (1984 - present). In many cases, the use of neural networks makes it possible to achieve greater accuracy of calculations than traditional methods of studying phenomena, with less labor. In recent years, a new computer vision approach based on neural networks has been successfully applied in solving various problems using satellite images. Thus, the use of convolutional neural networks (CNN) for processing color images of the earth's surface provided high accuracy in the recognition of various soil processes.

A prior analogue of the claimed technical solution is known, a method for remotely determining the state and use of agricultural land according to the Russian application for invention N ° 2018111761 of 04/02/2018, IPC: G01 W 1/00, G01W 1/08 published on 04.10.2019 The analog describes an information management tool and an automated system for monitoring agricultural land, which use unmanned aerial vehicles of an aircraft type with multispectral cameras placed on them, as well as ground equipment in the form of a multispectral camera and a thermal imager installed on a mobile tower on the ground, allowing to produce advance survey and allocate areas for detailed studies of areas of degraded land subject to water and wind erosion and waterlogging.

Despite the obvious differences in the technique used, a common feature with the proposed technical solution is the creation of maps to determine the location of eroded soils.

The disadvantages of this analogue are associated with the purpose of its development, mainly aimed at improving the effectiveness of the interpretation of the soil cover. The analogue characterizes the previous generation of methods for diagnosing soil with spectral survey and collection of thematic cartographic materials using aircraft. This analogue according to the Russian application for the invention Ne 2018111761 does not provide for the analysis of multi-time series of satellite imagery using neural networks.

An analogue of a new generation (prototype) is known, a method for pre-harvest desiccation of crops with a variable rate within one field according to the Russian patent for the invention M? equipment using satellite navigation for individual analysis of intra-field zones, previously designed on the cartographic contour of the field, and calculation of the dosage of the drug applied to improve the parameters of the field sown with crops, and for a single field, a task map previously developed by an agronomist based on the values of the vegetation index NDVI and based on the determination of the norm of the working solution during visual inspection, the fields are loaded into the on-board computer installed in the sprayer cabin, then using the on-board computer, geolocation is established using the GPS / GLONASS system, and the on-board computer in the process e sprayer operation reads the norms from the task card at a specific point in the field and in in the process of desiccation, adjusts the rate of consumption of the working solution of the drug in each section of the field according to the loaded map-task during desiccation as follows:

- based on three intrafield zones with adjustment of the norm of the working solution according to the scheme: - for the zone DO ° 1 - NDVI < 0.2 - the norm of the working solution is from 0 to 50%; - for the zone DO? 2 - NDVI from 0.2 to 0 ,45 - the rate of the working solution from 50 to 100%; - for the zone DO ° 3 - NDVI> 0.45 - the rate of the working solution from 100 to 125%;

- based on four intrafield zones with adjustment of the norm of the working solution, according to the scheme: - for the zone DO ° 1 - NDVI <0.2 - the norm of the working solution is from 0 to 50%; - for the zone DO? 2 - NDVI from 0.2 to 0.3 - working solution rate from 50 to 75%; - for zone N ° 3. NDVI from 0.3 to 0.45 - the norm of the working solution from 75 to 100%; - for the zone DO ° 4. NDVI>0.45 - the norm of the working solution is from 100 to 125%, and then the entire field is processed according to the developed task map with the established norms of the working solution.

The technical solution of the RF patent DO°2717933 has signs of similarity with the claimed technical solution, such as the creation of maps of the field inhomogeneity of the field using satellite information.

However, the relatively narrow range of application (desiccation), and the fact that the prototype method does not allow the analysis of electronic mapping in multi-time mode, reduces the effectiveness of the prototype method, while in the claimed solution, satellite imagery is selected over a multi-year period. The similarity is the study of the NDVI vegetation index, however, in the prototype, this index is calculated once, according to the analysis of one frame of a satellite imagery, and in the claimed technical solution, the NDVI index is calculated for each field for each image for a period of time since 1984. The frequency of occurrence of low NDVI in the proposed solution shows areas not with random, but with persistent annual problems in the development of plants, respectively, the ability to determine the location of soil degradation with the greatest accuracy is achieved.

The purpose of the proposed technical solution is based on the processing of multi-temporal satellite images to achieve the possibility of mapping areas of soil degradation within agricultural fields. The technical task is to create a reliable and as close as possible to reality cartographic representation of the results of the analysis of multi-temporal satellite images, for which the claimed method has been developed.

In the process of testing the method, degraded areas of arable land were indicated based on the selection of satellite images using deep machine learning methods and methods for calculating the average occurrence of low values of the vegetation index NDVI.

EFFECT: expanding the operational capabilities of mapping and automating the process of recognizing degraded soil cover areas inside agricultural fields using in practice the optimal parameters of the method for creating maps of soil cover degradation based on the results of processing multi-time satellite images. The operation of the training unit of the neural network recognition path, which provides the final classification solution for the studied field maps, increases the probability of correct classification of satellite information, and subsequent verification on the field allows to maximize the accuracy of the resulting maps. The achievement of this result is ensured by the features of the proposed method.

The essence of the proposed technical solution lies in the fact that the method of generating maps of soil degradation based on the processing of multi-temporal satellite images includes the creation of maps of soil cover degradation using satellite information, in which Landsat satellite images are downloaded for each studied point on the Earth's surface from different frames of satellite imagery for many years, then with the help of a trained neural network, images are selected that do not have defects that prevent the calculation of vegetation indices, the vegetation index NDVI is calculated for each suitable image, the average occurrence of low NDVI values (AOLNDVI ), make a binary map of the development of soil cover degradation within the agricultural field and conduct empirical field measurements of the thickness of the humus horizon to verify the results obtained.

The claimed method consistently applies the principles of binary logic and measurement of the frequency of occurrence of a binary attribute (low NDVI value) in a large number of satellite images. The manifestation of a binary feature of each specific satellite image is not an independent characteristic in this approach.

Indeed, low NDVI values in any particular year can be associated with soil degradation factors, weather fluctuations, deficiencies in agricultural technologies, properties of a particular crop, etc. A completely different situation arises when analyzing a set of binary maps of low NDVI values for 35 years or more for the same territory. If on dozens of Landsat satellite images over 35 years the pixel value most often (on more than 50% of scenes) fell into the zone of low NDVI value, then it can be assumed that soil fertility in this part of the field is reduced. In areas with low soil fertility, degradation of the soil cover can be assumed. This assumption was confirmed during the verification of the method by the method of field measurements of the thickness of the humus horizon.

The inventive method is illustrated by figures 1-4 and a table showing:

Figure 1 - Suitable and not suitable for calculating vegetation indices, images sorted by a neural network, where: 1 - suitable for calculations; 2 - not suitable due to cloudiness; 3 - unusable due to cloud shadows; 4 - parts of crops not suitable due to wetting; 5 - unsuitable due to the lack of crop vegetation; 6 - unsuitable due to snow cover; 7 - unsuitable due to coverage with crop residues; 8 - not suitable due to the presence of combustion and combustion products; 9 - 14 - not suitable due to the presence of errors in agricultural technology.

Fig.2 - Stages of calculating soil degradation maps, where: 1 - calculation of NDV1 maps for each field and each Landsat scene, selected by the neural network; 2 - highlighting the area of low NDV1 values for each NDVI map (one third of the map area); 3 - calculation of the map of the average occurrence of low NDVI values; 4 - construction of a binary map of the distribution of soil degradation, where more than 50% of the times in 35 years low NDVI values were encountered.

Fig.3 - Binary map of soil degradation, built on the basis of processing multi-temporal satellite images selected by the neural network, and the points of field measurements of the thickness of the humus horizon to verify the results obtained (the numbers indicate the points of field measurements, the results of processing, which presented in the table). Green - areas with a predicted absence of degradation, red - areas of potential degradation of the soil cover.

Fig.4 - Traditional soil map and locations of field measurements (the numbers indicate the points of field measurements, the results of processing, which are presented in the table). The numbers in circles on a traditional soil map mean: 2 - no degradation, 3 - weak wind degradation, 4 - weak water degradation, 6 - medium water degradation, 7 - strong water degradation, 9 - soils of another type.

An EXAMPLE of a specific implementation of the proposed method.

For a specific application of the proposed method, soil degradation zones were calculated based on the frequency of occurrence of low NDVI values from 1984 to 2021. Low NDVI values were calculated separately for each suitable satellite image fragment within each agricultural field. NDVI values one third of the field area and lower than the other two thirds were considered low (FIG. 2). An independent verification of the method was carried out on six agricultural fields on an area of 713.3 hectares. The content of humus and the thickness of the humus horizon were determined at 42 ground points (Fig. 3). During testing, the method gave 12.5% type I errors (false positives) and 3.8% type II errors (false negatives). The table shows the definition of the degree of soil degradation according to various criteria in the field measurements.

Table

Determination of the degree of soil degradation according to various criteria in the field measurements.

нейросети деградации*The number of the thickness of the degraded soil humus content containing the thickness of both the degradation map of the th horizon, the lower humus factor, and the factor built on the basis of the type of section _cmsaoa , humus horizon

degradation neural networks*

1 75 4.1 p

2 80 4.3 p

3 81 4.4 - - p

4 70 4.7 - - p

5 45 3.5 + - + + d

6 43 3.0 + + + e

7 78 3.1 + p

8 61 3.4 - - p

9 31 2.8 + + + + e

10 67 3.5 p 11 64 3.7 - - p 12 27 2.6 + + + + Е

*type of degradation: n - no degradation; e - wind erosion; d - water erosion.

In the areas of soil degradation identified by the proposed method, the probability of detecting actual degradation by field measurements was 87.5%. The probability of detecting soil degradation by field measurements outside the predicted areas of soil degradation was 3.8%. The results show that the use of a neural network is possible for the selection of satellite images suitable for the calculation of vegetation indices and the subsequent identification of degradation sites based on the processing of multi-temporal satellite images. This eliminates the need for intermediate filtering systems when selecting satellite images with difficulties in identifying clouds, cloud shadows, and open soil (Fig. 1). A direct selection by the neural network of Landsat satellite images suitable for calculations was made. An example of calculating the soil degradation map.

где For each NDVI calculation, the agricultural field was divided into three equal areas with NDVI values in the low, medium, and high ranges. Then the zones of medium and high NDVI values were assigned a value of "0", and the zone of low NDVI values was assigned a value of "1". A series of binary maps of the distribution of low NDVI values was obtained for each agricultural field. Then, the values of binary maps for each pixel were summed and divided by the number of Landsat satellite imagery scenes selected for the agricultural field:

where

• AOLNDVI - average occurrence of low NDVI values;

• LNDVIi - indicator of the lower NDVI zone for the i-th scene of Landsat;

• n - the number of Landsat scenes selected for calculations.

An AOLNDVI value > 0.5 was taken as a soil degradation zone. By dividing by a threshold value of 0.5, a binary soil degradation map was generated (FIG. 3).

Ground check example. Ground verification was carried out by classical methods of field research of soils. First, topographic maps and satellite images were analyzed. Then field routes and field survey sites were planned. At each point of the field survey, a soil section was laid, the soil profile was described, and soil samples were taken (Figs. 3, 4). The sampling coordinates and the location of the sections are fixed by GPS. The samples were then analyzed in the laboratory. Two indicators were measured - the thickness of the humus horizon and the content of humus in the plow horizon (shown in the table).

The accuracy of interpretation was determined by the percentage of coincidence of points of ground-based determination of the presence of soil degradation and soil degradation maps inside agricultural fields obtained by automated processing of multi-time satellite images (shown in the table). Field measurements show that the map of soil cover degradation, obtained on the basis of processing multi-temporal satellite images, is significantly more accurate than the traditional soil map. The probability of detecting non-eroded soils on the contours of a traditional soil map is 50%, while on binary maps of soil degradation obtained on the basis of processing multi-temporal satellite images, this accuracy is 87.5%. Based on the foregoing, it can be concluded that the proposed technical solution has developed a new method for generating soil degradation maps based on the processing of multi-temporal satellite images.

The proposed method for generating soil degradation maps can be used in automated mapping of degraded soils.

The proposed technical solution allows to solve the problem of clarification and optimization of cartographic information for agricultural purposes.

Such a combination of the versatility of the method with the relative ease of use to identify degraded areas of fields in the prototype is not achieved.

Based on the foregoing, we can conclude that the proposed technical solution meets the criteria of "novelty", "inventive step" and "industrial applicability".

Claims

FORM ULA AND INVENTION A method for generating maps of soil degradation A method for generating maps of soil cover degradation, including the creation of maps of soil cover degradation using satellite information, characterized in that Landsat satellite images are downloaded for each studied point on the Earth’s surface from various space imagery frames over many years, then images are selected using a trained neural network , which do not have defects that prevent the calculation of vegetation indices, calculate the NDVI vegetation index for each suitable image, calculate, using the principle of binary logic and measuring the frequency of occurrence of a binary attribute, the average occurrence of low NDVI values, sum the values of binary maps for each pixel and divide by the number of scenes Landsat selected for an agricultural field according to the formula:

AOLNDVI - average occurrence of low NDVI values; LNDVIi - the indicator of the zone of low NDVI values for the i-tooghot snapshot.); n is the number of images selected for calculations, where AOLNDVI > 0.5 values are taken as the soil degradation manifestation zone, and binary maps of the soil cover degradation manifestation are created using a threshold of 0.5, then empirical field measurements of the humus horizon thickness are carried out to verify the results obtained.