# Image Segmentation of the Sudd Wetlands in South Sudan for Environmental Analytics by GRASS GIS Scripts

Polina Lemenkova

► To cite this version:

Polina Lemenkova. Image Segmentation of the Sudd Wetlands in South Sudan for Environmental Analytics by GRASS GIS Scripts. Analytics, 2023, 2 (3), pp.745-780. 10.3390/analytics2030040. hal-04214295

HAL Id: hal-04214295

https://hal.science/hal-04214295v1

Submitted on 21 Sep 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
<figure><img src="image_2.png" /></figure>

Article

# Image Segmentation of the Sudd Wetlands in South Sudan for Environmental Analytics by GRASS GIS Scripts

Polina Lemenkova †, ‡

Laboratory of Image Synthesis and Analysis, École Polytechnique de Bruxelles, Université Libre de Bruxelles (ULB), Building L, Campus du Solbosch, ULB-LISA CP165/57, Avenue Franklin D. Roosevelt 50, 1050 Brussels, Belgium; polina.lemenkova@ulb.be or polina.lemenkova@plus.ac.at; Tel.: +32-471-86-04-59 or +43-677-617-327-72

† Current Affiliation: Department of Geoinformatics, Universität Salzburg, Schillerstraße 30, A-5020 Salzburg, Austria.

**Abstract:** This paper presents the object detection algorithms GRASS GIS applied for Landsat 8-9 OLI/TIRS data. The study area includes the Sudd wetlands located in South Sudan. This study describes a programming method for the automated processing of satellite images for environmental analytics, applying the scripting algorithms of GRASS GIS. This study documents how the land cover changed and developed over time in South Sudan with varying climate and environmental settings, indicating the variations in landscape patterns. A set of modules was used to process satellite images by scripting language. It streamlines the geospatial processing tasks. The functionality of the modules of GRASS GIS to image processing is called within scripts as subprocesses which automate operations. The cutting-edge tools of GRASS GIS present a cost-effective solution to remote sensing data modelling and analysis. This is based on the discrimination of the spectral reflectance of pixels on the raster scenes. Scripting algorithms of remote sensing data processing based on the GRASS GIS syntax are run from the terminal, enabling to pass commands to the module. This ensures the automation and high speed of image processing. The algorithm challenge is that landscape patterns differ substantially, and there are nonlinear dynamics in land cover types due to environmental factors and climate effects. Time series analysis of several multispectral images demonstrated changes in land cover types over the study area of the Sudd, South Sudan affected by environmental degradation of landscapes. The map is generated for each Landsat image from 2015 to 2023 using 481 maximum-likelihood discriminant analysis approaches of classification. The methodology includes image segmentation by 'i segment' module, image clustering and classification by 'i cluster' and 'i maxlike' modules, accuracy assessment by 'r kappa' module, and computing NDVI and cartographic mapping implemented using GRASS GIS. The benefits of object detection techniques for image analysis are demonstrated with the reported effects of various threshold levels of segmentation. The segmentation was performed 371 times with 90% of the threshold and minsize = 5; the process was converged in 37 to 41 iterations. The following segments are defined for images: 4515 for 2015, 4813 for 2016, 4114 for 2017, 5090 for 2018, 6021 for 2019, 3187 for 2020, 2445 for 2022, and 5181 for 2023. The percent convergence is 98% for the processed images. Detecting variations in land cover patterns is possible using spaceborne datasets and advanced applications of scripting algorithms. The implications of cartographic approach for environmental landscape analysis are discussed. The algorithm for image processing is based on a set of GRASS GIS wrapper functions for automated image classification.

**Keywords:** Africa; environmental monitoring; GRASS GIS; image analysis; image processing; image segmentation; remote sensing; satellite image; script; wetland

**PACS:** 91.10.Da; 91.10.Jf; 91.10.Sp; 91.10.Xa; 96.25.Vt; 91.10.Fc; 95.40.+s; 95.75.Qr; 95.75.Rs; 42.68.Wt

**MSC:** 86A30; 86-08; 86A99; 86A04

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** Lemenkova, P. Image Segmentation of the Sudd Wetlands in South Sudan for Environmental Analytics by GRASS GIS Scripts. *Analytics* **2023**, *2*, 745–780. https://doi.org/10.3390/analytics2030040

Academic Editor: Jose L. Salmeron

Received: 12 August 2023

Revised: 5 September 2023

Accepted: 19 September 2023

Published: 21 September 2023

<figure><img src="image_4.png" /></figure>

**Copyright:** © 2023 by the author. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
**JEL Classification:** Y91; Q20; Q24; Q23; Q3; Q01; R11; O44; O13; Q5; Q51; Q55; N57; C6; C61

# 1. Introduction

## 1.1. Background

Automatic segmentation is one of the important tasks required in environmental analytics and satellite image processing for pattern recognition. The Sudd is the largest wetland area in the world and the largest swamp area in the Nile Basin. The effects from hydrological, climate, and anthropogenic factors results in annual flooding of the Sudd, which yearly varies in extent and intensity. At the same time, the Sudd marshes provide key water and food resources for populations and habitats for species. A script-based framework of image processing is presented the GRASS GIS for monitoring the Sudd swamps using a time series of nine Landsat images from 2015 to 2023. The GRASS GIS techniques are successfully employed to detect changes in inundated areas of the Sudd wetlands during the last nine years to contribute to the environmental monitoring of East Africa, South Sudan (Figure 1).

<figure><img src="image_5.png" /><figcaption>**Figure 1.** Topographic map of South Sudan. Software: GMT v. 6.1.1. Data source: GEBCO. Rotated read square shows the study area of the Sudd wetlands. Map source: author.</figcaption></figure>
A satellite image is far from being a random configuration of pixels. Rather, it exhibits a high degree of organisation, e.g., reflected in its spatial and spectral properties, such as geometric shapes of land cover types, various levels of brightness, and texture of patterns. In this regard, one of the main tasks in satellite image processing is the detection of image structures, such as polygons and groups of pixels, which form a fundamental matrix of images [1]. The most widely accepted approach to this task is image segmentation, which partitions the image into several segments representing distinct regions based on the characteristics of the pixels constituting the image. The algorithms of image partition divide a whole image into multiple segments, with a general aim to discriminate relevant and important parts of the image from the entire array of pixels on the image [2–4].

## 1.2. Current Research Status

New scripting languages present more powerful approaches to remote sensing data processing, providing a cost-effective alternative to the state-of-the-art image processing software. In this work, we present a set of scripting algorithms through the Geographic Resources Analysis Support System (GRASS) GIS to satellite image processing. The biggest algorithmic challenge we face when putting the idea of script-based image processing into practice is pattern recognition, which is also the reason why the GRASS GIS concept was undertaken in this study. This paper builds on and extends our previous research on satellite image processing using scripting algorithms [5,6].

In contrast to the previously used Python and R tools, the GRASS GIS presents a more powerful toolset for image processing. Hence, after a first look at the data quality in the RGB image scenes, one may be tempted to not even try to interpret the land cover patterns due to the high similarity of the individual patches. However, in this paper, we show that an appropriate combination of the GRASS GIS modules, image enhancement, clustering, classification, and interpretation enables us to differentiate the land cover types through recognition of patches on the images taken at different time periods.

Satellite image segmentation have received significant attention in recent years. Models developed for it have been used in numerous applications, such as surficial materials mapping [7], machine-learning-based computer vision [8], change detection threshold techniques [9], contextual pattern recognition for object detection [10–12], statistical segmentation [13,14], fusion detection with spectral and thermal feature combination [15], and texture synthesis [16], among many others. The segmentation of a satellite image is based on probabilistic modelling, which is applicable to a wide range of image structures [17].

Nevertheless, much software for image processing is distributed on a commercial basis. In order to reduce the cost and increase the effectiveness of image processing, a lot of approaches are proposed as an alternative to proprietary software for image processing, some of them using programming and scripting methods and some of them using embedded algorithms and libraries. Alternative ways of image processing may be used as components of workflow in various approaches to image analysis, to enhance image quality, and to adjust bands of satellite channels for specific landscape types. At the same time, various software has different mechanisms and algorithms of image processing, controlling the techniques of raster data processing. Therefore, understanding the workflow and way of image processing, the effects of specific libraries, and their performance during scripting is of great importance.

Through detecting objects and boundaries, segmentation supplies essential information for detecting relevant landscape patches of Earth visible on spaceborne images at various scales [18]. The essential part of the segmentation algorithm consists in the partition of images using thresholds, which analyses the natural image and seeks to divide the image into the regions of interest and other parts through analysis of the properties of pixels characterising land cover types in the landscapes. Various modifications of threshold algorithms exist that aim to optimise the process. Examples include, for instance, multi-level thresholding [19–21], double threshold approaches for two-dimensional Otsu image
segmentation [22,23], weighted threshold algorithms using thread segmentation [24], and mean-shift clustering algorithms [25], to mention a few.

Specifically for remote sensing applications, segmentation recognises objects through grouping pixels with similar values of spectral reflectance identified as threshold values into unique segments on the image [26]. Threshold algorithms present an efficient method of image partition, which employs machine vision to derive the parameters of pixels discriminating the region of interest on the image from the background using properties of cells [27]. The extended approach of multilevel thresholding divides the image into multiple regions based on the level of colour intensity defined for each segment [28–30].

### 1.3. Examples of Tools and Software

Relevant examples of satellite image processing show that segmentation applied to environmental mapping gives rise to a semantically meaningful detection of vegetation assemblages, which are equivalent to habitats [31]. Selected previous works on satellite image segmentation include various developed algorithms, e.g., discriminating the regions against neighbours by semantic approach and normalisation using deep features in network convergence [32], contrasting land categories using diversity in pixels and smoothing shapes of the regions [33], iterative mean-shift clustering optimisation [34], layering images and segmenting through the R-Convolutional Neural Networks (CNNs) [35], evaluating the saliency in pixels using weighted dissimilarities in patches [36], and extracting contours by simplification [37,38].

Intuitively, using the patchy texture of images enables the detection of homogenous habitats for use in various image processing and computer vision applications in environmental analyses. The advantage of employing segmentation approaches in remote sensing data analysis is that the results are based on feature extraction independent of the choice of parametrisation of segments. A wide range of satellite image segmentations have also been reported, with case studies including the detection of shorelines [39], burned areas [40,41], or forest variables [42] and change detection in wetlands [43]. If a trackable parametrisation exists, similar to image classification, then it can be used directly with no loss of information in segmentation [44]. In such cases, the strategy of object detection in segmentation algorithms is based on the identification of the regions on the image which present an assembly of contiguous pixels that meet threshold criteria [45,46].

Conventional spectral clustering techniques have revealed critical links between polygonal approximation and the definition of the segments in image partition. This is achieved using the embedded segmentation algorithms [47–49]. More complex cases reduce the level of the fragmentation through contouring segment carcasses derived from the up-scaled colour texture features and adjusted to the level of fragmentation [50]. In this way, colour features perform better in the classification tasks, since the region is formed by an optimised size of clusters forming segments of pixels that depict their major contour and colours [51]. Such similarity between the segments and separated objects can further be used to convert the bitmap image into segments. Other algorithms include image filtering, which can be performed through the similarity of pixels analysed by Euclidean and Mahalanobis distances, as well as segmentation that splits the image into several clusters on the scene [52–54].

The success of image processing can be evaluated through measured characteristics of the land cover classes which indicate changes in landscapes over time. The traditional methods of evaluation land cover types include the classification of scenes and detecting diverse landscape types using image analysis. Such an approach provides direct information on variations in land cover types and landscape dynamics due to the relationship between spectral brightness and the properties of land cover types. As a general rule, the combination of various Landsat bands enables to detect various land cover types, along with the increased wetness or dryness of soil, which indicates the content of water and may point to the desertification of landscapes or, on the contrary, increases in swamp marshes.
## 1.4. Research Goals and Gaps

A major advantage of the machine learning approach to remote sensing applications and computer vision is that it allows optimised modelling through algorithms of automatic image processing [55-58]. Thus, specifically for image segmentation, machine learning proposes embedded techniques of image discrimination to find the contours of the objects [59-61]. However, the general workflow of image partition is still not fully automated for geospatial data and requires an optimised approach. This especially concerns such tasks as object recognition, image partition, and identification of image segments as separated objects for satellite imagery. Meanwhile, using segmentation techniques for remote sensing data processing suggests the benefits of automation for environmental monitoring through the extraction of spatial information [62].

Given the benefits of image segmentation algorithms, their geospatial application to satellite image partition promises to be an advantageous technique for environmental monitoring. Image processing techniques using only a classification approach are suitable for capturing categories of land patches, since they operate without any prior information or training samples. On the other hand, using segmentation prior to classification increases accuracy, since it helps extract features in an image using image partition, which improves the classification process. For example, landscape pattern recognition can be implemented using the partition of the bitmap satellite image by the optimisation technique of regrouping patches [63]. Moreover, the approaches of change detection based on image segmentation are often used for mapping based on remote sensing data [64].

In view of the discussed benefits of the advanced methods of image analysis, the GRASS GIS scripting framework was applied in this work for the environmental monitoring of East African landscapes with the case study of the Sudd wetlands, South Sudan. The strategy of scripting is successful in the case of automatic image processing and implementing the optimised workflow of image processing [65,66]. Automatic image processing using scripts enables avoiding the erroneous matching of pixels and misclassification while grouping cells into clusters due to finding correspondences among pixels with similar spectral reflectances using the machine-based algorithms of computer vision.

## 1.5. Motivation

The contribution of this research consists in the environmental monitoring of East Africa, South Sudan using advanced methods of programming applied to image processing. Using these methods, it was demonstrated how landscapes change over time. The innovation of this research consists in a novel developed workflow that includes several libraries of GRASS GIS for diversified steps of image processing: segmentation, clustering, classification, accuracy assessment, and mapping. It has been shown that segmentation serves as a useful seed for image classification and detection of land cover classes. Therefore, in this paper, an automated segmentation of the Landsat satellite images using a region growing and merging algorithm is presented. The employed approach includes a script-based framework by the Geographic Resources Analysis Support System (GRASS) Geographic Information System (GIS) [67]. The GRASS GIS was used due to its high computational functionality, cartographic functionality, logic, and flexibility of syntax [68]. Such advantages of the software enables one to improve the performance of image segmentation, clustering, and classification for the recognition of land cover classes in the Sudd wetlands of South Sudan, Eastern Africa (Figure 1).

The segmentation was compiled using nine Landsat images as a preprocessing step for image classification. The main idea behind segmentation is to use the collection of the raster scenes obtained from the archives of the United States Geological Survey (USGS) for the detection of landscape patches to map flooded areas of the Sudd wetlands which experience spatio-temporal changes over time [69-71]. To demonstrate the value of image segmentation techniques, we use them as priors in image classification and object detection as land cover classes. The classified series of images provides insights into the character-
istics of flooded marshes and surrounding landscapes in the Sudd region reflected in the images [72].

The problem of image segmentation is addressed by presenting the advanced scripting algorithm based on the GRASS GIS syntax [73–76]. In contrast to the existing image classification techniques which group pixels with similar spectral reflectance into classes [77–80], image segmentation is an object-based recognition techniques. It enables to identify contiguous region blocks on the images based on landscape categories. This presents a more advanced approach, which is useful both independently and linked to the next object-oriented classification process for noise reduction and to increase the effectiveness of image processing through increased accuracy and speed of image processing. Several modules of the GRASS GIS were applied to provide a new foundation for the automatic segmentation of the short-term time series of the satellite images. Of these, the most important module, ‘i segment’, was used to detect patches in wetlands, and ‘i.maxlik’ was used for image classification. Image analysis aimed at assessing the difference in the flooded areas by pixels assigned to segments as groups of the image processed as bitmap graphics. The example of satellite image segmentation and classification using GRASS GIS syntax is discussed to show how the general theory of image partition is applied to a particular case of East African wetlands.

## 2. Characterisation of Study Area

The Sudd is a large area of wetlands located in South Sudan (Figure 1). The total area of swamp varies between 30,000 to 40,000 km² according to the wet or dry seasons [81]. The origin of these wetlands is strongly related to the geologic evolution of the Nile River basin, which affected the development of nearby landscapes [82]. Thus, the Sudd wetlands were formed in the course of the geologic development of the Upper Nile. Specifically, the study area is situated around Lake No, near the Bahr al-Jabal section of the White Nile (Mountain Nile), a branch of the Greater Nile [83]. Major geologic units include the Quaternary outcrops with clayey sediments of the Cenozoic (QT) Nile floodplain, shown in Figure 2.

<figure><img src="image_6.png" /><figcaption>**Figure 2.** Surficial geologic units in South Sudan and surrounding area. Software: QGIS version 3.32. Data source: USGS. Map source: author.</figcaption></figure>
The Nile River basin forms a part of the Great Rift Valley, which originated from a system of rift and faults with correspondent geomorphic forms, such as hilly areas, valleys, and plains [84,85]. Possible movements along these faults are reflected in the recent relief and the structure of the Quaternary sediment complex around the Sudd. The topographic analysis shows that the Sudd region has a contrasting relief, with river meanders having northward-oriented general gradient [86]. Moreover, the topography of south Sudan is strongly connected to the hydrology of the Sudd swamps, which is reflected in morphological features on seasonally flooded grasslands and slopes. The effects from topography and hydrology of the Nile, together with climate factors (precipitation, atmospheric circulation, and temperature), determined the formation of the wetland ecosystems of the Sudd. Thus, the plain geomorphology of the Nile floodplain provided perfect conditions for a series of basins which serve as reservoirs and accumulate water in the Sudd marshes during wet periods [87].

