Next Article in Journal
SARDIMM: High-Speed Near-Memory Processing Architecture for Synthetic Aperture Radar Imaging
Next Article in Special Issue
Spatial Data Infrastructure and Mobile Big Data for Urban Planning Based on the Example of Mikolajki Town in Poland
Previous Article in Journal
Cost-Effective Localization of Mobile Robots Using Ultrasound Beacons and Differential Time-of-Flight Measurement
Previous Article in Special Issue
A Spatial Analysis of Urban Tree Canopy Using High-Resolution Land Cover Data for Chattanooga, Tennessee
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification of Pasture Degradation Using Remote Sensing Data and Machine Learning: A Case Study of Obichnik

1
Faculty of Electrical Engineering-Electronics and Automation, University of Ruse “Angel Kanchev”, 7017 Ruse, Bulgaria
2
Faculty of Agronomy, Agricultural University—Plovdiv, 4000 Plovdiv, Bulgaria
3
Faculty of Agriculture, Trakia University, 6000 Stara Zagora, Bulgaria
4
Faculty of Transport, Technical University of Sofia, 1756 Sofia, Bulgaria
*
Author to whom correspondence should be addressed.
Appl. Sci. 2024, 14(17), 7599; https://doi.org/10.3390/app14177599
Submission received: 8 July 2024 / Revised: 23 August 2024 / Accepted: 26 August 2024 / Published: 28 August 2024
(This article belongs to the Special Issue Geospatial Technology: Modern Applications and Their Impact)

Abstract

:
The degradation of pastures and meadows is a global problem with a wide range of impacts. It affects farmers in different ways, such as decreases in cattle production, milk yield, and forage quality. Still, it also has other side effects, such as loss of biodiversity, loss of resources, etc. In this study, the degradation of a semi-natural pasture near the village of Obichnik, Bulgaria, was evaluated using machine learning algorithms, and an unmanned aerial vehicle (UAV) obtained visual spectrum images. A high-quality (HQ) orthomosaic of the area was created and numerous regions of interest were manually marked for training and validation purposes. Three machine learning algorithms were used—Maximum likelihood, Random trees (RT), and Support Vector Machine (SVM). Furthermore, object-based and pixel-based approaches were utilized. The obtained results indicate that the object-based RT and SVM models provide significantly better accuracy, with their Cohen’s Kappa reaching 0.86 and 0.82, respectively. The performed classification showed that approximately 61% of the investigated pasture area is covered with grass, which indicates light-to-medium degradation.

1. Introduction

Pasture and meadow degradation is a global problem affecting not only cattle breeders who depend on healthy pastures for their livelihoods but also other people who suffer from the resulting hydrological disturbances, sandstorms, resource loss, and social consequences, including population displacement. Their health has a direct and indirect impact on biodiversity, as local flora and fauna are adapted to long-term evolutionary processes that have shaped these ecosystems [1].
Rangeland degradation has become an emerging topic in pasture management and environmental protection. Although many authors and international organizations call for the protection of natural grasslands, there is no consensus definition of what constitutes natural grassland degradation. Conceptual definitions for degradation are very broad, including declines in rangeland quality, decreases in vegetation cover, productivity, economic potential, service function, recovery capacity, and diversity [2,3,4]. Others focus only on reductions in forage production or grassland productivity, invasions by non-native plant species, dominance by unpalatable species, or functions related to biochemical and water cycles [5,6,7,8]. To sum up, the term ‘degradation’, referring to the condition of natural grasslands and pastures, is widely used in different circumstances by different stakeholders, from researchers and rural extension agents to policy makers, without a commonly agreed definition.
The causes of rangeland degradation are generally attributed to a combination of overgrazing, unscientific livestock management, historical-cultural barriers to adopting modern livestock management concepts, global climate change, and soil damage by small mammals [9]. From the perspective of farmers, pasture and meadow degradation reduces the economic performance and sustainability of livestock production, which may result in a reduction in the body weight and wool yield of sheep [10], forage quality [11], cattle production [12], milk production [13], etc.
Another important aspect is the criterion used for assessing the degradation of pastures and grasslands. The United Nations has proposed the changes in the land cover, productivity, and carbon reserves as sub-indicators for monitoring land degradation processes [14]. Other monitoring indicators that are often used are vegetation cover, vegetation height, aboveground biomass, grass yield, and net primary productivity, especially in remote-sensing scenarios [15]. Nevertheless, scientists classify the degradation of pasture/meadow/grassland into different categories, using different assessment criteria. In [16], the pastures in the state of Minas Gerais, Brazil, were categorized into only two classes, based on their chemical and physical characteristics: healthy pasture and degraded pasture. In another Brazilian study [17], three degradation levels of pastures are used, and the chosen criterion is the share of exposed soil. Level 1, no or low degradation (≤20% exposed soil); Level 2, medium degradation (20–45% exposed soil); and Level 3, high degradation (>45% soil exposed). Another study for Brazil assessed the degradation of humid tropic pastures in four classes, based on the erosion of the soil [18]: no visible degradation; light degradation (pastures with sheet erosion or shallow cattle tracks); moderate degradation (some vertical erosion is visible, such as gullies); and strong degradation (pastures that are partly interrupted and vertically broken with cattle tracks as well as large gullies). In a study dealing with Kazakhstani pastures [19], the vegetation cover was the main criterion for assessing their degradation. Furthermore, it applied the Drude scale [20], where the following classes are defined:
  • Soc.: plants grow close to each other, joining with their aboveground parts (90–100%);
  • Cop. 3: plants occur in very large numbers (75–90%);
  • Cop. 2: plants occur in large numbers (50–75%);
  • Cop. 1: plants occur in considerable numbers (25–50%);
  • Cop. 1-Sp.: the species are relatively abundant (15–25%);
  • Sp.: the species are abundant, but do not form a continuous cover (5–15%);
  • Sol.-Sp.: low abundance (1–5%);
  • Sol.: the species grow sparsely (<1%);
  • Un.: the species occur in single instances.
