Open Access
Multi-criteria evaluation of topographic correction methods
(Elsevier, 2016) Sola Torralba, Ion; González de Audícana Amenábar, María; Álvarez-Mozos, Jesús; Proyectos e Ingeniería Rural; Landa Ingeniaritza eta Proiektuak
In the last decades, several topographic correction methods (TOC) have been proposed, but there is not an agreement on the best method. Furthermore, different evaluation criteria have been used in the past, and there is not any simple and objective evaluation procedure to measure the quality of the correction. Consequently, a multicriteria analysis of widely used topographic correction methods is required that evaluates their performance over different sensors, terrain and temporal configurations. In this work, ten TOC methods were assessed using seven different evaluation strategies. The analysis was carried out for three SPOT5 images acquired over a mountainous area of northern Spain. The images had different acquisition dates and solar angles, so as to evaluate performance under varying illumination conditions. The results obtained showed that Statistic-Empiricalmethod, CCorrection and Sun-Canopy-Sensor+C performed the best, and differences were minor when favorable illumination conditions were considered. For the seven tested evaluation strategies, interquartile range reduction of land covers or the comparison of sunlit and shaded slopes gave very similar results, whereas there were greater contrasts among other criteria.
Open Access
Automatic detection of uprooted orchards based on orthophoto texture analysis
(MDPI, 2017) Ciriza Labiano, Raquel; Sola Torralba, Ion; Albizua, Lourdes; Álvarez-Mozos, Jesús; González de Audícana Amenábar, María; Proyectos e Ingeniería Rural; Landa Ingeniaritza eta Proiektuak
Permanent crops, such as olive groves, vineyards and fruit trees, are important in European agriculture because of their spatial and economic relevance. Agricultural geographical databases (AGDBs) are commonly used by public bodies to gain knowledge of the extension covered by these crops and to manage related agricultural subsidies and inspections. However, the updating of these databases is mostly based on photointerpretation, and thus keeping this information up-to-date is very costly in terms of time and money. This paper describes a methodology for automatic detection of uprooted orchards (parcels where fruit trees have been eliminated) based on the textural classification of orthophotos with a spatial resolution of 0.25 m. The textural features used for this classification were derived from the grey level co-occurrence matrix (GLCM) and wavelet transform, and were selected through principal components (PCA) and separability analyses. Next, a Discriminant Analysis classification algorithm was used to detect uprooted orchards. Entropy, contrast and correlation were found to be the most informative textural features obtained from the co-occurrence matrix. The minimum and standard deviation in plane 3 were the selected features based on wavelet transform. The classification based on these features achieved a true positive rate (TPR) of over 80% and an accuracy (A) of over 88%. As a result, this methodology enabled reducing the number of fields to photointerpret by 60–85%, depending on the membership threshold value selected. The proposed approach could be easily adopted by different stakeholders and could increase significantly the efficiency of agricultural database updating tasks.
Open Access
Monitoring rainfed alfalfa growth in semiarid agrosystems using Sentinel-2 imagery
(MDPI, 2021) Echeverría Obanos, Andrés; Urmeneta, Alejandro; González de Audícana Amenábar, María; González de Andrés, Ester; Zientziak; Ingeniaritza; Institute for Multidisciplinary Research in Applied Biology - IMAB; Institute on Innovation and Sustainable Development in Food Chain - ISFOOD; Ciencias; Ingeniería; Gobierno de Navarra / Nafarroako Gobernua
The aim of this study was to assess the utility of Sentinel-2 images in the monitoring of the fractional vegetation cover (FVC) of rainfed alfalfa in semiarid areas such as that of Bardenas Reales in Spain. FVC was sampled in situ using 1 m2 surfaces at 172 points inside 18 alfalfa fields from late spring to early summer in 2017 and 2018. Different vegetation indices derived from a series of Sentinel-2 images were calculated and were then correlated with the FVC measurements at the pixel and parcel levels using different types of equations. The results indicate that the normalized difference vegetation index (NDVI) and FVC were highly correlated at the parcel level (R 2 = 0.712), where as the correlation at the pixel level remained moderate across each of the years studied. Based on the findings, another 29 alfalfa plots (28 rainfed; 1 irrigated) were remotely monitored operationally for 3 years (2017–2019), revealing that location and weather conditions were strong determinants of alfalfa growth in Bardenas Reales. The results of this study indicate that Sentinel-2 imagery is a suitable tool for monitoring rainfed alfalfa pastures in semiarid areas, thus increasing the potential success of pasture management.
Open Access
Inter-comparison of atmospheric correction methods on Sentinel-2 images applied to croplands
(IEEE, 2018) Sola Torralba, Ion; Álvarez-Mozos, Jesús; González de Audícana Amenábar, María; Ingeniería; Ingeniaritza
Atmospheric correction of high resolution satellite scenery is a necessary preprocessing step for applications where bottom of atmosphere (BOA) reflectances are needed. The selection of the best atmospheric correction method to use on images acquired from new platforms, such as Sentinel-2, is essential to provide accurate BOA reflectances. In this work the performance of three atmospheric correction methods (6S, MAJA and SEN2COR) applied to Sentinel-2 scenes are compared by evaluating the resultant spectral signatures of six crop types on two specific dates, and their NDVI time series along a complete year. Although SEN2COR introduced greater corrections, especially in the infrared bands, the results suggest a varying performance of the methods depending on the land cover and the atmospheric conditions. Further research, particularly incorporating ground truth data, is recommended to rigorously validate the different atmospheric methods.
Open Access
Assessment of atmospheric correction methods for Sentinel-2 images in Mediterranean landscapes
(Elsevier, 2018) Sola Torralba, Ion; García-Martín, Alberto; Sandonís Pozo, Leire; Álvarez-Mozos, Jesús; González de Audícana Amenábar, María; Ingeniería; Ingeniaritza
Atmospheric correction of optical satellite imagery is an essential pre-processing for modelling biophysical variables, multi-temporal analysis, and digital classification processes. Sentinel-2 products available for users are distributed by the European Space Agency (ESA) as Top Of Atmosphere reflectance values in cartographic geometry (Level-1C product). In order to obtain Bottom Of Atmosphere reflectance images (Level-2A product) derived from this Level-1C products, ESA provides the SEN2COR module, which is implemented in the Sentinel Application Platform. Alternatively, ESA recently distributes Level-2A products processed by SEN2COR with a default configuration. On the other hand, the conversion from Level-1C to Level-2A product can be generated using alternative atmospheric correction methods, such as MAJA, 6S, or iCOR. In this context, this paper aims to evaluate the quality of Level-2A products obtained through different methods in Mediterranean shrub and grasslands by comparing data obtained from Sentinel-2 imagery with field spectrometry data. For that purpose, six plots with different land covers (asphalt, grass, shrub, pasture, and bare soil) were analyzed, by using synchronous imagery to fieldwork (from July to September 2016). The results suggest the suitability of the applied atmospheric corrections, with coefficients of determination higher than 0.90 and root mean square error lower than 0.04 achieving a relative error in bottom of atmosphere reflectance of only 2–3%. Nevertheless, minor differences were observed between the four tested methods, with slightly varying results depending on the spectral band and land cover.
Open Access
Identifying forest harvesting practices: clear-cutting and thinning in diverse tree species using dense Landsat time series
(Elsevier, 2025-02-15) Giambelluca, Ana Laura; Hermosilla, Txomin; Álvarez-Mozos, Jesús; González de Audícana Amenábar, María; Ingeniería; Ingeniaritza; Institute on Innovation and Sustainable Development in Food Chain - ISFOOD; Gobierno de Navarra / Nafarroako Gobernua
Forest monitoring plays a critical role in achieving sustainable forest management practices. The ability to identify ongoing harvesting activities is crucial for developing targeted strategies to maintain forest health. Traditional monitoring methods, which rely on field inventories, are often expensive and time-consuming. Remote sensing offers an interesting alternative, leveraging dense time series of satellite imagery and various algorithms for disturbance detection. This study presents and assesses a novel methodology for identifying forest harvesting practices (clear-cutting and thinning) using Continuous Change Detection and Classification (CCDC) algorithm, available in Google Earth Engine. The methodology comprises two steps. In the first step, performed at the pixel level, the CCDC algorithm was used to detect changes in the vegetation cover by considering Landsat 8 spectral bands, vegetation indices, and different combinations thereof. In the second step, two optimal thresholds were determined to identify forest harvesting practices based on the proportion of pixels flagged as change. This study was conducted in forest stands consisting of different conifer and broadleaf species. Accuracy was assessed using an independent set of photo-interpreted samples. The results indicated that the short-wave infrared 2 was the best individual band for forest harvesting practices identification, with an average F-score of 0.77 ± 0.06, overperforming vegetation indices. The combination of all spectral bands was the most effective to identify both clear-cuts and thinning (F-score = 0.85 ± 0.05). This combination was used to evaluate the accuracy of this approach for identifying harvesting practices over different tree species. Poplar (Populus sp.) had the highest identification rate (F-score = 0.99 ± 0.02), while black pine (Pinus nigra J.F. Arnold) stands had the lowest F-score (0.74 ± 0.05). These results highlight the ability to accurately identify forest harvesting practices even in heterogeneous forests with a high diversity of tree species using dense time series of Landsat imagery.
Open Access
Analysis of fire services coverage in Spain
(DYNA, 2018) Echeverría Iriarte, Francisco Javier; González de Audícana Amenábar, María; López Maestresalas, Ainara; Arazuri Garín, Silvia; Ciriza Labiano, Raquel; Jarén Ceballos, Carmen; Ingeniería; Ingeniaritza; Universidad Pública de Navarra / Nafarroako Unibertsitate Publikoa
Previous analysis of the locations of fire stations in Spain and the extent of the areas they cover revealed significant deficiencies with regard to the proportion of communities who would not receive fire service intervention within a reasonable time period. This article discusses and describes the use of Geographic Information Systems and related tools to determine the areas and population covered by existing fire services within a specific response time. This response time by road, is based on a survey of fire service interventions in other European countries. The analysis compares data from a statistical study with georeferenced ones and demonstrates that the areas and communities not covered within this response time is greater than previously believed. The article then describes an analysis an alternative solution to reinforce the current fire stations network with part-time firefighters to cover the areas not covered mainly in rural and remote locations.
Open Access
The added value of stratified topographic correction of multispectral images
(MDPI, 2016) Sola Torralba, Ion; González de Audícana Amenábar, María; Álvarez-Mozos, Jesús; Proyectos e Ingeniería Rural; Landa Ingeniaritza eta Proiektuak; Universidad Pública de Navarra / Nafarroako Unibertsitate Publikoa
Satellite images in mountainous areas are strongly affected by topography. Different studies demonstrated that the results of semi-empirical topographic correction algorithms improved when a stratification of land covers was carried out first. However, differences in the stratification strategies proposed and also in the evaluation of the results obtained make it unclear how to implement them. The objective of this study was to compare different stratification strategies with a non-stratified approach using several evaluation criteria. For that purpose, Statistic-Empirical and Sun-Canopy-Sensor + C algorithms were applied and six different stratification approaches, based on vegetation indices and land cover maps, were implemented and compared with the non-stratified traditional option. Overall, this study demonstrates that for this particular case study the six stratification approaches can give results similar to applying a traditional topographic correction with no previous stratification. Therefore, the non-stratified correction approach could potentially aid in removing the topographic effect, because it does not require any ancillary information and it is easier to implement in automatic image processing chains. The findings also suggest that the Statistic-Empirical method performs slightly better than the Sun-Canopy-Sensor + C correction, regardless of the stratification approach. In any case, further research is necessary to evaluate other stratification strategies and confirm these results.
Open Access
Evaluación multitemporal de métodos de corrección topográfica mediante el uso de imágenes sintéticas multiespectrales
(Asociación Española de Teledetección, 2014) Sola Torralba, Ion; Álvarez-Mozos, Jesús; González de Audícana Amenábar, María; Torres Escribano, José Luis; Proyectos e Ingeniería Rural; Landa Ingeniaritza eta Proiektuak
En este trabajo se presentan los resultados de la evaluación multitemporal de varios métodos de corrección topográfica (TOC), cuya bondad se determina de forma cuantitativa mediante el uso de imágenes sintéticas multiespectrales simuladas para diferentes fechas de adquisición a lo largo del año. Para cada fecha se generan dos imágenes sintéticas, una considerando el relieve real (imagen SR), y otra el relieve horizontal (imagen SH). Las imágenes SR se corrigen utilizando distintos TOC y estas imágenes corregidas se comparan con la corrección ideal (imagen SH) mediante el índice de similitud estructural (SSIM). Los valores de SSIM nos permiten evaluar la eficacia de cada corrección para distintas fechas, es decir, para distintos ángulos de elevación solar.
Open Access
Evaluation of 2D models for the prediction of surface depression storage using realistic reference values
(Wiley, 2016) Giménez Díaz, Rafael; Mezkiritz Barberena, Irantzu; Campo-Bescós, Miguel; Álvarez-Mozos, Jesús; González de Audícana Amenábar, María; Martínez de Aguirre Escobar, Alejandro; Casalí Sarasíbar, Javier; Proyectos e Ingeniería Rural; Landa Ingeniaritza eta Proiektuak
Depression storage (DS) is the maximum storage of precipitation and runoff in the soil surface at a given slope. The DS is determined by soil roughness that in agricultural soils is largely affected by tillage. The direct measurement of DS is not straightforward because of the natural permeability of the soil. Therefore, DS has generally been estimated from 2D/3D empirical relationships and numerical algorithms based on roughness indexes and height measurements of the soil surface, respectively. The objective of this work was to evaluate the performance of some 2D models for DS, using direct and reliable measurements of DS in an agricultural soil as reference values. The study was carried out in experimental microplots where DS was measured in six situations resulting from the combination of three types of tillage carried out parallel and perpendicular to the main slope. Those data were used as reference to evaluate four empirical models and a numerical method. Longitudinal altitudinal profiles of the relief were obtained by a laser profilometer. Infiltration measurements were carried out before and after tillage. The DS was largely affected by tillage and its direction. Highest values of DS are found on rougher surfaces mainly when macroforms cut off the dominant slope. The empirical models had a limited performance while the numerical method was the most effective, even so, with an important variability. In addition, a correct hydrological management should take into account that each type of soil tillage affects infiltration rate differently.