A Quantitative Assessment of Forest Cover Change in the Moulouya River Watershed (Morocco) by the Integration of a Subpixel-Based and Object-Based Analysis of Landsat Data

A quantitative assessment of forest cover change in the Moulouya River watershed (Morocco) was carried out by means of an innovative approach from atmospherically corrected reflectance Landsat images corresponding to 1984 (Landsat 5 Thematic Mapper) and 2013 (Landsat 8 Operational Land Imager). An object-based image analysis (OBIA) was undertaken to classify segmented objects as forested or non-forested within the 2013 Landsat orthomosaic. A Random Forest classifier was applied to a set of training data based on a features vector composed of different types of object features such as vegetation indices, mean spectral values and pixel-based fractional cover derived from probabilistic spectral mixture analysis). The very high spatial resolution image data of Google Earth 2013 were employed to train/validate the Random Forest classifier, ranking the NDVI vegetation index and the corresponding pixel-based percentages of photosynthetic vegetation and bare soil as the most statistically significant object features to extract forested and non-forested areas. Regarding classification accuracy, an overall accuracy of 92.34% was achieved. The previously developed classification scheme was applied to the 1984 Landsat data to extract the forest cover change between 1984 and 2013, showing a slight net increase of 5.3% (ca. 8800 ha) in forested areas for the whole region

Use of waveform lidar and hyperspectral sensors to assess selected spatial and structural patterns associated with recent and repeat disturbance and the abundance of sugar maple (Acer saccharum Marsh.) in a temperate mixed hardwood and conifer forest.

Abstract Waveform lidar imagery was acquired on September 26, 1999 over the Bartlett Experimental Forest (BEF) in New Hampshire (USA) using NASA\u27s Laser Vegetation Imaging Sensor (LVIS). This flight occurred 20 months after an ice storm damaged millions of hectares of forestland in northeastern North America. Lidar measurements of the amplitude and intensity of ground energy returns appeared to readily detect areas of moderate to severe ice storm damage associated with the worst damage. Southern through eastern aspects on side slopes were particularly susceptible to higher levels of damage, in large part overlapping tracts of forest that had suffered the highest levels of wind damage from the 1938 hurricane and containing the highest levels of sugar maple basal area and biomass. The levels of sugar maple abundance were determined through analysis of the 1997 Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) high resolution spectral imagery and inventory of USFS Northern Research Station field plots. We found a relationship between field measurements of stem volume losses and the LVIS metric of mean canopy height (r2 = 0.66; root mean square errors = 5.7 m3/ha, p \u3c 0.0001) in areas that had been subjected to moderate-to-severe ice storm damage, accurately documenting the short-term outcome of a single disturbance event

Nonlinear unmixing of hyperspectral images: Models and algorithms

When considering the problem of unmixing hyperspectral images, most of the literature in the geoscience and image processing areas relies on the widely used linear mixing model (LMM). However, the LMM may be not valid, and other nonlinear models need to be considered, for instance, when there are multiscattering effects or intimate interactions. Consequently, over the last few years, several significant contributions have been proposed to overcome the limitations inherent in the LMM. In this article, we present an overview of recent advances in nonlinear unmixing modeling

Robust Linear Spectral Unmixing using Anomaly Detection

This paper presents a Bayesian algorithm for linear spectral unmixing of hyperspectral images that accounts for anomalies present in the data. The model proposed assumes that the pixel reflectances are linear mixtures of unknown endmembers, corrupted by an additional nonlinear term modelling anomalies and additive Gaussian noise. A Markov random field is used for anomaly detection based on the spatial and spectral structures of the anomalies. This allows outliers to be identified in particular regions and wavelengths of the data cube. A Bayesian algorithm is proposed to estimate the parameters involved in the model yielding a joint linear unmixing and anomaly detection algorithm. Simulations conducted with synthetic and real hyperspectral images demonstrate the accuracy of the proposed unmixing and outlier detection strategy for the analysis of hyperspectral images

USING HYPERSPECTRAL IMAGERY TO ASSIST FEDERAL FOREST MONITORING AND RESTORATION PROJECTS IN THE SOUTHERN ROCKY MOUNTAINS, COLORADO

