Search

Upcoming satellite imaging spectroscopy missions will deliver spatiotemporal explicit data streams to be exploited for mapping vegetation properties, such as nitrogen (N) content. Within retrieval workflows for real-time mapping over agricultural regions, such crop-specific information products need to be derived precisely and rapidly. To allow fast processing, intelligent sampling schemes for training databases should be incorporated to establish efficient machine learning (ML) models. In this study, we implemented active learning (AL) heuristics using kernel ridge regression (KRR) to minimize and optimize a training database for variational heteroscedastic Gaussian processes regression (VHGPR) to estimate aboveground N content. Several uncertainty and diversity criteria were applied on a lookup table (LUT) composed of aboveground N content and corresponding hyperspectral reflectance simulated by the PROSAIL-PRO model. The best-performing AL criteria were Euclidian distance-based diversity (EBD) resulting in a reduction of the LUT training data set by 81% (50 initial samples plus 141 samples selected from a pool of 1000 samples). This reduced LUT was used for training VHGPR, which is not only a competitive algorithm but also provides uncertainty estimates. Validation against in situ N reference data provided excellent results with a root-mean-square error (RMSE) of 1.84 g/m2 and a coefficient of determination ($R^{2}$) of 0.92. Mapping aboveground N content over an agricultural region yielded reliable estimates and meaningful associated uncertainties. These promising results encourage the transfer of such hybrid workflows into space and time within the frame of future operational N monitoring from satellite imaging spectroscopy data. [ABSTRACT FROM AUTHOR]

The availability of Airborne Visible and Infrared Imaging Spectrometer-Next Generation (AVIRIS-NG) data has enormous possibilities for quantification of Leaf Chlorophyll Content (LCC). The present study used the AVIRIS-NG campaign site of Western India for generation and validation of new chlorophyll indices by denoising the AVIRIS-NG data. For validation, concurrent to AVIRIS-NG flight overpass, field samplings were performed. The acquired AVIRIS-NG was subjected to Spectral Angle Mapper (SAM) classifier for discriminating the crop types. Three smoothing techniques i.e., Fast-Fourier Transform (FFT), Mean and Savitzky-Golay filters were evaluated for their denoising capability. Raw and filtered data was used for developing new chlorophyll indices by optimizing AVIRIS-NG bands using VIs based on parametric regression algorithms. In total, 20 chlorophyll indices and corresponding 20 models were developed for mapping LCC in the area. SAM identified 17 crop types in the area, while FFT found to be the best for filtering. Performance of these models when checked based on Pearson correlation coefficient (${r}$) and Centered Root Mean Square Difference (CRMSD), indicated that LCC-CCI10 based on normalized difference type index formed through Near Infrared band and blue band is the best estimator of LCC $({r}_{\textit {cal}}=0.73,{r}_{\textit {val}}=0.66,CRMSD=4.97)$. The approach was also tested using AVIRIS-NG image of the year 2018, which also showed a promising correlation (${r} =0.704$ , CRSMD = 8.98, Bias = −0.5) between modeled and field LCC. [ABSTRACT FROM AUTHOR]

Physically based radiative transfer models (RTMs) are widely used in Earth observation to understand the radiation processes occurring on the Earth’s surface and their interactions with water, vegetation, and atmosphere. Through continuous improvements, RTMs have increased in accuracy and representativity of complex scenes at expenses of an increase in complexity and computation time, making them impractical in various remote sensing applications. To overcome this limitation, the common practice is to precompute large lookup tables (LUTs) 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 input–output points (nodes) to be included in an LUT. We present the gradient-based automatic LUT generator algorithm (GALGA), which relies on the notion of an acquisition function that incorporates: 1) the Jacobian evaluation of an RTM and 2) the information about the multivariate distribution of the current nodes. We illustrate the capabilities of GALGA in the automatic construction and optimization of MODTRAN-based LUTs of different dimensions of the input variables space. Our results indicate that when compared with a pseudorandom homogeneous distribution of the LUT nodes, GALGA reduces: 1) the LUT size by >24%; 2) the computation time by ~27%; and 3) the maximum interpolation relative errors by at least 10%. It is concluded that an automatic LUT design might benefit from the methodology proposed in GALGA to reduce interpolation errors and computation time in computationally expensive RTMs. [ABSTRACT FROM AUTHOR]

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 using interpolation and emulation: at canopy level, using PROSAIL; and at top-of-atmosphere level, using MODTRAN. Various interpolation (nearest-neighbor, inverse distance weighting, and piece-wice linear) and emulation [Gaussian process regression (GPR), kernel ridge regression, and neural networks] methods were evaluated against a dense reference LUT. In all experiments, the emulation methods clearly produced more accurate output spectra than classical interpolation methods. The GPR emulation performed up to ten times more accurately than the best performing interpolation method, and this with a speed that is competitive with the faster interpolation methods. It is concluded that emulation can function as a fast and more accurate alternative to commonly used interpolation methods for reconstructing RTM spectral data. [ABSTRACT FROM AUTHOR]