The Sudd wetlands are formed as a downstream of Lake Victoria and Lake Albert in the Nile Basin in the Sudd geologic province [88] (Figure 3). Recently detected diatoms proved the existence of the large Lake Sudd, which covered central and southern Sudan during the Holocene, when active tectonic structures significantly reduced their activity, acquiring the segmented character of the Sudd wetlands [89]. Other studies also reported the existing series of the interconnected basins along the Nile distributed over the territory of the modern Sudd province (Figure 3) during the Tertiary period [90]. Such palaeographic conditions contributed to further development of the current lacustrine environment in the Sudd. The dominated soil type in the Sudd wetlands is heavy clays and fine-grained sedimentary rocks [91]. Clayey soil creates favourable conditions for the formation of wetlands due to high impermeability and low porosity, which contribute to the accumulation of water [92]. As a consequence, a highly specific hydrogeological structure of the impermeable clays results in a very limited groundwater influence on the hydrology of the Sudd, where a top layer of vertisol is about 50 cm, and sands are distributed at depths of 30 m and below [93].

<figure><img src="image_7.png" /><figcaption>**Figure 3.** Geologic provinces in South Sudan and surroundings. Software: QGIS version 3.32. Data source: USGS. Map source: author.</figcaption></figure>

Climate factors affect the Sudd swamps through water balance [94] and evaporation [95]. The average temperature is 18.5 °C, while the precipitation is 320 mm. High
evaporation over the Sudd marshes results in strong effects on the regional water cycle of the Nile hydrology, which is amplified to a large extent by the Sudd wetland area: it is the largest wetland area existing in the world and the largest freshwater swamp region in the Nile Basin [96]. The environmental measures were undertaken to decrease evaporation from the Sudd by constructing Jonglei channels [97].

Climate effects on the Sudd wetlands are related to changes in precipitation and temperature: the increase in rainfall during the El Niño phases leads to warming and a rise in temperatures [98]. The dry season includes summer months, while the rainy season includes autumn–winter months. The transitional period is spring time. Climate factors threaten the hydrology and environmental sustainability of the Sudd wetlands [99]. Other climatic issues are related to the rise of Lake Victoria in 1960s, which triggered water losses in the Sudd [100]. The integrated effects from all these factors result in the highly unstable dynamics of the Sudd’s hydrology. Thus, the intensity of annual flooding differs significantly by years and affects the extent of the wetlands [101]. During the wet season, the Sudd increases in extent by almost twice due to the excess of water, which results in an extended area of floodplains affected by recurring inundation.

The Sudd wetlands play a strategic role in livelihoods, environmental sustainability, biodiversity balance, and the maintenance of water resources in South Sudan [102]. The value of the water resources in the Sudd relies on its economic and environmental services, high biodiversity impact, and fishery and food resources [103], which are necessary for social development and the existence of the local population [104]. At the same time, the ecosystems of the Sudd form a part of the global tropical wetland system, which is an important source of biodiversity and carbon storage in soils and vegetation [105], contributing to biogeochemical cycles and climate regulation [106]. The Sudd wetlands are known for hierarchical and complex food webs with diverse types of aquatic plants, animals, and microbial communities. The dominating vegetation types in the Sudd include papyrus, herbaceous plants, water hyacinth, marsh sedges, and grasslands [107]. Their distribution differs by habitat in the open-water areas with floating and submerged plants, as well as seasonally flooded grasslands occupied by the adapted plants.

Climate–hydrological fluctuations have cumulative environmental effects on the sustainability of the Sudd ecosystems. Thus, during the flood period, large areas of grasslands in the permanent Sudd swamps are inundated, which triggers fish migration into other sections of the floodplain [108]. Human-related factors affecting the Sudd ecosystems include overexploitation of the natural resources, increased pollution [109], landscape fragmentation [110], and habitat changes [111], which are reflected in the recent dynamics of land cover types in the Sudd region. At the same time, small grassland patches are hotspots of biodiversity in the fragmented landscapes and should be conserved for environmental sustainability [112].

### 3. Materials and Methods

In the following section, we present the scripting algorithm of GRASS GIS to process a series of satellite images through the powerful functionality of its tuned modules that are adjusted to diverse tasks of image processing and modelling. Our procedure works by running the script via the terminal command console shell, which makes it flexible and easier to handle and interpret the spaceborne data. Our pipeline then creates the set of commands used for image processing via the modules that include the multiline code snippets integrated in a script for image processing (Figure 4).

#### 3.1. Software and Tools

A general scheme summarising the approach in this study is visualised in Figure 4. The data include topographic DEM, geologic layers, and remote sensing data processed with GRASS GIS, GMT, and QGIS. While the GRASS GIS tool was used for image processing, other software was used as an auxiliary tool for creating the methodology flowchart and topographical and geological mapping.
<figure><img src="image_8.png" /><figcaption>**Figure 4.** Methodological flowchart scheme. Software: R version 4.3.1. Graph source: author.</figcaption></figure>

## 3.2. *Data Collection and Import*

The full dataset included in the framework is available at the United States Geological Survey (USGS) (Figure 5).

The Landsat images from 8 and 9 OLI/TIRS sensors were downloaded from the EarthExplorer repository: URL https://earthexplorer.usgs.gov/ (accessed on 30 July 2023). The EarthExplorer aims at Earth observation data, collected for the multiscale monitoring of Earth-related processes using remote sensing data. It is supported by the USGS, which coordinates and promotes storage of the datasets in digital format for queries. It enables downloading the satellite images and provides cartographic information for data. Besides Landsat, it also supports other remote sensing products, such as radar data, aerial imagery, Digital Elevation Models (DEM), Advanced Very High Resolution Radiometer (AVHRR), etc. The reason for the nine satellite images consists in demonstrating the dynamics of the land cover types within a comparable period. Another reason consists in the available images with low cloudiness, which are below 10% for the scenes taken in the appropriate seasons.

The images cover 9 years, from 2015 to 2023, Figure 5. These years were chosen for the two following factors: the availability of the cloud-free images or those with minimised cloudiness (below 10%) and the gap between the years enabling the comparison and analysis of environmental dynamics. The choice of image data is explained by the two criteria: (1) the cloudiness of the image is below 10% for all the scenes; (2) the images are captured in the dry period to objectively assess the postflood scenario. The wet period of the Sudd region is accompanied by heavy rainfalls that last throughout the summer period. The exact period differs significantly by year and, according to various information sources, may last between April and September [113] or from March to October, according to the Climate Change Knowledge Portal. The data import was performed using the 'r.import' command in the GRASS GIS module. The scripts are shown in Appendix A of this study with the code for data import.

The greatest concentration of rainfalls in the Sudd is recorded between June and September [114]. Therefore, all the images were taken during the dry period to avoid the months from June to September. Since flooding and seasonal dynamics vary yearly, each
scene was inspected to evaluate the quality and distinguished contours of the diverse land cover types. Each Landsat OLI/TIRS image included eleven spectral bands in visual, panchromatic, and near-infrared channels. Besides the satellite Landsat images, this study also included auxiliary data, such as topographic data (GEBCO/SRTM grid), geologic USGS data, and descriptive information from textual sources regarding the social and economic activities in the Sudd region, as well as statistical and descriptive environmental reports on South Sudan available online), which were used for environmental analysis.

<figure><img src="image_9.png" /><figcaption>**Figure 5.** The Sudd wetlands in the Landsat images in natural colours for 8 recent years (2015–2023). The acquisition period of each image is as follows: (a) 8 January 2015, (b) 12 February 2016; (c) 31 December 2017; (d) 1 February 2018; (e) 8 March 2019; (f) 26 March 2020; (g) 29 March 2021; (h) 19 January 2022; (i) 14 May 2023.</figcaption></figure>

### 3.3. *Data Preprocessing*

The overview topographic map of the study area shown in Figure 1 was mapped using the Generic Mapping Tools (GMT) software version 6.1.1. [115]. The applied GMT
scripting technique was derived from the existing works [116,117]. The geological maps were plotted using QGIS. The remaining workflow was performed using the GRASS GIS using diverse modules, following existing similar works using GRASS GIS in environmental applications [118–121]. A folder with the uploaded Landsat imagery was stored in the 'Location/Mapset' working directory of the GRASS GIS with relevant subdirectories. Here, all the map layers were located and hierarchy-supported for imported satellite images following the standard GRASS GIS workflow [122]. The codes were written using the Xcode and run from the GRASS GIS console. The files were imported in TIFF format using Listing A1 and stored in the WGS84 coordinate system.

The atmospheric correction method included the algorithms of the embedded GRASS GIS module 'i.landsat.toar'. This module enables to convert the DN pixel values to reflectance values using DOS1 from Digital Number (DN) to reflectance. Before creating an RGB composite, it is important to perform the atmospheric correction and thus convert the digital number data (DN) to reflectance or radiance. Otherwise, the colours of a natural RGB composite do not look convincing but rather hazy. This conversion is performed using the metadata file which is included in the dataset with i.landsat.toar, which calculates the top-of-atmosphere radiance or reflectance and temperature for the Landsat MSS/TM/ETM+/OLI.

The general outline of the data processing includes data capture of the Landsat 8-9 OLI/TIRS images, data preprocessing, data conversion, segmentation, classification, validation, computing the land cover classes by covered area, and mapping. The workflow aimed at detecting variations in the wetland areas and identifying the extent of the flooded area in the Sudd marshes over the nine-year period. For a comparison of gradual changes, a one-year interval was selected between each pair of images. To ensure that technical requirements of code quality are met, several tests were carried out for Landsat images on various years using different parameters of segmentation threshold. This aimed at analysing the behaviour of the algorithm for various levels of image fragmentation using different threshold levels. The obtained image samples were stored in a separate folder of the GRASS GIS with a path to the working folder of the repository.

### 3.4. Metadata and Extent

A dataset of the nine Landsat 8-9 OLI/TIRS images containing TIFF raster files was analysed for metadata, with the parameters summarised in Table 1. The Landsat images were taken on the following periods: 8 January 2015, 12 February 2016, 31 December 2017, 1 February 2018, 8 March 2019, 26 March 2020, 29 March 2021, 19 January 2022, 14 May 2023.

<table><caption>**Table 1.** Metadata for Landsat 8-9 OLI/TIRS images.</caption><thead><tr><th>Proj.</th><th>Zone</th><th>Dat.</th><th>Ellips.</th><th>N</th><th>S</th><th>W</th><th>E</th><th>Nsres</th><th>Ewres</th><th>Rows</th><th>Col</th><th>Cells</th></tr></thead><tbody><tr><td>UTM</td><td>36</td><td>WGS84</td><td>WGS84</td><td>915,615</td><td>682,785</td><td>190,785</td><td>419,115</td><td>30</td><td>30</td><td>7761</td><td>7611</td><td>59,068,971</td></tr></tbody></table>

Abbreviations in Table 1: Proj—projection; Dat.—datum; Ellips—ellipsoid; N—north; S—south; W—west; E—east; nsres—resolution in north–south direction; ewres—resolution in east–west direction; cols—columns.

Visual bands of the original Landsat images (channels 1 to 7) were uploaded in TIFF format, processed and converted into several segments. The region extent and groups of the bands were defined on the Landsat images using the parameters in the scene by the snippet of code presented in Listing A2. The region borders for the Landsat scenes covering the Sudd area are as follows: north: n = 915,615; south: s = 683,085; west: w = 185,985; east: e = 414,315.

### 3.5. Defining Segments

The approach of the GRASS GIS to image segmentation is based on the use of the module 'i segment'. The algorithm groups similar pixels on the satellite image into unique segments. The thresholding algorithm assigns pixels on the image based on the similarity between two neighbour segments and detects the segments. This enables the detection of
the flooded area and automatically recognises changes in the landscapes. The segmentation was performed with a 90% threshold and minsize = 5; the process was converged in 41 iterations. The process was completed in 50 min for the Landsat satellite image for 2023. The region IDs were assigned to all the regions, including the remaining single-cell regions. Overall, 8464 segments were created for one Landsat image (2023).

Afterwards, the minimal size was changed to 100 to check the effects from the modified parameters on the results of the segmentation process. For modified parameters, the seeds were used to optimise the procedure and to provide the basis information for image classification. These included random segments which were selected automatically using previous segmentation rounds and used to start the segmentation process anew, using the code shown in Listing A3.

### 3.6. Threshold Algorithm

The threshold algorithm searches for the bounds of each segment on the image and plots the image generated using threshold parameters according to the similarity level below the input threshold for a coarse analysis. The rise of threshold level increases the fragmentation of the segments accordingly; see Figure 6. In turn, if the similarity distance is smaller, the pixels are assigned to other neighbour segment. Afterwards, the algorithm sets a start-end position, and the process is repeated iteratively until no more merges are possible for the segments of landscape patches during a complete pass of image segmentation. The segmented image is then visualised using the code in Listing A4 and saved as standard TIFF output format in full-resolution mode using GRASS GIS. The information on the Landsat scene is retrieved from the file, and the segments are visible in the visualised image.

<figure><img src="image_10.png" /><figcaption>**Figure 6.** Segmentation of the Landsat 8-9 OLI/TIRS image of the Sudd area for 2023: (a) segmentation parameters include minsize = 5 and threshold = 0.90, (b) minsize = 100 and threshold = 0.05 and included seeds from the previous segmentation.</figcaption></figure>

### 3.7. Image Segmentation

The images resulted from the segmentation at the two levels of threshold are shown in Figure 6. The code for processing the image with script is based on defining the objects through image segmentation that contains patches of regions with similar parameters of the pixels located inside. Here, the similarity between the current segment and each of its neighbours is computed using a search algorithm, which includes the given distance formula for a target segment. The processing of these segments is possible at both high and low threshold modes and defined similarity parameters. It displays merged segments if they meet technical criteria and analyses the coverage of the valid segments using the
algorithm function on an image. Then, it computes the position for each pair of the
segments, which shows the best mutual similarity in a target region of the image using
iteration for each next region.

Accordingly, the search for the closest segment is based on a similarity between the
segments and objects. In such a way, the algorithm is propagated along the image, searching
for every consecutive object iteratively to determine which objects are merged. The values
of the smaller distance between the objects were evaluated to indicate a closer match within
each iteration on the image. Thus, a similarity score of zero is assigned for identical pixels
which are assigned to an identical segment. In case of the lower threshold of segments
on the images, the similarity between the two segments is lower than a given threshold
value. In such a case, the combination of these region is performed using the minimal size
parameter. According to this principle, all the segments with a lesser number of pixels
are merged with their similar neighbour. Such an approach enables the optimisation of
segments through the distribution of the array of pixels into the segments.

3.8. *Parameter Estimation*

The estimation of segmentation parameters was based on the tested variants of the threshold value of the segments in a relative number, which is always between 0.0 and 1.0. A changed degree of segment fragmentation is visible in trial cases (Figure 5). The tested segment size has a thickness varying from a threshold = 0.90 to a threshold = 0.05, and a seed minimal size is defined at 100 pixels. The repetitive iterations described above divided the image into several segments, indicating land cover classes. The resolution of 30 pixels for the Landsat image is used as the optimal parameter for a given landscape patch, allowing one to indicate small segments on an image. A lower threshold allows only large groups of vegetation to be merged using valued pixels with similar spectral reflectance values. In contrast, a higher threshold (close to one) allows neighbouring land cover classes to be merged. Thus, the threshold level of image segmentation is scaled to the actual data range.

To reduce the noise effect and to optimise data processing, a minimal size greater than
one was added as an additional step of image processing. During this step, the threshold
of segmentation is ignored to avoid too fragmented images. Thus, for segments smaller
then the defined size, this parameter merges tiny patches with their most appropriate
neighbours. In such a way, the original Landsat scene is partitioned by the complete image
according to the values of the pixels' colour and intensity. The process of thresholding
was based on the analysis and separation of the pixels compared against the value of the
minimum segment size. It aims to discriminate the meaningful part of the image containing
landscape patches from the noise pixels.

3.9. *Clustering*

The performance of the segmentation was a slow process and required further post-
processing for identification of the land cover classes within the Sudd area. To address
this, a three-step process was implemented: First, the segmentation was executed as dis-
cussed in the previous subsection. This produces segments that include identity grid
pixels, in which we expect to landscape patches. The grid cells were then generalised,
and the classification process was performed using the ‘i.maxlik’ module, which uses the
maximum-likelihood discriminant analysis classifier. Third, the accuracy assessment is
performed using the GRASS GIS module ‘r.kappa’, which computes error matrices and
kappa parameters for accuracy assessment of classification for each Landsat scene. The
code for the maximum-likelihood classification is presented in Listing A5.

The clustering algorithm approach by the 'i.cluster' module is a powerful tool for
image partition prior to classification. It operates on an image using defined parameters,
such as initial number of clusters, minimum distance between them, coherence between
loops, and minimum area for each cluster. The robustness of clustering is tuned by the
correspondence between the iterations and the defined maximum number of loops. The
initial cluster means for each band are defined by values of the first cluster as a band mean corrected for standard deviation, and all other clusters are distributed equally between the first and last clusters, as implemented in GRASS GIS.

The flexibility of clustering is that all clusters less than the defined minimum are merged smoothly, which outperforms similar algorithms of image partition. Hence, the clusters are regrouped accordingly in the iterative way. Here, each pixel is assigned according to the closest distance to a given cluster using the algorithm of Euclidean distance, and the results are saved in a signature file, which is then used for classification by the 'i.maxlik' module of GRASS GIS. The clustering report is generated automatically for an image, with an example presented in Appendix C. Thus, clustering presents the advanced object-based detection method. which creates signatures for the next step of image classification, which computes the distance between pixels using the similarity method.

### 3.10. Classification

Afterwards, the maximum-likelihood discriminant analysis classifies the generated clusters, segments and covariance matrices, computed previously. These are used to define categories of each evaluated cell. The 30-resolution remote sensing imagery corresponds to the following land cover types in the Sudd region [123]: (1) Cropland; (2) Herbaceous coverage; (3) Forest; (4) Mosaic tree canopies; (5) Shrubland; (6) Grassland; (7) Flooded and inundated areas; (8) Bare areas; (9) Built-up areas; (10) and Water areas. These land cover types were used for landscape analysis. Hence, the pixels are assigned to the categories of the land cover classes according to the calculated highest probability of belonging to the given class on the image based on their spectral reflectance. The assignment of pixels into classes of land cover types is based on the signature file ("signaturefile"), which contains the cluster and covariance matrices calculated by the module 'i.cluster' and is shown in Listing A5. In such a way, the maximum-likelihood classifier partitions the total number of pixels on the whole Landsat scene using the segmentation and clustering results as preprocessing steps for image classification that sequentially examine all current segments in the raster map.

