In May 2019, Collection 2 of the Copernicus Global Land Cover layers was released. Next to a global discrete land cover map at 100 m resolution, a set of cover fraction layers is provided depicting the percentual cover of the main land cover types in a pixel. This additional continuous classification scheme represents areas of heterogeneous land cover better than the standard discrete classification scheme. Overall, 20 layers are provided which allow customization of land cover maps to specific user needs or applications (e.g., forest monitoring, crop monitoring, biodiversity and conservation, climate modeling, etc.). However, Collection 2 was not just a global up-scaling, but also includes major improvements in the map quality, reaching around 80% or more overall accuracy. The processing system went into operational status allowing annual updates on a global scale with an additional implemented training and validation data collection system. In this paper, we provide an overview of the major changes in the production of the land cover maps, that have led to this increased accuracy, including aligning with the Sentinel 2 satellite system in the grid and coordinate system, improving the metric extraction, adding better auxiliary data, improving the biome delineations, as well as enhancing the expert rules. An independent validation exercise confirmed the improved classification results. In addition to the methodological improvements, this paper also provides an overview of where the different resources can be found, including access channels to the product layer as well as the detailed peer-review product documentation.

Airborne remote sensing with a CASI‐550 sensor has been used to map the benthic coverage and the bottom topography of the Pulau Nukaha coral reef located in the Tanimbar Archipelago (Southeast Moluccas, Eastern Indonesia). The image classification method adopted was performed in three steps. Firstly, five geomorphological reef components were identified using a supervised spectral angle mapping algorithm in combination with data collected during the field survey, i.e. benthic cover type, percentage cover and depth. Secondly, benthic cover mapping was performed for each of the five geomorphological components separately using an unsupervised hierarchical clustering algorithm followed by class aggregation using both spectral and spatial information. Finally, 16 benthic cover classes could be labelled using the benthic cover data collected during the field survey. The overall classification accuracy, calculated on the biological diverse fore reef, was 73% with a kappa coefficient of 0.63. A reliable bathymetric model (up to a depth of 15 m) of the Pulau Nukaha reef was also obtained using a semi‐analytical radiative transfer model. When compared with independent in‐situ depth measurements, the result proved relatively accurate (mean residual error: −0.9 m) and was consistent with the seabed topography (Pearson correlation coefficient: 86%).

Atmospheric correction plays a crucial role among the processing steps applied to remotely sensed hyperspectral data. Atmospheric correction comprises a group of procedures needed to remove atmospheric effects from observed spectra, i.e. the transformation from at-sensor radiances to at-surface radiances or reflectances. In this paper we present the different steps in the atmospheric correction process for APEX hyperspectral data as applied by the Central Data Processing Center (CDPC) at the Flemish Institute for Technological Research (VITO, Mol, Belgium). The MODerate resolution atmospheric TRANsmission program (MODTRAN) is used to determine the source of radiation and for applying the actual atmospheric correction. As part of the overall correction process, supporting algorithms are provided in order to derive MODTRAN configuration parameters and to account for specific effects, e.g. correction for adjacency effects, haze and shadow correction, and topographic BRDF correction. The methods and theory underlying these corrections and an example of an application are presented

Water body detection worldwide using spaceborne remote sensing is a challenging task. A global scale multi-temporal and multi-spectral image analysis method for water body detection was developed. The PROBA-V microsatellite has been fully operational since December 2013 and delivers daily near-global synthesis with a spatial resolution of 1 km and 333 m. The Red, Near-InfRared (NIR) and Short Wave InfRared (SWIR) bands of the atmospherically corrected 10-day synthesis images are first Hue, Saturation and Value (HSV) color transformed and subsequently used in a decision tree classification for water body detection. To minimize commission errors four additional data layers are used: the Normalized Difference Vegetation Index (NDVI), Water Body Potential Mask (WBPM), Permanent Glacier Mask (PGM) and Volcanic Soil Mask (VSM). Threshold values on the hue and value bands, expressed by a parabolic function, are used to detect the water bodies. Beside the water bodies layer, a quality layer, based on the water bodies occurrences, is available in the output product. The performance of the Water Bodies Detection Algorithm (WBDA) was assessed using Landsat 8 scenes over 15 regions selected worldwide. A mean Commission Error (CE) of 1.5% was obtained while a mean Omission Error (OE) of 15.4% was obtained for minimum Water Surface Ratio (WSR) = 0.5 and drops to 9.8% for minimum WSR = 0.6. Here, WSR is defined as the fraction of the PROBA-V pixel covered by water as derived from high spatial resolution images, e.g., Landsat 8. Both the CE = 1.5% and OE = 9.8% (WSR = 0.6) fall within the user requirements of 15%. The WBDA is fully operational in the Copernicus Global Land Service and products are freely available.