In [21], the grassland degradation was assessed in Northern Tibet. The suggested classification was based on the grassland degradation index (GDI):
G D I = i = 1 5 D i × A i A ,
where D i is the grading score of grassland degradation, A i is the distribution area of grade i , and A is the total area of the investigated grassland. The degradation was then categorized into five classes: undegraded (GDI ≤ 1); light degradation (GDI ≤ 2); moderate degradation (GDI ≤ 3); serious degradation (GDI ≤ 4); and extremely serious degradation (GDI > 4).
Another approach is the classification of degradation, based on the coverage of a certain plant or the canopy coverage. In [22], the grassland degradation in China was assessed based on the coverage of Stellera chamaejasme (SC). Five degradation classes were defined as follows: degradation stage 1 (SC ≤ 19%); degradation stage 2 (SC ≤ 39%); degradation stage 3 (SC ≤ 59%); degradation stage 4 (SC ≤ 79%); and degradation stage 5 (SC ≤ 95%). In another study aimed at assessing the grassland degradation in Inner Mongolia, China, the canopy coverage (CC) was selected as a criterion [23], and three degradation classes were defined based on it: dense grassland (CC > 50%); moderately dense grassland (CC > 20%); sparse grassland (CC > 5%). Similarly, in [24] the vegetation cover (VC) was used as the criterion for assessing degradation: none (VC > 80%); light to moderate (VC > 60%); moderate to severe (VC > 30%); and extreme (VC < 30%).
From the performed analysis, it can be seen that numerous criteria and classifications exist, which could be used to assess the degradation of pastures, meadows and grasslands. In each individual study, the criterion and the scale/categories of degradation were chosen individually, depending on the local specifics of the investigated terrain: the observed problems, the species composition, the biodiversity, their impact on the farmers’ production, etc.
With the advancement of unmanned aerial vehicles (UAV) and satellite imaging technologies, the interest in management of grassland and pastures, their degradation, desertification, etc., has increased [25,26,27]. High-spatial-resolution drone and satellite-obtained data are becoming increasingly cheap and globally available, meaning that the use of these datasets for rangeland mapping, monitoring, and tracking is affordable [28,29]. The results shown in [30] confirm the possibility of UAV imagery being applied as a tool for precision grassland surveys. The authors reached R2 = 0.62 for their best-performing biomass estimation model and concluded that such an approach can be an alternative to conventional field methods. Different approaches are used when assessing the quality and efficiency of pastures with remote sensing, which are based mostly on vegetation indices and machine learning (ML). According to [31], the first Scopus publication on the application of ML in pasture evaluation is from 2013, and by 2020, the number of articles on the topic increased to 456, which shows that the application of ML algorithms with other metrics gives promising results. The study also stated that the most commonly used ML algorithms for this purpose are random forest (RF) and partial least squares regression (PLSR).
Different studies have used ML regression models for the evaluation of aboveground biomass and dry matter of grasslands and pastures. In [32], machine learning-based predictive models were trained for the estimation of aboveground biomass using satellite and UAV-obtained data. For large-scale evaluation, the Huber model returned the best results with R2 = 0.60, while for short-scale the best-performing models were K-nearest neighbours with R2 = 0.76 and Extra trees with R2 = 0.75. The grassland biomass was also estimated in [33], using a combination of Sentinel-2 imagery with machine learning algorithms. Even though the authors did not test their approach on large-scale grasslands, their best-performing model achieved R2 = 0.6, which they consider enough to support farmers in the optimization of the biomass management of their pastures. In [34], a UAV-mounted multispectral camera was used, in order to estimate the grass stand of pastures. Using a combination of statistical models and machine learning allowed predicting the pre-grazing pasture biomass production with an error of less than 20%. In [35], computer vision algorithms were used to extract the biophysical parameters from images, which were then used to train machine learning algorithms that predict grassland height and biomass. The study achieved up to R2 = 0.656 for pasture height estimation and R2 = 0.496 for biomass estimation.
Other studies have used ML classification algorithms for categorization of the types of grasslands/meadows based on different factors. In [36], the authors used support vector machine (SVM), RF, neural network, 3D fully convolutional network, 3D convolutional neural network, and the systematic classification method (SCM) with satellite-obtained multispectral images of grassland, in order to identify the types of vegetation. A large amount of ground-based reference data were used for training and validation, which allowed reaching accuracies between 43.42% to 99% for the different species with the SCM method. In [37], the authors tried to identify pastures in Brazil using Landsat 8 data and the RF algorithm. The study created a new pasture map for Brazil and reached an overall accuracy of 87%. In [38], an object-based image analysis was applied to Sentinel-2 data with seven vegetation indices in order to classify different types of pastures (herbaceous and shrubby). The authors used the simple linear iterative clustering (SLIC) algorithm for image segmentation and the random forest for image classification, which allowed them to achieve 87% accuracy. In another study, a pasture was classified into grazed and ungrazed categories, using ground-based camera images and the Modified Green Red Vegetation Index [39].
Different studies dealing specifically with the assessment of grasslands and pasture degradation exist, which are based on appropriate vegetation indices. In [40], shrub degradation of grasslands in Kyrgyzstan was identified based on Google Earth data. The study results showed that the shrubs’ normalized difference vegetation index (NDVI) is higher, by approximately 0.20, than the grass NDVI during the summer months. In another study, the grassland degradation index (GDI) and satellite-obtained images were used to classify the degradation of grasslands [41]. The vegetable pasture cover was derived from NDVI, and with its help the GDI was estimated. NDVI and satellite-obtained data were also used in [42], in order to assess the degradation of pastures in Brazil over a five-level scale. The obtained degradation map was validated using ground samples, and 84.2% overall accuracy was achieved. In [16], pixel-based analysis was also applied, in order to estimate the degradation of pastures. LANDSAT-8 and Google Earth data were analyzed over a three-year period in order to obtain four vegetation index signatures of healthy and degraded pastures (NDVI and three RGB indices: total brightness quotient of blue, green and blue). Similarly, in [19], a degraded pasture land in Kazakhstan was investigated with the use of satellite imaging from Google Earth, Landsat 5 and 8 as well as other reference data. Thereafter, isolinear maps were created by dividing the area of herbaceous or woody vegetation by the total area. In [24], NDVI and the modified soil-adjusted vegetation index 2 (MSAVI2) were used on WorldView-2 satellite data to assess grassland degradation in the Georgian Caucasus. Thereafter, the random forest regression model was used to fit the results.
Other studies have used machine learning to evaluate grassland degradation. In [22], the authors tried to assess it based on the percent of Stellera chamaejasme coverage. Data from ground-based cameras and deep neural networks allowed a success rate between 84% and 96% to be achieved for the different classes. In another study, multispectral and multitemporal data from two satellites (RapidEye Ortho 3A and Landsat 5 TM L1T) were used to identify the degradation of a pasture [18]. The images were analyzed using the random forest classification algorithm, where a total of 210 regions of interest from the four degradation classes were identified. Furthermore, a total of 460 polygons were used for training and validation purposes. The obtained overall accuracy for the pasture degradation data reached 61% using the Landsat data, 75% using the RapidEye data and 77.5% using the combined data. In [43], the decision tree method was used to analyze NDVI data, based on Landsat images, aimed at predicting land degradation in Brazil. The validation of the decision tree achieved Kappa 0.82.
As can be seen from the performed analysis, there are many studies dealing with management of pastures and meadows. They are aimed at assessing the pasture biomass and dry matter, identifying the vegetation species, monitoring the regenerative grazing, large-scale categorization of pasture lands, etc. There are also many studies specifically aimed at assessing and monitoring degradation, though most of them rely on vegetation indices and multispectral data. So far, very few authors have used the machine learning approach in a combination with visual spectrum images for remote assessment of pasture degradation. Considering also the broad understanding of the term degradation by the different stakeholders, a visible research gap can be observed.
The main objective of this work is to investigate the efficacy of different supervised pixel-based and object-based machine learning algorithms for remote assessment of the degradation of a pasture near the village of Obichnik, Bulgaria. To achieve this, appropriate degradation classes should be defined, corresponding to the specifics and problems of the investigated location. This study relies on UAV-obtained data from visual spectrum sensors, though satellite-obtained data are also applicable. The obtained results should provide new knowledge about the application of remote sensing and machine learning for evaluation of pasture and meadow degradation.
The rest of the paper is organized as follows. Section 2 describes the object and means of the investigation, as well as the methodology for data collection and data analysis. Section 3 presents the obtained row data and their processing, application, and evaluation by three machine learning algorithms—maximum likelihood (MxL), random trees (RT), and SVM. Furthermore, the optimal models are compared in terms of Kappa, precision, and recall in order to gain an in-depth understanding of their performance. The obtained accuracies are also compared with those obtained in other studies, in order to evaluate the effectiveness of the proposed approach. Finally, Section 4 presents the conclusions from the study.

