ALBERT

Description: Land-cover monitoring is one of the core applications of remote sensing. Monitoring and mapping changes in the distribution of agricultural land covers provide a reliable source of information that helps environmental sustainability and supports agricultural policies. Synthetic Aperture Radar (SAR) can contribute considerably to this monitoring effort. The first objective of this research is to extend the use of time series of polarimetric data for land-cover classification using a decision tree classification algorithm. With this aim, RADARSAT-2 (quad-pol) and Sentinel-1 (dual-pol) data were acquired over an area of 600 km2 in central Spain. Ten polarimetric observables were derived from both datasets and seven scenarios were created with different sets of observables to evaluate a multitemporal parcel-based approach for classifying eleven land-cover types, most of which were agricultural crops. The study demonstrates that good overall accuracies, greater than 83%, were achieved for all of the different proposed scenarios and the scenario with all RADARSAT-2 polarimetric observables was the best option (89.1%). Very high accuracies were also obtained when dual-pol data from RADARSAT-2 or Sentinel-1 were used to classify the data, with overall accuracies of 87.1% and 86%, respectively. In terms of individual crop accuracy, rapeseed achieved at least 95% of a producer’s accuracy for all scenarios and that was followed by the spring cereals (wheat and barley), which achieved high producer’s accuracies (79.9%-95.3%) and user’s accuracies (85.5% and 93.7%).

Description: This paper proposes a photogrammetric procedure able to determine out-of-plane movements experienced by a masonry structure subjected to a quasi-static cyclic test. The method tracks the movement of circular targets by means of a coarse-to-fine strategy. These targets were captured by means of a photogrammetric network, made up of four cameras optimized following the precepts of a zero-, first-, and second-order design. The centroid of each circular target was accurately detected for each image using the Hough transform, a sub-pixel edge detector based on the partial area effect, and a non-linear square optimization strategy. The three-dimensional (3D) coordinates of these targets were then computed through a photogrammetric bundle adjustment considering a self-calibration model of the camera. To validate the photogrammetric method, measurements were carried out in parallel to an ongoing test on a full-scale two-story unreinforced masonry structure (5.4 × 5.2 × 5.4-m) monitored with more than 200 contact sensors. The results provided by the contact sensors during one of the load phases were compared with those obtained by the proposed approach. According to this accuracy assessment, the method was able to determine the out-of-plane displacement during the quasi-static cyclic test with a sub-pixel accuracy of 0.58.

Description: Over the last two decades, a wide range of applications have been developed from Light Detection and Ranging (LiDAR) point clouds. Most LiDAR-derived products require the distinction between ground and non-ground points. Because of this, ground filtering its being one of the most studied topics in the literature and robust methods are nowadays available. However, these methods have been designed to work with offline data and they are generally not well suited for real-time scenarios. Aiming to address this issue, this paper proposes an efficient method for ground filtering of airborne LiDAR data based on scan-line processing. In our proposal, an iterative 1-D spline interpolation is performed in each scan line sequentially. The final spline knots of a scan line are taken into account for the next scan line, so that valuable 2-D information is also considered without compromising computational efficiency. Points are labelled into ground and non-ground by analysing their residuals to the final spline. When tested against synthetic ground truth, the method yields a mean kappa value of 88.59% and a mean total error of 0.50%. Experiments with real data also show satisfactory results under visual inspection. Performance tests on a workstation show that the method can process up to 1 million points per second. The original implementation was ported into a low-cost development board to demonstrate its feasibility to run in embedded systems, where throughput was improved by using programmable logic hardware acceleration. Analysis shows that real-time filtering is possible in a high-end board prototype, as it can process the amount of points per second that current lightweight scanners acquire with low-energy consumption.

Description: In viticulture, deficit irrigation strategies are often implemented to control vine canopy growth and to impose stress at critical stages of vine growth to improve wine grape quality. To support deficit irrigation scheduling, remote sensing technologies can be employed in the mapping of evapotranspiration (ET) at the field to sub-field scales, quantifying time-varying vineyard water requirements and actual water use. In the current study, we investigate the utility of ET maps derived from thermal infrared satellite imagery over a vineyard in the Central Valley of California equipped with a variable rate drip irrigation (VRDI) system which enables differential water applications at the 30 × 30 m scale. To support irrigation management at that scale, we utilized a thermal-based multi-sensor data fusion approach to generate weekly total actual ET (ETa) estimates at 30 m spatial resolution, coinciding with the resolution of the Landsat reflectance bands. Crop water requirements (ETc) were defined with a vegetative index (VI)-based approach. To test capacity to capture stress signals, the vineyard was sub-divided into four blocks with different irrigation management strategies and goals, inducing varying degrees of stress during the growing season. Results indicate derived weekly total ET from the thermal-based data fusion approach match well with observations. The thermal-based method was also able to capture the spatial heterogeneity in ET over the vineyard due to a water stress event imposed on two of the four vineyard blocks. This transient stress event was not reflected in the VI-based ETc estimate, highlighting the value of thermal band imaging. While the data fusion system provided valuable information, latency in current satellite data availability, particularly from Landsat, impacts operational applications over the course of a growing season.