The pixels were grouped into segments representing land cover types as separate segment objects, and the map is generated for each Landsat image from 2015 to 2023 using the maximum-likelihood discriminant analysis approach of classification. The description is created for the segments of each land cover type of the Sudd in the relevant images. The segments of land cover classes which have a sharp transient from neighbouring class are considered as another category. Afterwards, valid segments in connected landscapes are regrouped using the trial tests for various threshold parameters. The segments from the Sudd landscapes are combined into maps and compared for various years from 2015 to 2023. The distinct classes are detected using the GRASS GIS algorithms iteratively for each segment using image analysis.

### 3.11. *Calculating the NDVI*

To evaluate the distribution of healthy vegetation over the study area in the postflood scenario, the NDVI was computed using the script in Listing A6. For atmospheric corrections, the values of the Digital Number (DN) of the pixel were converted to reflectance values using DOS1 to avoid a hazy background. This conversion is performed using the metadata file which is included in the dataset with 'i.landsat.toar' and information about the sun elevation for each scene. Using the 'i.landsat.toar' module, the top-of-atmosphere radiance and temperature were corrected for the Landsat OLI sensor. Afterwards, the NDVI was computed using the existing combination of the Red and NIR bands of the image, and the images were visualised using the script presented in Listing A6.

### 3.12. Accuracy Assessment

The accuracy assessment was performed in two ways. First, the error matrix and kappa parameters were computed by the 'r.kappa' module of GRASX GIS. Second, the
rejection probability classes were calculated to estimate the pixel classified according to
confidence levels based on the classification of the satellite images. This was implemented
using the 'i.maxlik' module and resulted in plotted maps showing the rejected threshold
results, as shown above in Listing A5. The confusion matrix was estimated by kappa with
computed possible misclassification cases and derived kappa index of agreement. This was
performed using the code in Listing A7.

Image processing by GRASS GIS also included the removal of the noise signals from
the images through the adjusted segmentation threshold and image partitioning into
segments for monitoring inundated areas and classification. The results of the kappa
calculation are presented as confusion matrices in a tabular format (kappa.csv). These
tables were computed for each classified Landsat image and present the calculation results,
reporting data for every category of land cover classes, as summarised in the Appendix B.

4. Results

4.1. *Remote Sensing Data Analysis*

The algorithms of the GRASS GIS described above and summarised in scripts were applied to process the Landsat satellite images, with the results of the segmentation shown in Figure 6. Scenes were segmented for each year with the visualised maps. Remote sensing data organisation and management were performed using the GRASS GIS software. It presents a multifunctional GIS as well as workspace and editing system for remote sensing and cartographic data storage and processing [122]. Its effective functionality enables one to perform various steps of image processing and spatial data processing: storage and organising, navigating and visual inspection, projecting, formatting and converting, image analysis, handling metadata adding annotations on maps, visualising and analysing diverse features related to Earth observation data, and mapping. Other advantages include open-source availability and double-mode functionality: using either scripts or a Graphical User Interface (GUI). Moreover, GRASS GIS enables one to operate large datasets in vector and raster formats. Here, each folder of Landsat images contained 800–900 MB, which resulted in the processing of 9 Gb for nine satellite images, effectively processed by the GRASS GIS.

4.2. *Detection of Segmented Areas*

The results of Landsat 8-9 OLI/TIRS image segmentation are presented in Figure 7.

The computed areas by land cover classes and their changes by years are summarised in Table 2. The classes signify the following land cover types: (1) Cropland; (2) Herbaceous coverage; 473 areas (3) Forest; (4) Mosaic tree canopies; (5) Shrubland; (6) Grassland; (7) Flooded and inundated; 474 areas; (8) Bare areas; (9) Built-up areas; (10) Water areas. These land cover types were used 475 times for landscape analysis.

A multiscale time series analysis demonstrated changes in flooded areas of the Sudd region, revealed in the maps by gradual changes in the landscapes, which are visible on the segmented Landsat images taken yearly from 2015 to 2023. For the Landsat 8-9 OLI/TIRS scenes, the results of the image segmentation performed using the identical parameters defined for all the images are summarised in Table 3. Here, the number of iterations (passes) depends on the level of fragmentation of the image and varies by years. Segmentation and processing of the Landsat 8-9 OLI/TIRS datasets by GRASS GIS served for detecting flooded areas and monitoring the inundated areas in the Sudd marshes. The detected large seasonal changes in wetlands result from variation in the flow of Nile tributaries and Victoria Lake, as well as climate-hydrological pulsing. The approach is based on the analysis of landscape patches, grouping segments on the images, and hierarchical clustering (subgroups of landscapes using various threshold levels) for classification following image segmentation (Figure 7).
<figure><img src="image_11.png" /><figcaption>**Figure 7.** Segmentation maps of the satellite images of the Sudd wetlands, South Sudan based on the time series of the Landsat 8-9 OLI/TIRS images (2015–2023).</figcaption></figure>

Detecting and recognising the segments on the Landsat images implies identifying the fields of regions which correspond to the 10 major land cover types in the Sudd wetlands of South Sudan. These regions are grouped into semantic categories (landscape patches and land cover types). The hierarchical level of the geometric objects with regard to their scale (small-, middle-, and large-size) was applied, and the threshold was optimised. The detection of segments on the images was based on the mean shift image segmentation, including filtering and clustering (Figure 7). Since both algorithms are embedded in the 'isegment' module, the images were segmented without the priori landscape-dependent information, which ensured independent interpretation.
<table><caption>**Table 2.** Computed areas of land cover classes derived from the classified Landsat 8-9 OLI/TIRS images for the region of the Sudd, South Sudan in the period from 2015 to 2023.</caption><thead><tr><th>Year</th><th>Class 1</th><th>Class 2</th><th>Class 3</th><th>Class 4</th><th>Class 5</th><th>Class 6</th><th>Class 7</th><th>Class 8</th><th>Class 9</th><th>Class 10</th></tr></thead><tbody><tr><td>2015</td><td>1,501,537.8</td><td>243,574.1</td><td>317,231.5</td><td>17,077.0</td><td>241,763.7</td><td>734,523.3</td><td>447,070.3</td><td>528,704.9</td><td>466,546.4</td><td>281,024.0</td></tr><tr><td>2016</td><td>606,130.3</td><td>594,056.0</td><td>251,200.3</td><td>478,775.3</td><td>108,958.6</td><td>202,940.2</td><td>558,836.7</td><td>132,573.3</td><td>471,370.5</td><td>386,764.2</td></tr><tr><td>2017</td><td>344,062</td><td>395,409.0</td><td>458,500.0</td><td>458,986.4</td><td>526,807.8</td><td>540,912.0</td><td>761,353.5</td><td>468,330.1</td><td>664,728.9</td><td>873,149.5</td></tr><tr><td>2018</td><td>427,690.5</td><td>1,054,373.0</td><td>1,083,009.3</td><td>881,093.1</td><td>881,190.4</td><td>705,596.0</td><td>526,685.1</td><td>315,125.8</td><td>1,128,926.5</td><td>495,258.0</td></tr><tr><td>2019</td><td>376,759.4</td><td>332,136.8</td><td>251,513.7</td><td>982,144.2</td><td>124,956.3</td><td>757,819.1</td><td>1,090,100.0</td><td>625,188.6</td><td>109,114.2</td><td>689,025.5</td></tr><tr><td>2020</td><td>272,113.4</td><td>180,816.8</td><td>419,207.8</td><td>402,059.3</td><td>493,457.7</td><td>1,156,585.4</td><td>937,007.5</td><td>584,898.4</td><td>1,567,576.9</td><td>752,311.4</td></tr><tr><td>2021</td><td>416,413.7</td><td>234,023.4</td><td>189,669.0</td><td>414,165.1</td><td>490,691.7</td><td>578,060.7</td><td>118,864.8</td><td>467,398.9</td><td>283,889.0</td><td>355,065.7</td></tr><tr><td>2022</td><td>307,658.6</td><td>1,168,006.6</td><td>147,964.0</td><td>423,542.0</td><td>617,279.5</td><td>442,903.9</td><td>383,683.5</td><td>788,245.4</td><td>770,970.0</td><td>267,592.8</td></tr><tr><td>2023</td><td>0.0</td><td>29,203.0</td><td>1,599,665.7</td><td>434,129.5</td><td>1,012,475.1</td><td>1,032,878.3</td><td>934,255.0</td><td>1,214,757.6</td><td>167,164.5</td><td>687,965.9</td></tr></tbody></table>

<table><caption>**Table 3.** Results of the segmentation procedure for Landsat 8-9 images with No of created segments.</caption><thead><tr><th>Year</th><th>Scene ID</th><th>Iterations</th><th>Segments</th></tr></thead><tbody><tr><td>8 January 2015</td><td>LC08_L1TP_173055_20150108_20200910_02_T1</td><td>37</td><td>4515</td></tr><tr><td>12 February 2016</td><td>LC08_L1TP_173055_20160212_20200907_02_T1</td><td>37</td><td>4813</td></tr><tr><td>31 December 2017</td><td>LC08_L1TP_173055_20171231_20200902_02_T1</td><td>38</td><td>4114</td></tr><tr><td>1 February 2018</td><td>LC08_L1TP_173055_20180201_20200902_02_T1</td><td>36</td><td>5090</td></tr><tr><td>8 March 2019</td><td>LC08_L1TP_173055_20190308_20200829_02_T1</td><td>34</td><td>6021</td></tr><tr><td>26 March 2020</td><td>LC08_L1TP_173055_20200326_20200822_02_T1</td><td>39</td><td>3187</td></tr><tr><td>29 March 2021</td><td>LC08_L1TP_173055_20210329_20210408_02_T1</td><td>35</td><td>2445</td></tr><tr><td>19 January 2022</td><td>LC09_L1TP_173055_20220119_20230501_02_T1</td><td>35</td><td>4413</td></tr><tr><td>14 May 2023</td><td>LC09_L1TP_173055_20230514_20230514_02_T1</td><td>41</td><td>5181</td></tr></tbody></table>

Notation for Table 3: The Landsat images were selected with cloudiness below 10% to achieve maximal distinguishability of the contours of the images.

The grouping decision is made using the features of pixels that match target classes ('segment'/'not segment'), colour, distance to the threshold in pixels, and the spectral reflectance of the pixel. The shape of the segment is defined through the boundary constraints, which limits the adjacency of pixels and segments on a satellite image. In such a way, the image is represented as a vector geometric structure, recognised and identified by a computer vision approach. The assessment is based on the connectivity of pixels constituting the segment, except for the threshold fitness and the difference between several segments, which break the image scene into a mosaic of patches.

Defining the segments on an image series enabled us to detect inundated areas after flood disasters and compare changes in complex channel and lagoon systems within the area of the Sudd marshes (Figure 8). The areas covered by water are distinct from the neighbour regions and only include pixels with corresponding spectral reflectance for water, which shows a contrast between the red and near-infrared (NIR) areas. The segments show the location of the Bahr al Jabal flow and its floodplain with distributed small tributaries separated based on the values of pixels, which are different from the forest land cover types or those covered by agricultural vegetation. The consequences of the flood events are visible on the corresponding image scenes by an analysis of land cover changes. Thus, the changed pattern of the postflood landscapes as a consequence of severe floods in July 2014 triggered by seasonal rainfalls is visible in the image in Figure 8a. Social consequences of such disasters include worsened living conditions and the displacement of 68% of 1.3 M people [124].

Winter months after the annual flood peak are characterised by the saturated soils with the highest moisture level, which cannot retain additional water due to the reached saturation level. Therefore, even minimal rainfall triggers further catastrophic flooding and
results in chained consequences. Thus, occasional rainfalls contribute to the high level of inundated lands (Figure 8b-d). Flash floods and heavy rains between June and November 2018 affected over 142,000 people, with damaged households and livelihoods [125]. The consequences of these disasters with flooded grassland ecosystems and inundated areas (bright green areas) are visible in the image taken in early 2019, shown in Figure 8e. Abnormally heavy seasonal flooding in South Sudan in July 2019 devastated large areas of the Sudd and surrounding areas of the White Nile tributaries, swamps, and lakes [126]. The postflood image in early 2020 shows inundated areas and landscapes (Figure 8f).

<figure><img src="image_12.png" /><figcaption>**Figure 8.** Classification maps of the satellite images of the Sudd wetlands, South Sudan based on the time series of the Landsat 8-9 OLI/TIRS images (2015–2023). The classification shows the categories of land cover classes based on the maximum-likelihood algorithm using GRASS GIS.</figcaption></figure>

The increased in flooded areas (bright green areas; Figure 8g) resulted from the worst floods in Sudan in August 2020, which were reported by the United Nations Office for the Coordination of Humanitarian Affairs (UN OCHA) [127]. The social consequences of such events are around 600,000 people being affected by disasters along the White Nile since July
2020, and widespread flooding continued in until autumn 2020. The largest affected state is Ayod, with 150,000 affected and displaced people due to the flood's effects. The affected land cover changes are reflected in the postflood image of March 2021, shown in Figure 8g.

The NDVI maps for a time series of the images show the distribution of vegetation in the postflooded scenario, shown in Figure 9. The floods in August 2021 represent the disaster in the Sudd due to the heavy showers in southwestern South Sudan, resulting in floods and river overflow. The flooding in 2021 was particularly dire in central regions of the country along the Upper Nile region, including the Jonglei state where the Sudd area extends [128].

<figure><img src="image_13.png" /><figcaption>**Figure 9.** Normalised difference vegetation index (NDVI) computed based on the satellite images of the Sudd wetlands, South Sudan: Landsat 8-9 OLI/TIRS images (2015–2023).</figcaption></figure>

The social consequences of floods for the local population are reported accordingly [129]. It is also noted by the World Health Organization (WHO) [130] that severe floods in July 2022 resulted in evacuations and damage to one million people, and 7380 people were displaced in South Sudan. The consequences for infrastructure include blocked routes by increased flood waters, which disabled humanitarian actions. The southern part of
the country, including Bentiu and Jongleis state, were the most affected. The increased inundated areas as a postflood consequence in 2023 are visible as bright-green-coloured areas, shown in Figure 9i. To evaluate the rejection probability classes, maps of pixels classified according to confidence levels were made based on the classification of nine satellite images, shown in Figure 10.

<figure><img src="image_14.png" /><figcaption>**Figure 10.** Rejection probability classes with pixels classified according to confidence levels based on the classification of the satellite images of the Sudd wetlands, South Sudan: Landsat 8-9 OLI/TIRS images (2015–2023).</figcaption></figure>

Spectral information relevant to land cover classes is obtained from the landscape analysis of the Sudd area, while the information on flooding detected in each image is retrieved through segmentation. The land cover classes are identified during classification using the distinct colours of segments and constructing the morphological shapes of the represented objects (e.g., the flow of the Nile River) in each of the images. All the pixels encompassed inside each segment are identified iteratively by the machine and interpreted
accordingly. The information from landscape patches is used to find the variations by years and propagate inundated areas using comparative analysis. Afterwards, the segments are identified as land cover classes for each segment of the image, and the metadata of the segments are updated for each scene accordingly.

# 5. Discussion

## 5.1. *Advantages of the Tools*

This paper presents an alternative algorithm of scripts-based remote sensing data analysis for environmental modelling. A multiple-methodology approach was applied to various steps of research. Image segmentation, classification, NDVI calculation, computing land cover classes, and validation of the results with accuracy assessment were performed using the GRASS GIS using a time series of the Landsat satellite images, while the QGIS was applied for geological mapping, and the topographic map was plotted using GMT. The presented series of maps supports the evaluation of the environmental setting in the Sudd region, South Sudan. Thus, based on the satellite computations of the flooded area in the Sudd, the obtained maps estimated the extent of 10 major land cover classes and flooded areas, with notable extent especially in the years 2016, 2018, and 2020. Additionally, the maximum flood extent occurs early, with an earlier peak in the flooded extent of inundated areas, which is related to the Sudd's hydrodynamics. Thus, the expanded flooding in the wetlands of the Sudd might have caused backwater effects that affected the wetland's extent and behaviour.

In this way, the current paper contributed to monitoring the Sudd wetlands through image analysis, including segmentation and classification, with the aim to determine landscape changes over the past nine years. The important deliverables of this work include image segmentation and classification to identify the diverse land cover classes for analysis of vegetation and flooded areas based on image analysis techniques. In contrast to the former reported results [131] where low-resolution SAR imagery was used and processed using ENVISAT software for the period of 2007–2011, this study presented a scripts-based approach using an open-access high-resolution dataset, which enables the repeatability and continuation of this study. Previous works also reported an increase in peak flood area, which ranges considerably from 2007 to 2009, and classified permanent flooding, seasonal flooding, and intermittent flooding areas. Likewise, this study reports an increase in the inundated areas of the Sudd for the years 2016, 2018, and 2020.

## 5.2. *Key Deliverables*

The wetland systems of the Sudd are one of the most important ecosystems of South Sudan and are included in the list of the Ramsar Convention on Wetlands of International Importance Especially as Waterfowl Habitat due to their hydrological importance in the Nile Basin, as well as the high number of endangered and vulnerable species therein. The Sudd plays a crucial role in regulating the balance of floodwater and accumulating sediments from the Mountain Nile. Moreover, since over half of the water is evaporated in the Sudd, it serves as an important mechanism of hydrological stability in the Nile River. Therefore, the disturbed flooding system will necessarily affect the Nile Basin and involve negative environmental consequences. The decrease in wetland area and the changed hydrological regime of marshes can directly and indirectly affect the Nile Basin and thus increase the negative effects of climate change. In this regard, conservation actions focused on the Sudd wetlands support the regulation of the following climate-environmental issues: increase in temperatures, unstable precipitation patterns, unbalanced crop planting, deforestation, carbon emissions related to regional agriculture sector, etc.