2. Materials and Methods

2.1. The Study Area

The researched area (Figure 1) is a pasture near the village of Obichnik, Momchilgrad municipality, Kardzhali region, Bulgaria (coordinates 41.494062, 25.478968). The pasture is located in the mountainous region of the Eastern Rhodopes and is characterized as semi-natural dry pastures. The altitude is about 546 m and also falls within the Natura 2000 site under Directive 92/43/EEC (Habitats Directive).
The observed pasture territory is freely used by cattle, and no grazing system is applied. Therefore, in the long term, the species composition of the pasture is greatly disturbed, which also leads to a change in the ratio between cereal grasses, leguminous grasses, and other species, as well as lower-quality grasses. As a result, exposed areas are also observed, due to a decrease in species that are involved in the formation of healthy turf. Previously conducted field studies show that the observed territory is characterized by a rich species composition and diversity of the biological and economic qualities of the species. In the composition of the vegetation, 73 species of higher plants from 22 botanical families are found.
Three types of degradation could be identified on the investigated pasture:
  • Shrubs and trees (Figure 2a)—since no grazing system is applied on the pasture, shrubs and trees are steadily increasing their share;
  • Stones and rocks (Figure 2b)—the pasture is surrounded by abandoned old houses and other constructions, which are the source of this type of degradation;
  • Exposed soil and cattle tracks (Figure 2c)—they are created mostly by the grazing animals themselves. One of the reasons for this is the slope of the terrain, which is up to 15%.

2.2. Data Acquisition and Data Analysis

Data collection was carried out using a DJI Mavic 2 Enterprise Advanced UAV by SZ DJI Technology Co., Ltd. (Shenzhen, China). The characteristics regarding the applicability of this type of drone are related to the range of visibility that is provided, the possible speeds of movement, the capabilities of the camera, and the possibilities of use in windy weather. DJI Mavic 2 Enterprise Advanced is characterized by a speed of up to 50 kph in positioning mode, and the vertical and horizontal accuracy ranges are ±0.5 m and ±1.5 m, respectively. Its max flight time is between 24 and 31 min and it has resistance to wind speeds of up to 10 m/s. The gimbal has 3-axis (tilt, roll, and pan) stabilization, and the max control speed on a tilt is 120°/s with an angular vibration range of ±0.005°. The most important element for this study is the camera of the UAV. It uses a ½” CMOS 48 Mpx visual spectrum sensor with a resolution of 8000 × 6000 px. For the needs of this study, the UAV has been set up to take photos with ISO 100, an 80% overlap, and the relative flying altitude and flying speed were 40 m and 30 km/h, respectively.
The image processing and analysis have been implemented using the software ArcGIS Pro v. 2.8.6, developed by Esri Inc. (California, CA 92373, USA). This tool provides a wide range of options, such as combining overlapping images into a single orthomosaic map, generation of classification maps based on vegetation indices, machine learning, deep learning, etc. Furthermore, ArcGIS has numerous integrated tools for raster- and object-based analysis.

2.3. Methodology of the Study

Considering the available types of degradation on the investigated pasture, the following classes are defined as follows:
  • Grass—this class corresponds to areas with no or insignificant degradation, which are covered by grass vegetation;
  • Shrubs—this class corresponds to areas which are covered with bushes and in the long term with trees;
  • Soil—this class represents areas with exposed soil and cattle tracks;
  • Stones—this class represents areas which are covered with stones and rocks.