Description: Transport networks need periodic inspections to increase their safety and improve their management. In the last few years, LiDAR (light detection and ranging) technology has become a tool for helping to create a precise database of almost any type of infrastructure. Mobile laser scanning (MLS) systems use a laser beam to collect dense three dimensional (3D) point clouds, which include geometric and radiometric data of the environment in which they are placed. In the context of this paper, a methodology for automatically inspecting the clearance gauge and the deflection of the aerial contact line in railway tunnels is presented. The main objective is to compare results and verify their compliance with the Spanish norm. The 3D data are provided by a LYNX Mobile Mapper System (MMS). First, the area is surveyed and then the obtained (3D) point cloud is classified into contact wire, suspension wire, and remaining points. Finally, the inspection of the railway’s power line is performed. The validation of the proposed methodology has been carried out in three different tunnel point clouds, obtaining both qualitative and quantitative results for points’ classification, together with the results of the measures performed.

Description: It is well recognized that exposure to fine particulate matter (PM2.5) affects health adversely, yet few studies from South America have documented such associations due to the sparsity of PM2.5 measurements. Lima’s topography and aging vehicular fleet results in severe air pollution with limited amounts of monitors to effectively quantify PM2.5 levels for epidemiologic studies. We developed an advanced machine learning model to estimate daily PM2.5 concentrations at a 1 km2 spatial resolution in Lima, Peru from 2010 to 2016. We combined aerosol optical depth (AOD), meteorological fields from the European Centre for Medium-Range Weather Forecasts (ECMWF), parameters from the Weather Research and Forecasting model coupled with Chemistry (WRF-Chem), and land use variables to fit a random forest model against ground measurements from 16 monitoring stations. Overall cross-validation R2 (and root mean square prediction error, RMSE) for the random forest model was 0.70 (5.97 μg/m3). Mean PM2.5 for ground measurements was 24.7 μg/m3 while mean estimated PM2.5 was 24.9 μg/m3 in the cross-validation dataset. The mean difference between ground and predicted measurements was −0.09 μg/m3 (Std.Dev. = 5.97 μg/m3), with 94.5% of observations falling within 2 standard deviations of the difference indicating good agreement between ground measurements and predicted estimates. Surface downwards solar radiation, temperature, relative humidity, and AOD were the most important predictors, while percent urbanization, albedo, and cloud fraction were the least important predictors. Comparison of monthly mean measurements between ground and predicted PM2.5 shows good precision and accuracy from our model. Furthermore, mean annual maps of PM2.5 show consistent lower concentrations in the coast and higher concentrations in the mountains, resulting from prevailing coastal winds blown from the Pacific Ocean in the west. Our model allows for construction of long-term historical daily PM2.5 measurements at 1 km2 spatial resolution to support future epidemiological studies.

Description: Convective hydrothermal systems have been extensively studied using electrical and electromagnetic methods given the strong correlation between low conductivity anomalies associated with hydrothermal brines and high temperature areas. However, studies addressing the application of similar geophysical methods to hot dry rock geothermal systems are very limited in the literature. The Timanfaya volcanic area, located on Lanzarote Island (Canary Islands), comprises one of these hot dry rock systems, where ground temperatures ranging from 250 to 605 °C have been recorded in pyroclastic deposits at shallow (〈70 m) depths. With the aim of characterizing the geophysical signature of the high ground temperature areas, three different geophysical techniques (ground penetrating radar, electromagnetic induction and magnetic prospecting) were applied in a well-known geothermal area located inside Timanfaya National Park. The area with the highest ground temperatures was correlated with the location that exhibited strong ground penetrating radar reflections, high resistivity values and low magnetic anomalies. Moreover, the high ground temperature imaging results depicted a shallow, bowl-shaped body that narrowed and deepened vertically to a depth greater than 45 m. The ground penetrating radar survey was repeated three years later and exhibited subtle variations of the signal reflection patterns, or signatures, suggesting a certain temporal variation of the ground temperature. By identifying similar areas with the same geophysical signature, up to four additional geothermal areas were revealed. We conclude that the combined use of ground penetrating radar, electromagnetic induction and magnetic methods constitutes a valuable tool to locate and study both the geometry at depth and seasonal variability of geothermal areas associated with hot dry rock systems.

