The unprecedented availability of optical satellite data in cloud-based computing platforms, such as Google Earth Engine (GEE), opens new possibilities to develop crop trait retrieval models from the local to the planetary scale. Hybrid retrieval models are of interest to run in these platforms as they combine the advantages of physically- based radiative transfer models (RTM) with the flexibility of machine learning regression algorithms. Previous research with GEE primarily relied on processing bottom-of-atmosphere (BOA) reflectance data, which requires atmospheric correction. In the present study, we implemented hybrid models directly into GEE for processing Sentinel-2 (S2) Level-1C (L1C) top-of-atmosphere (TOA) reflectance data into crop traits. To achieve this, a training dataset was generated using the leaf-canopy RTM PROSAIL in combination with the atmospheric model 6SV. Gaussian process regression (GPR) retrieval models were then established for eight essential crop traits namely leaf chlorophyll content, leaf water content, leaf dry matter content, fractional vegetation cover, leaf area index (LAI), and upscaled leaf variables (i.e., canopy chlorophyll content, canopy water content and canopy dry matter content). An important pre-requisite for implementation into GEE is that the models are sufficiently light in order to facilitate efficient and fast processing. Successful reduction of the training dataset by 78% was achieved using the active learning technique Euclidean distance-based diversity (EBD). With the EBD-GPR models, highly accurate validation results of LAI and upscaled leaf variables were obtained against in situ field data from the validation study site Munich-North-Isar (MNI), with normalized root mean square errors (NRMSE) from 6% to 13%. Using an independent validation dataset of similar crop types (Italian Grosseto test site), the retrieval models showed moderate to good performances for canopy-level variables, with NRMSE ranging from 14% to 50%, but failed for the leaf-level estimates. Obtained maps over the MNI site were further compared against Sentinel-2 Level 2 Prototype Processor (SL2P) vegetation estimates generated from the ESA Sentinels’ Application Platform (SNAP) Biophysical Processor, proving high consistency of both retrievals (R² from 0.80 to 0.94). Finally, thanks to the seamless GEE processing capability, the TOA-based mapping was applied over the entirety of Germany at 20 m spatial resolution including information about prediction uncertainty. The obtained maps provided confidence of the developed EBD-GPR retrieval models for integration in the GEE framework and national scale mapping from S2-L1C imagery. In summary, the proposed retrieval workflow demonstrates the possibility of routine processing of S2 TOA data into crop traits maps at any place on Earth as required for operational agricultural applications.

In this letter, a coherence-based technique for atmospheric artifact removal in ground-based (GB) zero-baseline synthetic aperture radar (SAR) acquisitions is proposed. For this purpose, polarimetric measurements acquired using the GB-SAR sensor developed at the Universitat Politecnica de Catalunya are employed. The heterogeneous environment of Collserola Park in the outskirts of Barcelona, Spain, was selected as the test area. Data sets were acquired at X-band during one week in June 2005. The effects of the atmosphere variations between successive zero-baseline SAR polarimetric acquisitions are treated here in detail. The need to compensate for the resulting phase-difference errors when retrieving interferometric information is put forward. A compensation technique is then proposed and evaluated using the control points placed inside the observed scene.

The integration of rooftop greenhouses (RTGs) in urban buildings is a practice that is becoming increasingly important in the world for their contribution to food security and sustainable development. However, the supply of tools and procedures to facilitate their implementation at the city scale is limited and laborious. This work aims to develop a specific and automated methodology for identifying the feasibility of implementation of rooftop greenhouses in non-residential urban areas, using airborne sensors. The use of Light Detection and Ranging (LIDAR) and Long Wave Infrared (LWIR) data and the Leica ALS50-II and TASI-600 sensors allow for the identification of some building roof parameters (area, slope, materials, and solar radiation) to determine the potential for constructing a RTG. This development represents an improvement in time and accuracy with respect to previous methodology, where all the relevant information must be acquired manually. The methodology has been applied and validated in a case study corresponding to a non-residential urban area in the industrial municipality of Rubí, Barcelona (Spain). Based on this practical application, an area of 36,312m2 out of a total area of 1,243,540m2 of roofs with ideal characteristics for the construction of RTGs was identified. This area can produce approximately 600tons of tomatoes per year, which represents the average yearly consumption for about 50% of Rubí total population. The use of this methodology also facilitates the decision making process in urban agriculture, allowing a quick identification of optimal surfaces for the future implementation of urban agriculture in housing. It also opens new avenues for the use of airborne technology in environmental topics in cities.

For the last decade, Gaussian process regression (GPR) proved to be a competitive machine learning regression algorithm for Earth observation applications, with attractive unique properties such as band relevance ranking and uncertainty estimates. More recently, GPR also proved to be a proficient time series processor to fill up gaps in optical imagery, typically due to cloud cover. This makes GPR perfectly suited for large-scale spatiotemporal processing of satellite imageries into cloud-free products of biophysical variables. With the advent of the Google Earth Engine (GEE) cloud platform, new opportunities emerged to process local-to-planetary scale satellite data using advanced machine learning techniques and convert them into gap-filled vegetation properties products. However, GPR is not yet part of the GEE ecosystem. To circumvent this limitation, this work proposes a general adaptation of GPR formulation to parallel processing framework and its integration into GEE. To demonstrate the functioning and utility of the developed workflow, a GPR model predicting green leaf area index (LAI${}_G$) from Sentinel-2 imagery was imported. Although by running this GPR model into GEE any corner of the world can be mapped into LAI${}_G$ at a resolution of 20 m, here we show some demonstration cases over western Europe with zoom-ins over Spain. Thanks to the computational power of GEE, the mapping takes place on-the-fly. Additionally, a GPR-based gap filling strategy based on pre-optimized kernel hyperparameters is also put forward for the generation of multi-orbit cloud-free LAI${}_G$ maps with an unprecedented level of detail, and the extraction of regularly-sampled LAI${}_G$ time series at a pixel level. The ability to plugin a locally-trained GPR model into the GEE framework and its instant processing opens up a new paradigm of remote sensing image processing.