The methodology of this study is summarized in Figure 3. It includes the following steps:
  • Step 1. A UAV is used to make overlapping images of the investigated pasture;
  • Step 2. The overlapping images are combined into one big high-quality (HQ) map, representing the investigated pasture area;
  • Step 3. The HQ map is used to select numerous polygons (regions of interest) from each class, to be used for training and validation purposes;
  • Step 4. The HQ map is used to generate a segmentation map, by manipulating three parameters of the pixels: spectral detail, spatial detail, and minimum segment size in pixels;
  • Step 5. Using the combined HQ map, the training data, and the segmentation map, a classification map is created, using one of the available algorithms;
  • Step 6. The accuracy of the classification map is assessed by using randomly selected pixels from the chosen polygons in Step 3 and comparing them with the reference class. Steps 4, 5, and 6 are repeated numerous times until the best-performing models are obtained;
  • Step 7. Once the optimal classification models are obtained, they are used to evaluate and analyze the degradation of the pasture.
Figure 3. Overview of the study methodology.
Figure 3. Overview of the study methodology.
Applsci 14 07599 g003
In this study, two types of classification approaches are used and compared: pixel-based and object-based. Furthermore, for each one of them the following algorithms are applied as follows:
-
Maximum likelihood—the idea behind the maximum likelihood is to estimate the parameters of a probabilistic model by maximizing the likelihood function, which measures how well the model explains the observed data. By iteratively adjusting the model parameters, MxL finds the parameter values that make the observed data most probable under the model [44];
-
Random trees—a very straightforward and easy-to-implement algorithm, based on the idea of using a random subset of features and samples from the training data. The randomness improves the models’ accuracy and the overfitting, leading to a robust and more generalized model [45];
-
Support Vector Machine—it is especially useful for models that work well with a small dataset or to separate the data into clearly distinguished classes. The goal of SVM is to find a hyperplane in multidimensional space that separates different classes in classification tasks. A hyperplane is a geometric representation of the dividing line in two-dimensional space or the plane in three-dimensional space, for example. The margin is the distance between the nearest points of the two classes and the hyperplane. The goal of SVM is to maximize this margin, because a wider margin usually leads to better generalization and fewer errors on new data. The support vectors are the data points that are closest to the hyperplane. They are key to defining the hyperplane and pattern, as the margin is calculated relative to them [46].

3. Results and Discussion

An experimental study at the investigated pasture was conducted on 21 April 2023. The UAV made a total of 39 overlapping images of parts of the pasture according to the trajectory shown in Figure 4a. Next, using the ArgGIS Pro v. 2.8.6 software, a combined HQ map of the pasture was generated (Figure 4b).
According to Step 3 of the developed methodology, the HQ map was used to select a total of 280 regions of interest (ROI) for training and validation purposes (Figure 4c), covering a total of 4,004,502 pixels. Their distribution among the classes is summarized in Table 1. It can be seen that a high imbalance exists between the different classes. The used ROI from the “Soil” and “Stones” classes are 0.6% and 3.3% of the total number of pixels, respectively, and this mostly corresponds to the actual situation on the investigated pasture. Therefore, this characteristic of the training data should be considered when analyzing the obtained results.
Next, steps 4, 5, and 6 of the methodology were repeated numerous times, in order to estimate the optimal segmentation parameters, which returned the highest accuracy results for the object-based classification. The highest results were obtained with the following segmentation characteristics: spectral detail—14, spatial detail—2, and minimum segment size in pixels—50. The generated optimal classification maps using the three object-based algorithms are presented in Figure 5. Similarly, the generated maps using the pixel-based algorithms are shown in Figure 6.
The accuracy of the six classification maps has been assessed using 10,000 stratified random points, which are randomly distributed within each class, proportionally to their relative area. The overall accuracy of each classification method is presented in Table 2 and is expressed in terms of Cohen’s Kappa. The Kappa coefficient is a statistical measure, used for evaluating the performance of the trained models. It could be interpreted as follows [47]: k ≤ 0—no agreement; 0.01 ≤ k ≤ 0.20—none to slight agreement; 0.21 ≤ k ≤ 0.40—fair agreement; 0.41 ≤ k ≤ 0.60—moderate agreement; 0.61 ≤ k ≤ 0.80—substantial agreement; 0.81 ≤ k ≤ 1.00—almost perfect agreement. The obtained results show that the object-based RT classification is the best-performing model with Kappa 0.86, followed by the object-based SVM classification with Kappa 0.82. These values indicate that there is almost perfect agreement between the validation data and the predicted data. All other algorithms (pixel-based and object-based) returned worse results, which is indicated by their lower Kappa, ranging between 0.41 and 0.54, i.e., between “moderate agreement” and “substantial agreement”.
The obtained results can be compared with those obtained in similar studies, in order to evaluate the efficacy of the proposed approach. In [43], a decision tree model with validation Kappa 0.82 was obtained. These results are practically identical to the optimal SVM model in this study and slightly lower than the optimal RT model. Similar results were also received in [18], where for a RF model the obtained Cohen’s Kappa varies between 0.85 and 0.88, depending on the satellite map used. Similarly, in [24], a coefficient of determination between the predicted and observed values of R2 = 0.79 was obtained for the pixel-based analysis of grassland degradation, which is slightly lower compared to our results. Nevertheless, when performing such analyses, one should consider that the other studies used different classification categories, so no direct comparison could be made.
Considering the high imbalance in the training data between the different classes, it is worth going deeper into the confusion matrices of the object-based RT and SVM classifications, as the best-performing algorithms (Table 3 and Table 4, respectively). If compared, it can be seen that their precision and recall when identifying grasses and shrubs are almost identical. There is some difference when it comes to the identification of stones. It can be seen that the RT classification has 0.568 precision, compared to 0.401 of the SVM one, which indicates that the latter misclassified more than half of the stones as shrubs. On the other hand, the SVM algorithm returns fewer false positives, as its recall is 0.812, compared to 0.787 of the RT one. Once again, the false positives come mostly from misclassifying stones as shrubs.
Significantly different is the situation when it comes to the identification of soil. The SVM algorithm performs significantly better, and its precision reaches 0.929, compared to 0.346 for RT. The reason for this is that the RT algorithm identifies many shrubs as soil. On the other hand, SVM returns more false positives and the recall reaches 0.684, compared to 0.825 for the RT algorithm.
Examples of the misclassification of stones and shrubs are presented in Figure 7. In the upper part of the image, a large area of different shrubs (Shrubs 1) is selected, some of which were falsely classified as stones by the SVM algorithm. False-positive stones are also classified with the RT algorithm, though the areas are smaller. If we zoom in, it can be seen that the main reason for this is the color of the shrubs, as some of them are peeling. This allows us to speculate that if the investigation of the pasture degradation is made one month later when all the bushes and trees have completely turned green, the accuracy of both algorithms will increase.
In the bottom part of Figure 7 are shown samples for misclassifying stones as shrubs. The first marked stones area (Stones 1) is characterized by small white stones with a black background, and the second area (Stones 2) is characterized by small white stones surrounded by green background (grass and/or shrubs). It can be seen that both algorithms fail to correctly identify large parts of the area “Stones 1”. For area “Stones 2”, the SVM algorithm performs slightly better than the RT one. Nevertheless, this particular misclassification is not that crucial, as such mixed areas are hard to classify even with a human intellect.
The top part of Figure 8 presents a large marked area of bare soil. It can be seen that both algorithms fail to identify some parts of it, though the RT classification generally performs significantly better in this particular case. In all cases, the misclassified soil areas are falsely classified as shrubs. Furthermore, two areas of shrubs were marked. It can be seen that the RT algorithm has some misclassified grass areas, which are considered to be shrubs. On the other hand, the SVM algorithm is misclassifying some of the shrub areas as grass.
The bottom part of Figure 8 shows another example that represents a large grass area with numerous tiny defects, such as small stones and small soil areas. Both algorithms are able to identify them relatively well; however, RT classifies all small soil areas as shrubs, while SVM classifies them correctly in most cases.
Finally, the degradation of the pasture was assessed using the two best-performing algorithms and is summarized in Table 5. It can be seen that there is no significant difference between the classification models, based on RT and SVM. The relative area for classes “Grass” and “Shrubs” is almost identical, with less than 1% difference. The difference in the results of RT and SVM for the “Stones” class is slightly more than 1%, though this corresponds to an almost 100% relative difference. Similarly, for the “Soil” class, the difference between the two algorithms is less than 2%, though this corresponds to almost 200% relative difference. Nevertheless, from the point of view of a farmer, such differences can be considered insignificant, as the results for the “Grass” class are practically identical, and this is enough to evaluate the degree of degradation. The reason for this is that any further actions that are going to be taken to recover the degraded pasture will require an in-situ evaluation and micro-management.
If the degradation scale from [24] is used, the investigated pasture can be classified with “Light to moderate” degradation, which is characterized by vegetation coverage between 60% and 80%.