Description: The prediction of growing stock volume is one of the commonest applications of remote sensing to support the sustainable management of forest ecosystems. In this study, we used data from the 4th Spanish National Forest Inventory (SNFI-4) and from the 1st nationwide Airborne Laser Scanning (ALS) survey to develop predictive yield models for the three major commercial tree forest species (Eucalyptus globulus, Pinus pinaster and Pinus radiata) grown in north-western Spain. Integration of both types of data required prior harmonization because of differences in timing of data acquisition and difficulties in accurately geolocating the SNFI plots. The harmonised data from 477 E. globulus, 760 P. pinaster and 191 P. radiata plots were used to develop predictive models for total over bark volume, mean volume increment and total aboveground biomass by relating SNFI stand variables to metrics derived from the ALS data. The multiple linear regression methods and several machine learning techniques (k-nearest neighbour, random trees, random forest and the ensemble method) were compared. The study findings confirmed that multiple linear regression is outperformed by machine learning techniques. More specifically, the findings suggest that the random forest and the ensemble method slightly outperform the other techniques. The resulting stand level relative RMSEs for predicting total over bark volume, annual increase in total volume and total aboveground biomass ranged from 30.8–38.3%, 34.2–41.9% and 31.7–38.3% respectively. Although the predictions can be considered accurate, more precise geolocation of the SNFI plots and coincide temporarily with the ALS data would have enabled use of a much larger and robust field database to improve the overall accuracy of estimation.

Description: This paper presents a new methodological process for detecting the instantaneous land-water border at sub-pixel level from mid-resolution satellite images (30 m/pixel) that are freely available worldwide. The new method is based on using an iterative procedure to compute Laplacian roots of a polynomial surface that represents the radiometric response of a set of pixels. The method uses a first approximation of the shoreline at pixel level (initial pixels) and selects a set of neighbouring pixels to be part of the analysis window. This adaptive window collects those stencils in which the maximum radiometric variations are found by using the information given by divided differences. Therefore, the land-water surface is computed by a piecewise interpolating polynomial that models the strong radiometric changes between both interfaces. The assessment is tested on two coastal areas to analyse how their inherent differences may affect the method. A total of 17 Landsat 7 and 8 images (L7 and L8) were used to extract the shorelines and compare them against other highly accurate lines that act as references. Accurate quantitative coastal data from the satellite images is obtained with a mean horizontal error of 4.38 ± 5.66 m and 1.79 ± 2.78 m, respectively, for L7 and L8. Prior methodologies to reach the sub-pixel shoreline are analysed and the results verify the solvency of the one proposed.

Description: This work describes a new procedure aimed to semi-automatically identify clusters of active persistent scatterers and preliminarily associate them with different potential types of deformational processes over wide areas. This procedure consists of three main modules: (i) ADAfinder, aimed at the detection of Active Deformation Areas (ADA) using Persistent Scatterer Interferometry (PSI) data; (ii) LOS2HV, focused on the decomposition of Line Of Sight (LOS) displacements from ascending and descending PSI datasets into vertical and east-west components; iii) ADAclassifier, that semi-automatically categorizes each ADA into potential deformational processes using the outputs derived from (i) and (ii), as well as ancillary external information. The proposed procedure enables infrastructures management authorities to identify, classify, monitor and categorize the most critical deformations measured by PSI techniques in order to provide the capacity for implementing prevention and mitigation actions over wide areas against geological threats. Zeri, Campiglia Marittima–Suvereto and Abbadia San Salvatore (Tuscany, central Italy) are used as case studies for illustrating the developed methodology. Three PSI datasets derived from the Sentinel-1 constellation have been used, jointly with the geological map of Italy (scale 1:50,000), the updated Italian landslide and land subsidence maps (scale 1:25,000), a 25 m grid Digital Elevation Model, and a cadastral vector map (scale 1:5000). The application to these cases of the proposed workflow demonstrates its capability to quickly process wide areas in very short times and a high compatibility with Geographical Information System (GIS) environments for data visualization and representation. The derived products are of key interest for infrastructures and land management as well as decision-making at a regional scale.