At present, accessing and processing Earth Observation (EO) data on different cloud platforms requires users to exercise distinct communication strategies as each backend platform is designed differently. The openEO API (Application Programming Interface) standardises EO-related contracts between local clients (R, Python, and JavaScript) and cloud service providers regarding data access and processing, simplifying their direct comparability. Independent of the providers’ data storage system, the API mimics the functionalities of a virtual EO raster data cube. This article introduces the communication strategy and aspects of the data cube model applied by the openEO API. Two test cases show the potential and current limitations of processing similar workflows on different cloud platforms and a comparison of the result of a locally running workflow and its openEO-dependent cloud equivalent. The outcomes demonstrate the flexibility of the openEO API in enabling complex scientific analysis of EO data collections on cloud platforms in a homogenised way.

A data assimilation method to produce complete temporal sequences of synthetic medium-resolution images is presented. The method implements a Kalman filter recursive algorithm that integrates medium and moderate resolution imagery. To demonstrate the approach, time series of 30-m spatial resolution NDVI images at 16-day time steps were generated using Landsat NDVI images and MODIS NDVI products at four sites with different ecosystems and land cover-land use dynamics. The results show that the time series of synthetic NDVI images captured seasonal land surface dynamics and maintained the spatial structure of the landscape at higher spatial resolution. The time series of synthetic medium-resolution NDVI images were validated within a Monte Carlo simulation framework. Normalized residuals decreased as the number of available observations increased, ranging from 0.2 to below 0.1. Residuals were also significantly lower for time series of synthetic NDVI images generated at combined recursion (smoothing) than individually at forward and backward recursions (filtering). Conversely, the uncertainties of the synthetic images also decreased when the number of available observations increased and combined recursions were implemented.

Motivated by the increasing availability of open and free Earth observation data through the Copernicus Sentinel missions, this study investigates the capacity of advanced computational models to automatically generate thematic layers, which in turn contribute to and facilitate the creation of land cover products. In concrete terms, we assess the practical and computational aspects of multi-class Sentinel-2 image segmentation based on a convolutional neural network and random forest approaches. The annotated learning set derives from data that is made available as result of the implementation of European Union’s INSPIRE Directive. Since this network of data sets remains incomplete in regard to some geographic areas, another objective of this work was to provide consistent and reproducible ways for machine-driven mapping of these gaps and a potential update of the existing ones. Finally, the performance analysis identifies the most important hyper-parameters, and provides hints on the models’ deployment and their transferability.

Processing large amounts of image data such as the Sentinel-2 archive is a computationally demanding task. However, for most applications, many of the images in the archive are redundant and do not contribute to the quality of the final result. An optimization scheme is presented here that selects a subset of the Sentinel-2 archive in order to reduce the amount of processing, while retaining the quality of the resulting output. As a case study, we focused on the creation of a cloud-free composite, covering the global land mass and based on all the images acquired from January 2016 until September 2017. The total amount of available images was 2,128,556. The selection of the optimal subset was based on quicklooks, which correspond to a spatial and spectral subset of the original Sentinel-2 products and are lossy compressed. The selected subset contained 94,093 image tiles in total, reducing the amount of images to be processed to 4.42% of the full set.

Appropriate management of (semi-)natural areas requires detailed knowledge of the ecosystems present and their status. Remote sensing can provide a systematic, synoptic view at regular time intervals, and is therefore often suggested as a powerful tool to assist with the mapping and monitoring of protected habitats and vegetation. In this study, we present a multi-step mapping framework that enables detailed NATURA 2000 (N2000) heathland habitat patch mapping and the assessment of their conservation status at patch level. The method comprises three consecutive steps: (1) a hierarchical land/vegetation type (LVT) classification using airborne AHS imaging spectroscopy and field reference data; (2) a spatial re-classification to convert the LVT map to a patch map based on life forms; and (3) identification of the N2000 habitat type and conservation status parameters for each of the patches. Based on a multivariate analysis of 1325 vegetation reference plots acquired in 2006–2007, 24 LVT classes were identified that were considered relevant for the assessment of heathland conservation status. These labelled data were then used as ground reference for the supervised classification of the AHS image data to an LVT classification map, using Linear Discriminant Analysis in combination with Sequential-Floating-Forward-Search feature selection. Overall classification accuracies for the LVT mapping varied from 83% to 92% (Kappa ≈ 0.82–0.91), depending on the level of detail in the hierarchical classification. After converting the LVT map to a N2000 habitat type patch map, an overall accuracy of 89% was obtained. By combining the N2000 habitat type patch map with the LVT map, two important conservation status parameters were directly deduced per patch: tree and shrub cover, and grass cover, showing a strong similarity to an independent dataset with estimates made in the field in 2009. The results of this study indicate the potential of imaging spectroscopy for detailed heathland habitat characterization of N2000 sites in a way that matches the current field-based workflows of the user.