4. Conclusions

This study investigates the effectiveness of several machine learning methods for the identification of different types of pasture degradation using UAV-obtained visual spectrum images. The data were obtained from a pasture, located in the Rhodopa mountain near the village of Obichnik, Kardzhali region, Bulgaria. Unlike previous studies, where the degradation was estimated in terms of severity, in this paper it was assessed by classifying the pasture’s surface into four categories (grass, stones, soil, and shrubs) and by obtaining their shares.
Three pixel-based and three object-based models were trained, and their performance was assessed using Cohen’s Kappa. The obtained results indicate that the object-based RT and SVM models perform better, which is indicated by their Kappa coefficients—0.86 and 0.82, respectively. The comparison of their confusion matrices showed that the two models returned very accurate results for the “Grass” and “Shrubs” classes, with F1 scores ranging between 0.926 and 0.943. The RT model gave better results when identifying “Stones” (F1 = 0.660); however, SVM performed better with the identification of “Soil” (F1 = 0.788). In general, both algorithms showed that approximately 61% of the investigated pasture’s surface can be classified as “grass” and approximately 35% can be classified as “shrubs”, which indicates that, currently, it has light to medium degradation. It is important to continue the monitoring of the pasture’s degradation and identify appropriate measures for improving its state.
The obtained results also prove that the application of machine learning with visual spectrum data for the identification of pasture degradation can be as effective as other approaches, based on multispectral data and vegetation indices, with or without machine learning. This makes it applicable to small and medium-sized farmers using low-cost UAVs. Such an approach would allow them to optimize pastures and grasslands management and increase their financial results.
It is important to verify if the accuracies of the ML algorithms would be influenced by the annual time of UAV inspections and compare them with the results in this study. It is also of interest to investigate the application of satellite-obtained RGB images with machine learning for identification of pasture degradation, in order to compare the results with the UAV-obtained ones. Such comparisons were not performed in this study and are an object for future research.

Author Contributions