The floods presented in simulated models are also related to the drought periods, which are reflected in the periods of minimum flood extent in the years 2015, 2019, and 2021. This also supports the findings from previous works on climate-environmental relationships and report on the impact of the Sudd wetlands on atmospheric moisture fluxes in South Sudan through evaporation [132]. The presented geospatial visualisation
supported sustainable monitoring of the Sudd wetlands and the mapping of the inundated areas using segmentation and classification as comparative analysis and NDVI calculation to detect vegetation areas. Furthermore, topographic and geologic data were used for the analysis of geomorphic structures, main geologic units, and provinces to show the distribution of Quaternary sediments characterised by clayey soils, which create perfect conditions for the extension of swamps and wetlands.

### 5.3. Reliability of Methods

The demonstrated results are based on Landsat 8-9 OLI/TIRS products, which continues existing studies using Landsat products for the Sudd area [133]. The images were processed and tested using GRASS GIS to contribute to the initiatives on digital environmental monitoring of African wetlands. In relevant studies [95], spatial constraints in flooded areas were exploited using data from MODIS imagery and reported a connectivity between swamps in the wetlands of the Sudd. In this regards, this study supports the analysis of changes in flooded areas of the Sudd through presenting a cartographic mapping based on the open-source remote sensing data and advanced techniques of image processing. Landscape analysis included the detection of segments corresponding to the flooded land areas from pixel-based data extraction. Creating a novel series of maps based on the segmentation and classification of the remote sensing data aims at the environmental monitoring of South Sudan and, specifically, detecting the inundated areas of the Sudd wetlands. It is furthermore intended to present novel information accessible to ecologists and environmental modellers as an information source for conservation actions and detecting vulnerable regions prone to inundation during flood periods with links to the ecology in the Sudd wetlands.

The study aimed at monitoring flooded areas of the Sudd and affected land cover types for analysis of the climate-environmental effects on the sustainability of wetlands. To this end, a series of numerical experiments and cartographic data processing using remote sensing data were performed to evaluate the changes in the Sudd wetlands in regard to the Nile's environmental setting in South Sudan. The detected variation in segments by years indicated difference in peaks of flooding, which were visualised using processing Landsat 8-9 OLI/TIRS images. The fluctuated flooded areas of the Sudd wetlands in South Sudan (dry land or filled by water) were recognised during the period of 9 years. The recognition of the pixels was based on the discrimination of spectral reflectance properties based on threshold criteria, which is an essential part of the segmentation algorithm. Thus, pixels were assigned to segments of the images that were distinct from the others, which identified the following land cover classes of the Sudd wetlands: sand, marsh, flooded areas, bare land, etc.

From a cartographic perspective, the demonstrated application and functionality of the GRASS GIS also contributes to the continuation of the environmental research focused in the Sudd ecosystems due to the open-source availability of the used tools. Thus, the GRASS GIS software was used for processing the remote sensing data and image analysis, which can be continued in similar studies using presented scripts. The demonstrated and explained cartographic tasks included the conversion of raster satellite images into the maps of segmented patches and classification of the land cover types. Technically, the algorithms for region growing and merging were employed for discrimination of various land cover classes using unique IDs in segmentation by the 'i' module. A collection of contiguous pixels that meet these criteria is merged and assigned to segments as objects. The classification was based on the 'i.maxlike' module. The input dataset included 9 satellite images in a raster TIFF format obtained from the USGS. The algorithms of the GRASS GIS were discussed in detail with comments provided on scripts, demonstrating the efficiency of this software for the processing and segmentation of satellite images.

## 6. Conclusions

This research developed links between the technical approach of cartographic data processing by GRASS GIS and environmental analysis of the Sudd area using image processing. In this way, it presents the first data-driven approach that can make its own decisions on
the variations in the flooded areas of the unique Sudd wetland system in South Sudan. The cartographic interpretation of the vegetation and inundated areas was performed using data collected in a sequence of nine years (from 2015 to 2023) for a retrospective analysis of changes in the Sudd marshes. In this respect, the research performed monitoring and mapping of the extent of floods in South Sudan using a comparison of images as a short-term time series of satellite images. The analysis of the remote sensing data and supplementary information supported the detection of the areas prone to flooding. Future similar studies may also consider the overlay of the presented maps with additional cartographic materials, as well as the use of biogeochemical and environmental data as additional information for extended research.

The use of additional data for the presented research would enable one to extend the environmental analysis and monitoring in further directions. The methods of image processing by GRASS GIS are applicable to other research areas, since the Landsat images have comparable technical characteristics standardised for the satellite products of Landsat. Furthermore, this study can be continued by using new Landsat images covering other periods. Since access to the USGS EarthExplorer repositories with Landsat data are open and freely available, the use of images for various periods can support long-term environmental monitoring of South Sudan. As a continuation of this work, these data can be reused for recommendations regarding preventive measures to avoid the risks of flooding using data on wetland boundaries, as well as the inflow and outflow of water affecting vegetation communities in the Sudd marshes. Furthermore, the communication and results dissemination to involved parties can assist with decision making regarding the extent of the inundated areas. This especially concerns fishery communities or farmers depending on the flood periods in the Sudd area.

**Funding:** This research received no external funding.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** GRASS GIS scripts for image segmentation, clustering and classification, clustering results, computed kappa parameters, and errors matrices are available in the GitHub repository: https://github.com/paulinelemenkova/Sudd_South_Sudan_Image_Analysis (accessed on 10 August 2023).

**Acknowledgments:** The author thanks the reviewers for reading and reviewing this manuscript.

**Conflicts of Interest:** The author declares no conflict of interest.

## Abbreviations

The following abbreviations are used in this manuscript:

<table><tr><td>AVHRR</td><td>Advanced Very High-Resolution Radiometer</td></tr><tr><td>CNNs</td><td>Convolutional Neural Networks</td></tr><tr><td>DCW</td><td>Digital Chart of the World</td></tr><tr><td>DEM</td><td>Digital Elevation Model</td></tr><tr><td>FAO UN</td><td>Food and Agriculture Organization of the United Nations</td></tr><tr><td>GEBCO</td><td>General Bathymetric Chart of the Oceans</td></tr><tr><td>GMT</td><td>Generic Mapping Tools</td></tr><tr><td>GRASS</td><td>Geographic Resources Analysis Support System</td></tr><tr><td>GIS</td><td>Geographic Information System</td></tr><tr><td>Landsat OLI/TIRS</td><td>Landsat Operational Land Imager and Thermal Infrared Sensor</td></tr><tr><td>NDVI</td><td>Normalized Difference Vegetation Index</td></tr><tr><td>TIFF</td><td>Tag Image File Format</td></tr><tr><td>UN OCHA</td><td>United Nations Office for the Coordination of Humanitarian Affairs</td></tr><tr><td>USGS</td><td>United States Geological Survey</td></tr><tr><td>WHO</td><td>World Health Organization</td></tr></table>
# Appendix A. GRASS GIS Scripts for Image Processing, Segmentation and Classification

## Listing A1. GRASS GIS code for importing data for the Landsat OLI/TIRS bands.

1 r.import input=/Users/polinalemenkova/grassdata/SSudan/
LC09_L2SP_173055_20230514_20230516_02_T1_SR_B1.TIF output=L9_2023_01 resample=
bilinear extent=region resolution=region --overwrite

2 r.import input=/Users/polinalemenkova/grassdata/SSudan/
LC09_L2SP_173055_20230514_20230516_02_T1_SR_B2.TIF output=L9_2023_02 extent=region
resolution=region

3 r.import input=/Users/polinalemenkova/grassdata/SSudan/
LC09_L2SP_173055_20230514_20230516_02_T1_SR_B3.TIF output=L9_2023_03 extent=region
resolution=region

4 # repeated for the rest of bands

5 r.import input=/Users/polinalemenkova/grassdata/SSudan/
LC09_L2SP_173055_20230514_20230516_02_T1_SR_B7.TIF output=L9_2023_07 extent=region
resolution=region

## Listing A2. GRASS GIS code for creating semantic labels for the Landsat OLI/TIRS.

1 # listing the available raster bands in GRASS GIS mapset

2 g.list rast

3 # obtaining information on raster metadata:

4 r.info -r L9_2023_07

5 # grouping data by i.group and set up a computational region to match the scene

6 g.region raster=L9_2023_01 -p

7 # store VIZ, NIR, MIR into group/subgroup and leaving out TIR as redundant

8 i.group group=L9_2023 subgroup30m \

9 input=L9_2023_01,L9_2023_02,L9_2023_03,L9_2023_04,L9_2023_05,L9_2023_06,L9_2023_07

10 # Set the region test area with the resolution taken from the input Landsat bands

11 g.region -p raster=L9_2023_01

## Listing A3. GRASS GIS code for segmentation for image tested with 2 levels of threshold.

1 # Threshold tested with values between > 0 and < 1

2 isegment group=L9_2023 output=segs_L9 threshold=0.90 similarity=euclidean method=
region_growing

3 minsize=100 --overwrite

4 isegment group=L9_2023 output=segs_L9_2 threshold=0.05 similarity=euclidean method=
region_growing

5 seeds=segs_L9 minsize=100 iterations=10

## Listing A4. GRASS GIS code for mapping the segmented raster image Landsat 9 OLI/TIRS.

1 d.mon wx0

2 g.region raster=segs_L9 -p

3 rcolors segs_L9 color=roygbiv -e

4 d.rast segs_L9

5 d.legend raster=segs_L9 title="2023" titlefontsize=12 font="Helvetica" fontsize=10
bgcolor=white border_color=white

6 d.out.file output=segs_L9 format=jpg --overwrite

## Listing A5. GRASS GIS code for classification of the Sudd region based on the segmented raster image Landsat 9 OLI/TIRS.

1 # Importing bands, here for Band 1, repeated for all the Landsat bands likewise

2 r.import input=/Users/polinalemenkova/grassdata/SSudan/
LC08_L2SP_173055_20150108_20200910_02_T1_SR_B1.TIF output=L8_2015_01 resample=
bilinear extent=region resolution=region --overwrite

3 # check up the imported files

4 g.list rast

5 # Defining computational region to match the scene's extent

6 g.region raster=L8_2015_01 -p

7 # creating groups from the VIZ, NIR, MIR Bands into groups/subgroups:

8 i.group group=L8_2015 subgroup=1res_30m \

9 input=L8_2015_01,L8_2015_02,L8_2015_03,L8_2015_04,L8_2015_05,L8_2015_06,L8_2015_07

10 # Clustering by k-means algorithm and generating signature file and clustering report ( presented in Appendix B)

11 i.cluster group=L8_2015 subgroup=1res_30m \

12 signaturefile=cluster_L8_2015 \

13 classes=10 reportfile=rep_clust_L8_2015.txt --overwrite

14 # Running classification

15 i.maxlik group=L8_2015 subgroup=1res_30m signaturefile=cluster_L8_2015 \

16 output=L8_2015_cluster_classes reject=L8_2015_cluster_reject
**Listing A6.** GRASS GIS code for computing the NDVI for assessment of vegetation coverage over Sudd (example for 2015).

1 i.landsat.toar input=lsat8_2015. output=lsat8_2015_toar. sensor=oli8 \
2 method=dos1 date=2015-01-08 sun_elevation=50.07334117 \
3 product_date=2015-01-08 gain=HHHLHLHHL
4 # Calculation of NDVI
5 g.region raster=lsat8_2015_toar.4 -p
6 i.vi red=lsat8_2015_toar.4 nir=lsat8_2015_toar.5 viname=ndvi \
7 output=lsat8_2015.ndvi --overwrite
8 r.colorlsat8_2015.ndvi color=ndvi
9 # displaying the map
10 d.mon wx0
11 g.region raster=lsat8_2015_toar.4 -p
12 d.rast lsat8_2015.ndvi
13 d.legend raster=lsat8_2015.ndvi range=-1,1 title="NDVI" title_fontsize=14 font=Helvetica
14 d.out.file output=SSudan_NDVI_2015 format=jpg --overwrite

**Listing A7.** GRASS GIS code for computing the error matrix and kappa parameters for accuracy assessment of Landsat classification.

1 # r.kappa - Calculates error matrix and kappa parameter for accuracy assessment of classification result.
2 g.region raster=L8_2015_cluster_classes -p
3 r.kappa -w classification=L8_2015_cluster_classes reference=training_classes_Sudd
4 # export Kappa matrix as CSV file "kappa.csv"
5 r.kappa classification=L8_2015_cluster_classes reference=training_classes_Sudd output=kappa.csv -m -h --overwrite

## Appendix B. Accuracy Assessment: Calculated Error Matrices and Kappa Parameters

<table><caption>**Table A1.** Calculated error matrix and kappa parameter for accuracy assessment of the classification results for Landsat 8 image on 2015 using 'r.kappa' module of GRASS GIS.</caption><thead><tr><th>cat#</th><th>Class 1</th><th>Class 2</th><th>Class 3</th><th>Class 4</th><th>Class 5</th><th>Class 6</th><th>Class 7</th><th>Class 8</th><th>Class 9</th><th>Class 10</th><th>RowSum</th></tr></thead><tbody><tr><td>Class 1</td><td>733,990</td><td>71,250</td><td>51,480</td><td>7839</td><td>46,757</td><td>700,709</td><td>113,227</td><td>113,324</td><td>9887</td><td>5118</td><td>1,853,581</td></tr><tr><td>Class 2</td><td>26,060</td><td>6727</td><td>35,538</td><td>26,438</td><td>364,431</td><td>2,554,250</td><td>860,056</td><td>672,721</td><td>31,768</td><td>28,371</td><td>4,606,360</td></tr><tr><td>Class 3</td><td>126,808</td><td>1,031,318</td><td>265,035</td><td>526,075</td><td>315,900</td><td>19,612</td><td>264,441</td><td>41,069</td><td>252,749</td><td>7260</td><td>2,850,267</td></tr><tr><td>Class 4</td><td>49,698</td><td>1,019,916</td><td>414,261</td><td>1,502,475</td><td>729,561</td><td>151,428</td><td>544,346</td><td>56,599</td><td>699,290</td><td>36,435</td><td>5,204,009</td></tr><tr><td>Class 5</td><td>55,812</td><td>124,601</td><td>65,706</td><td>587,382</td><td>1,069,099</td><td>2,306,196</td><td>1,927,058</td><td>744,408</td><td>55,621</td><td>67,236</td><td>7,003,119</td></tr><tr><td>Class 6</td><td>52,480</td><td>54,546</td><td>5866</td><td>311,905</td><td>743,280</td><td>2,588,755</td><td>1,866,894</td><td>982,650</td><td>119,871</td><td>76,889</td><td>6,803,136</td></tr><tr><td>Class 7</td><td>24,987</td><td>458,562</td><td>337,404</td><td>1,432,717</td><td>276,379</td><td>72,536</td><td>202,169</td><td>59,283</td><td>1,294,871</td><td>6729</td><td>4,165,637</td></tr><tr><td>Class 8</td><td>85,044</td><td>497,573</td><td>2,216,937</td><td>118,633</td><td>30,004</td><td>449</td><td>10,108</td><td>1598</td><td>90,003</td><td>1007</td><td>3,051,356</td></tr><tr><td>Class 9</td><td>5751</td><td>10,876</td><td>1737</td><td>36,475</td><td>89,652</td><td>424,421</td><td>705,550</td><td>1,286,583</td><td>70,834</td><td>17,102</td><td>2,648,981</td></tr><tr><td>Class 10</td><td>19,907</td><td>13,236</td><td>51,275</td><td>160,668</td><td>121,074</td><td>40,570</td><td>76,701</td><td>87,681</td><td>861,525</td><td>291</td><td>1,432,928</td></tr><tr><td>ColSum</td><td>1,180,537</td><td>3,288,605</td><td>3,445,239</td><td>4,710,607</td><td>3,786,137</td><td>8,858,926</td><td>6,570,550</td><td>4,045,916</td><td>3,486,419</td><td>246,438</td><td>39,619,374</td></tr></tbody></table>