Hyperspectral imagery and the corresponding ability to conduct analysis below the pixel level have tremendous potential to aid in landcover monitoring. During large ecosystem restoration projects, being able to monitor specific aspects of the recovery over large and often inaccessible areas under constrained finances are major challenges. The Civil Air Patrol's Airborne Real-time Cueing Hyperspectral Enhanced Reconnaissance (ARCHER) can provide hyperspectral data in most parts of the United States at relatively low cost. Although designed specifically for use in locating downed aircraft, the imagery holds the potential to identify specific aspects of landcover at far greater fidelity than traditional multispectral means. The goals of this research were to improve the use of ARCHER hyperspectral imagery to classify sub-canopy and open-area vegetation in coniferous forests located in the Southern Rockies and to determine how much fidelity might be lost from a baseline of 1 meter spatial resolution resampled to 2 and 5 meter pixel size to simulate higher altitude collection. Based on analysis comparing linear spectral unmixing with a traditional supervised classification, the linear spectral unmixing proved to be statistically superior. More importantly, however, linear spectral unmixing provided additional sub-pixel information that was unavailable using other techniques. The second goal of determining fidelity loss based on spatial resolution was more difficult to determine due to how the data are represented. Furthermore,the 2 and 5 meter imagery were obtained by resampling the 1 meter imagery and therefore may not be representative of the quality of actual 2 or 5 meter imagery. Ultimately, the information derived from this research may be useful in better utilizing hyperspectral imagery to conduct forest monitoring and assessment

Collaborative sparse regression using spatially correlated supports - Application to hyperspectral unmixing

This paper presents a new Bayesian collaborative sparse regression method for linear unmixing of hyperspectral images. Our contribution is twofold; first, we propose a new Bayesian model for structured sparse regression in which the supports of the sparse abundance vectors are a priori spatially correlated across pixels (i.e., materials are spatially organised rather than randomly distributed at a pixel level). This prior information is encoded in the model through a truncated multivariate Ising Markov random field, which also takes into consideration the facts that pixels cannot be empty (i.e, there is at least one material present in each pixel), and that different materials may exhibit different degrees of spatial regularity. Secondly, we propose an advanced Markov chain Monte Carlo algorithm to estimate the posterior probabilities that materials are present or absent in each pixel, and, conditionally to the maximum marginal a posteriori configuration of the support, compute the MMSE estimates of the abundance vectors. A remarkable property of this algorithm is that it self-adjusts the values of the parameters of the Markov random field, thus relieving practitioners from setting regularisation parameters by cross-validation. The performance of the proposed methodology is finally demonstrated through a series of experiments with synthetic and real data and comparisons with other algorithms from the literature

Gradients in urban material composition: A new concept to map cities with spaceborne imaging spectroscopy data

To understand processes in urban environments, such as urban energy fluxes or surface temperature patterns, it is important to map urban surface materials. Airborne imaging spectroscopy data have been successfully used to identify urban surface materials mainly based on unmixing algorithms. Upcoming spaceborne Imaging Spectrometers (IS), such as the Environmental Mapping and Analysis Program (EnMAP), will reduce the time and cost-critical limitations of airborne systems for Earth Observation (EO). However, the spatial resolution of all operated and planned IS in space will not be higher than 20 to 30 m and, thus, the detection of pure Endmember (EM) candidates in urban areas, a requirement for spectral unmixing, is very limited. Gradient analysis could be an alternative method for retrieving urban surface material compositions in pixels from spaceborne IS. The gradient concept is well known in ecology to identify plant species assemblages formed by similar environmental conditions but has never been tested for urban materials. However, urban areas also contain neighbourhoods with similar physical, compositional and structural characteristics. Based on this assumption, this study investigated (1) whether cover fractions of surface materials change gradually in urban areas and (2) whether these gradients can be adequately mapped and interpreted using imaging spectroscopy data (e.g. EnMAP) with 30 m spatial resolution. Similarities of material compositions were analysed on the basis of 153 systematically distributed samples on a detailed surface material map using Detrended Correspondence Analysis (DCA). Determined gradient scores for the first two gradients were regressed against the corresponding mean reflectance of simulated EnMAP spectra using Partial Least Square regression models. Results show strong correlations with R2 = 0.85 and R2 = 0.71 and an RMSE of 0.24 and 0.21 for the first and second axis, respectively. The subsequent mapping of the first gradient reveals patterns that correspond to the transition from predominantly vegetation classes to the dominance of artificial materials. Patterns resulting from the second gradient are associated with surface material compositions that are related to finer structural differences in urban structures. The composite gradient map shows patterns of common surface material compositions that can be related to urban land use classes such as Urban Structure Types (UST). By linking the knowledge of typical material compositions with urban structures, gradient analysis seems to be a powerful tool to map characteristic material compositions in 30 m imaging spectroscopy data of urban areas