Conceptualization, B.E., A.S., T.Z. and D.S.; methodology, B.E., I.V., T.K. and N.V.; software, T.K.; validation, G.S., G.K., M.G., M.T., N.G. and I.D.; formal analysis, B.E., A.S. and T.Z.; investigation, D.S., I.D., T.Z., M.G., M.T., N.G., N.V., G.S. and G.K.; resources, A.S., D.S., T.Z. and B.E.; data curation, A.S., T.Z. and B.E.; writing—original draft preparation, B.E., I.V., T.K., A.S. and D.S.; writing—review and editing, B.E., A.S., T.Z. and D.S.; visualization, B.E.; supervision, A.S., B.E., T.Z. and D.S.; project administration, A.S.; funding acquisition, B.E. and A.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministry of Education and Science of Bulgaria under the National Research Program “Intelligent Animal Husbandry”, grant number Д01-62/18.03.2021.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The datasets used in this study are published under the CC BY 4.0 license and can be found at https://doi.org/10.6084/m9.figshare.26210276.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Gantulga, N.; Iimaa, T.; Batmunkh, M.; Surenjav, U.; Tserennadmin, E.; Turmunkh, T.; Denchingungaa, D.; Dorjsuren, B. Impacts of natural and anthropogenic factors on soil erosion. Proc. Mong. Acad. Sci. 2023, 63, 3–18. [Google Scholar] [CrossRef]
  2. Behmanesh, B.; Barani, H.; Abedi Sarvestani, A.; Shahraki, M.R.; Sharafatmandrad, M. Rangeland degradation assessment: A new strategy based on indigenous ecological knowledge of pastoralists. Solid Earth Discuss 2015, 7, 2999–3019. [Google Scholar] [CrossRef]
  3. Mganga, K.Z.; Nyariki, D.M.; Musimba, N.K.; Amwata, D.A. Determinants and rates of land degradation: Application of stationary time-series model to data from a semi-arid environment in Kenya. J. Arid Land 2018, 10, 1–11. [Google Scholar] [CrossRef]
  4. Nasiyev, B.; Shibaikin, V.; Bekkaliyev, A.; Zhanatalapov, N.Z.; Bekkaliyeva, A. Changes in the Quality of Vegetation Cover and Soil of Pastures in Semi-Deserts of West Kazakhstan, Depending on the Grazing Methods. J. Ecol. Eng. 2022, 23, 50–60. [Google Scholar] [CrossRef]
  5. Lu, X.; Kelsey, K.C.; Yan, Y.; Sun, J.; Wang, X.; Cheng, G.; Neff, J.C. Effects of grazing on ecosystem structure and function of alpine grasslands in Qinghai-Tibetan Plateau: A synthesis. Ecosphere 2017, 8, e01656. [Google Scholar] [CrossRef]
  6. Bekele, N.; Kebede, G. Rangeland Degradation and Restoration in Semi-arid Areas of Southern Ethiopia: The Case of Borana Rangeland. Int. J. Environ. Sci. 2014, 3, 94–103. [Google Scholar]
  7. Shamsutdinov, Z. Ecological restoration of biodiversity and forage productivity of degraded pasture ecosystems in the Central Asian Desert. BIO Web Conf. 2022, 43, 01025. [Google Scholar] [CrossRef]
  8. Ragimov, A.; Mazirov, M.; Nikolaev, V.; Shitikova, A.; Malakhova, S. Impact of Different Type of Cattle Grazing on the Processes of Agrochemical Degradation and Digression of Soil Cover. E3S Web Conf. 2020, 220, 01002. [Google Scholar] [CrossRef]
  9. Cao, J.; Adamowski, J.F.; Deo, R.C.; Xu, X.; Gong, Y.; Feng, Q. Grassland Degradation on the Qinghai-Tibetan Plateau: Reevaluation of Causative Factors. Rangel. Ecol. Manag. 2019, 72, 988–995. [Google Scholar] [CrossRef]
  10. Suresh, S.; Gupta, D.C.; Mann, J.S. Degradation of Common Pastures: An Economics Perspective of its Impact on Livestock Farming and Coping Strategies. Agric. Econ. Res. Rev. 2010, 23, 47–56. [Google Scholar]
  11. Stonecipher, C. Mitigation of Medusahead (Teaniatherum Caput-Medusae) through Grazing and Revegetation on the Channeled Scablands of Eastern Washington. Ph.D. Dissertation, Utah State University, Logan, UT, USA, 2015. [Google Scholar] [CrossRef]
  12. Feltran-Barbieri, R.; Féres, J.G. Degraded pastures in Brazil: Improving livestock production and forest restoration. R. Soc. Open Sci. 2021, 8, 201854. [Google Scholar] [CrossRef]
  13. Yamamoto, W.; Dewi, I.A.; Ibrahim, M. Effects of silvopastoral areas on milk production at dual-purpose cattle farms at the semi-humid old agricultural frontier in central Nicaragua. Agric. Syst. 2007, 94, 368–375. [Google Scholar] [CrossRef]
  14. United Nations. SDG Indicators. Available online: https://unstats.un.org/sdgs/metadata/ (accessed on 22 May 2024).
  15. Zhang, N.; Li, Z.; Feng, Y.; Li, X.; Tang, J. Development and application of a vegetation degradation classification approach for the temperate grasslands of northern China. Ecol. Indic. 2023, 154, 110857. [Google Scholar] [CrossRef]
  16. Quinaia, T.L.; do Valle Junior, R.F.; de Miranda Coelho, V.P.; da Cunha, R.C.; Valera, C.A.; Fernandes, L.F.S.; Pacheco, F.A.L. Application of an improved vegetation index based on the visible spectrum in the diagnosis of degraded pastures: Implications for development. Land Degrad Dev. 2021, 32, 4693–4707. [Google Scholar] [CrossRef]
  17. Galdino, S. Classification of Pasture Degradation Levels in Terms of Hydric Erosion Risk in Quartzipsamments Areas at Alto Taquari Watershed (MS/MT, Brazil). Geografia 2013, 38, 95–107. [Google Scholar]
  18. de Torres, F.N.; Richter, R.; Vohland, M. A multisensoral approach for high-resolution land cover and pasture degradation mapping in the humid tropics: A case study of the fragmented landscape of Rio de Janeiro. Int. J. Appl. Earth Obs. Geoinf. 2019, 78, 189–201. [Google Scholar] [CrossRef]
  19. Yesmagulova, B.Z.; Assetova, A.Y.; Tassanova, Z.B.; Zhildikbaeva, A.N.; Molzhigitova, D.K. Determination of the Degradation Degree of Pasture Lands in the West Kazakhstan Region Based on Monitoring Using Geoinformation Technologies. J. Ecol. Eng. 2023, 24, 179–187. [Google Scholar] [CrossRef]
  20. Drude, O. Die Ökologie der Pflanzen; F. Vieweg & Sohn: Braunschweig, Germany, 1913. (In German) [Google Scholar]
  21. Gao, Q.; Wan, Y.; Li, Y.; Lin, E. Grassland degradation in Northern Tibet based on remote sensing data. J. Geogr. Sci. 2006, 16, 165–173. [Google Scholar] [CrossRef]
  22. Yan, X.; Jiang, Y.; Chen, S.; He, Z.; Li, C.; Xia, S.-T.; Dai, T.; Dong, S.; Zheng, F. Automatic Grassland Degradation Estimation Using Deep Learning. In Proceedings of the Twenty-Eighth International Joint Conference on Artificial Intelligence AI for Improving Human Well-Being, Macao, China, 10–16 August 2019; pp. 6028–6034. [Google Scholar] [CrossRef]
  23. Hu, Y.; Nacun, B. An Analysis of Land-Use Change and Grassland Degradation from a Policy Perspective in Inner Mongolia, China, 1990–2015. Sustainability 2018, 10, 4048. [Google Scholar] [CrossRef]
  24. Wiesmair, M.; Feilhauer, H.; Magiera, A.; Otte, A.; Waldhardt, R. Estimating vegetation cover from high- resolution satellite data to assess grassland degradation in the Georgian Caucasus. Res. Dev. 2016, 36, 56–65. [Google Scholar] [CrossRef]
  25. Ogungbuyi, M.G.; Mohammed, C.; Ara, I.; Fischer, A.M.; Harrison, M.T. Advancing Skyborne Technologies and High-Resolution Satellites for Pasture Monitoring and Improved Management: A Review. Remote Sens. 2023, 15, 4866. [Google Scholar] [CrossRef]
  26. Xu, X.; Liu, L.; Han, P.; Gong, X.; Zhang, Q. Accuracy of Vegetation Indices in Assessing Different Grades of Grassland Desertification from UAV. Int. J. Environ. Res. Public Health 2022, 19, 16793. [Google Scholar] [CrossRef]
  27. Jin, E.; Du, J.; Bi, Y.; Wang, S.; Gao, X. Research on Classification of Grassland Degeneration Indicator Objects Based on UAV Hyperspectral Remote Sensing and 3D_RNet-O Model. Sensors 2024, 24, 1114. [Google Scholar] [CrossRef] [PubMed]
  28. Akumu, C.E.; Amadi, E.O.; Dennis, S. Application of drone and WorldView-4 satellite data in mapping and monitoring grazing land cover and pasture quality: Pre-and post-flooding. Land 2021, 10, 321. [Google Scholar] [CrossRef]
  29. Barnetson, J.; Phinn, S.; Scarth, P. Estimating plant pasture biomass and quality from UAV imaging across Queensland’s Rangelands. AgriEngineering 2020, 2, 523–543. [Google Scholar] [CrossRef]
  30. Michez, A.; Lejeune, P.; Bauwens, S.; Herinaina, A.A.L.; Blaise, Y.; Muñoz, E.C.; Lebeau, F.; Bindelle, J. Mapping and monitoring of biomass and grazing in pasture with an unmanned aerial system. Remote Sens. 2019, 11, 473. [Google Scholar] [CrossRef]
  31. Morais, T.G.; Teixeira, R.F.; Figueiredo, M.; Domingos, T. The use of machine learning methods to estimate aboveground biomass of grasslands: A review. Ecol. Indic. 2021, 130, 108081. [Google Scholar] [CrossRef]
  32. Alvarez-Mendoza, C.I.; Guzman, D.; Casas, J.; Bastidas, M.; Polanco, J.; Valencia-Ortiz, M.; Montenegro, F.; Arango, J.; Ishitani, M.; Selvaraj, M.G. Predictive Modeling of Above-Ground Biomass in Brachiaria Pastures from Satellite and UAV Imagery Using Machine Learning Approaches. Remote Sens. 2022, 14, 5870. [Google Scholar] [CrossRef]
  33. Chen, Y.; Guerschman, J.; Shendryk, Y.; Henry, D.; Harrison, M.T. Estimating Pasture Biomass Using Sentinel-2 Imagery and Machine Learning. Remote Sens. 2021, 13, 603. [Google Scholar] [CrossRef]
  34. De Rosa, D.; Basso, B.; Fasiolo, M.; Friedl, J.; Fulkerson, B.; Grace, P.R.; Rowlings, D.W. Predicting pasture biomass using a statistical model and machine learning algorithm implemented with remotely sensed imagery. Comput. Electron. Agric. 2021, 180, 105880. [Google Scholar] [CrossRef]
  35. Franco, V.R.; Hott, M.C.; Andrade, R.G.; Goliatt, L. Hybrid machine learning methods combined with computer vision approaches to estimate biophysical parameters of pastures. Evol. Intell. 2023, 16, 1271–1284. [Google Scholar] [CrossRef]
  36. Wei, D.; Liu, K.; Xiao, C.; Sun, W.; Liu, W.; Liu, L.; Huang, X.; Feng, C. A Systematic Classification Method for Grassland Community Division Using China’s ZY1-02D Hyperspectral Observations. Remote Sens. 2022, 14, 3751. [Google Scholar] [CrossRef]
  37. Parente, L.; Ferreira, L.; Faria, A.; Nogueira, S.; Araújo, F.; Teixeira, L.; Hagen, S. Monitoring the brazilian pasturelands: A new mapping approach based on the landsat 8 spectral and temporal domains. Int. J. Appl. Earth Obs. Geoinf. 2017, 62, 135–143. [Google Scholar] [CrossRef]
  38. Girolamo-Neto, C.D.; Sato, L.Y.; Sanches, I.D.; Silva, I.C.O.; Rocha, J.C.S.; Almeida, C.A. Object Based Image Analysis and Texture Features for Pasture Classification in Brazilian Savannah. ISPRS Ann. Photogramm. Remote Sens. Spat. Inf. Sci. 2020, V-3-2020, 453–460. [Google Scholar] [CrossRef]
  39. Evstatiev, B.; Mladenova, T.; Valov, N.; Zhelyazkova, T.; Gerdzhikova, M.; Todorova, M.; Grozeva, N.; Sevov, A.; Stanchev, G. Fast Pasture Classification Method using Groundbased Camera and the Modified Green Red Vegetation Index (MGRVI). Int. J. Adv. Comput. Sci. Appl. 2023, 14, 45–51. [Google Scholar] [CrossRef]
  40. Khazieva, E.; Verburg, P.H.; Pazúr, R. Grassland degradation by shrub encroachment: Mapping patterns and drivers of encroachment in Kyrgyzstan. J. Arid Environ. 2022, 207, 104849. [Google Scholar] [CrossRef]
  41. Motta, J.J.S.; Encina, C.C.; Guaraldo, E.; Gonçalves, A.; Gamarra, R.M.; Filho, A.C.P. The Analysis of the Degree of Grassland Degradation Using Remote Sensing. The paths of the Geography. Life 2021, 22, 201–219. [Google Scholar] [CrossRef]
  42. do Valle Júnior, R.F.; Siqueira, H.E.; Valera, C.A.; Oliveira, C.F.; Fernandes, L.F.S.; Moura, J.P.; Pacheco, F.A.L. Diagnosis of degraded pastures using an improved NDVI-based remote sensing approach: An application to the Environmental Protection Area of Uberaba River Basin (Minas Gerais, Brazil). Remote Sens. Appl. Soc. Environ. 2019, 14, 20–33. [Google Scholar] [CrossRef]
  43. Vieira, R.M.D.S.P.; Tomasella, J.; Barbosa, A.A.; Polizel, S.P.; Ometto, J.P.H.B.; Santos, F.C.; da Cruz Ferreira, Y.; de Toledo, P.M. Land degradation mapping in the MATOPIBA region (Brazil) using remote sensing data and decision-tree analysis. Sci. Total Environ. 2021, 782, 146900. [Google Scholar] [CrossRef]
  44. Le Cam, L. Maximum Likelihood: An Introduction. Int. Stat. Rev./Rev. Int. De Stat. 1990, 58, 153–171. [Google Scholar] [CrossRef]
  45. Le Gall, J.-F. Random trees and applications. Probab. Surv. 2005, 2, 245–311. [Google Scholar] [CrossRef]
  46. Pisner, D.A.; Schnyer, D.M. Chapter 6—Support vector machine. In Machine Learning; Academic Press: Cambridge, MA, USA, 2020; pp. 101–121. [Google Scholar] [CrossRef]
  47. McHugh, M.L. Interrater reliability: The kappa statistic. Biochem. Med. 2012, 22, 276–282. [Google Scholar] [CrossRef] [PubMed] [PubMed Central]
