Open Access
Identifying forest harvesting practices: clear-cutting and thinning in diverse tree species using dense Landsat time series
(Elsevier, 2024-12-07) 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
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.
Open Access
Evaluation of R tools for downloading MODIS images and their use in urban growth analysis of the city of Tarija (Bolivia)
(MDPI, 2022) Campero Taboada, Milton J.; Luquin Oroz, Eduardo Adrián; Montesino San Martín, Manuel; González de Audícana Amenábar, María; Campo-Bescós, Miguel; Estadística, Informática y Matemáticas; Estatistika, Informatika eta Matematika; Ingeniería; Ingeniaritza; Institute on Innovation and Sustainable Development in Food Chain - ISFOOD; Universidad Pública de Navarra / Nafarroako Unibertsitate Publikoa
The aim of this study was to compare the available tools in R for downloading and processing Moderate Resolution Imaging Spectroradiometer (MODIS) data, specifically the Enhanced Vegetation Index (EVI) product. The R tools evaluated were the MODIS package, RGISTools, MODISTools, R Google Earth Engine (RGEE) package, MODIStsp, and the Application for Extracting and Exploring Analysis Ready Samples (AppEEARS) application. Each tool was used to download the same product (EVI) corresponding to the same day (3 December 2015), and downloaded data were used to analyze the urban growth of Tarija (Bolivia) as an interesting application. The following features were analyzed: download time and memory used during the download, additional postprocessing time, local memory occupied on the computer, and downloaded file formats. Results showed that the most efficient R tools were those that work directly in the “cloud” or use text queries (RGEE and AppEEARS, respectively) and provide, as a final product, a cropped.tif image according to the area of interest.
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
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
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
Validation of a simplified model to generate multispectral synthetic images
(MDPI, 2015) 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
A new procedure to assess the quality of topographic correction (TOC) algorithms applied to remote sensing imagery was previously proposed by the authors. This procedure was based on a model that simulated synthetic scenes, representing the radiance an optical sensor would receive from an area under some specific conditions. TOC algorithms were then applied to synthetic scenes and the resulting corrected scenes were compared with a horizontal synthetic scene free of topographic effect. This comparison enabled an objective and quantitative evaluation of TOC algorithms. This approach showed promising results but had some shortcomings that are addressed herein. First, the model, originally built to simulate only broadband panchromatic scenes, is extended to multispectral scenes in the visible, near infrared (NIR), and short wave infrared (SWIR) bands. Next, the model is validated by comparing synthetic scenes with four Satellite pour l'Observation de la Terre 5 (SPOT5) real scenes acquired on different dates and different test areas along the Pyrenees mountain range (Spain). The results obtained show a successful simulation of all the spectral bands. Therefore, the model is deemed accurate enough for its purpose of evaluating TOC algorithms.
Open Access
Synthetic images for evaluating topographic correction algorithm
(IEEE, 2013) Sola Torralba, Ion; González de Audícana Amenábar, María; Álvarez-Mozos, Jesús; Torres Escribano, José Luis; Proyectos e Ingeniería Rural; Landa Ingeniaritza eta Proiektuak; Gobierno de Navarra / Nafarroako Gobernua
In the last years, many topographic correction (TOC) methods have been proposed to correct the illumination differences between the areas observed by optical remote sensors. Although the available number of TOC methods is high, the evaluation of their performance generally relies on the existence of precise land-cover information, and a standardized and objective evaluation procedure has not been proposed yet. In this paper, we propose an objective procedure to assess the accuracy of these TOC methods on the basis of simulated scenes, i.e., synthetically generated images. These images represent the radiance an optical sensor would receive under specific geometric and temporal acquisition conditions and assuming a certain land-cover type. A simplified method for creating synthetic images using the stateof- the-art irradiance models is proposed, both considering the real topography of a certain area [synthetic real (SR) image] or considering the relief of this area as being completely flat [synthetic horizontal image (SH)]. The comparison between the corrected image obtained by applying a TOC method to the SR and SH images of the same area, allows assessing the performance of each TOC algorithm. This comparison is quantitatively carried out using the structural similarity index. The proposed TOC evaluation procedure is applied to a specific case study in northern Spain to explain its implementation and demonstrate its potential. The procedure proposed in this paper could be also used to assess the behavior of TOC methods operating under different scenarios considering diverse topographic, geometrical, and temporal acquisition configurations.
Open Access
Multitemporal evaluation of topographic correction algorithms using synthetic images
(SPIE, 2012) 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
Land cover classification and quantitative analysis of multispectral data in mountainous regions is considerably hampered by the influence of topography on the spectral response pattern. In the last years, different topographic correction (TOC) algorithms have been proposed to correct illumination differences between sunny and shaded areas observed by optical remote sensors. Although the available number of TOC methods is high, the evaluation of their performance usually relies on the existence of precise land cover information, and a standardised and objective evaluation procedure has not been proposed yet. Besides, previous TOC assessment studies only considered a limited set of illumination conditions, normally assuming favourable illumination conditions. This paper presents a multitemporal evaluation of TOC methods based on synthetically generated images in order to evaluate the influence of solar angles on the performance of TOC methods. These synthetic images represent the radiance an optical sensor would receive under specific geometric and temporal acquisition conditions and assuming a certain land-cover type. A method for creating synthetic images using state-of-the-art irradiance models has been tested for different periods of the year, which entails a variety of solar angles. Considering the real topography of a specific area a Synthetic Real image (SR) is obtained, and considering the relief of this area as being completely flat a Synthetic Horizontal image (SH) is obtained. The comparison between corrected image obtained applying a TOC method to SR image and SH image of the same area, i.e. considered the ideal correction, allows assessing the performance of each TOC algorithm.
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.