A Method for Comparing Real and Simulated Airborne Hyperspectral Imagery

This thesis proposes a method for comparing real and simulated hyperspectral imagery by examining the characteristics of simulated imagery in comparison to real imagery acquired with multiple sensors hosted on an airborne platform. The dataset includes aerial multi- and hyperspectral imagery with spatial resolutions of one meter or less. The multispectral imagery includes data from an airborne sensor with three-band visible color and calibrated radiance imagery in the long-, mid-, and short-wave infrared. The airborne hyperspectral imagery includes 360 bands of calibrated radiance and reflectance data spanning 400 to 2450 nm in wavelength. Collected in September 2012, the imagery is of a park in Avon, NY, and includes a dirt track and areas of grass, gravel, forest, and agricultural fields. A number of artificial targets were deployed in the scene prior to collection for purposes of target detection, subpixel detection, spectral unmixing, and 3D object recognition. A synthetic reconstruction of the collection site was created in DIRSIG, an image generation and modeling tool developed by the Rochester Institute of Technology, based on ground-measured reflectance data, ground photography, and previous airborne imagery. Simulated airborne images were generated using the scene model, time of observation, atmospheric conditions, and knowledge of the sensor characteristics. The thesis provides a comparison between the empirical and simulated images, including a comparison of achieved performance for classification, detection and unmixing applications. It was found that several differences exist due to the way the image is generated, including finite sampling and incomplete knowledge of atmospheric conditions and sensor characteristics. The lessons learned from these differences can be used to refine the modeling tool and its use as part of ongoing efforts to increase the realism of the simulated data

Remote Sensing Methods and Applications for Detecting Change in Forest Ecosystems

Forest ecosystems are being altered by climate change, invasive species, and additional stressors. Our ability to detect these changes and quantify their impacts relies on detailed data across spatial and temporal scales. This dissertation expands the ecological utility of long-term satellite imagery by developing high quality forest mapping products and examining spatiotemporal changes in tree species abundance and phenology across the northeastern United States (US; the ‘Northeast’). Species/genus-level forest composition maps were developed by integrating field data and Landsat images to model abundance at a sub-pixel scale. These abundance maps were then used to 1) produce a more detailed, accurate forest classification compared to similar products and 2) construct a 30-year time-series of abundance for eight common species/genera. Analyzing the time-series data revealed significant abundance trends in notable species, including increases in American beech (Fagus grandifolia) at the expense of sugar maple (Acer saccharum). Climate was the dominant predictor of abundance trends, indicating climate change may be altering competitive relationships. Spatiotemporal trends in deciduous forest phenology – start and end of the growing season (SOS/EOS) – were examined based on MODIS imagery from 2001-2015. SOS exhibited a slight advancing trend across the Northeast, but with a distinct spatial pattern: eastern ecoregions showed advance and western ecoregions delay. EOS trended substantially later almost everywhere. SOS trends were linked to winter-spring temperature and precipitation trends; areas with higher elevation and fall precipitation anomalies had negative associations with EOS trends. Together, this work demonstrates the value of remote sensing in furthering our understanding of long-term forest responses to changing environmental conditions. By highlighting potential changes in forest composition and function, the research presented here can be used to develop forest conservation and management strategies in the Northeast