Figure 1. Geographic location of the investigated pasture.
Figure 1. Geographic location of the investigated pasture.
Applsci 14 07599 g001
Figure 2. Different types of degradation on the investigated pasture: shrubs and trees (a); stones and rocks (b); exposed soil and cattle tracks (c).
Figure 2. Different types of degradation on the investigated pasture: shrubs and trees (a); stones and rocks (b); exposed soil and cattle tracks (c).
Applsci 14 07599 g002
Figure 4. The process of data acquisition and processing: collection of UAV images (a); generation of a HQ map of the investigated pasture (b); labeling of the training objects (c); and generation of a segmented map (d).
Figure 4. The process of data acquisition and processing: collection of UAV images (a); generation of a HQ map of the investigated pasture (b); labeling of the training objects (c); and generation of a segmented map (d).
Applsci 14 07599 g004
Figure 5. Optimal object-based classification maps: MxL (a), RT (b); SVM (c).
Figure 5. Optimal object-based classification maps: MxL (a), RT (b); SVM (c).
Applsci 14 07599 g005
Figure 6. Pixel-based classification maps: MxL (a), RT (b), SVM (c).
Figure 6. Pixel-based classification maps: MxL (a), RT (b), SVM (c).
Applsci 14 07599 g006
Figure 7. Examples of false-positive classification of shrubs and stones.
Figure 7. Examples of false-positive classification of shrubs and stones.
Applsci 14 07599 g007
Figure 8. Examples of false-positive classification of grass and soil.
Figure 8. Examples of false-positive classification of grass and soil.
Applsci 14 07599 g008
Table 1. ROI and pixels for each class, that are used for training and validation purposes.
Table 1. ROI and pixels for each class, that are used for training and validation purposes.
ClassNumber of ROINumber of PixelsRelative Area
Grass941,220,99530.5%
Shrubs922,626,64265.6%
Soil3224,4520.6%
Stones62132,4133.3%
Total2804,004,502100.0%
Table 2. Comparison of the different algorithms’ overall accuracy.
Table 2. Comparison of the different algorithms’ overall accuracy.
No.Classification AlgorithmCohen’s Kappa
1Optimal object-based MxL0.53
2Optimal object-based RT0.86
3Optimal object-based SVM0.82
4Pixel-based MxL0.54
5Pixel-based RT0.41
6Pixel-based SVM0.43
Table 3. Confusion matrix for object-based RT classification.
Table 3. Confusion matrix for object-based RT classification.
ReferencePrecisionRecallF-Score
GrassStonesSoilShrubs
ResultsGrass30520103430.8960.9930.942
Stones025101910.5680.7870.660
Soil0047890.3460.8250.487
Shrubs2168059220.9850.9050.943
Table 4. Confusion matrix for object-based SVM classification.
Table 4. Confusion matrix for object-based SVM classification.
ReferencePrecisionRecallF-Score
GrassStonesSoilShrubs
ResultsGrass30208104130.8750.9830.926
Stones025903870.4010.8120.537
Soil033900.9290.6840.788
Shrubs5349857440.9810.8780.927
Table 5. Analysis of the pasture’s degradation using the selected optimal models.
Table 5. Analysis of the pasture’s degradation using the selected optimal models.
ClassArea, m2Relative Area, %
Optimal object-based RT classification
Grass606561.50
Stones1021.03
Soil2902.94
Shrubs340534.53
Optimal object-based SVM classification
Grass603761.21
Stones2212.24
Soil1111.13
Shrubs349335.42
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Evstatiev, B.; Valova, I.; Kaneva, T.; Valov, N.; Sevov, A.; Stanchev, G.; Komitov, G.; Zhelyazkova, T.; Gerdzhikova, M.; Todorova, M.; et al. Identification of Pasture Degradation Using Remote Sensing Data and Machine Learning: A Case Study of Obichnik. Appl. Sci. 2024, 14, 7599. https://doi.org/10.3390/app14177599

