Assessing the performance of machine learning algorithms in Google Earth Engine for land use and land cover analysis: A case study of Muğla province, Türkiye

Regions with high tourism density are very sensitive to human activities. Ensuring sustainability by preserving the cultural characteristics and natural structure of these regions is of critical importance in order to transfer these assets to the future world heritage. Detecting and mapping changes in land use and land cover (LULC) using innovative methods within short time intervals are of great importance for both monitoring the regional change and making administrative planning by taking necessary measures in a timely manner. In this context, this study focuses on the creation of a 4-class LULC map of Muğla province over the Google Earth Engine (GEE) platform by utilizing three different machine learning algorithms, namely, Support Vector Machines (SVM), Random Forest (RF), and Classification and Regression Tree (CART), and on comparison of their accuracy assessments. For improved classification accuracy, as well with the Sentinel-2 and Landsat-8 satellite images, the Normalized Difference Vegetation Index (NDVI) and Normalized Difference Water Index (NDWI) are also derived and used in classification of the major land use classes, which are ‘built-up area & barren land’, ‘dense vegetation’, ‘water surface’, and ‘shrub, grassland & sparse vegetation’. Experimental results show that the most relevant algorithm is RF with 0.97 overall accuracy and 0.96 Kappa value, followed by SVM and CART algorithms, respectively. These results indicate that the RF classifier outperforms both SVM and CART classifiers in terms of accuracy. Moreover, based on the results of the RF classifier, 19% (2,429 km2) of the study region is classified as built-up area & barren land, 48% (6,135 km2) as dense vegetation, 2% (301 km2) as water surface and 30% (3,832 km2) as shrub, grassland & sparse vegetation class

The Influence of the Time Equation on Remote Sensing Data Interpretation

The interpretation of optical Earth observation data (remote sensing data from satellites) requires knowledge of the exact geographic position of each pixel as well as the exact local acquisition time. But these parameters are not available in each case. If a satellite has a sun-synchronous orbit, equator crossing time (ECT) can be used to determine the local crossing time (LCT) and its corresponding solar zenith distance. Relation between local equator crossing time (LECT) and LCT is given by orbit geometry. The calculation is based on ECT of satellite. The method of actual ECT determination for different satellites on basis of the two-line-elements (TLE), available for their full lifetime period and with help of orbit prediction package is well known. For land surface temperature (LST) studies mean solar conditions are commonly used in the relation between ECT given in Coordinated Universal Time (UTC) and LECT given in hours, thus neglecting the difference between mean and real Sun time (MST, RST). Its difference is described by the equation of time (ET). Of particular importance is the variation of LECT during the year within about ±15 minutes. This is in each case the variation of LECT of a satellite, including satellites with stable orbit as LANDSAT (L8 around 10:05 a.m.) or ENVISAT (around 10:00 a.m.). In case of NOAA satellites the variation of LECT is overlaid by a long-term orbital drift. Ignatov et al. (2004) developed a method to describe the drift-based variation of LECT that can be viewed as a formal mathematical approximation of a periodic function with one or two Fourier terms. But, nevertheless, ET is not included in actual studies of LST. Our paper aims to demonstrate the possible influence of equation of time on simple examples of data interpretation, e.g. NDVI

Spatial Analysis of Post-Hurricane Katrina Thermal Pattern and Intensity in Greater New Orleans: Implications for Urban Heat Island Research

In 2005, Hurricane Katrina’s diverse impacts on the Greater New Orleans area included damaged and destroyed trees, and other despoiled vegetation, which also increased the exposure of artificial and bare surfaces, known factors that contribute to the climatic phenomenon known as the urban heat island (UHI). This is an investigation of UHI in the aftermath of Hurricane Katrina, which entails the analysis of pre and post-hurricane Katrina thermal imagery of the study area, including changes to surface heat patterns and vegetative cover. Imagery from Landsat TM was used to show changes to the pattern and intensity of the UHI effect, caused by an extreme weather event. Using remote sensing visualization methods, field data, and local knowledge, the author found there was a measurable change in the pattern and intensity of the New Orleans UHI effect, as well as concomitant changes to vegetative land cover. This finding may be relevant for urban planners and citizens, especially in the context of recovery from a large-scale disaster of a coastal city, regarding future weather events, and other natural and human impacts

Développement des algorithmes pour l’automatisation de la classification des données utilisant les réseaux de neurones probabilistes (PNN). Application à l’analyse, la catégorisation et la cartographie des images de télédétection.