<table><caption>**Table A2.** Calculated error matrix and kappa parameter for accuracy assessment of the classification results for Landsat 8 image on 2016 using 'r.kappa' module of GRASS GIS.</caption><thead><tr><th>cat#</th><th>Class 1</th><th>Class 2</th><th>Class 3</th><th>Class 4</th><th>Class 5</th><th>Class 6</th><th>Class 7</th><th>Class 8</th><th>Class 9</th><th>Class 10</th><th>RowSum</th></tr></thead><tbody><tr><td>Class 1</td><td>780,231</td><td>220,046</td><td>154,297</td><td>83,530</td><td>30,508</td><td>5503</td><td>20,356</td><td>2895</td><td>47,941</td><td>1275</td><td>1,346,582</td></tr><tr><td>Class 2</td><td>99,741</td><td>1,580,160</td><td>638,455</td><td>893,432</td><td>63,464</td><td>1502</td><td>16,019</td><td>1711</td><td>347,409</td><td>2430</td><td>3,644,323</td></tr><tr><td>Class 3</td><td>45,287</td><td>105,436</td><td>146,747</td><td>352,943</td><td>857,966</td><td>1,395,889</td><td>842,830</td><td>283,690</td><td>379,210</td><td>41,515</td><td>4,451,513</td></tr><tr><td>Class 4</td><td>35,525</td><td>99,491</td><td>19,632</td><td>174,655</td><td>357,086</td><td>3,414,838</td><td>689,202</td><td>854,448</td><td>255,453</td><td>45,549</td><td>5,945,879</td></tr><tr><td>Class 5</td><td>62,103</td><td>104,822</td><td>101,515</td><td>593,389</td><td>1,081,257</td><td>1,214,029</td><td>1,864,351</td><td>467,083</td><td>179,258</td><td>81,250</td><td>5,749,057</td></tr><tr><td>Class 6</td><td>78,202</td><td>734,566</td><td>2,116,828</td><td>415,473</td><td>20,173</td><td>166</td><td>6818</td><td>564</td><td>101,522</td><td>2036</td><td>3,476,348</td></tr><tr><td>Class 7</td><td>34,081</td><td>325,250</td><td>183,816</td><td>1,522,681</td><td>609,621</td><td>148,043</td><td>895,476</td><td>76,289</td><td>838,978</td><td>29,955</td><td>4,664,190</td></tr><tr><td>Class 8</td><td>20,087</td><td>100,341</td><td>47,817</td><td>255,331</td><td>336,311</td><td>2,069,073</td><td>1,207,148</td><td>1,054,234</td><td>176,795</td><td>22,133</td><td>5,289,270</td></tr><tr><td>Class 9</td><td>26,399</td><td>62,204</td><td>48,815</td><td>358,068</td><td>358,608</td><td>529,045</td><td>819,897</td><td>518,414</td><td>845,407</td><td>14,161</td><td>3,581,018</td></tr><tr><td>Class 10</td><td>4591</td><td>43,312</td><td>37,546</td><td>143,710</td><td>82,955</td><td>93,770</td><td>258,290</td><td>809,648</td><td>410,504</td><td>6326</td><td>1,890,652</td></tr><tr><td>ColSum</td><td>1,186,247</td><td>3,375,628</td><td>3,495,468</td><td>4,793,212</td><td>3,797,949</td><td>8,871,858</td><td>6,620,387</td><td>4,068,976</td><td>3,582,477</td><td>246,630</td><td>40,038,832</td></tr></tbody></table>
<table><caption>**Table A3.** Calculated error matrix and kappa parameter for accuracy assessment of the classification results for Landsat 8 image on 2017 using 'r.kappa' module of GRASS GIS.</caption><thead><tr><th>cat#</th><th>Class 1</th><th>Class 2</th><th>Class 3</th><th>Class 4</th><th>Class 5</th><th>Class 6</th><th>Class 7</th><th>Class 8</th><th>Class 9</th><th>Class 10</th><th>RowSum</th></tr></thead><tbody><tr><td>Class 1</td><td>714,689</td><td>117,487</td><td>52,980</td><td>29,136</td><td>14,787</td><td>67,746</td><td>17,701</td><td>9756</td><td>14,785</td><td>723</td><td>1,039,790</td></tr><tr><td>Class 2</td><td>29,352</td><td>89,841</td><td>42,227</td><td>189,497</td><td>504,534</td><td>2,227,043</td><td>827,928</td><td>594,172</td><td>32,439</td><td>44,234</td><td>4,581,267</td></tr><tr><td>Class 3</td><td>132,797</td><td>1,146,691</td><td>330,041</td><td>822,072</td><td>207,013</td><td>12,648</td><td>202,651</td><td>11,771</td><td>1,125,908</td><td>7090</td><td>3,998,682</td></tr><tr><td>Class 4</td><td>41,460</td><td>802,806</td><td>164,039</td><td>1,157,656</td><td>1,081,506</td><td>581,552</td><td>1,024,952</td><td>173,449</td><td>519,751</td><td>25,839</td><td>5,573,010</td></tr><tr><td>Class 5</td><td>31,282</td><td>101,832</td><td>47,418</td><td>331,621</td><td>491,478</td><td>1,274,293</td><td>1,413,071</td><td>693,162</td><td>42,470</td><td>39,397</td><td>4,466,024</td></tr><tr><td>Class 6</td><td>94,598</td><td>767,242</td><td>1,889,767</td><td>803,189</td><td>138,593</td><td>5507</td><td>58,560</td><td>6832</td><td>831,596</td><td>3347</td><td>4,599,231</td></tr><tr><td>Class 7</td><td>60,993</td><td>201,077</td><td>890,179</td><td>350,265</td><td>126,975</td><td>21,032</td><td>61,723</td><td>47,984</td><td>720,003</td><td>432</td><td>2,480,663</td></tr><tr><td>Class 8</td><td>56,736</td><td>177,652</td><td>99,278</td><td>980,998</td><td>910,525</td><td>1,832,830</td><td>1,627,262</td><td>653,205</td><td>200,345</td><td>67,755</td><td>6,606,586</td></tr><tr><td>Class 9</td><td>28,932</td><td>33,432</td><td>11,456</td><td>142,960</td><td>309,440</td><td>2,717,290</td><td>1,130,225</td><td>1,143,763</td><td>76,522</td><td>52,290</td><td>5,646,310</td></tr><tr><td>Class 10</td><td>4903</td><td>29,081</td><td>7521</td><td>51,253</td><td>49,527</td><td>129,146</td><td>315,224</td><td>775,406</td><td>94,041</td><td>6616</td><td>1,462,718</td></tr><tr><td>ColSum</td><td>1,195,742</td><td>3,467,141</td><td>3,534,906</td><td>4,858,647</td><td>3,834,378</td><td>8,869,087</td><td>6,679,297</td><td>4,109,500</td><td>3,657,860</td><td>247,723</td><td>4,045,4281</td></tr></tbody></table>

<table><caption>**Table A4.** Calculated error matrix and kappa parameter for accuracy assessment of the classification results for Landsat 8 image on 2018 using 'r.kappa' module of GRASS GIS.</caption><thead><tr><th>cat#</th><th>Class 1</th><th>Class 2</th><th>Class 3</th><th>Class 4</th><th>Class 5</th><th>Class 6</th><th>Class 7</th><th>Class 8</th><th>Class 9</th><th>Class 10</th><th>RowSum</th></tr></thead><tbody><tr><td>Class 1</td><td>732,406</td><td>145,516</td><td>55,377</td><td>18,756</td><td>2042</td><td>567</td><td>2028</td><td>265</td><td>3889</td><td>540</td><td>961,386</td></tr><tr><td>Class 2</td><td>42,867</td><td>83,461</td><td>39,794</td><td>120,546</td><td>646,654</td><td>1,668,545</td><td>60,5985</td><td>381,538</td><td>67,729</td><td>44,928</td><td>3,702,047</td></tr><tr><td>Class 3</td><td>50,806</td><td>72,438</td><td>243,34</td><td>192,753</td><td>550,591</td><td>3,041,550</td><td>1,180,439</td><td>877,680</td><td>45,519</td><td>67,934</td><td>6,104,044</td></tr><tr><td>Class 4</td><td>39,524</td><td>549,507</td><td>130,403</td><td>1,162,081</td><td>956,566</td><td>163,501</td><td>1,005,861</td><td>76,938</td><td>592,355</td><td>40,140</td><td>4,716,876</td></tr><tr><td>Class 5</td><td>19,153</td><td>74,871</td><td>29,427</td><td>184,968</td><td>292,231</td><td>1,738,221</td><td>1,323,091</td><td>903,817</td><td>42,208</td><td>24,979</td><td>4,632,966</td></tr><tr><td>Class 6</td><td>103,865</td><td>1,380,744</td><td>464,199</td><td>991,808</td><td>155,944</td><td>2218</td><td>115,542</td><td>2619</td><td>675,785</td><td>4254</td><td>3,896,978</td></tr><tr><td>Class 7</td><td>108,450</td><td>727,966</td><td>2,543,054</td><td>684,347</td><td>22,111</td><td>767</td><td>13,922</td><td>755</td><td>586,134</td><td>1725</td><td>4,689,231</td></tr><tr><td>Class 8</td><td>40,667</td><td>144,063</td><td>127,960</td><td>1,069,080</td><td>804,237</td><td>591,309</td><td>1,066,923</td><td>232,496</td><td>725,319</td><td>28,954</td><td>4,831,008</td></tr><tr><td>Class 9</td><td>19,387</td><td>60,859</td><td>25,457</td><td>154,727</td><td>238,558</td><td>1,584,571</td><td>993,358</td><td>960,600</td><td>88,933</td><td>25,916</td><td>4,152,366</td></tr><tr><td>Class 10</td><td>24,735</td><td>80,861</td><td>21,160</td><td>165,180</td><td>115,558</td><td>65,516</td><td>258,729</td><td>595,195</td><td>693,860</td><td>6799</td><td>2,027,593</td></tr><tr><td>ColSum</td><td>1,181,860</td><td>3,320,286</td><td>3,461,165</td><td>4,744,246</td><td>3,784,492</td><td>8,856,765</td><td>6,565,878</td><td>4,031,903</td><td>3,521,731</td><td>246,169</td><td>39,714,495</td></tr></tbody></table>

<table><caption>**Table A5.** Calculated error matrix and kappa parameter for accuracy assessment of the classification results for Landsat 8 image on 2019 using 'r.kappa' module of GRASS GIS.</caption><thead><tr><th>cat#</th><th>Class 1</th><th>Class 2</th><th>Class 3</th><th>Class 4</th><th>Class 5</th><th>Class 6</th><th>Class 7</th><th>Class 8</th><th>Class 9</th><th>Class 10</th><th>RowSum</th></tr></thead><tbody><tr><td>Class 1</td><td>718,640</td><td>203302</td><td>172,239</td><td>93,576</td><td>43,600</td><td>4256</td><td>18,785</td><td>1009</td><td>154,123</td><td>961</td><td>1,410,491</td></tr><tr><td>Class 2</td><td>125,598</td><td>1,493,295</td><td>662,957</td><td>875,909</td><td>62,791</td><td>233</td><td>12,810</td><td>1003</td><td>199,881</td><td>2890</td><td>3,437,367</td></tr><tr><td>Class 3</td><td>54,049</td><td>77,429</td><td>71,198</td><td>263,136</td><td>915,884</td><td>1,601,573</td><td>838,669</td><td>292,464</td><td>542,080</td><td>38,583</td><td>4,695,065</td></tr><tr><td>Class 4</td><td>54,698</td><td>420,455</td><td>224,468</td><td>1,593,385</td><td>860,119</td><td>197,570</td><td>1,208,920</td><td>145,740</td><td>498,281</td><td>68,378</td><td>5,272,014</td></tr><tr><td>Class 5</td><td>100,175</td><td>681,102</td><td>2,082,822</td><td>353,990</td><td>9774</td><td>12</td><td>1168</td><td>33</td><td>55,344</td><td>1488</td><td>3,285,908</td></tr><tr><td>Class 6</td><td>62,037</td><td>74,841</td><td>42,628</td><td>229,455</td><td>708,692</td><td>3,498,641</td><td>1,282,209</td><td>997,764</td><td>207,553</td><td>73,493</td><td>7,177,313</td></tr><tr><td>Class 7</td><td>27,615</td><td>96,635</td><td>47,774</td><td>265,935</td><td>444,419</td><td>2,280,163</td><td>1,435,551</td><td>1,089,070</td><td>182,786</td><td>27,049</td><td>5,896,997</td></tr><tr><td>Class 8</td><td>22,403</td><td>103,576</td><td>69,604</td><td>564,151</td><td>491,174</td><td>491,554</td><td>1,015,617</td><td>231,168</td><td>782,133</td><td>20,739</td><td>3,792,119</td></tr><tr><td>Class 9</td><td>7734</td><td>84,827</td><td>34,720</td><td>226,934</td><td>165,732</td><td>728,783</td><td>621,082</td><td>968,785</td><td>172,290</td><td>5761</td><td>3,016,648</td></tr><tr><td>Class 10</td><td>8316</td><td>75,416</td><td>46,918</td><td>267,452</td><td>80,567</td><td>52,556</td><td>124,414</td><td>300,513</td><td>716,496</td><td>6747</td><td>1,679,395</td></tr><tr><td>ColSum</td><td>1,181,265</td><td>3,310,878</td><td>3,455,328</td><td>4,733,923</td><td>3,782,752</td><td>8,855,341</td><td>6,559,225</td><td>4,027,549</td><td>3,510,967</td><td>246,089</td><td>39,663,317</td></tr></tbody></table>
<table><caption>**Table A6.** Calculated error matrix and kappa parameter for accuracy assessment of the classification results for Landsat 8 image on 2020 using 'r.kappa' module of GRASS GIS.</caption><thead><tr><th>cat#</th><th>Class 1</th><th>Class 2</th><th>Class 3</th><th>Class 4</th><th>Class 5</th><th>Class 6</th><th>Class 7</th><th>Class 8</th><th>Class 9</th><th>Class 10</th><th>RowSum</th></tr></thead><tbody><tr><td>Class 1</td><td>683,544</td><td>309,023</td><td>99,294</td><td>155,266</td><td>51,772</td><td>166,487</td><td>46,734</td><td>33,393</td><td>90,153</td><td>1763</td><td>1,637,429</td></tr><tr><td>Class 2</td><td>140,068</td><td>1,405,103</td><td>502,303</td><td>960,693</td><td>301,650</td><td>1969</td><td>244,355</td><td>7645</td><td>463,950</td><td>15,581</td><td>4,043,317</td></tr><tr><td>Class 3</td><td>34,985</td><td>29,270</td><td>6265</td><td>56,130</td><td>239,734</td><td>1,980,369</td><td>393,702</td><td>256,511</td><td>33,095</td><td>39,034</td><td>3,069,095</td></tr><tr><td>Class 4</td><td>37,825</td><td>13,154</td><td>4829</td><td>51,406</td><td>404,025</td><td>3,157,656</td><td>857,426</td><td>692,151</td><td>43,581</td><td>53,723</td><td>5,315,776</td></tr><tr><td>Class 5</td><td>138,477</td><td>671,024</td><td>2,455,571</td><td>568,894</td><td>34,995</td><td>1184</td><td>15,184</td><td>5799</td><td>152,695</td><td>3956</td><td>4,047,779</td></tr><tr><td>Class 6</td><td>51,528</td><td>717,269</td><td>267,346</td><td>1,543,446</td><td>1,005,959</td><td>256,180</td><td>637,860</td><td>126,776</td><td>660,357</td><td>42,351</td><td>5,309,072</td></tr><tr><td>Class 7</td><td>18,831</td><td>12,591</td><td>2094</td><td>90,961</td><td>386,704</td><td>1,777,537</td><td>1,583,559</td><td>956,054</td><td>54,353</td><td>31,397</td><td>4,914,081</td></tr><tr><td>Class 8</td><td>44,343</td><td>134,707</td><td>147,092</td><td>1,003,947</td><td>1,035,822</td><td>832,105</td><td>1,440,973</td><td>404,912</td><td>853,895</td><td>42,981</td><td>5,940,777</td></tr><tr><td>Class 9</td><td>19,763</td><td>43,411</td><td>10,566</td><td>249,256</td><td>272,111</td><td>597,942</td><td>1,194,978</td><td>989,956</td><td>638,501</td><td>9637</td><td>4,026,121</td></tr><tr><td>Class 10</td><td>18,842</td><td>61,878</td><td>10,717</td><td>128,166</td><td>71,966</td><td>106,674</td><td>222,352</td><td>605,939</td><td>610,029</td><td>6357</td><td>1,842,920</td></tr><tr><td>ColSum</td><td>1,188,206</td><td>3,397,430</td><td>3,506,077</td><td>4,808,165</td><td>3,804,738</td><td>8,878,103</td><td>6,637,123</td><td>4,079,136</td><td>3,600,609</td><td>246,780</td><td>40,146,367</td></tr></tbody></table>

<table><caption>**Table A7.** Calculated error matrix and kappa parameter for accuracy assessment of the classification results for Landsat 8 image on 2021 using 'r.kappa' module of GRASS GIS.</caption><thead><tr><th>cat#</th><th>Class 1</th><th>Class 2</th><th>Class 3</th><th>Class 4</th><th>Class 5</th><th>Class 6</th><th>Class 7</th><th>Class 8</th><th>Class 9</th><th>Class 10</th><th>RowSum</th></tr></thead><tbody><tr><td>Class 1</td><td>702,916</td><td>170,016</td><td>144,023</td><td>121,055</td><td>160,622</td><td>146,539</td><td>319,847</td><td>201,335</td><td>17,125</td><td>11,327</td><td>1,994,805</td></tr><tr><td>Class 2</td><td>137,245</td><td>772,838</td><td>405,525</td><td>861,919</td><td>482,726</td><td>297,935</td><td>637,744</td><td>319,340</td><td>133,967</td><td>40,635</td><td>4,089,874</td></tr><tr><td>Class 3</td><td>100,084</td><td>553,381</td><td>1,979,310</td><td>375,399</td><td>32,820</td><td>15,876</td><td>19,359</td><td>17,344</td><td>62,178</td><td>772</td><td>3,156,523</td></tr><tr><td>Class 4</td><td>52,869</td><td>1,280,901</td><td>736,804</td><td>1,341,217</td><td>227,239</td><td>203,122</td><td>179,091</td><td>173,844</td><td>332,915</td><td>21,818</td><td>4,549,820</td></tr><tr><td>Class 5</td><td>31,298</td><td>348,878</td><td>82,337</td><td>988,002</td><td>669,137</td><td>367,343</td><td>590,315</td><td>179,041</td><td>732,481</td><td>23,781</td><td>4,012,613</td></tr><tr><td>Class 6</td><td>82,156</td><td>20,454</td><td>2383</td><td>73,286</td><td>843,448</td><td>2,928,722</td><td>1,185,373</td><td>683,802</td><td>52,542</td><td>112,618</td><td>5,984,784</td></tr><tr><td>Class 7</td><td>14,080</td><td>2677</td><td>12</td><td>10,134</td><td>199,861</td><td>3,664,407</td><td>1,106,320</td><td>869,961</td><td>27,095</td><td>16,303</td><td>5,910,850</td></tr><tr><td>Class 8</td><td>11,133</td><td>3988</td><td>2954</td><td>140,173</td><td>514,342</td><td>987,416</td><td>1,809,103</td><td>776,155</td><td>79,990</td><td>9808</td><td>4,335,062</td></tr><tr><td>Class 9</td><td>53,852</td><td>247,468</td><td>154,669</td><td>881,230</td><td>602,744</td><td>79,115</td><td>361,528</td><td>69,245</td><td>1,863,006</td><td>6719</td><td>4,319,576</td></tr><tr><td>Class 10</td><td>2991</td><td>1500</td><td>291</td><td>19,321</td><td>73,820</td><td>188,976</td><td>432,371</td><td>791,227</td><td>303,449</td><td>3076</td><td>1,817,022</td></tr><tr><td>ColSum</td><td>1,188,624</td><td>3,402,101</td><td>3,508,308</td><td>4,811,736</td><td>3,806,759</td><td>8,879,451</td><td>6,641,051</td><td>4,081,294</td><td>3,604,748</td><td>246,857</td><td>40,170,929</td></tr></tbody></table>