AMA Style

Evstatiev B, Valova I, Kaneva T, Valov N, Sevov A, Stanchev G, Komitov G, Zhelyazkova T, Gerdzhikova M, Todorova M, et al. Identification of Pasture Degradation Using Remote Sensing Data and Machine Learning: A Case Study of Obichnik. Applied Sciences. 2024; 14(17):7599. https://doi.org/10.3390/app14177599

Chicago/Turabian Style

Evstatiev, Boris, Irena Valova, Tsvetelina Kaneva, Nikolay Valov, Atanas Sevov, Georgi Stanchev, Georgi Komitov, Tsenka Zhelyazkova, Mariya Gerdzhikova, Mima Todorova, and et al. 2024. "Identification of Pasture Degradation Using Remote Sensing Data and Machine Learning: A Case Study of Obichnik" Applied Sciences 14, no. 17: 7599. https://doi.org/10.3390/app14177599

APA Style

Evstatiev, B., Valova, I., Kaneva, T., Valov, N., Sevov, A., Stanchev, G., Komitov, G., Zhelyazkova, T., Gerdzhikova, M., Todorova, M., Grozeva, N., Saliev, D., & Damyanov, I. (2024). Identification of Pasture Degradation Using Remote Sensing Data and Machine Learning: A Case Study of Obichnik. Applied Sciences, 14(17), 7599. https://doi.org/10.3390/app14177599

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop