FLEX End-to-End Mission Performance Simulator
The FLuorescence EXplorer (FLEX) mission, selected as the European Space Agency's eighth Earth Explorer, aims to globally measure the sun-induced-chlorophyll-fluorescence spectral emission from terrestrial vegetation. In the frame of the FLEX mission, several industrial and scientific studies have analyzed the instrument design, image processing algorithms, or modeling aspects. At the same time, a common tool is needed to address the overall FLEX mission performance by combining all these features. For this reason, an end-to-end mission performance simulator has been developed for the FLEX mission (FLEX-E). This paper describes the FLEX-E software design, which combines the generation of co…
Global Sensitivity Analysis of Leaf-Canopy-Atmosphere RTMs: Implications for Biophysical Variables Retrieval from Top-of-Atmosphere Radiance Data.
Knowledge of key variables driving the top of the atmosphere (TOA) radiance over a vegetated surface is an important step to derive biophysical variables from TOA radiance data, e.g., as observed by an optical satellite. Coupled leaf-canopy-atmosphere Radiative Transfer Models (RTMs) allow linking vegetation variables directly to the at-sensor TOA radiance measured. Global Sensitivity Analysis (GSA) of RTMs enables the computation of the total contribution of each input variable to the output variance. We determined the impacts of the leaf-canopy-atmosphere variables into TOA radiance using the GSA to gain insights into retrievable variables. The leaf and canopy RTM PROSAIL was coupled with…
HICO L1 and L2 data processing: Radiometric recalibration, atmospheric correction and retrieval of water quality parameters
The Hyperspectral Imager for the Coastal Ocean (HICO) is an imaging spectrometer designed with a very high signal-to-noise ratio to monitor coastal ocean and inland waters. The processing of Top-Of-Atmosphere radiance data down to surface reflectance is fundamental for the retrieval of water quality products. However, the current HICO processing chain does not provide atmospheric corrected data nor higher-level water quality products. This paper describes the algorithms implemented within an HICO data processing chain that includes image pre-processing, atmospheric correction and the retrieval of water quality parameters. The implemented algorithms have been validated over a set of HICO ima…
Automatic emulator and optimized look-up table generation for radiative transfer models
This paper introduces an automatic methodology to construct emulators for costly radiative transfer models (RTMs). The proposed method is sequential and adaptive, and it is based on the notion of the acquisition function by which instead of optimizing the unknown RTM underlying function we propose to achieve accurate approximations. The Automatic Gaussian Process Emulator (AGAPE) methodology combines the interpolation capabilities of Gaussian processes (GPs) with the accurate design of an acquisition function that favors sampling in low density regions and flatness of the interpolation function. We illustrate the good capabilities of the method in toy examples and for the construction of an…
HICO level-2 data processing toolbox for the atmospheric correction and the retrieval of water quality parameters
The Hyperspectral Imager for the Coastal Ocean (HICO) is an imaging spectrometer specifically designed to monitor the coastal ocean. The processing of Top-Of-Atmosphere (TOA) radiance data down to surface reflectance is fundamental for the retrieval of water quality products. However, the current HICO processing chain does not provide atmospheric corrected data nor higher-level water quality products. This work describes a toolbox for the atmospheric correction of HICO data and the retrieval of water quality products. The HICO toolbox, consisting on three main modules (image pre-processing, atmospheric correction and retrieval of water quality products), has been used over a set of HICO ima…
Oxygen transmittance correction for solar-induced chlorophyll fluorescence measured on proximal sensing: Application to the NASA-GSFC fusion tower
Since oxygen (O 2 ) absorption of light becomes more pronounced at higher pressure levels, even a few meters distance between the target and the sensor can strongly affect canopy-leaving Solar-Induced chlorophyll Fluorescence (SIF) retrievals. This study was conducted to quantify the consequent error propagation and the impact of ignoring oxygen absorption effects on proximal sensing SIF measurements based on the O 2 -A absorption band with field-acquired and simulated data. It was demonstrated that the uncorrected oxygen transmittance between target and sensor distance of 10 m can lead to SIF relative errors ranging from 66% to higher than 100% when using a Spectral Fitting (SF) technique …
Synthetic scene simulator for hyperspectral spaceborne passive optical sensors. Application to ESA's FLEX/sentinel-3 tandem mission
The simulation of synthetic images serve scientists and engineers to study the instrument configuration as well as to develop image processing and retrieval strategies for a sensor in development. Despite synthetic scene simulators have been developed in the past in the frame of satellite missions, their functionality and flexibility to create a user-defined scene is limited by their architecture, design and implementation. This paper introduces the design of a generic scene simulator with the flexibility to generate realistic synthetic scenes by configuration of the surface and atmosphere. Following this generic design, a scene simulator is being developed for the ESA's Earth Explorer 8th …
Automatic Emulation by Adaptive Relevance Vector Machines
This paper introduces an automatic methodology to construct emulators for costly radiative transfer models (RTMs). The proposed method is sequential and adaptive, and it is based on the notion of the acquisition function by which instead of optimizing the unknown RTM underlying function we propose to achieve accurate approximations. The proposed methodology combines the interpolation capabilities of a modified Relevance Vector Machine (RVM) with the accurate design of an acquisition function that favors sampling in low density regions and flatness of the interpolation function. The proposed Relevance Vector Machine Automatic Emulator (RAE) is illustrated in toy examples and for the construc…
Gradient-based Automatic Look-Up Table Generator for Atmospheric Radiative Transfer Models
Atmospheric correction of Earth Observation data is one of the most critical steps in the data processing chain of a satellite mission for successful remote sensing applications. Atmospheric Radiative Transfer Models (RTM) inversion methods are typically preferred due to their high accuracy. However, the execution of RTMs on a pixel-per-pixel basis is impractical due to their high computation time, thus large multi-dimensional look-up tables (LUTs) are precomputed for their later interpolation. To further reduce the RTM computation burden and the error in LUT interpolation, we have developed a method to automatically select the minimum and optimal set of nodes to be included in a LUT. We pr…
Atmospheric and Instrumental Effects on the Fluorescence Remote Sensing Retrieval
Accurately disentangling the tiny Solar–Induced Chlorophyll Fluorescence (SIF) from canopy reflected solar irradiance by using passive remote sensing techniques is always challenging. Regardless the scale at which SIF is measured, i.e., proximal sensing, airborne or satellite level; instrumental and atmospheric effects must be accounted for and compensated as part of the SIF retrieval strategy. Regarding the instrumental effects, the use of very high spectral resolution spectrometers makes mandatory an accurate characterization of the Instrument Spectral Response Function (ISRF); and – in the case of imager spectrometers – an accurate characterization of the full instrument response in the …
FLEX/S3 Tandem Mission Performance Assessment: Evolution of the End-to-End Simulator Flex-E
An End-to-end simulator (E2ES) is a tool to evaluate the performance of a satellite mission. Once a mission is approved for operation, E2ES evolves during Phase C/D to become a supporting tool for the development and validation of the ground data processor, as well as for simulating data sets to test the Prototype and Operational Processors. FLEX-E is the E2ES of the FLEX/Sentinel-3 tandem mission, which was selected in 2015 as ESA's eighth Earth Explorer. The FLEX-E evolution implies the consolidation of all the retrieval algorithms (e.g. fluorescence, reflectance, biophysical variables), the implementation of new scientific developments, as well the improvement of the co-registration proc…
Misión FLEX (Fluorescence Explorer): Observación de la fluorescencia por teledetección como nueva técnica de estudio del estado de la vegetación terrestre a escala global
[EN] FLEX (Fluorescence EXplorer) is a candidate for the 8th ESA’s Earth Explorer mission. Is the first space mission specifically designed for the estimation of vegetation fluorescence on a global scale. The mission is proposed to fly in tandem with the future ESA´s Sentinel-3 satellite. It is foreseen that the information obtained by Sentinel-3 will be supplemented with that provided by FLORIS (Fluorescence Imaging Spectrometer) onboard FLEX. FLORIS will measure the radiance between 500 and 800 nm with a bandwidth between 0.1 nm and 2 nm, providing images with a 150 km swath and 300 m pixel size. This information will allow a detailed monitoring of vegetation dynamics, by improving the me…
Improving the analysis of biogeochemical patterns associated with internal waves in the strait of Gibraltar using remote sensing images
High Amplitude Internal Waves (HAIWs) are physical processes observed in the Strait of Gibraltar (the narrow channel between the Atlantic Ocean and the Mediterranean Sea). These internal waves are generated over the Camarinal Sill (western side of the strait) during the tidal outflow (toward the Atlantic Ocean) when critical hydraulic conditions are established. HAIWs remain over the sill for up to 4 h until the outflow slackens, being then released (mostly) towards the Mediterranean Sea. These have been previously observed using Synthetic Aperture Radar (SAR), which captures variations in surface water roughness. However, in this work we use high resolution optical remote sensing, with the…
Gaussian Processes Retrieval of LAI from Sentinel-2 Top-of-Atmosphere Radiance Data
Abstract Retrieval of vegetation properties from satellite and airborne optical data usually takes place after atmospheric correction, yet it is also possible to develop retrieval algorithms directly from top-of-atmosphere (TOA) radiance data. One of the key vegetation variables that can be retrieved from at-sensor TOA radiance data is leaf area index (LAI) if algorithms account for variability in atmosphere. We demonstrate the feasibility of LAI retrieval from Sentinel-2 (S2) TOA radiance data (L1C product) in a hybrid machine learning framework. To achieve this, the coupled leaf-canopy-atmosphere radiative transfer models PROSAIL-6SV were used to simulate a look-up table (LUT) of TOA radi…
Statistical Learning for End-to-End Simulations
End-to-end mission performance simulators (E2ES) are suitable tools to accelerate satellite mission development from concet to deployment. One core element of these E2ES is the generation of synthetic scenes that are observed by the various instruments of an Earth Observation mission. The generation of these scenes rely on Radiative Transfer Models (RTM) for the simulation of light interaction with the Earth surface and atmosphere. However, the execution of advanced RTMs is impractical due to their large computation burden. Classical interpolation and statistical emulation methods of pre-computed Look-Up Tables (LUT) are therefore common practice to generate synthetic scenes in a reasonable…
Top-of-Atmosphere Retrieval of Multiple Crop Traits Using Variational Heteroscedastic Gaussian Processes within a Hybrid Workflow.
In support of cropland monitoring, operational Copernicus Sentinel-2 (S2) data became available globally and can be explored for the retrieval of important crop traits. Based on a hybrid workflow, retrieval models for six essential biochemical and biophysical crop traits were developed for both S2 bottom-of-atmosphere (BOA) L2A and S2 top-of-atmosphere (TOA) L1C data. A variational heteroscedastic Gaussian process regression (VHGPR) algorithm was trained with simulations generated by the combined leaf-canopy reflectance model PROSAILat the BOA scale and further combined with the Second Simulation of a Satellite Signal in the Solar Spectrum (6SV) atmosphere model at the TOA scale. Establishe…
Multioutput Automatic Emulator for Radiative Transfer Models
This paper introduces a methodology to construct emulators of costly radiative transfer models (RTMs). The proposed methodology is sequential and adaptive, and it is based on the notion of acquisition functions in Bayesian optimization. Here, instead of optimizing the unknown underlying RTM function, one aims to achieve accurate approximations. The Automatic Multi-Output Gaussian Process Emulator (AMO-GAPE) methodology combines the interpolation capabilities of Gaussian processes (GPs) with the accurate design of an acquisition function that favors sampling in low density regions and flatness of the interpolation function. We illustrate the promising capabilities of the method for the const…
A sun-induced vegetation fluorescence retrieval method from top of atmosphere radiance for the FLEX/Sentinel-3 TanDEM mission
A new fluorescence retrieval method is proposed to support ESA's 8th Earth Explorer FLuorescence EXplorer/Sentinel-3 (FLEX-S3) candidate tandem mission. FLEX is the first mission specially dedicated to measure the Sun-Induced vegetation chlorophyll fluorescence (SIF) strongly related with the vegetation photosynthetic activity. Most hyperspectral fluorescence retrieval algorithms available in the literature are very sensitive to true reflectance modelization and/or they assume the atmospheric status as known. The proposed algorithm delivers the retrieval of full fluorescence spectrum at canopy level by using only Top Of Atmosphere (TOA) radiances from S3 and FLEX as input. Once the spatial …
A fluorescence retrieval method for the flex sentinel-3 tandem mission
A new fluorescence retrieval method is proposed to support ESA's 8th Earth Explorer Fluorescence EXplorer (FLEX) candidate mission. Most hyperspectral fluorescence retrieval algorithms available in the literature are very sensitive to true reflectance modelization and/or they assume the atmospheric status as known. The proposed algorithm delivers the retrieval of full fluorescence spectrum at canopy level by using only Top Of Atmosphere (TOA) radiances as input. The proposed method starts with (1) the atmospheric correction of TOA radiances, characterizing the state of the atmosphere without assuming any a-priori classification on aerosols models, (2) performing a first estimation of fluore…
The flex end-to-end simulator: From concept phase (A/B1) to ground segment and operations (C/D)
ESA's FLEX/Sentinel-3 tandem mission aims at mapping Sun-induced fluorescence (SIF) as a proxy to quantify photosynthetic activity of terrestrial vegetation. Due to the complexity of the mission concept and stringent requirements for the data processing algorithms, ESA developed a Phase A/B1 End-to-End Mission Performance Simulator (E2ES) tool to reproduce the expected mission performance and check the mission and instrument concepts. In the current Phase C/D, the E2ES concept must evolve to consolidate the whole data processing chain, providing an accurate figures of the whole mission error budget and serving as a roadmap for the future development of FLEX Ground Segment. This paper gives …
Compensation of Oxygen Transmittance Effects for Proximal Sensing Retrieval of Canopy–Leaving Sun–Induced Chlorophyll Fluorescence
Estimates of Sun–Induced vegetation chlorophyll Fluorescence (SIF) using remote sensing techniques are commonly determined by exploiting solar and/or telluric absorption features. When SIF is retrieved in the strong oxygen (O 2 ) absorption features, atmospheric effects must always be compensated. Whereas correction of atmospheric effects is a standard airborne or satellite data processing step, there is no consensus regarding whether it is required for SIF proximal–sensing measurements nor what is the best strategy to be followed. Thus, by using simulated data, this work provides a comprehensive analysis about how atmospheric effects impact SIF estimations on proximal sensing, regarding: (…
Retrieval of sun-induced fluorescence using advanced spectral fitting methods
Abstract The FLuorescence EXplorer (FLEX) satellite mission, candidate of ESA's 8th Earth Explorer program, is explicitly optimized for detecting the sun-induced fluorescence emitted by plants. It will allow consistent measurements around the O2-B (687 nm) and O2-A (760 nm) bands, related to the red and far-red fluorescence emission peaks respectively, the photochemical reflectance index, and the structural-chemical state variables of the canopy. The sun-induced fluorescence signal, overlapped to the surface reflected radiance, can be accurately retrieved by employing the powerful spectral fitting technique. In this framework, a set of fluorescence retrieval algorithms optimized for FLEX ar…
Approximating Empirical Surface Reflectance Data through Emulation: Opportunities for Synthetic Scene Generation
Collection of spectroradiometric measurements with associated biophysical variables is an essential part of the development and validation of optical remote sensing vegetation products. However, their quality can only be assessed in the subsequent analysis, and often there is a need for collecting extra data, e.g., to fill in gaps. To generate empirical-like surface reflectance data of vegetated surfaces, we propose to exploit emulation, i.e., reconstruction of spectral measurements through statistical learning. We evaluated emulation against classical interpolation methods using an empirical field dataset with associated hyperspectral spaceborne CHRIS and airborne HyMap reflectance spectra…
Emulation as an Accurate Alternative to Interpolation in Sampling Radiative Transfer Codes
Computationally expensive radiative transfer models (RTMs) are widely used to realistically reproduce the light interaction with the earth surface and atmosphere. Because these models take long processing time, the common practice is to first generate a sparse look-up table (LUT) and then make use of interpolation methods to sample the multidimensional LUT input variable space. However, the question arise whether common interpolation methodsperform most accurate. As an alternative to interpolation, this paper proposes to use emulation, i.e., approximating the RTM output by means of the statistical learning. Two experiments were conducted to assess the accuracy in delivering spectral outputs…
Adaptive Sequential Interpolator Using Active Learning for Efficient Emulation of Complex Systems
Many fields of science and engineering require the use of complex and computationally expensive models to understand the involved processes in the system of interest. Nevertheless, due to the high cost involved, the required study becomes a cumbersome process. This paper introduces an interpolation procedure which belongs to the family of active learning algorithms, in order to construct cheap surrogate models of such costly complex systems. The proposed technique is sequential and adaptive, and is based on the optimization of a suitable acquisition function. We illustrate its efficiency in a toy example and for the construction of an emulator of an atmosphere modeling system.
Alg: a Toolbox for the Generation of Look-Up tables Based on Atmospheric Radiative Transfer Models
Atmospheric radiative transfer models (RTMs) are software tools describing the radiation processes occurring on the Earth’s atmosphere. While the evolution of these tools have achieved better representations of the light-atmosphere interactions, the increase of complexity, interpretability and computation time bears implications towards practical applications in Earth observation. Despite of existing RTM-specific graphical user interfaces, none of these tools allow common streamlining model setup for a wide variety of atmospheric RTMs. In addition, the automatic generation of atmospheric look-up tables (LUTs) can hardly be done with the use of these graphical tools. This paper presents the …