<table><caption>**Table A8.** Calculated error matrix and kappa parameter for accuracy assessment of the classification results for Landsat 8 image on 2022 using 'r.kappa' module of GRASS GIS.</caption><thead><tr><th>cat#</th><th>Class 1</th><th>Class 2</th><th>Class 3</th><th>Class 4</th><th>Class 5</th><th>Class 6</th><th>Class 7</th><th>Class 8</th><th>Class 9</th><th>Class 10</th><th>RowSum</th></tr></thead><tbody><tr><td>Class 1</td><td>731,146</td><td>290,319</td><td>237,623</td><td>291,223</td><td>132,078</td><td>263,431</td><td>319,057</td><td>187,293</td><td>56,788</td><td>22,701</td><td>2,531,659</td></tr><tr><td>Class 2</td><td>35,991</td><td>11,448</td><td>5617</td><td>16,586</td><td>180,690</td><td>2,141,375</td><td>456,161</td><td>458,916</td><td>27,192</td><td>49,483</td><td>3,383,459</td></tr><tr><td>Class 3</td><td>18,712</td><td>5274</td><td>1374</td><td>21,109</td><td>225,586</td><td>2,424,597</td><td>624,656</td><td>524,088</td><td>18,643</td><td>21,329</td><td>3,885,368</td></tr><tr><td>Class 4</td><td>130,449</td><td>1,176,764</td><td>483,567</td><td>1,149,433</td><td>546,547</td><td>352,960</td><td>717,444</td><td>317,804</td><td>486,570</td><td>39,343</td><td>5,400,881</td></tr><tr><td>Class 5</td><td>30,884</td><td>35,859</td><td>16,397</td><td>213,410</td><td>535,419</td><td>1,358,447</td><td>1,173,149</td><td>582,758</td><td>40,234</td><td>37,240</td><td>4,023,797</td></tr><tr><td>Class 6</td><td>26,411</td><td>18,092</td><td>8246</td><td>150,852</td><td>504,812</td><td>1,178,353</td><td>1,640,165</td><td>765,099</td><td>49,838</td><td>32,619</td><td>4,374,487</td></tr><tr><td>Class 7</td><td>53,164</td><td>1,017,166</td><td>636,794</td><td>1,598,491</td><td>736,672</td><td>375,514</td><td>548,660</td><td>219,776</td><td>849,537</td><td>30,208</td><td>6,065,982</td></tr><tr><td>Class 8</td><td>108,133</td><td>594,918</td><td>1,948,383</td><td>348,009</td><td>32,300</td><td>59,628</td><td>18,387</td><td>25,360</td><td>138,083</td><td>1559</td><td>3,274,760</td></tr><tr><td>Class 9</td><td>24,740</td><td>168,328</td><td>93,599</td><td>521,826</td><td>575,816</td><td>606,032</td><td>927,642</td><td>694,886</td><td>707,640</td><td>9358</td><td>4,329,867</td></tr><tr><td>Class 10</td><td>32,749</td><td>110,613</td><td>91,007</td><td>522,968</td><td>352,275</td><td>131,057</td><td>246,218</td><td>329,246</td><td>1,254,761</td><td>3620</td><td>3,074,514</td></tr><tr><td>ColSum</td><td>1,192,379</td><td>3,428,781</td><td>3,522,607</td><td>4,833,907</td><td>3,822,195</td><td>8,891,394</td><td>6,671,539</td><td>4,105,226</td><td>3,629,286</td><td>247,460</td><td>40,344,774</td></tr></tbody></table>
<table><caption>**Table A9.** Calculated error matrix and kappa parameter for accuracy assessment of the classification results for Landsat 8 image on 2023 using 'r.kappa' module of GRASS GIS.</caption><thead><tr><th>cat#</th><th>Class 1</th><th>Class 2</th><th>Class 3</th><th>Class 4</th><th>Class 5</th><th>Class 6</th><th>Class 7</th><th>Class 8</th><th>Class 9</th><th>Class 10</th><th>RowSum</th></tr></thead><tbody><tr><td>Class 1</td><td>696,181</td><td>265,690</td><td>210,677</td><td>147,949</td><td>47,848</td><td>2351</td><td>28,532</td><td>11,048</td><td>318,677</td><td>916</td><td>1,729,869</td></tr><tr><td>Class 2</td><td>149,373</td><td>875,176</td><td>345,372</td><td>679,954</td><td>298,325</td><td>17,823</td><td>276,754</td><td>36,126</td><td>491,199</td><td>13,422</td><td>3,183,524</td></tr><tr><td>Class 3</td><td>68,987</td><td>1,103,574</td><td>616,225</td><td>1,449,488</td><td>384,491</td><td>27,522</td><td>383,103</td><td>58,688</td><td>592,878</td><td>31,005</td><td>4,715,961</td></tr><tr><td>Class 4</td><td>113,808</td><td>718,554</td><td>1,928,287</td><td>399,929</td><td>82,193</td><td>29,429</td><td>53,107</td><td>31,572</td><td>253,169</td><td>1602</td><td>3,611,650</td></tr><tr><td>Class 5</td><td>36,422</td><td>248,019</td><td>189,972</td><td>890,162</td><td>567,087</td><td>292,147</td><td>1,021,289</td><td>370,809</td><td>523,327</td><td>29,114</td><td>4,168,348</td></tr><tr><td>Class 6</td><td>44,185</td><td>40,446</td><td>42,802</td><td>309,544</td><td>969,941</td><td>1,256,540</td><td>1,076,897</td><td>398,925</td><td>377,458</td><td>54,559</td><td>4,571,297</td></tr><tr><td>Class 7</td><td>51,650</td><td>18,914</td><td>13,974</td><td>115,112</td><td>323,824</td><td>4,387,107</td><td>966,427</td><td>899,834</td><td>55,327</td><td>72,308</td><td>6,904,477</td></tr><tr><td>Class 8</td><td>17,292</td><td>54,481</td><td>80,843</td><td>400,775</td><td>627,113</td><td>1,343,379</td><td>1,726,289</td><td>675,519</td><td>195,048</td><td>23,051</td><td>5,143,790</td></tr><tr><td>Class 9</td><td>9240</td><td>5475</td><td>4241</td><td>86,789</td><td>210,760</td><td>1,163,138</td><td>817,121</td><td>1,116,546</td><td>54,160</td><td>12,431</td><td>3,479,901</td></tr><tr><td>Class 10</td><td>7212</td><td>129,323</td><td>102,017</td><td>367,368</td><td>311,372</td><td>314,712</td><td>311,368</td><td>501,372</td><td>792,399</td><td>8292</td><td>2,845,435</td></tr><tr><td>ColSum</td><td>1,194,350</td><td>3,459,652</td><td>3,534,410</td><td>4,847,070</td><td>3,822,954</td><td>8,834,148</td><td>6,660,887</td><td>4,100,439</td><td>3,653,642</td><td>246,700</td><td>40,354,252</td></tr></tbody></table>

## Appendix C. Clustering Report of GRASS GIS: Calculated for the Landsat Image

The example is given for the scene on 2016. All the other reports are provided in the author's GitHub repository along with GRASS GIS scripts: https://github.com/paulinelemenkova/Sudd_South_Sudan_Image_Analysis (accessed on 10 August 2023).

########## CLUSTER (Sun Jul 2 13:35:38 2023) ##########

Location: SSudan

Mapset: PERMANENT

Group: L8_2016

Subgroup: res_30m

L8_2016_01@PERMANENT

L8_2016_02@PERMANENT

L8_2016_03@PERMANENT

L8_2016_04@PERMANENT

L8_2016_05@PERMANENT

L8_2016_06@PERMANENT

L8_2016_07@PERMANENT

Result signature file: cluster_L8_2016

Region

North: 915615.00 East: 416415.00

South: 682785.00 West: 189555.00

Res: 30.00 Res: 30.00

Rows: 7761 Cells: 58688682

Mask: no

Cluster parameters

Nombre de classes initiales: 10

Minimum class size: 17

Minimum class separation: 0.000000

Percent convergence: 98.000000

Maximum number of iterations: 30

Row sampling interval: 77

Col sampling interval: 75

Sample size: 7018 points

means and standard deviations for 7 bands

moyennes 8341.81 8839.87 9796.06 10347 13521 14268 12593.9

écart-type 333.559 397.401 538.123 866.886 1982.91 1927.26 1630.37

initial means for each band
<table><tr><td>classe 1</td><td>8008.25</td><td>8442.47</td><td>9257.94</td><td>9480.1</td><td>11538.1</td><td>12340.7</td><td>10963.6</td></tr><tr><td>classe 2</td><td>8082.38</td><td>8530.78</td><td>9377.52</td><td>9672.74</td><td>11978.8</td><td>12769</td><td>11325.9</td></tr><tr><td>classe 3</td><td>8156.5</td><td>8619.09</td><td>9497.1</td><td>9865.39</td><td>12419.4</td><td>13197.3</td><td>11688.2</td></tr><tr><td>classe 4</td><td>8230.63</td><td>8707.41</td><td>9616.69</td><td>10058</td><td>12860.1</td><td>13625.5</td><td>12050.5</td></tr><tr><td>classe 5</td><td>8304.75</td><td>8795.72</td><td>9736.27</td><td>10250.7</td><td>13300.7</td><td>14053.8</td><td>12412.8</td></tr><tr><td>classe 6</td><td>8378.87</td><td>8884.03</td><td>9855.85</td><td>10443.3</td><td>13741.3</td><td>14482.1</td><td>12775.1</td></tr><tr><td>classe 7</td><td>8453</td><td>8972.34</td><td>9975.43</td><td>10636</td><td>14182</td><td>14910.4</td><td>13137.4</td></tr><tr><td>classe 8</td><td>8527.12</td><td>9060.65</td><td>10095</td><td>10828.6</td><td>14622.6</td><td>15338.7</td><td>13499.7</td></tr><tr><td>classe 9</td><td>8601.25</td><td>9148.96</td><td>10214.6</td><td>11021.2</td><td>15063.3</td><td>15766.9</td><td>13862</td></tr><tr><td>classe 10</td><td>8675.37</td><td>9237.27</td><td>10334.2</td><td>11213.9</td><td>15503.9</td><td>16195.2</td><td>14224.3</td></tr></table>

class means/stddev for each band

class 1 (742)

moyennes 7951.81 8339.41 9061.5 9178.93 11076.1 10852.5 10158.2
écart-type 257 262.216 327.764 454.974 1467.52 1336.4 1392.45

class 2 (402)

moyennes 8135.75 8577.49 9368.62 9690.99 12024.7 12474.3 11626.2
écart-type 206.535 190.022 178.166 289.903 1227.6 363.851 1008.55

class 3 (548)

moyennes 8233.58 8669.82 9470.17 9846.65 12343.6 13017.8 12042.3
écart-type 245.584 210.629 176.849 316.624 1339.17 373.941 1071.32

class 4 (767)

moyennes 8279.24 8738.97 9588.8 10030.2 12750 13501.6 12361
écart-type 242.165 238.415 231.391 373.139 1383.88 409.177 1106.05

class 5 (973)

moyennes 8313.09 8784.31 9689.54 10170 13314.3 14001.5 12545.2
écart-type 239.434 244.552 210.543 463.114 1560.73 423.344 1161.71

class 6 (1048)

moyennes 8315.49 8800.87 9775.46 10268.6 14087 14416.4 12578.4
écart-type 241.751 268.06 233.945 562.31 1841.81 557.805 1274.18

class 7 (810)

moyennes 8395.39 8911.3 9925.59 10555 14249.5 14988.4 12993.8
écart-type 226.309 252.265 238.656 532.931 1667.51 541.879 1186.87

class 8 (589)

moyennes 8451.32 9018.86 10096.9 10879.7 14397.7 15615.7 13388.1
écart-type 244.289 296.493 347.211 510.795 1292.85 513.351 1030.24

class 9 (383)

moyennes 8542.06 9143.82 10280.9 11175.9 14651.7 16187.5 13737.7
écart-type 226.471 215.073 232.586 443.034 1002.04 370.811 894.101

class 10 (756)

moyennes 8805.39 9451.85 10737.7 11805 15797.2 17600.5 14593.2
écart-type 355.91 426.665 563.559 744.099 1255.16 1086.85 1233.36

Distribution des classes

<table><tr><td>742</td><td>402</td><td>548</td><td>767</td><td>973</td></tr><tr><td>1048</td><td>810</td><td>589</td><td>383</td><td>756</td></tr></table>

###### iteration 1

10 classes, 63.02% points stable

Distribution des classes

<table><tr><td>494</td><td>665</td><td>533</td><td>840</td><td>908</td></tr><tr><td>1068</td><td>608</td><td>721</td><td>664</td><td>517</td></tr></table>

###### iteration 2

10 classes, 75.24% points stable

Distribution des classes

<table><tr><td>369</td><td>624</td><td>667</td><td>799</td><td>988</td></tr></table>
1017 765 661 709 419

#### iteration 3
10 classes, 86.09% points stable
Distribution des classes
293 598 833 757 927
833 944 761 720 352

#### iteration 4
10 classes, 91.58% points stable
Distribution des classes
249 599 869 818 947
716 943 818 747 312

#### iteration 5
10 classes, 94.69% points stable
Distribution des classes
229 622 824 874 1009
648 896 865 751 300

#### iteration 6
10 classes, 96.21% points stable
Distribution des classes
217 640 795 921 1023
604 851 930 742 295

#### iteration 7
10 classes, 97.08% points stable
Distribution des classes
210 649 770 972 1018
582 807 984 735 291

#### iteration 8
10 classes, 98.10% points stable
Distribution des classes
205 647 756 1014 994
574 779 1025 733 291

#### final results
10 classes (convergence=98.1%)

class separability matrix

<table><thead><tr><th></th><th>1</th><th>2</th><th>3</th><th>4</th><th>5</th><th>6</th><th>7</th><th>8</th><th>9</th><th>10</th></tr></thead><tbody><tr><th>1</th><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><th>2</th><td>1.3</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><th>3</th><td>1.6</td><td>1.0</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><th>4</th><td>2.6</td><td>1.8</td><td>1.1</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><th>5</th><td>2.6</td><td>1.4</td><td>1.3</td><td>0.8</td><td>0</td><td></td><td></td><td></td><td></td><td></td></tr><tr><th>6</th><td>1.9</td><td>0.8</td><td>1.6</td><td>2.3</td><td>1.8</td><td>0</td><td></td><td></td><td></td><td></td></tr><tr><th>7</th><td>2.5</td><td>1.3</td><td>1.5</td><td>1.3</td><td>0.7</td><td>1.2</td><td>0</td><td></td><td></td><td></td></tr><tr><th>8</th><td>3.2</td><td>2.2</td><td>2.1</td><td>1.2</td><td>1.0</td><td>2.3</td><td>1.1</td><td>0</td><td></td><td></td></tr><tr><th>9</th><td>3.2</td><td>2.2</td><td>2.3</td><td>1.8</td><td>1.4</td><td>2.0</td><td>1.1</td><td>0.8</td><td>0</td><td></td></tr><tr><th>10</th><td>3.6</td><td>2.8</td><td>2.9</td><td>2.4</td><td>2.2</td><td>2.7</td><td>2.0</td><td>1.5</td><td>1.0</td><td>0</td></tr></tbody></table>

class means/stddev for each band

class 1 (205)
moyennes 7792.09 8125.77 8749.8 8685.25 9720.39 9033 8554.86
écart-type 269.203 248.508 325.667 367.886 1192.41 1032.1 923.673

class 2 (647)
moyennes 7958.2 8386.62 9338.01 9492.04 13584.7 12120 10320.6
écart-type 184.172 217.176 342.491 469.583 889.819 826.243 578.126
class 3 (756)
moyennes 8203 8635.6 9330.98 9704.04 11117.4 12400.5 11990.4
écart-type 180.79 155.652 189.472 264.785 613.41 635.13 637.803

class 4 (1014)
moyennes 8474.28 8940.79 9705.86 10252.6 11769.3 13891.5 13495.2
écart-type 245.105 171.081 169.501 243.028 399.04 444.049 495.903

class 5 (994)
moyennes 8293.53 8811.29 9742.85 10406.3 13008.9 14232 12696
écart-type 152.479 155.328 213.418 325.007 427.795 473.935 481.045

class 6 (574)
moyennes 8070.27 8430.14 9509.54 9368.75 17090.2 13486.3 10600.1
écart-type 167.24 165.95 253.195 345.468 1170.35 700.286 496.695

class 7 (779)
moyennes 8270.21 8760.55 9784.77 10342.1 14674.7 14922.2 12197.4
écart-type 140.611 149.126 220.25 390.38 674.5 660.764 550.25

class 8 (1025)
moyennes 8556.69 9135.1 10148 11001.2 13308.2 15441.2 14126.3
écart-type 230.6 188.688 211.968 304.661 564.956 535.661 698.561