The main topic of this thesis is modeling and classification for analyzing and processing digital data. Our objective is the realization of a set of algorithms to automate data classification using probabilistic neural networks (PNN). The considered data are remote sensing images.We propose a novel procedure for automatic classification based on PNN. We used unsupervised methods to search for classes in the learning phase; we built a function to validate classes inspired from fuzzy clustering techniques and we adapted our procedure to be applied on multidimensional data such as satellite images. Before applying this approach toremote sensing images, we conducted a series of tests on various types of synthetic and real data. These tests have produced very convincing results in comparison to usual unsupervised methods.We applied our algorithm in order to build land cover maps from satellite images. So, we have to analyze high spatial resolution images representing the study area for identifying all existing vegetation patterns. The classification is applied on sequence of NDVI time series data derived from satellite images. The aim is to test all of the developed algorithms on realcases to highlight their performance. These tests have demonstrated once again that the proposed approaches are entirely capable to produce successful classification. In semi-arid regions such as the region of Marrakech Tensift El Haouz, mapping land cover contributes extremely in evapotranspiration flow estimation used for water management.Le thème principal de cette thèse est la modélisation et la classification pour l’analyse et le traitement de l’information contenue dans les données numériques. Notre contribution est la réalisation d’un ensemble d’algorithmes pour l’automatisation de la classification des données en utilisant les réseaux de neurones probabilistes (PNN). Les données considérées sont les images de télédétection. Cette thèse s’adresse non seulement aux numériciens et spécialistes du traitement des images, mais aussi aux chercheurs et praticiens dans plusieurs domaines tels que la télédétection qui utilisent la classification des données en général et l’analyse de l’information pour la modélisation en particulier. Nous proposons une nouvelle procédure de classification automatique fondée sur les PNN. Nous avons utilisé des méthodes non supervisées pour la recherche des classes à la phase de l’apprentissage. Nous avons construit une fonction pour la validité des classes en s’inspirant des techniques de classification automatique floue. Nous avons aussi adapté notre procédure pour l’application sur des données multidimensionnelles telles que les images satellitaires. Avant d’appliquer cette approche sur des images de télédétection, nous avons mené une série de tests sur plusieurs types de données synthétiques et réelles. Ces tests ont abouti à des résultats très convaincants en comparaison avec les méthodes non superviséesusuelles, ce qui a conduit à valider la performance de nos algorithmes. Dans la partie application à la télédétection, l’objectif est d’élaborer des cartes d’occupation du sol à partir des images satellitaires. Nous avons analysé les images à haute résolution spatiale représentant la région étudiée pour identifier tous les profils de végétation existants. La classification est menée en se basant sur les images de l’indice de végétation NDVI extraites des images satellitaires SPOT. L’objet à classifier est une série temporelle de sept scènes NDVI. Le but est de tester l’ensemble des algorithmes développés sur des cas réels pour mettre en évidence leur performance. Ces tests ont démontré encore une fois de plus que les approches proposées sont tout à fait aptes à produire une classificationperformante.Nous avons classifié et analysé les images satellitaires d’une région semi-aride de Marrakech Tensift El Haouz. Le résultat obtenu est une cartographie précise de l’occupation du sol. Ce résultat contribuera d’une façon importante dans le dressage des cartes du flux d’évapotranspiration pour établir un bilan hydrique de la région

Spectroscopy-supported digital soil mapping

