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

: 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 reﬂectance 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 classiﬁcation. The methodology includes image segmentation by ‘i.segment’ module, image clustering and classiﬁcation by ’i.cluster’ and ’i.maxlike’ modules, accuracy assessment by ‘r.kappa’ module, and computing NDVI and cartographic mapping implemented using GRASS GIS. The beneﬁts 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 deﬁned 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 classiﬁcation.


Introduction 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).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][3][4].

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][11][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 images 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, multilevel thresholding [19][20][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][29][30].

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][48][49].More complex cases reduce the level of the fragmentation through contouring segment carcasses derived from the upscaled 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][53][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.

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][56][57][58].Thus, specifically for image segmentation, machine learning proposes embedded techniques of image discrimination to find the contours of the objects [59][60][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.

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][70][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][74][75][76].In contrast to the existing image classification techniques which group pixels with similar spectral reflectance into classes [77][78][79][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 objectoriented 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.

Characterisation of Study Area
The Sudd is a large area of wetlands located in South Sudan (Figure 1).The total area of swamp varyies between 30,000 to 40,000 km 2 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. 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].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].

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

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.

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 takes 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 evaluated 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.

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][119][120][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.

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

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.

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.

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.

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.

Clustering
The performance of the segmentation was a slow process and required further postprocessing 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 discussed 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.

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] (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.

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.

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.

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.
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).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 'i.segment' module, the images were segmented without the priori landscape-dependent information, which ensured independent interpretation.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).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].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.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.

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.

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.

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 opensource 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 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.segment' 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 efficiently of this software for the processing and segmentation of satellite images.

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.

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.Listing A2.GRASS GIS code for creating semantic labels for the Landsat OLI/TIRS.
# listing the available raster bands in GRASS GIS mapset g .list rast # obtaining information on raster metadata : r .info -r L9_2023_07 # grouping data by i .group and set up a computational region to match the scene g .region raster = L9_2023_01 -p # store VIZ , NIR , MIR into group / subgroup and leaving out TIR as redundant i .group group = L9_2023 subgroup = res_30m \ input = L9_2023_01 , L9_2023_02 , L9_2023_03 , L9_2023_04 , L9_2023_05 , L9_2023_06 , L9_2023_07 # Set the region test area with the resolution taken from the input Landsat bands g .region -p raster = L9_2023_01 Listing A3.GRASS GIS code for segmentation for image tested with 2 levels of threshold.

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

Figure 1 .
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.

Figure 7 .
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).

Figure 8 .
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.

Figure 10 .
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).

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.

Table 3 .
Results of the segmentation procedure for Landsat 8-9 images with No of created segments.

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.

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.

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.

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.

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.

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.

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.

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.

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.