class 9 (733)
moyennes 8568.91 9185.04 10382 11327.4 15261.8 16677 13644.6
écart-type 254.495 334.43 432.144 536.085 844.786 635.139 713.034

class 10 (291)
moyennes 9044.36 9738.61 11135.6 12383.7 16390.7 18607.7 15523
écart-type 348.47 409.735 566.885 667.853 1254.96 936.286 1038.31

CLASSES

10 classes, 98.10% points stable

CLUSTER END (Sun Jul 2 13:35:38 2023) CLASSES

References

1. Solomon, C.; Breckon, T. Image Segmentation. In Fundamentals of Digital Image Processing; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 2010; Chapter 10, pp. 263–290. [CrossRef]

2. Li, W.; Zhao, W.; Yu, J.; Zheng, J.; He, C.; Fu, H.; Lin, D. Joint semantic–geometric learning for polygonal building segmentation from high-resolution remote sensing images. *ISPRS J. Photogramm. Remote Sens.* **2023**, *201*, 26–37. [CrossRef]

3. Dong, X.; Zhang, C.; Fang, L.; Yan, Y. A deep learning based framework for remote sensing image ground object segmentation. *Appl. Soft Comput.* **2022**, *130*, 109695. [CrossRef]

4. Wang, J.; Feng, Z.; Jiang, Y.; Yang, S.; Meng, H. Orientation Attention Network for semantic segmentation of remote sensing images. *Knowl.-Based Syst.* **2023**, *267*, 110415. [CrossRef]

5. Lemenkova, P.; Debeir, O. Satellite Image Processing by Python and R Using Landsat 9 OLI/TIRS and SRTM DEM Data on Côte d'Ivoire, West Africa. *J. Imaging* **2022**, *8*, 317. [CrossRef] [PubMed]

6. Lemenkova, P.; Debeir, O. Multispectral Satellite Image Analysis for Computing Vegetation Indices by R in the Khartoum Region of Sudan, Northeast Africa. *J. Imaging* **2023**, *9*, 98. [CrossRef]

7. Li, F.; Wong, A.; Clausi, D.A. Comparison of unsupervised segmentation methods for surficial materials mapping in Nunavut, Canada using RADARSAT-2 polarimetric, Landsat-7, and DEM data. In Proceedings of the 2014 IEEE Geoscience and Remote Sensing Symposium, Quebec City, QC, Canada, 13–18 July 2014; pp. 2727–2730. [CrossRef]

8. Bona, D.S.; Murni, A.; Mursanto, P. Semantic Segmentation And Segmentation Refinement Using Machine Learning Case Study: Water Turbidity Segmentation. In Proceedings of the 2019 IEEE International Conference on Aerospace Electronics and Remote Sensing Technology (ICARES), Yogyakarta, Indonesia, 17–18 October 2019; pp. 1–5. [CrossRef]

9. Raši, R.; Kissiyar, O.; Vollmar, M. Land cover change detection thresholds for Landsat data samples. In Proceedings of the 2011 6th International Workshop on the Analysis of Multi-Temporal Remote Sensing Images (Multi-Temp), Trento, Italy, 12–14 July 2011; pp. 205–208. [CrossRef]
10. Herlawati, H.; Handayanto, R.T.; Atika, P.D.; Sugiyatno, S.; Rasim, R.; Mugiarso, M.; Hendharsetiawan, A.A.; Jaja, J.; Purwanti, S. Semantic Segmentation of Landsat Satellite Imagery. In Proceedings of the 2022 Seventh International Conference on Informatics and Computing (ICIC), Denpasar, Bali, Indonesia, 8–9 December 2022; pp. 1–6. [CrossRef]

11. Yang, D.; Wu, D.; Ding, H. Study on land use change detection based on Landsat data with object-oriented method. In Proceedings of the 2021 International Conference on Computer Information Science and Artificial Intelligence (CISAI), Kunming, China, 17–19 September 2021; pp. 268–272. [CrossRef]

12. Tunay, M.; Marangoz, M.A.; Karakis, S.; Atesoglu, A. Detecting Urban Vegetation from Different Images Using an Object-Based Approach in Bartın, Turkey. In Proceedings of the 2007 3rd International Conference on Recent Advances in Space Technologies, Istanbul, Turkey, 14–16 June 2007; pp. 636–640. [CrossRef]

13. Xiong, Y.; Chen, Y.; Han, W.; Tong, L. A new aerosol retrieval algorithm based on statistical segmentation using Landsat-8 OLI data. In Proceedings of the 2016 IEEE International Geoscience and Remote Sensing Symposium (IGARSS), Beijing, China, 10–15 July 2016; pp. 4059–4062. [CrossRef]

14. Kochher, R.; Sharma, A. Improved principle component analysis based gray stretch algorithm for landsat image segmentation. In Proceedings of the 2016 2nd International Conference on Next Generation Computing Technologies (NGCT), Dehradun, India, 14–16 October 2016; pp. 765–771. [CrossRef]

15. Liu, Y.; Yao, L.; Xiong, W.; Zhou, Z. Fusion detection of ship targets in low resolution multi-spectral images. In Proceedings of the 2016 IEEE International Geoscience and Remote Sensing Symposium (IGARSS), Beijing, China, 10–15 July 2016; pp. 6545–6548. [CrossRef]

16. Lohmann, G. Co-occurrence-based analysis and synthesis of textures. In Proceedings of the 12th International Conference on Pattern Recognition, Jerusalem, Israel, 9–13 October 1994; Volume 1, pp. 449–453. [CrossRef]

17. Buscombe, D.; Goldstein, E.B. A Reproducible and Reusable Pipeline for Segmentation of Geoscientific Imagery. *Earth Space Sci.* **2022**, 9, e2022EA002332. [CrossRef]

18. Tzotsos, A.; Karantzalos, K.; Argialas, D. Multiscale Segmentation and Classification of Remote Sensing Imagery with Advanced Edge and Scale-Space Features. In *Scale Issues in Remote Sensing*; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 2014; Chapter 9, pp. 170–196. [CrossRef]

19. Turajlic, E.; Buza, E.; Akagic, A. Honey Badger Algorithm and Chef-based Optimization Algorithm for Multilevel Thresholding Image Segmentation. In Proceedings of the 2022 30th Telecommunications Forum (TELFOR), Belgrade, Serbia, 5–16 November 2022; pp. 1–4. [CrossRef]

20. Chen, H.; Deng, X.; Yan, L.; Ye, Z. Multilevel thresholding selection based on the fireworks algorithm for image segmentation. In Proceedings of the 2017 International Conference on Security, Pattern Analysis, and Cybernetics (SPAC), Shenzhen, China, 15–17 December 2017; pp. 175–180. [CrossRef]

21. Liu, W.; Shi, H.; Pan, S.; Huang, Y.; Wang, Y. An Improved Otsu Multi-Threshold Image Segmentation Algorithm Based on Pigeon-Inspired Optimization. In Proceedings of the 2018 11th International Congress on Image and Signal Processing, BioMedical Engineering and Informatics (CISP-BMEI), Beijing, China, 13–15 October 2018; pp. 1–5. [CrossRef]

22. Chao, J.; Xiaoxiao, Y.; Xiaohai, W. Algorithm of Double Threshold Image Segmentation Combined QGA with Two-Dimensional Otsu. In Proceedings of the 2020 5th International Conference on Mechanical, Control and Computer Engineering (ICMCCE), Harbin, China, 25–27 December 2020; pp. 2219–2223. [CrossRef]

23. Tang, Z.; Wu, Y. One image segmentation method based on Otsu and fuzzy theory seeking image segment threshold. In Proceedings of the 2011 International Conference on Electronics, Communications and Control (ICECC), Ningbo, China, 9–11 September 2011; pp. 2170–2173. [CrossRef]

24. Zhao, N.; Sui, S.K.; Kuang, P. Research on image segmentation method based on weighted threshold algorithm. In Proceedings of the 2015 12th International Computer Conference on Wavelet Active Media Technology and Information Processing (ICCWAMTIP), Chengdu, China, 18–20 December 2015; pp. 307–310. [CrossRef]

25. Barbato, M.P.; Napoletano, P.; Piccoli, F.; Schettini, R. Unsupervised segmentation of hyperspectral remote sensing images with superpixels. *Remote Sens. Appl. Soc. Environ.* **2022**, 28, 100823. [CrossRef]

26. Pal, R.; Mukhopadhyay, S.; Chakraborty, D.; Sughanthan, P.N. Very high-resolution satellite image segmentation using variable-length multi-objective genetic clustering for multi-class change detection. *J. King Saud Univ. Comput. Inf. Sci.* **2022**, 34, 9964–9976. [CrossRef]

27. Du, W.; Tian, X.; Sun, Y. A dynamic threshold edge-preserving smoothing segmentation algorithm for anterior chamber OCT images based on modified histogram. In Proceedings of the 2011 4th International Congress on Image and Signal Processing, Shanghai, China, 15–17 October 2011; Volume 2, pp. 1123–1126. [CrossRef]

28. Choi, J.; Choi, H.H.S.; Chen, M. Multi-Level Thresholding Grayscale Image Segmentation Implemented with Genetic Algorithm. In Proceedings of the 2018 IEEE MIT Undergraduate Research Technology Conference (URTC), Cambridge, MA, USA, 5–7 October 2018; pp. 1–5. [CrossRef]

29. Kaihua, W.; Tao, B. Optimal Threshold Image Segmentation Method Based on Genetic Algorithm in Wheel Set Online Measurement. In Proceedings of the 2011 Third International Conference on Measuring Technology and Mechatronics Automation, Shanghai, China, 6–7 January 2011; Volume 2, pp. 799–802. [CrossRef]
30. Chaofu, Z.; Li-Ni, M.; Lu-Na, J. Threshold infrared image segmentation based on improved genetic algorithm. In Proceedings of the IET International Conference on Information Science and Control Engineering 2012 (ICISCE 2012), Shenzhen, China, 7–9 December 2012; pp. 1–4. [CrossRef]

31. Munyati, C.; Ratshibvumo, T.; Ogola, J. Landsat TM image segmentation for delineating geological zone correlated vegetation stratification in the Kruger National Park, South Africa. *Phys. Chem. Earth Parts A/B/C* **2013**, *55–57*, 1–10. [CrossRef]

32. Wang, X.; Jing, S.; Dai, H.; Shi, A. High-resolution remote sensing images semantic segmentation using improved UNet and SegNet. *Comput. Electr. Eng.* **2023**, *108*, 108734. [CrossRef]

33. Maurya, A.; Akashdeep, Mittal, P.; Kumar, R. A modified U-net-based architecture for segmentation of satellite images on a novel dataset. *Ecol. Inform.* **2023**, *75*, 102078. [CrossRef]

34. Banerjee, B.; Varma G., S.; Buddhiraju, K.M. Satellite image segmentation: A novel adaptive mean-shift clustering based approach. In Proceedings of the 2012 IEEE International Geoscience and Remote Sensing Symposium, Munich, Germany, 22–27 July 2012; pp. 4319–4322. [CrossRef]

35. He, Y.; Sun, X.; Gao, L.; Zhang, B. Ship Detection Without Sea-Land Segmentation for Large-Scale High-Resolution Optical Satellite Images. In Proceedings of the IGARSS 2018—2018 IEEE International Geoscience and Remote Sensing Symposium, Valencia, Spain, 22–27 July 2018; pp. 717–720. [CrossRef]

36. Zhao, J.; Chen, S.; Zhao, D.; Zhu, H.; Chen, X. Unsupervised saliency detection and a-contrario based segmentation for satellite images. In Proceedings of the 2013 Seventh International Conference on Sensing Technology (ICST), Wellington, New Zealand, 3–5 December 2013; pp. 678–681. [CrossRef]

37. Wu, Q.; Gan, Y.; Lin, B.; Zhang, Q.; Chang, H. An active contour model based on fused texture features for image segmentation. *Neurocomputing* **2015**, *151*, 1133–1141. [CrossRef]

38. Ratajczak, R.; Crispim-Junior, C.F.; Faure, E.; Fervers, B.; Tougne, L. Automatic Land Cover Reconstruction From Historical Aerial Images: An Evaluation of Features Extraction and Classification Algorithms. *IEEE Trans. Image Process.* **2019**, *28*, 3357–3371. [CrossRef] [PubMed]

39. Erdem, F.; Bayram, B.; Bakirman, T.; Bayrak, O.C.; Akpinar, B. An ensemble deep learning based shoreline segmentation approach (WaterNet) from Landsat 8 OLI images. *Adv. Space Res.* **2021**, *67*, 964–974. [CrossRef]

40. Kotaridis, I.; Lazaridou, M. Integrating image segmentation in the delineation of burned areas on Sentinel-2 and Landsat 8 data. *Remote Sens. Appl. Soc. Environ.* **2023**, *30*, 100944. [CrossRef]

41. Toulouse, T.; Rossi, L.; Akhloufi, M.; Celik, T.; Maldague, X. Benchmarking of wildland fire colour segmentation algorithms. *IET Image Process.* **2015**, *9*, 1064–1072. [CrossRef]

42. Mäkelä, H.; Pekkarinen, A. Estimation of timber volume at the sample plot level by means of image segmentation and Landsat TM imagery. *Remote Sens. Environ.* **2001**, *77*, 66–75. [CrossRef]

43. Zhang, M.; Zhang, H.; Yao, B.; Lin, H.; An, X.; Liu, Y. Spatiotemporal changes of wetlands in China during 2000–2015 using Landsat imagery. *J. Hydrol.* **2023**, *621*, 129590. [CrossRef]

44. Kharma, N.; Mazhurin, A.; Saigol, K.; Sabahi, F. Adaptable image segmentation via simple pixel classification. *Comput. Intell.* **2018**, *34*, 734–762. [CrossRef]

45. Aalan Babu, A.; Mary Anita Rajam, V. Water-body segmentation from satellite images using Kapur’s entropy-based thresholding method. *Comput. Intell.* **2020**, *36*, 1242–1260. [CrossRef]

46. Awad, M.M.; Chehdi, K. Satellite image segmentation using hybrid variable genetic algorithm. *Int. J. Imaging Syst. Technol.* **2009**, *19*, 199–207. [CrossRef]

47. Saha, S.; Maulik, U. A new line symmetry distance based automatic clustering technique: Application to image segmentation. *Int. J. Imaging Syst. Technol.* **2011**, *21*, 86–100. [CrossRef]

48. A, P.; Kumar, L.S. Automatic cloud segmentation from INSAT-3D satellite image via IKM and IFCM clustering. *IET Image Process.* **2020**, *14*, 1273–1280. [CrossRef]

49. Cong, L.; Ding, S.; Wang, L.; Zhang, A.; Jia, W. Image segmentation algorithm based on superpixel clustering. *IET Image Process.* **2018**, *12*, 2030–2035. [CrossRef]

50. Vansteenkiste, E.; Gautama, S.; Philips, W. Analysing multispectral textures in very high resolution satellite images. In Proceedings of the IGARSS 2004. 2004 IEEE International Geoscience and Remote Sensing Symposium, Anchorage, AK, USA, 20–24 September 2004; Volume 5, pp. 3062–3065. [CrossRef]

51. Vansteenkiste, E.; Schoutteet, A.; Gautama, S.; Philips, W. Comparing color and textural information in very high resolution satellite image classification. In Proceedings of the 2004 International Conference on Image Processing, ICIP '04, Singapore, 24–27 October 2004; Volume 5, pp. 3351–3354. [CrossRef]

52. Zhang, J.; Cui, Y.; Lu, S.; Xiao, L. Multilayer image segmentation based on Gaussian weighted Euclidean distance and nonlinear interpolation. In Proceedings of the 2017 10th International Congress on Image and Signal Processing, BioMedical Engineering and Informatics (CISP-BMEI), Shanghai, China, 14–16 October 2017; pp. 1–5. [CrossRef]

53. Bayram, E.; Nabiyev, V. Image segmentation by using K-means clustering algorithm in Euclidean and Mahalanobis distance calculation in camouflage images. In Proceedings of the 2020 28th Signal Processing and Communications Applications Conference (SIU), Istanbul, Turkey, 5–7 October 2020; pp. 1–4. [CrossRef]
54. Selvarasu, N.; Nachiappan, A.; Nandhitha, N.M. Abnormality Detection from Medical Thermographs in Human Using Euclidean Distance Based Color Image Segmentation. In Proceedings of the 2010 International Conference on Signal Acquisition and Processing, Bangalore, India, 9–10 February 2010; pp. 73–75. [CrossRef]

55. Yamashita, H.; Sonobe, R.; Hirono, Y.; Ikka, T. Dissection of hyperspectral reflectance to estimate nitrogen and chlorophyll contents in tea leaves based on machine learning algorithms. *Sci. Rep.* **2020**, *10*, 17360. [CrossRef]

56. Oliveira, R.A.; Näsi, R.; Niemeläinen, O.; Nyholm, L.; Alhonoja, K.; Kaivosoja, J.; Jauhiainen, L.; Viljanen, N.; Nezami, S.; Markelin, L.; et al. Machine learning estimators for the quantity and quality of grass swards used for silage production using drone-based imaging spectrometry and photogrammetry. *Remote Sens. Environ.* **2020**, *246*, 111830. [CrossRef]

57. Bauer, A.; Bostrom, A.G.; Ball, J.; Applegate, C.; Cheng, T.; Laycock, S.; Rojas, S.M.; Kirwan, J.; Zhou, J. Combining computer vision and deep learning to enable ultra-scale aerial phenotyping and precision agriculture: A case study of lettuce production. *Hortic. Res.* **2019**, *6*, 70. [CrossRef]

58. Onishi, M.; Ise, T. Explainable identification and mapping of trees using UAV RGB image and deep learning. *Sci. Rep.* **2021**, *11*, 903. [CrossRef]