During the design phase of a satellite mission, end-to-end mission performance simulator (E2ES) tools allow scientists and engineers evaluating the mission concept, consolidating system technical requirements and analyzing the suitability of the implemented technical solutions and data processing algorithms. The generation of synthetic scenes is one of the core parts of an E2ES, providing scenes (ground truth) as would be observed by satellite instruments and used as reference against simulated retrieved mission products. An appropriate generation of the scene also allows assessing the performance of the ground data processing chain replacing real instrument data before the mission is in orbit, for which the fidelity of the scene generation is critical. This paper describes the design of a generic scene generator (GSG) with capabilities to generate complex 3-D synthetic scenes that combine the effects of surface, heterogeneity, topography and atmosphere, and viewing/illumination geometry. The proposed design allows generating consistent high spatial and spectral resolution top-of-atmosphere radiance scenes for multiple instruments based on the use of thematic maps, radiative transfer models, and reflectance databases. The described GSG was implemented within the FLuorescence EXplorer (FLEX) E2ES software tool and showed its capabilities to generate compatible scenes for the fluorescence imaging spectrometer, ocean and land color instrument, and sea and land surface temperature radiometer instruments of ESA's FLEX/Sentinel-3 tandem mission and to validate the fulfillment of the FLEX mission requirements. [ABSTRACT FROM AUTHOR]

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 complex synthetic scenes with an advanced modeling of the instrument behavior and the full processing scheme up to the final fluorescence product. The results derived from FLEX-E simulations indicate that the instrument and developed image processing algorithms are able to retrieve the sun-induced fluorescence with an accuracy below the 0.2 \textmW\cdot\textm^-2\cdot\textsr^-1\cdot \textnm^-1 mission requirement. It is expected that FLEX-E will not only optimize the FLEX retrieval algorithms and technical requirements, but also serve as the baseline for the ground processing implementation and testing of calibration/validation procedures. [ABSTRACT FROM AUTHOR]

Biophysical parameters such as leaf chlorophyll content (LCC) and leaf area index (LAI) are standard vegetation products that can be retrieved from Earth observation imagery. This paper introduces a new machine learning regression algorithms (MLRAs) toolbox into the scientific Automated Radiative Transfer Models Operator (ARTMO) software package. ARTMO facilitates retrieval of biophysical parameters from remote observations in a MATLAB graphical user interface (GUI) environment. The MLRA toolbox enables analyzing the predictive power of various MLRAs in a semiautomatic and systematic manner, and applying a selected MLRA to multispectral or hyperspectral imagery for mapping applications. It contains both linear and nonlinear state-of-the-art regression algorithms, in particular linear feature extraction via principal component regression (PCR), partial least squares regression (PLSR), decision trees (DTs), neural networks (NNs), kernel ridge regression (KRR), and Gaussian processes regression (GPR). The performance of multiple implemented regression strategies has been evaluated against the SPARC dataset (Barrax, Spain) and simulated Sentinel-2 (8 bands), CHRIS (62 bands) and HyMap (125 bands) observations. In general, nonlinear regression algorithms (NN, KRR, and GPR) outperformed linear techniques (PCR and PLSR) in terms of accuracy, bias, and robustness. Most robust results along gradients of training/validation partitioning and noise variance were obtained by KRR while GPR delivered most accurate estimations. We applied a GPR model to a hyperspectral HyMap flightline to map LCC and LAI. We exploited the associated uncertainty intervals to gain insight in the per-pixel performance of the model. [ABSTRACT FROM PUBLISHER]

An accurate estimation of biophysical variables is the key to monitor our Planet. Leaf chlorophyll content helps in interpreting the chlorophyll fluorescence signal from space, whereas oceanic chlorophyll concentration allows us to quantify the healthiness of the oceans. Recently, the family of Bayesian nonparametric methods has provided excellent results in these situations. A particularly useful method in this framework is the Gaussian process regression (GPR). However, standard GPR assumes that the variance of the noise process is independent of the signal, which does not hold in most of the problems. In this letter, we propose a nonstandard variational approximation that allows accurate inference in signal-dependent noise scenarios. We show that the so-called variational heteroscedastic GPR (VHGPR) is an excellent alternative to standard GPR in two relevant Earth observation examples, namely, Chl vegetation retrieval from hyperspectral images and oceanic Chl concentration estimation from in situ measured reflectances. The proposed VHGPR outperforms the tested empirical approaches, as well as statistical linear regression (both least squares and least absolute shrinkage and selection operator), neural nets, and kernel ridge regression, and the homoscedastic GPR, in terms of accuracy and bias, and proves more robust when a low number of examples is available. [ABSTRACT FROM PUBLISHER]

