Google Earth Engine cloud computing platform for remote sensing big data applications: a comprehensive review

Remote sensing (RS) systems have been collecting massive volumes of datasets for decades, managing and analyzing of which are not practical using common software packages and desktop computing resources. In this regard, Google has developed a cloud computing platform, called Google Earth Engine (GEE), to effectively address the challenges of big data analysis. In particular, this platformfacilitates processing big geo data over large areas and monitoring the environment for long periods of time. Although this platformwas launched in 2010 and has proved its high potential for different applications, it has not been fully investigated and utilized for RS applications until recent years. Therefore, this study aims to comprehensively explore different aspects of the GEE platform, including its datasets, functions, advantages/limitations, and various applications. For this purpose, 450 journal articles published in 150 journals between January 2010 andMay 2020 were studied. It was observed that Landsat and Sentinel datasets were extensively utilized by GEE users. Moreover, supervised machine learning algorithms, such as Random Forest, were more widely applied to image classification tasks. GEE has also been employed in a broad range of applications, such as Land Cover/land Use classification, hydrology, urban planning, natural disaster, climate analyses, and image processing. It was generally observed that the number of GEE publications have significantly increased during the past few years, and it is expected that GEE will be utilized by more users from different fields to resolve their big data processing challenges.Peer ReviewedPostprint (published version

The third generation of pan-canadian wetland map at 10 m resolution using multisource earth observation data on cloud computing platform

Development of the Canadian Wetland Inventory Map (CWIM) has thus far proceeded over two generations, reporting the extent and location of bog, fen, swamp, marsh, and water wetlands across the country with increasing accuracy. Each generation of this training inventory has improved the previous results by including additional reference wetland data and focusing on processing at the scale of ecozone, which represent ecologically distinct regions of Canada. The first and second generations attained relatively highly accurate results with an average approaching 86% though some overestimated wetland extents, particularly of the swamp class. The current research represents a third refinement of the inventory map. It was designed to improve the overall accuracy (OA) and reduce wetlands overestimation by modifying test and train data and integrating additional environmental and remote sensing datasets, including countrywide coverage of L-band ALOS PALSAR-2, SRTM, and Arctic digital elevation model, nighttime light, temperature, and precipitation data. Using a random forest classification within Google Earth Engine, the average OA obtained for the CWIM3 is 90.53%, an improvement of 4.77% over previous results. All ecozones experienced an OA increase of 2% or greater and individual ecozone OA results range between 94% at the highest to 84% at the lowest. Visual inspection of the classification products demonstrates a reduction of wetland area overestimation compared to previous inventory generations. In this study, several classification scenarios were defined to assess the effect of preprocessing and the benefits of incorporating multisource data for large-scale wetland mapping. In addition, the development of a confidence map helps visualize where current results are most and least reliable given the amount of wetland test and train data and the extent of recent landscape disturbance (e.g., fire). The resulting OAs and wetland areal extent reveal the importance of multisource data and adequate test and train data for wetland classification at a countrywide scale

Wetland Mapping in Great Lakes Using Sentinel-1/2 Time-Series Imagery and DEM Data in Google Earth Engine

The Great Lakes (GL) wetlands support a variety of rare and endangered animal and plant species. Thus, wetlands in this region should be mapped and monitored using advanced and reliable techniques. In this study, a wetland map of the GL was produced using Sentinel-1/2 datasets within the Google Earth Engine (GEE) cloud computing platform. To this end, an object-based supervised machine learning (ML) classification workflow is proposed. The proposed method contains two main classification steps. In the first step, several non-wetland classes (e.g., Barren, Cropland, and Open Water), which are more distinguishable using radar and optical Remote Sensing (RS) observations, were identified and masked using a trained Random Forest (RF) model. In the second step, wetland classes, including Fen, Bog, Swamp, and Marsh, along with two non-wetland classes of Forest and Grassland/Shrubland were identified. Using the proposed method, the GL were classified with an overall accuracy of 93.6% and a Kappa coefficient of 0.90. Additionally, the results showed that the proposed method was able to classify the wetland classes with an overall accuracy of 87% and a Kappa coefficient of 0.91. Non-wetland classes were also identified more accurately than wetlands (overall accuracy = 96.62% and Kappa coefficient = 0.95)

Wetland Mapping in Great Lakes Using Sentinel-1/2 Time-Series Imagery and DEM Data in Google Earth Engine

The Great Lakes (GL) wetlands support a variety of rare and endangered animal and plant species. Thus, wetlands in this region should be mapped and monitored using advanced and reliable techniques. In this study, a wetland map of the GL was produced using Sentinel-1/2 datasets within the Google Earth Engine (GEE) cloud computing platform. To this end, an object-based supervised machine learning (ML) classification workflow is proposed. The proposed method contains two main classification steps. In the first step, several non-wetland classes (e.g., Barren, Cropland, and Open Water), which are more distinguishable using radar and optical Remote Sensing (RS) observations, were identified and masked using a trained Random Forest (RF) model. In the second step, wetland classes, including Fen, Bog, Swamp, and Marsh, along with two non-wetland classes of Forest and Grassland/Shrubland were identified. Using the proposed method, the GL were classified with an overall accuracy of 93.6% and a Kappa coefficient of 0.90. Additionally, the results showed that the proposed method was able to classify the wetland classes with an overall accuracy of 87% and a Kappa coefficient of 0.91. Non-wetland classes were also identified more accurately than wetlands (overall accuracy = 96.62% and Kappa coefficient = 0.95)

Advanced machine learning algorithms for Canadian wetland mapping using polarimetric synthetic aperture radar (PolSAR) and optical imagery

Wetlands are complex land cover ecosystems that represent a wide range of biophysical conditions. They are one of the most productive ecosystems and provide several important environmental functionalities. As such, wetland mapping and monitoring using cost- and time-efficient approaches are of great interest for sustainable management and resource assessment. In this regard, satellite remote sensing data are greatly beneficial, as they capture a synoptic and multi-temporal view of landscapes. The ability to extract useful information from satellite imagery greatly affects the accuracy and reliability of the final products. This is of particular concern for mapping complex land cover ecosystems, such as wetlands, where complex, heterogeneous, and fragmented landscape results in similar backscatter/spectral signatures of land cover classes in satellite images. Accordingly, the overarching purpose of this thesis is to contribute to existing methodologies of wetland classification by proposing and developing several new techniques based on advanced remote sensing tools and optical and Synthetic Aperture Radar (SAR) imagery. Specifically, the importance of employing an efficient speckle reduction method for polarimetric SAR (PolSAR) image processing is discussed and a new speckle reduction technique is proposed. Two novel techniques are also introduced for improving the accuracy of wetland classification. In particular, a new hierarchical classification algorithm using multi-frequency SAR data is proposed that discriminates wetland classes in three steps depending on their complexity and similarity. The experimental results reveal that the proposed method is advantageous for mapping complex land cover ecosystems compared to single stream classification approaches, which have been extensively used in the literature. Furthermore, a new feature weighting approach is proposed based on the statistical and physical characteristics of PolSAR data to improve the discrimination capability of input features prior to incorporating them into the classification scheme. This study also demonstrates the transferability of existing classification algorithms, which have been developed based on RADARSAT-2 imagery, to compact polarimetry SAR data that will be collected by the upcoming RADARSAT Constellation Mission (RCM). The capability of several well-known deep Convolutional Neural Network (CNN) architectures currently employed in computer vision is first introduced in this thesis for classification of wetland complexes using multispectral remote sensing data. Finally, this research results in the first provincial-scale wetland inventory maps of Newfoundland and Labrador using the Google Earth Engine (GEE) cloud computing resources and open access Earth Observation (EO) collected by the Copernicus Sentinel missions. Overall, the methodologies proposed in this thesis address fundamental limitations/challenges of wetland mapping using remote sensing data, which have been ignored in the literature. These challenges include the backscattering/spectrally similar signature of wetland classes, insufficient classification accuracy of wetland classes, and limitations of wetland mapping on large scales. In addition to the capabilities of the proposed methods for mapping wetland complexes, the use of these developed techniques for classifying other complex land cover types beyond wetlands, such as sea ice and crop ecosystems, offers a potential avenue for further research

Remote sensing of wetlands in the Lake Whangape catchment, Waikato, New Zealand.

Wetlands are among the world's most valuable ecosystems. They provide numerous ecological and socio-economic benefits. However, wetlands continue to disappear due to the increasing demand for wetland resources. In New Zealand, more than 90% of the original extent of wetlands has been lost since the mid-eighteenth century. Therefore, legislation has been identified for the protection of wetlands as a matter of national importance. Geographic Information System (GIS) and Remote Sensing (RS) techniques have proven helpful for mapping and monitoring wetland resources. This study aims to understand how RStechniques can classify wetlands in the Lake Whangape catchment, Waikato. The parameters that can be extracted from available data and their effectiveness in the classification process are also studied. Four types of input data are collectively employed in the study. The data types are optical RS data, Synthetic Aperture Radar (SAR) data, a Digital Elevation Model (DEM), and wetland polygons provided by the Waikato Regional Council (WRC). All the steps including, accessing satellite scenes and data processing were performed within Google Earth Engine (GEE) computing platform using JavaScript language. The classification process for this study includes feature extraction, feature selection, model training, classification, and validation. Finally, the accuracy of the classification results is checked visually and statistically. The classification was carried out in two stages. In Stage one, open water, wetland, and non-wetland areas are classified (simple classification). The combined wetlands class is separated into marsh and swamp in the second stage (detailed classification). Based on the results, the Topographic Position Index (TPI) is the most influential parameter in identifying wetlands, while the Modified Normalized Water Index (MNDWI) successfully identifies open water. The overall accuracy reached 91% at the simple classification stage. However, the detailed classification results received comparatively low classification accuracies (the overall accuracy is 76%)

Status and trends of wetland studies in Canada using remote sensing technology with a focus on wetland classification: a bibliographic analysis

A large portion of Canada is covered by wetlands; mapping and monitoring them is of great importance for various applications. In this regard, Remote Sensing (RS) technology has been widely employed for wetland studies in Canada over the past 45 years. This study evaluates meta-data to investigate the status and trends of wetland studies in Canada using RS technology by reviewing the scientific papers published between 1976 and the end of 2020 (300 papers in total). Initially, a meta-analysis was conducted to analyze the status of RS-based wetland studies in terms of the wetland classification systems, methods, classes, RS data usage, publication details (e.g., authors, keywords, citations, and publications time), geographic information, and level of classification accuracies. The deep systematic review of 128 peer-reviewed articles illustrated the rising trend in using multi-source RS datasets along with advanced machine learning algorithms for wetland mapping in Canada. It was also observed that most of the studies were implemented over the province of Ontario. Pixel-based supervised classifiers were the most popular wetland classification algorithms. This review summarizes different RS systems and methodologies for wetland mapping in Canada to outline how RS has been utilized for the generation of wetland inventories. The results of this review paper provide the current state-of-the-art methods and datasets for wetland studies in Canada and will provide direction for future wetland mapping research.Peer ReviewedPostprint (published version

Automatic mapping of burned areas using Landsat 8 time-series images in Google Earth engine: a case study from Iran

Due to the natural conditions and inappropriate management responses, large part of plains and forests in Iran have been burned in recent years. Given the increasing availability of open-access satellite images and open-source software packages, we developed a fast and cost-effective remote sensing methodology for characterizing burned areas for the entire country of Iran. We mapped the fire-affected areas using a post-classification supervised method and Landsat 8 time-series images. To this end, the Google Earth Engine (GEE) and Google Colab computing services were used to facilitate the downloading and processing of images as well as allowing for effective implementation of the algorithms. In total, 13 spectral indices were calculated using Landsat 8 images and were added to the nine original bands of Landsat 8. The training polygons of the burned and unburned areas were accurately distinguished based on the information acquired from the Iranian Space Agency (ISA), Sentinel-2 images, and Fire Information for Resource Management System (FIRMS) products. A combination of Genetic Algorithm (GA) and Neural Network (NN) approaches was then implemented to specify 19 optimal features out of the 22 bands. The 19 optimal bands were subsequently applied to two classifiers of NN and Random Forest (RF) in the timespans of 1 January 2019 to 30 December 2020 and of 1 January 2021 to 30 September 2021. The overall classification accuracies of 94% and 96% were obtained for these two classifiers, respectively. The omission and commission errors of both classifiers were also less than 10%, indicating the promising capability of the proposed methodology in detecting the burned areas. To detect the burned areas caused by the wildfire in 2021, the image differencing method was used as well. The resultant models were finally compared to the MODIS fire products over 10 sampled polygons of the burned areas. Overall, the models had a high accuracy in detecting the burned areas in terms of shape and perimeter, which can be further implicated for potential prevention strategies of endangered biodiversity.Peer ReviewedPostprint (published version

Multitemporal optical and radar metrics for wetland mapping at national level in Albania

Wetlands are highly dynamic, with many natural and anthropogenic drivers causing seasonal, periodic or permanent changes in their structure and composition. Thus, it is necessary to use time series of images for accurate classifications and monitoring. We used all available Sentinel-1 and Sentinel-2 images to produce a national wetlands map for Albania. We derived different indices and temporal metrics and investigated their impacts and synergies in terms of mapping accuracy. Best results were achieved when combining Sentinel-1 with Sentinel-2 and its derived indices. We reduced systematic errors and increased the thematic resolution using morphometric characteristics and knowledge-based rules, achieving an overall accuracy of 82%. Results were also validated against field inventories. This methodology can be reproducible to other countries and can be made operational for an integrated planning that considers the food, water, and energy nexus

Mapping intra- and inter-annual dynamics in wetlands with multispectral, thermal and SAR time series

Kartierung der intra- und interannuellen Dynamik von Feuchtgebieten mit multispektralen, thermischen und SAR-Zeitreihen Die Analyse der aktuellen räumlichen Verbreitung und der zeitlichen Entwicklung von Feuchtgebieten stellt eine äußerst komplexe Aufgabe dar, welche durch die Saisonalität, die schwierige Zugänglichkeit und die besonderen Eigenschaften als Ökoton bedingt ist. Erdbeobachtungssysteme sind somit das am besten geeignete Werkzeug, um zeitliche und räumliche Muster von Feuchtgebieten auf globaler Ebene zu beobachten (saisonale Veränderungen und Langzeit-Trends) und um den Einfluss der menschlichen Aktivitäten auf ihre physischen und biologischen Eigenschaften zu untersuchen. Zur Kartierung von raum-zeitlichen Mustern wurden Zeitreihen von Radar- (Sentinel-1), Multispektral- (Sentinel-2) und Thermal-Satellitendaten (MODIS) in fünf Untersuchungsgebieten, mit für Feuchtgebiete unterschiedlichen typischen Charakteristika, untersucht. In Kapitel 1 werden die Problematik in Bezug auf die Definition von Feuchtgebieten erläutert und allgemeine Degradations-Trends beschrieben. Die Kapitel 2 und 3 behandeln einen Algorithmus, der Veränderungen mithilfe von SAR-Zeitreihen feststellt, sowie die Vorteile des Cloud-Computings für das operationelle Monitoring saisonaler Muster und die Erkennung kurzfristig auftretender Veränderungen. In den Kapiteln 4 und 5 werden die zwei Hauptursachen für den Verlust von Feuchtgebieten betrachtet: der Staudammbau und die Ausdehnung landwirtschaftlicher Flächen. In Kapitel 4 werden dichte Zeitreihen multispektraler (Sentinel-2) und SAR-Daten (Sentinel-1) verwendet, um die Feuchtgebiete Albaniens – eines Landes in dem konträre Pläne zum Ausbau seines Wasserkraftpotentials und dem Schutz intakter Flussökosysteme zu Spannungen führen – landesweit zu kartieren. Die synergetischen Vorteile, die sich durch die Fusionierung von multispektralen und SAR-Daten für die Klassifikation ergeben, werden dabei herausgestellt. Kapitel 5 veranschaulicht, dass die Kilombero-Überschwemmungsebene in Tansania ein großes und bedeutendes Feuchtgebiet ist, das in den vergangenen Jahren infolge der weitgehend unkontrollierten Ausbreitung landwirtschaftlicher Flächen in seiner Ausdehnung und seiner Ökologie stark beeinträchtigt wurde. Um die Auswirkungen der Landnutzungsänderungen des Feuchtgebietes während der vergangenen 18 Jahre zu analysieren, wurden eine Zeitreihe (2000 bis 2017) thermaler Daten (MODIS) analysiert. Die drei für die Zeitreihenanalyse angewandten Modelle zeigen, wie landwirtschaftliche Praktiken die Landoberflächentemperatur in den landwirtschaftlich genutzten Gebieten sowie in den angrenzenden natürlichen Feuchtgebieten erhöht haben.Due to wetlands’ seasonality, their difficult access and ecotone character, determining their actual extension and trends over time is a complex task. Earth Observation systems are the most appropriate tool to monitor their spatio-temporal patterns (seasonal changes and long term trends) at global scales, and to study the effects that human activities have in their physical and biological properties. In this work I use time series of radar (Sentinel-1), multispectral (Sentinel-2) and thermal (MODIS) imagery to map the spatio-temporal patterns in 5 wetlands of different characteristics. First, I introduce in chapter 1 the problematic of wetlands’ definitions and their degradation trends. I continue with a brief introduction on remote sensing, time series analysis, and their applications on wetlands’ research and management. In chapters 2 and 3 I implement an algorithm for change detection of time series of Sentinel-1 images and demonstrate the advantages of cloud computation for operational monitoring. In chapters 4 and 5 I address two of the main causes of wetland degradation: dam building and agricultural expansion. In chapter 4 I use dense time series of Sentinel-1 and Sentinel-2 images map all the wetlands of Albania; a country struggling between developing its large hydropower potential or preserving its intact and valuable river ecosystems. I evaluate the synergic advantages of fusing multispectral and radar imagery in combination with knowledge-based rules to produce classification of higher thematic and spatial resolutions. In chapter 5 I present how the Kilombero Floodplain, in Tanzania, has been degraded during the last years due to uncontrolled farmland expansion. I use a time series of thermal imagery (MODIS) from 2000 until 2017 to analyze the effect of land use changes on the wetland. I compare three models for time series analysis and reveal how farming practices have increased the surface temperature of the farmed area, as well as in adjacent natural wetlands.Mapeo de las dinámicas inter- e intra-anuales en humedales con series temporales de imágenes multiespectrales, termales y de radar Debido a la estacionalidad de los humedales, su difícil acceso y sus características de ecotono, determinar su actual extensión y sus tendencias a lo largo del tiempo es una tarea compleja. Los sistemas de observación terrestres son la herramienta más apropiada para monitorear sus patrones espacio-temporales (estacionalidad y tendencias a largo plazo) a escalas globales, y para estudiar los efectos que las actividades humanas causan en sus propiedades físicas y biológicas. En esta tesis uso series temporales de imágenes radar (Sentinel-1), multiespectrales (Sentinel-2) y termales (MODIS) para mapear los patrones espacio-temporales de 5 humedales de diferentes características. En el capítulo 1 describo los retos que derivan de las diferentes definiciones que existen de los humedales. También presento las tendencias globales de degradación que la mayoría de los humedales continúan experimentando en los últimos años. Continúo con una breve introducción de los sistemas de teledetección remota, análisis de series temporales, y sus aplicaciones a la investigación y gestión de los humedales. En los capítulos 2 y 3 implemento un algoritmo de detección de cambios para series temporales de imágenes radar, y muestro las ventajas de usar sistemas de computación en la nube para monitorear cambios en la cobertura del suelo a corto plazo. En los capítulos 4 y 5 trato con dos de las causas más comunes de degradación de humedales: la construcción de presas y la expansión de la agricultura. En el capítulo 4 uso series temporales de imágenes multiespectrales (Sentinel-2) y radar (Sentinel-1) para mapear todos los humedales Albania; un país que se debate entre desarrollar su potencial hidroenergético o preservar sus valiosos e intactos ecosistemas de rivera. Mediante la fusión de imágenes radar y multiespectrales y el uso de reglas de decisión genero un mapa de suficiente resolución espacial y temática para que pueda ser usado por sectores interesados y gestores. En el capítulo 5 presento como las llanuras inundables de Kilombero, en Tanzania, han sido degradadas durante los últimos años debido a la expansión incontrolada de la agricultura. Usando series temporales de imágenes termales (MODIS) desde 2000 hasta 2017 y mapas de cambios de usos del suelo, determino los efectos que estos cambios han tenido en el humedal. Comparo 3 modelos diferentes de análisis de series temporales y muestro cómo la expansión de la agricultura ha incrementado la temperatura superficial terrestre, no solo de la zona cultivada, sino también de zonas adyacentes aún naturales