59. Pavani, V.; Divya, K.; Likhitha, V.V.; Mounika, G.S.; Harshitha, K.S. Image Segmentation based Imperative Feature Subset Model for Detection of Vehicle Number Plate using K Nearest Neighbor Model. In Proceedings of the 2023 Third International Conference on Artificial Intelligence and Smart Energy (ICAIS), Coimbatore, India, 2–4 February 2023; pp. 704–709. [CrossRef]

60. Zhang, J.H.; Chen, Y.J.; Kuo, Y.F.; Chen, C.Y. Fast automatic segmentation of cells and nucleuses in large-scale liquid-based monolayer smear images. In Proceedings of the 2017 International Conference on Image and Vision Computing New Zealand (IVCNZ), Christchurch, New Zealand, 4–6 December 2017; pp. 1–6. [CrossRef]

61. Mohamed, C.; Nsiri, B.; Abdelmajid, S.; Abdelghani, E.M.; Brahim, B. Deep Convolutional Networks for Image Segmentation: Application to Optic Disc detection. In Proceedings of the 2020 International Conference on Electrical and Information Technologies (ICEIT), Rabat, Morocco, 4–7 March 2020; pp. 1–3. [CrossRef]

62. Colwell, R. Remote Sensing and Spatial Information. *Nature* **1981**, *293*, 364. [CrossRef]

63. Li, C.; Balla-Arabe, S.; Yang-Song, F. Embedded Implementation of VHR Satellite Image Segmentation. In *Architecture-Aware Optimization Strategies in Real-Time Image Processing*; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 2017; Chapter 6, pp. 113–139. [CrossRef]

64. Lei, T.; Nandi, A.K. Image Segmentation for Remote Sensing Analysis. In *Image Segmentation*; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 2022; Chapter 10, pp. 229–262. [CrossRef]

65. Borcard, D.; Gillet, F.; Legendre, P. *Numerical Ecology with R*, 1st ed.; Use R!; Springer: New York, NY, USA, 2011. [CrossRef]

66. Lemenkova, P. Mapping Climate Parameters over the Territory of Botswana Using GMT and Gridded Surface Data from TerraClimate. *ISPRS Int. J. Geo-Inf.* **2022**, *11*, 473. [CrossRef]

67. Neteler, M.; Beaudette, D.E.; Cavallini, P.; Lami, L.; Cepicky, J. GRASS GIS. In *Open Source Approaches in Spatial Data Handling*; Springer: Berlin/Heidelberg, Geramny, 2008; pp. 171–199. [CrossRef]

68. Neteler, M.; Mitasova, H. *Open Source GIS—A GRASS GIS Approach*, 3rd ed.; Springer: New York, NY, USA, 2008.

69. Petersen, G.; Fohrer, N. Flooding and drying mechanisms of the seasonal Sudd flood plains along the Bahr el Jebel in southern Sudan. *Hydrol. Sci. J.* **2010**, *55*, 4–16. [CrossRef]

70. Sosnowski, A.; Ghoneim, E.; Burke, J.J.; Hines, E.; Halls, J. Remote regions, remote data: A spatial investigation of precipitation, dynamic land covers, and conflict in the Sudd wetland of South Sudan. *Appl. Geogr.* **2016**, *69*, 51–64. [CrossRef]

71. Mulatu, D.W.; Ahmed, J.; Semereab, E.; Arega, T.; Yohannes, T.; Akwany, L.O. Stakeholders, Institutional Challenges and the Valuation of Wetland Ecosystem Services in South Sudan: The Case of Machar Marshes and Sudd Wetlands. *Environ. Manag.* **2022**, *69*, 666–683. [CrossRef] [PubMed]

72. Chen, L.; Jin, Z.; Michishita, R.; Cai, J.; Yue, T.; Chen, B.; Xu, B. Dynamic monitoring of wetland cover changes using time-series remote sensing imagery. *Ecol. Inform.* **2014**, *24*, 17–26. [CrossRef]

73. Lemenkova, P. Dataset compilation by GRASS GIS for thematic mapping of Antarctica: Topographic surface, ice thickness, subglacial bed elevation and sediment thickness. *Czech Polar Rep.* **2021**, *11*, 67–85. [CrossRef]

74. Lemenkova, P.; Debeir, O. Computing Vegetation Indices from the Satellite Images Using GRASS GIS Scripts for Monitoring Mangrove Forests in the Coastal Landscapes of Niger Delta, Nigeria. *J. Mar. Sci. Eng.* **2023**, *11*, 871. [CrossRef]

75. Hofierka, J.; Mitášová, H.; Neteler, M. Chapter 17 Geomorphometry in GRASS GIS. In *Geomorphometry; Developments in Soil Science*; Hengl, T., Reuter, H.I., Eds.; Elsevier: Amsterdam, The Netherlands, 2009; Volume 33, pp. 387–410. [CrossRef]

76. Lemenkova, P.; Debeir, O. GDAL and PROJ Libraries Integrated with GRASS GIS for Terrain Modelling of the Georeferenced Raster Image. *Technologies* **2023**, *11*, 46. [CrossRef]

77. Di Vittorio, C.A.; Georgakakos, A.P. Land cover classification and wetland inundation mapping using MODIS. *Remote Sens. Environ.* **2018**, *204*, 1–17. [CrossRef]

78. Lemenkova, P.; Debeir, O. Recognizing the Wadi Fluvial Structure and Stream Network in the Qena Bend of the Nile River, Egypt, on Landsat 8-9 OLI Images. *Information* **2023**, *14*, 249. [CrossRef]

79. Campos-Taberner, M.; García-Haro, F.; Martínez, B.; Izquierdo-Verdiguier, E.; Atzberger, C.; Camps-Valls, G.; Gilabert, M.A. Understanding deep learning in land use classification based on Sentinel-2 time series. *Sci. Rep.* **2020**, *10*, 17188. [CrossRef]

80. Lemenkova, P.; Debeir, O. R Libraries for Remote Sensing Data Classification by k-means Clustering and NDVI Computation in Congo River Basin, DRC. *Appl. Sci.* **2022**, *12*, 12554. [CrossRef]
81. Mohamed, Y.; Savenije, H. Impact of climate variability on the hydrology of the Sudd wetland: Signals derived from long term (1900–2000) water balance computations. *Wetl. Ecol. Manag.* **2014**, *22*, 191–198. [CrossRef]

82. Adamson, D.; Gasse, F.; Street, F.; Williams, M.A.J. Late Quaternary history of the Nile. *Nature* **1980**, *288*, 50–55. [CrossRef]

83. Broun, A.F. Some Notes on the "Sudd"-Formation of the Upper Nile. *J. Linn. Soc. Lond. Bot.* **1905**, *37*, 51–58. [CrossRef]

84. Chorowicz, J. The East African rift system. *J. Afr. Earth Sci.* **2005**, *43*, 379–410. [CrossRef]

85. Lemenkova, P. Tanzania Craton, Serengeti Plain and Eastern Rift Valley: Mapping of geospatial data by scripting techniques. *Est. J. Earth Sci.* **2022**, *71*, 61–79. [CrossRef]

86. Petersen, G.; Fohrer, N. Two-dimensional numerical assessment of the hydrodynamics of the Nile swamps in southern Sudan. *Hydrol. Sci. J.* **2010**, *55*, 17–26. [CrossRef]

87. Sutcliffe, J.V. A Hydrological Study of the Southern Sudd Region of the Upper Nile. *Hydrol. Sci. Bull.* **1974**, *19*, 237–255. [CrossRef]

88. Berry, L.; Whiteman, A.J. The Nile in the Sudan. *Geogr. J.* **1968**, *134*, 1–33. [CrossRef]

89. El Shafie, A.G.A.; Elsayed Zeinelabdein, K.A.; Eisawi, A.A. Paleogeographic evolution and paleoenvironmental reconstruction of the Sudd area during the Early-Mid Holocene, Sudan. *J. Afr. Earth Sci.* **2011**, *60*, 13–18. [CrossRef]

90. Salama, R.B. The evolution of the River Nile. The buried saline rift lakes in Sudan—I. Bahr El Arab Rift, the Sudd buried saline lake. *J. Afr. Earth Sci.* (1983) **1987**, *6*, 899–913. [CrossRef]

91. Wolman, M.G.; Le Meur, C.; Giegengack, R.F. The Nile River: Geology, Hydrology, Hydraulic Society. In *Large Rivers*; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 2022; Chapter 24, pp. 704–736. [CrossRef]

92. Lindh, P.; Lemenkova, P. Permeability, compressive strength and Proctor parameters of silts stabilised by Portland cement and ground granulated blast furnace slag (GGBFS). *Arch. Mech. Eng.* **2022**, *69*, 667–692. [CrossRef]

93. Whiteman, A.J. *Geology of the Sudan Republic*; Cambridge University Press: Cambridge, UK, 1971.

94. Sutcliffe, J.V.; Parks, Y.P. Comparative water balances of selected African wetlands. *Hydrol. Sci. J.* **1989**, *34*, 49–62. [CrossRef]

95. Di Vittorio, C.A.; Georgakakos, A.P. Hydrologic Modeling of the Sudd Wetland using Satellite-based Data. *J. Hydrol. Reg. Stud.* **2021**, *37*, 100922. [CrossRef]

96. Woodward, J.C.; Macklin, M.G.; Krom, M.D.; Williams, M.A. The River Nile: Evolution and Environment. In *Large Rivers*; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 2022; Chapter 14, pp. 388–432. [CrossRef]

97. Sutcliffe, J.V.; Parks, Y.P. Hydrological modelling of the Sudd and Jonglei Canal. *Hydrol. Sci. J.* **1987**, *32*, 143–159. [CrossRef]

98. Birkett, C.; Murtugudde, R.; Allan, T. Indian Ocean Climate event brings floods to East Africa's lakes and the Sudd Marsh. *Geophys. Res. Lett.* **1999**, *26*, 1031–1034. [CrossRef]

99. Mitchell, S.A. The status of wetlands, threats and the predicted effect of global climate change: The situation in Sub-Saharan Africa. *Aquat. Sci.* **2013**, *75*, 95–112. [CrossRef]

100. Sutcliffe, J.; Brown, E. Water losses from the Sudd. *Hydrol. Sci. J.* **2018**, *63*, 527–541. [CrossRef]

101. Gabr, S.; El Bastawesy, M. The Implications of the Topographic, Hydrologic and Tectonic Settings On the Development of Bahr El-Ghazal Catchment, South Sudan. *Int. J. Water Resour. Arid. Environ.* **2013**, *2*, 90–101.

102. Benansio, J.S.; Funk, S.M.; Lino, J.L.; Balli, J.J.; Dante, J.O.; Dendi, D.; Fa, J.E.; Luiselli, L. Perceptions and attitudes towards climate change in fishing communities of the Sudd Wetlands, South Sudan. *Reg. Environ. Chang.* **2022**, *22*, 78. [CrossRef]

103. Bailey, R.G. An appraisal of the fisheries of the Sudd wetlands, River Nile, southern Sudan. *Aquac. Res.* **1989**, *20*, 79–89. [CrossRef]

104. Assessment, M.E. Ecosystems and Human Well-Being: Wetlands and Water; World Resources Institute: Washington, DC, USA, 2005.

105. Thompson, J.R.; Polet, G. Hydrology and land use in a sahelian floodplain wetland. *Wetlands* **2000**, *20*, 639–659. [CrossRef]

106. Fynn, R.W.S.; Murray-Hudson, M.; Dhliwayo, M.; Scholte, P. African wetlands and their seasonal use by wild and domestic herbivores. *Wetl. Ecol. Manag.* **2015**, *23*, 559–581. [CrossRef]

107. Pacini, N.; Hesslerová, P.; Pokorný, J.; Mwinami, T.; Morrison, E.H.; Cook, A.A.; Zhang, S.; Harper, D.M. Papyrus as an ecohydrological tool for restoring ecosystem services in Afrotropical wetlands. *Ecohydrol. Hydrobiol.* **2018**, *18*, 142–154.

108. Hickley, P.; Bailey, R.G. Fish communities in the eastern, seasonal-floodplain of the Sudd, Southern Sudan. *Hydrobiologia* **1987**, *144*, 243–250. [CrossRef]

109. Löw, F.; Stieglitz, K.; Diemar, O. Terrestrial oil spill mapping using satellite earth observation and machine learning: A case study in South Sudan. *J. Environ. Manag.* **2021**, *298*, 113424. [CrossRef]

110. Collins, C.D.; Banks-Leite, C.; Brudvig, L.A.; Foster, B.L.; Cook, W.M.; Damschen, E.I.; Andrade, A.; Austin, M.; Camargo, J.L.; Driscoll, D.A.; et al. Fragmentation affects plant community composition over time. *Ecography* **2017**, *40*, 119–130. [CrossRef]

111. Nagendra, H.; Lucas, R.; Honrado, J.P.; Jongman, R.H.; Tarantino, C.; Adamo, M.; Mairota, P. Remote sensing for conservation monitoring: Assessing protected areas, habitat extent, habitat condition, species diversity, and threats. *Ecol. Indic.* **2013**, *33*, 45–59. [CrossRef]

112. Yan, Y.; Jarvie, S.; Zhang, Q.; Zhang, S.; Han, P.; Liu, Q.; Liu, P. Small patches are hotspots for biodiversity conservation in fragmented landscapes. *Ecol. Indic.* **2021**, *130*, 108086. [CrossRef]

113. Martin, E.; Burgess, N. Sudd Flooded Grasslands. Online, OneEarth. 2023. Available online: https://www.oneearth.org/ ecoregions/sudd-flooded-grasslands/ (accessed on 10 August 2023).

114. Climatic Research Unit (CRU) of University of East Anglia. Climate Change Knowledge Portal. 2023. Available online: https://climateknowledgeportal.worldbank.org/country/sudan/climate-data-historical (accessed on 8 August 2023).
115. Wessel, P.; Luis, J.F.; Uieda, L.; Scharroo, R.; Wobbe, F.; Smith, W.H.F.; Tian, D. The Generic Mapping Tools Version 6. *Geochem. Geophys. Geosystems* **2019**, *20*, 5556–5564. [CrossRef]

116. Lemenkova, P. Console-Based Mapping of Mongolia Using GMT Cartographic Scripting Toolset for Processing TerraClimate Data. *Geosciences* **2022**, *12*, 140. [CrossRef]

117. Lemenkova, P. Handling Dataset with Geophysical and Geological Variables on the Bolivian Andes by the GMT Scripts. *Data* **2022**, *7*, 74. [CrossRef]

118. Hofierka, J.; Lacko, M.; Zubal, S. Parallelization of interpolation, solar radiation and water flow simulation modules in GRASS GIS using OpenMP. *Comput. Geosci.* **2017**, *107*, 20–27. [CrossRef]

119. Jasiewicz, J.; Metz, M. A new GRASS GIS toolkit for Hortonian analysis of drainage networks. *Comput. Geosci.* **2011**, *37*, 1162–1173. [CrossRef]

120. Jasiewicz, J. A new GRASS GIS fuzzy inference system for massive data analysis. *Comput. Geosci.* **2011**, *37*, 1525–1531. [CrossRef]

121. Sorokine, A. Implementation of a parallel high-performance visualization technique in GRASS GIS. *Comput. Geosci.* **2007**, *33*, 685–695. [CrossRef]

122. Neteler, M.; Bowman, M.H.; Landa, M.; Metz, M. GRASS GIS: A multi-purpose open source GIS. *Environ. Model. Softw.* **2012**, *31*, 124–130. [CrossRef]

123. Food and Agriculture Organization of the United Nations (FAO UN). *Land Cover Atlas of the Republic of South Sudan*; FAO: Rome, Italy, 2023. [CrossRef]

124. ReliefWeb. South Sudan: Floods—August 2014. 2014. Available online: https://mreliefweb.int/disaster/14337/fl-2014-000123-ssd?lang=fr (accessed on 9 August 2023).

125. ReliefWeb. Sudan: Floods—July 2018. 2018. Available online: https://reliefweb.int/disaster/fl-2018-000128-sdn (accessed on 11 August 2023).

126. United Nations Office for the Coordination of Humanitarian Affairs (UN OCHA). South Sudan: Floods Emergency Response Strategy and Funding Requirements as of 14 November 2019. 2019. Available online: https://reliefweb.int/report/south-sudan/south-sudan-floods-emergency-response-strategy-and-funding-requirements-14 (accessed on 11 August 2023).

127. United Nations Office for the Coordination of Humanitarian Affairs (UN OCHA). South Sudan Flooding Snapshot. 3 September 2020. Available online: https://www.unocha.org/ (accessed on 10 August 2023).

128. South Sudan Crisis Group. Floods, Displacement and Violence in South Sudan. 2021. Available online: https://southsudan.crisisgroup.org/ (accessed on 10 August 2023).

129. ReliefWeb. South Sudan: Floods 2021–2022. 2022. Available online: https://reliefweb.int/disaster/fl-2021-000108-ssd (accessed on 10 August 2023).

130. World Health Organization (WHO). Weekly Bulletin on Outbreaks and Other Emergencies. 2022. Available online: https://www.afro.who.int/health-topics/disease-outbreaks/outbreaks-and-other-emergencies-updates (accessed on 10 August 2023).

131. Wilusz, D.C.; Zaitchik, B.F.; Anderson, M.C.; Hain, C.R.; Yilmaz, M.T.; Mladenova, I.E. Monthly flooded area classification using low resolution SAR imagery in the Sudd wetland from 2007 to 2011. *Remote Sens. Environ.* **2017**, *194*, 205–218. [CrossRef]

132. Mohamed, Y.A.; van den Hurk, B.J.J.M.; Savenije, H.H.G.; Bastiaanssen, W.G.M. Impact of the Sudd wetland on the Nile hydroclimatology. *Water Resour. Res.* **2005**, *41*, 1–14. [CrossRef]

133. Petersen, G.; Sutcliffe, J.V.; Fohrer, N. Morphological analysis of the Sudd region using land survey and remote sensing data. *Earth Surf. Process. Landforms* **2008**, *33*, 1709–1720. [CrossRef]

**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.