Inversion of radiative transfer models (RTM) using a lookup-table (LUT) approach against satellite reflectance data can lead to concurrent retrievals of biophysical parameters such as leaf chlorophyll content (Chl) and leaf area index (LAI), but optimization strategies are not consolidated yet. ESA's upcoming satellites Sentinel-2 (S2) and Sentinel-3 (S3) aim to ensure continuity of old generation satellite sensors by providing superspectral images of high spatial and temporal resolution. This unprecedented data availability leads to an urgent need for developing robust, accurate, and operational retrieval methods. For three simulated Sentinel settings (S2-10 m: 4 bands, S2-20 m: 8 bands and S3-OLCI: 19 bands) various optimization strategies in LUT-based RTM inversion have been evaluated, being the role of i) added noise, ii) multiple best solutions, iii) combined parameters (Chl × LAI), and iv) applied cost functions. By inverting the PROSAIL model and using data from the ESA-led field campaign SPARC (Barrax, Spain), it was demonstrated that introducing noise and opting for multiple best solutions in the inversion considerably improved retrievals. However, the widely used RMSE was not the best performing cost function. Three families of alternative cost functions were applied here: information measures, minimum contrast, and M-estimates. We found that so-called "Power divergence measure", "Trigonometric", and spectral measure with "Contrast function K(x) = -log(x) + x", yielded more accurate results, although this also depended on the biophysical parameter. Particularly, when simultaneous retrieval of multiple biophysical parameters is the objective then "Contrast function K(x) = -log(x) + x" provided most consistent optimized estimates of leaf Chl, LAI and canopy Chl across the different Sentinel configurations (relative RMSE: 24-29%). [ABSTRACT FROM AUTHOR]

Precise and spatially-explicit knowledge of leaf chlorophyll content (Chl) is crucial to adequately interpret the chlorophyll fluorescence (ChF) signal from space. Accompanying information about the reliability of the Chl estimation becomes more important than ever. Recently, a new statistical method was proposed within the family of nonparametric Bayesian statistics, namely Gaussian Processes regression (GPR). GPR is simpler and more robust than their machine learning family members while maintaining very good numerical performance and stability. Other features include: i) GPR requires a relatively small training data set and can adopt very flexible kernels, ii) GPR identifies the relevant bands and observations in establishing relationships with a variable, and finally iii) along with pixelwise estimations GPR provides accompanying confidence intervals. We used GPR to retrieve Chl from hyperspectral reflectance data and evaluated the portability of the regression model to other images. Based on field Chl measurements from the SPARC dataset and corresponding spaceborne CHRIS spectra (acquired in 2003, Barrax, Spain), GPR developed a regression model that was excellently validated (r^2: 0.96, RMSE: 3.82 \mug/cm^2). The SPARC-trained GPR model was subsequently applied to CHRIS images (Barrax, 2003, 2009) and airborne CASI flightlines (Barrax 2009) to generate Chl maps. The accompanying confidence maps provided insight in the robustness of the retrievals. Similar confidences were achieved by both sensors, which is encouraging for upscaling Chl estimates from field to landscape scale. Because of its robustness and ability to deliver confidence intervals, GPR is evaluated as a promising candidate for implementation into ChF processing chains. [ABSTRACT FROM PUBLISHER]

This letter proposes a multioutput support vector regression (M-SVR) method for the simultaneous estimation of different biophysical parameters from remote sensing images. General retrieval problems require multioutput (and potentially nonlinear) regression methods. M-SVR extends the single-output SVR to multiple outputs maintaining the advantages of a sparse and compact solution by using an \varepsilon-insensitive cost function. The proposed M-SVR is evaluated in the estimation of chlorophyll content, leaf area index and fractional vegetation cover from a hyperspectral compact high-resolution imaging spectrometer images. The achieved improvement with respect to the single-output regression approach suggests that M-SVR can be considered a convenient alternative for nonparametric biophysical parameter estimation and model inversion. [ABSTRACT FROM PUBLISHER]

The Compact High Resolution Imaging Spectrometer (CHRIS) mounted onboard the Project for Onboard Autonomy (PROBA) spacecraft is capable of sampling reflected radiation at five viewing angles over the visible and near-infrared regions of the solar spectrum with high spatial resolution. We combined the spectral domain with the angular domain of CHRIS data in order to map the surface heterogeneity of an Alpine coniferous forest during winter. In the spectral domain, linear spectral unmixing of the nadir image resulted in a canopy cover map. In the angular domain, pixelwise inversion of the Rahman-Pinty-Verstraete (RPV) model at a single wavelength at the red edge (722 nm) yielded a map of the Minnaert-k parameter that provided information on surface heterogeneity at a subpixel scale. However, the interpretation of the Minnaert- k parameter is not always straightforward because fully vegetated targets typically produce the same type of reflectance anisotropy as non-vegetated targets. Merging both maps resulted in a forest cover heterogeneity map, which contains more detailed information on canopy heterogeneity at the CHRIS subpixel scale than is possible to realize from a single-source optical data set. [ABSTRACT FROM AUTHOR]