For maintaining the tidal waterways in the Scheldt basin, including the rivers Rupel and Durme and a large part of the Nete catchment, and for ecological monitoring of the mud flats, salt marshes and riverbank vegetation, the Flemish government needs detailed maps of these rivers and their bank structures. These maps indicate not only vegetation types, plant associations and sediment types but also hard structures, such as quays, locks, sluices and roads. Different remote sensing techniques were used to collect the data necessary to produce the required detailed maps. During the months of July and August 2007 an airborne flight campaign took place to collect hyperspectral and LiDAR data of the Scheldt basin and the Nete catchments. These rivers have a total length of about 240 km. The Airborne Imaging Spectrometer for Applications (AISA) Eagle sensor acquired hyperspectral data in 32 spectral bands covering the visible/near-infrared (VIS/NIR) part of the electromagnetic spectrum with a ground resolution of 1 m. A multiple binary classification algorithm based on Fisher's linear discriminant analysis (LDA) was used to map the salt marshes and riverbank vegetation. Ground truth information, that is vegetation and sediment types, together with their geographical locations collected around the time of the flight campaign, was used to train the classifier in the later classification step. Laser scanning was performed using the Riegl LMS-Q560. The LiDAR dataset obtained had a resolution of at least 1 point per m² and was used to produce a digital elevation model (DEM) that contains all elements of the terrain. From this DEM a digital terrain model (DTM) was derived by applying appropriate filtering techniques. The elevation models were used primarily to derive information on the height, slope and aspect of the banks and dikes, but they also served as expert knowledge in the classification of the mud flats and bank vegetation. Overall, this work illustrates how airborne hyperspectral and LiDAR data can be used to derive highly detailed maps of the sediments, vegetation and hard structures along tidal rivers in large river basins. It also shows how large datasets can be handled in an expert system, in combination with different classification techniques, to produce the required result and accuracy.

PROBA-V (PRoject for On-Board Autonomy–Vegetation) was launched in May-2013 as an operational continuation to the vegetation (VGT) instruments on-board the Système Pour l’Observation de la Terre (SPOT)-4 and -5 satellites. The first reprocessing campaign of the PROBA-V archive from Collection 0 (C0) to Collection 1 (C1) aims at harmonizing the time series, thanks to improved radiometric and geometric calibration and cloud detection. The evaluation of PROBA-V C1 focuses on (i) qualitative and quantitative assessment of the new cloud detection scheme; (ii) quantification of the effect of the reprocessing by comparing C1 to C0; and (iii) evaluation of the spatio-temporal stability of the combined SPOT/VGT and PROBA-V archive through comparison to METOP/advanced very high resolution radiometer (AVHRR). The PROBA-V C1 cloud detection algorithm yields an overall accuracy of 89.0%. Clouds are detected with very few omission errors, but there is an overdetection of clouds over bright surfaces. Stepwise updates to the visible and near infrared (VNIR) absolute calibration in C0 and the application of degradation models to the SWIR calibration in C1 result in sudden changes between C0 and C1 Blue, Red, and NIR TOC reflectance in the first year, and more gradual differences for short-wave infrared (SWIR). Other changes result in some bias between C0 and C1, although the root mean squared difference (RMSD) remains well below 1% for top-of-canopy (TOC) reflectance and below 0.02 for the normalized difference vegetation index (NDVI). Comparison to METOP/AVHRR shows that the recent reprocessing campaigns on SPOT/VGT and PROBA-V have resulted in a more stable combined time series.