Global environmental changes have resulted in changes in key ecosystem services that soils provide. It is necessary to have up to date soil information on regional and global scales to ensure that these services continue to be provided. As a result, Digital Soil Mapping (DSM) research priorities are among others, advancing methods for data collection and analyses tailored towards large-scale mapping of soil properties. Scientifically, this thesis contributed to the development of methodologies, which aim to optimally use remote and proximal sensing (RS and PS) for DSM to facilitate regional soil mapping. The main contributions of this work with respect to the latter are (I) the critical evaluation of recent research achievements and identification of knowledge gaps for large-scale DSM using RS and PS data, (II) the development of a sparse RS-based sampling approach to represent major soil variability at regional scale, (III) the evaluation and development of different state-of-the-art methods to retrieve soil mineral information from PS, (IV) the improvement of spatially explicit soil prediction models and (V) the integration of RS and PS methods with geostatistical and DSM methods. A review on existing literature about the use of RS and PS for soil and terrain mapping was presented in Chapter 2. Recent work indicated the large potential of using RS and PS methods for DSM. However, for large-scale mapping, current methods will need to be extended beyond the plot. Improvements may be expected in the fields of developing more quantitative methods, enhanced geostatistical analysis and improved transferability to other areas. From these findings, three major research interests were selected: (I) soil sampling strategies, (II) retrieval of soil information from PS and (III) spatially continuous mapping of soil properties at larger scales using RS. Budgetary constraints, limited time and available soil legacy data restricted the soil data acquisition, presented in Chapter 3. A 15.000 km2 area located in Northern Morocco served as test case. Here, a sample was collected using constrained Latin Hypercube Sampling (cLHS) of RS and elevation data. The RS data served as proxy for soil variability, as alternative for the required soil legacy data supporting the sampling strategy. The sampling aim was to optimally sample the variability in the RS data while minimizing the acquisition efforts. This sample resulted in a dataset representing major soil variability. The cLHS sample failed to express spatial correlation; constraining the LHS by a distance criterion favoured large spatial variability over short distances. The absence of spatial correlation in the sampled soil variability precludes the use of additional geostatistical analyses to spatially predict soil properties. Predicting soil properties using the cLHS sample is thus restricted to a modelled statistical relation between the sample and exhaustive predictor variables. For this, the RS data provided the necessary spatial information because of the strong spatial correlation while the spectral information provided the variability of the environment (Chapter 3 and 6). Concluding, the RS-based cLHS approach is considered a time and cost efficient method for acquiring information on soil resources over extended areas. This sample was further used for developing methods to derive soil mineral information from PS, and to characterize regional soil mineralogy using RS. In Chapter 4, the influences of complex scattering within the mixture and overlapping absorption features were investigated. This was done by comparing the success of PRISM’s MICA in determining mineralogy of natural samples and modelled spectra. The modelled spectra were developed by a linearly forward model of reflectance spectra, using the fraction of known constituents within the sample. The modelled spectra accounted for the co-occurrence of absorption features but eluded the complex interaction between the components. It was found that more minerals could be determined with higher accuracy using modelled reflectance. The absorption features in the natural samples were less distinct or even absent, which hampered the classification routine. Nevertheless, grouping the individual minerals into mineral categories significantly improved the classification accuracy. These mineral categories are particularly useful for regional scale studies, as key soil property for parent material characterization and soil formation. Characterizing regional soil mineralogy by mineral categories was further described in Chapter 6. Retrieval of refined information from natural samples, such as mineral abundances, is more complex; estimating abundances requires a method that accounts for the interaction between minerals within the intimate mixture. This can be done by addressing the interaction with a non-linear model (Chapter 5). Chapter 5 showed that mineral abundances in complex mixtures could be estimated using absorption features in the 2.1–2.4 µm wavelength region. First, the absorption behaviour of mineral mixtures was parameterized by exponential Gaussian optimization (EGO). Next, mineral abundances were successfully predicted by regression tree analysis, using these parameters as inputs. Estimating mineral abundances using prepared mixes of calcite, kaolinite, montmorillonite and dioctahedral mica or field samples proved the validity of the proposed method. Estimating mineral abundances of field samples showed the necessity to deconvolve spectra by EGO. Due to the nature of the field samples, the simple representation of the complex scattering behaviour by a few Gaussian bands required the parameters asymmetry and saturation to accurately deconvolve the spectra. Also, asymmetry of the EGO profiles showed to be an important parameter for estimating the abundances of the field samples. The robustness of the method in handling the omission of minerals during the training phase was tested by replacing part of the quartz with chlorite. It was found that the accuracy of the predicted mineral content was hardly affected. Concluding, the proposed method allowed for estimating more than two minerals within a mixture. This approach advances existing PS methods and has the potential to quantify a wider set of soil properties. With this method the soil science community was provided an improved inference method to derive and quantify soil properties The final challenge of this thesis was to spatially explicit model regional soil mineralogy using the sparse sample from Chapter 3. Prediction models have especially difficulties relating predictor variables to sampled properties having high spatial correlation. Chapter 6 presented a methodology that improved prediction models by using scale-dependent spatial variability observed in RS data. Mineral predictions were made using the abundances from X-ray diffraction analysis and mineral categories determined by PRISM. The models indicated that using the original RS data resulted in lower model performance than those models using scaled RS data. Key to the improved predictions was representing the variability of the RS data at the same scale as the sampled soil variability. This was realized by considering the medium and long-range spatial variability in the RS data. Using Fixed Rank Kriging allowed smoothing the massive RS datasets to these ranges. The resulting images resembled more closely the regional spatial variability of soil and environmental properties. Further improvements resulted from using multi-scale soil-landscape relationships to predict mineralogy. The maps of predicted mineralogy showed agreement between the mineral categories and abundances. Using a geostatistical approach in combination with a small sample, substantially improves the feasibility to quantitatively map regional mineralogy. Moreover, the spectroscopic method appeared sufficiently detailed to map major mineral variability. Finally, this approach has the potential for modelling various natural resources and thereby enhances the perspective of a global system for inventorying and monitoring the earth’s soil resources. With this thesis it is demonstrated that RS and PS methods are an important but also an essential source for regional-scale DSM. Following the main findings from this thesis, it can be concluded that: Improvements in regional-scale DSM result from the integrated use of RS and PS with geostatistical methods. In every step of the soil mapping process, spectroscopy can play a key role and can deliver data in a time and cost efficient manner. Nevertheless, there are issues that need to be resolved in the near future. Research priorities involve the development of operational tools to quantify soil properties, sensor integration, spatiotemporal modelling and the use of geostatistical methods that allow working with massive RS datasets. This will allow us in the near future to deliver more accurate and comprehensive information about soils, soil resources and ecosystem services provided by soils at regional and, ultimately, global scale.</p