Spectral mixture analysis to assess post-fire vegetation regeneration using Landsat Thematic Mapper imagery: Accounting for soil brightness variation

Post-fire vegetation cover is a crucial parameter in rangeland management. This study aims to assess the post-fire vegetation recovery 3 years after the large 2007 Peloponnese (Greece) wildfires. Post-fire recovery landscapes typically are mixed vegetation-substrate environments which makes spectral mixture analysis (SMA) a very effective tool to derive fractional vegetation cover maps. Using a combination of field and simulation techniques this study aimed to account for the impact of background brightness variability on SMA model performance. The field data consisted out of a spectral library of in situ measured reflectance signals of vegetation and substrate and 78 line transect plots. In addition, a Landsat Thematic Mapper (TM) scene was employed in the study. A simple SMA, in which each constituting terrain feature is represented by its mean spectral signature, a multiple endmember SMA (MESMA) and a segmented SMA, which accounts for soil brightness variations by forcing the substrate endmember choice based on ancillary data (lithological map), were applied. In the study area two main spectrally different lithological units were present: relatively bright limestone and relatively dark flysch (sand-siltstone). Although the simple SMA model resulted in reasonable regression fits for the flysch and limestones subsets separately (coefficient of determination R² of respectively 0.67 and 0.72 between field and TM data), the performance of the regression model on the pooled dataset was considerably weaker (R² = 0.65). Moreover, the regression lines significantly diverged among the different subsets leading to systematic over-or underestimations of the vegetative fraction depending on the substrate type. MESMA did not solve the endmember variability issue. The MESMA model did not manage to select the proper substrate spectrum on a reliable basis due to the lack of shape differences between the flysch and limestone spectra,. The segmented SMA model which accounts for soil brightness variations minimized the variability problems. Compared to the simple SMA and MESMA models, the segmented SMA resulted in a higher overall correlation (R² = 0.70), its regression slope and intercept were more similar among the different substrate types and its resulting regression lines more closely resembled the expected one-one line. This paper demonstrates the improvement of a segmented approach in accounting for soil brightness variations in estimating vegetative cover using SMA. However, further research is required to evaluate the model's performance for other soil types, with other image data and at different post-fire timings.     

Prior-knowledge-based spectral mixture analysis for impervious surface mapping

In this study, we developed a prior-knowledge-based spectral mixture analysis (PKSMA) to map impervious surfaces by using endmembers derived separately for high- and low-density urban regions. First, an urban area was categorized into high- and low-density urban areas, using a multi-step classification method. Next, in high-density urban areas that were assumed to have only vegetation and impervious surfaces (ISs), the vegetation–impervious model (V–I) was used in a spectral mixture analysis (SMA) with three endmembers: vegetation, high albedo, and low albedo. In low-density urban areas, the vegetation–impervious–soil model (V–I–S) was used in an SMA analysis with four endmembers: high albedo, low albedo, soil, and vegetation. The fraction of IS with high and low albedo in each pixel was combined to produce the final IS map. The root mean-square error (RMSE) of the IS map produced using PKSMA was about 11.0%, compared to 14.52% only using four-endmember SMA. Particularly in high-density urban areas, PKSMA (RMSE = 6.47%) showed better performance than four-endmember (15.91%). The results indicate that PKSMA can improve IS mapping compared to traditional SMA by using appropriately selected endmembers and is particularly strong in high-density urban areas.     

Estimating fractional land cover in semi-arid central Kalahari: the impact of mapping method (spectral unmixing vs. object-based image analysis) and vegetation morphology

Focusing on the central Kalahari, this study utilized fractional cover of photosynthetic vegetation (f_PV), non-photosynthetic vegetation (f_NPV) and bare soil (f_BS), derived in situ and estimated from GeoEye-1 imagery using Multiple Endmember Spectral Mixture Analysis (MESMA) and object-based image analysis (OBIA) to determine superior method for fractional cover estimation and the impact of vegetation morphology on the estimation accuracy. MESMA mapped fractional cover by testing endmember models of varying complexity. Based on OBIA, image was segmented at five segmentation scales followed by classification. MESMA provided more accurate fractional cover estimates than OBIA. The increasing segmentation scale in OBIA resulted in a consistent increase in error. Different vegetation morphology types showed varied responses to the changing segmentation scale, reflecting their unique ecology and physiognomy. While areas under woody cover produced lower error even at coarse segmentation scales, those with herbaceous cover provided low error only at the fine segmentation scale.     

Lithological mapping from hyperspectral data by improved use of spectral angle mapper

The spectral angle mapper (SAM), as a spectral matching method, has been widely used in lithological type identification and mapping using hyperspectral data. The SAM quantifies the spectral similarity between an image pixel spectrum and a reference spectrum with known components. In most existing studies a mean reflectance spectrum has been used as the reference spectrum for a specific lithological class. However, this conventional use of SAM does not take into account the spectral variability, which is an inherent property of many rocks and is further magnified in remote sensing data acquisition process. In this study, two methods of determining reference spectra used in SAM are proposed for the improved lithological mapping. In first method the mean of spectral derivatives was combined with the mean of original spectra, i.e., the mean spectrum and the mean spectral derivative were jointly used in SAM classification, to improve the class separability. The second method is the use of multiple reference spectra in SAM to accommodate the spectral variability. The proposed methods were evaluated in lithological mapping using EO-1 Hyperion hyperspectral data of two arid areas. The spectral variability and separability of the rock types under investigation were also examined and compared using spectral data alone and using both spectral data and first derivatives. The experimental results indicated that spectral variability significantly affected the identification of lithological classes with the conventional SAM method using a mean reference spectrum. The proposed methods achieved significant improvement in the accuracy of lithological mapping, outperforming the conventional use of SAM with a mean spectrum as the reference spectrum, and the matching filtering, a widely used spectral mapping method.     

Stratified spectral mixture analysis of medium resolution imagery for impervious surface mapping

Linear spectral mixture analysis (LSMA) is widely employed in impervious surface estimation, especially for estimating impervious surface abundance in medium spatial resolution images. However, it suffers from a difficulty in endmember selection due to within-class spectral variability and the variation in the number and the type of endmember classes contained from pixel to pixel, which may lead to over or under estimation of impervious surface. Stratification is considered as a promising process to address the problem. This paper presents a stratified spectral mixture analysis in spectral domain (Sp_SSMA) for impervious surface mapping. It categorizes the entire data into three groups based on the Combinational Build-up Index (CBI), the intensity component in the color space and the Normalized Difference Vegetation Index (NDVI) values. A suitable endmember model is developed for each group to accommodate the spectral variation from group to group. The unmixing into the associated subset (or full set) of endmembers in each group can make the unmixing adaptive to the types of endmember classes that each pixel actually contains. Results indicate that the Sp_SSMA method achieves a better performance than full-set-endmember SMA and prior-knowledge-based spectral mixture analysis (PKSMA) in terms of R, RMSE and SE.     

Evaluation of MIVIS Hyperspectral Data for Mapping Soil Degradation in an Upland Ecosystem of Southern Italy

Abstract

This paper describes the first stage of an experiment aiming to evaluate the potential and limitations of MIVIS data for mapping the degradational state of soils in a sub‐scene of a southern Apennines study area (Italy). After radiometric rectification of the image data and the collection of a field/laboratory spectral library, linear spectral mixture modelling (SMA) was used to decompose image spectra into fractions of spectrally distinct mixing components. Spectral endmember selection was based upon a principal component analysis (PCA) applied to a set of soil spectra, collected from the spectral library. The resulting abundance estimates (fractions) trough SMA were then analysed to identify soil conditions and to obtain an improved measure of dry and green vegetation cover. A map of soil conditions and dry‐green vegetation abundance, based upon MIVIS data was then derived from normalised fractions of soil‐vegetation endmembers obtained from SMA.     

Enhancing endmember selection in multiple endmember spectral mixture analysis (MESMA) for urban impervious surface area mapping using spectral angle and spectral distance parameters

Successful retrieval of urban impervious surface area is achieved with remote sensing data using the multiple endmember spectral mixture analysis (MESMA). MESMA is well suited for studying the urban impervious surface area because it allows the number and types of the endmembers to vary on a per-pixel basis, thereby, allowing the control of the large spectral variability. However, MESMA must calculate all potential endmember combinations of each pixel to determine the best-fit one. Therefore, it is a time-consuming and inefficient unmixing technology, especially for hyperspectral images because these images have more complicated endmember categories. Hence, in this paper, we design an improved MESMA (SASD-MESMA: spectral angle and spectral distance MESMA) to enhance the computational efficiency of conventional MESMA, and we validate this new method by analyzing the Hyperion image (Jan-2011) and the field-spectra data of Guangzhou (China). In SASD-MESMA, the parameters of spectral angle (SA) and spectral distance (SD) are used to evaluate the similarity degree between library spectra and image spectra in order to identify the most representative endmember combination for each pixel. Results demonstrate that the SA and SD parameters are useful to reduce misjudgment in selecting candidate endmembers and effective for determining the appropriate endmembers in one pixel. Meanwhile, this research indicates that the proposed SASD-MESMA performs very well in retrieving impervious surface area, forest, grass and soil distributions on the sub-pixel level (the overall root mean square error (RMSE) is 0.15 and the correlation coefficient of determination (R²) is 0.68).     

Hyperspectral analysis of soil polluted with four types of hydrocarbons

In this study, a high spectral resolution GER-2600 spectroradiometer was used to obtain the spectral data of soil samples that were polluted with four different types of petroleum–hydrocarbons products: Diesel, Gasoline, Crude Oil and Fuel Oil. The polluted soil samples were prepared in the laboratory at five concentrations levels: unpolluted soil, 2500, 100,000, 250,000 ppm and pure pollutant. Spectral data were pre-processed and then analysed with various approaches: Principal Components Transformation and ANOVA, Spectral Angle Mapper (SAM), Hydrocarbon Index (HI) and Spectral Mixture Analysis (SMA). The results showed that it was possible to determine the different spectral response between clean soil and some of the polluted soils: crude oil at concentrations higher than 100,000 ppm were the easiest to recognize; while samples polluted with gasoline at concentrations below 250,000 ppm were the most difficult to distinguish from non-polluted samples.     

流域尺度的不透水面遥感提取

一个地区的不透水面覆盖度不仅是该地区城镇化程度重要指示因子,也是该地区生态环境状况的重要指示因子.现有的不透水面遥感提取方法,多集中在城区尺度上.而流域尺度上快速、准确的不透水面遥感提取方法在国内外还鲜有研究.本研究以覆盖海河流域同一季节的Landsat影像为数据源,利用已有土地利用数据集中的道路、城市、农村和工业用地...     

混合像元分解技术及其进展

混合像元分解模型是定量遥感研究的重要组成部分,为各种地学应用提供了更精细的亚像元级地物信息,这一领域受到国内外学者们广泛关注。本文围绕混合像元分解研究的4个核心问题——光谱混合模型、端元提取、模型反演方法以及解混精度评估,总结了近20年来混合像元分解的重要研究进展,分析和介绍了典型算法模型的原理和思路。进一步阐述了现有研究在一些关键问题上存在的不足,如目前仍缺乏公认的线性和非线性模型的选择判据、已有的混合像元分解模型无法抑制由端元光谱相关造成的共线性问题。最后总结了混合像元分解未来的发展趋势和值得探索的研究方向。如结合辐射传输模型和地面试验,定量分析多次散射的影响机制,以及结合克服共线性的统计回归模型。     

Spectral mixture analysis to monitor defoliation in mixed-aged Eucalyptus globulus Labill plantations in southern Australia using Landsat5-TM and EO-1Hyperion data

Defoliation is a key parameter of forest health and is associated with reduced productivity and tree mortality. Assessing the health of forests requires regular observations over large areas. Satellite remote sensing provides a cost-effective alternative to traditional ground-based assessment of forest health, but assessing defoliation can be difficult due to mixed pixels where vegetation cover is low or fragmented. In this study we apply a novel spectral unmixing technique, referred to as weighted Multiple Endmember Spectral Mixture Analysis (wMESMA), to Landsat 5-TM and EO-1 Hyperion data acquired over a Eucalyptus globulus (Labill.) plantation in southern Australia. This technique combines an iterative mixture analysis cycle allowing endmembers to vary on a per pixel basis (MESMA) and a weighting algorithm that prioritizes wavebands based on their robustness against endmember variability. Spectral mixture analysis provides an estimate of the physically interpretable canopy cover, which is not necessarily correlated with defoliation in mixed-aged plantations due to natural variation in canopy cover as stands age. There is considerable variability in the degree of defoliation as well as in stand age among sites and in this study we found that results were significantly improved by the inclusion of an age correction algorithm for both the multi-spectral (R²_{no age correction} = 0.55 vs R²_{age correction} = 0.73 for Landsat) and hyperspectral (R²_{no age correction} = 0.12 vs R²_{age correction} = 0.50 for Hyperion) image data. The improved accuracy obtained from Landsat compared to the Hyperion data illustrates the potential of applying SMA techniques for analysis of multi-spectral datasets such as MODIS and SPOT-VEGETATION.     

Soil Organic Carbon mapping of partially vegetated agricultural fields with imaging spectroscopy

Soil Organic Carbon (SOC) is one of the key soil properties, but the large spatial variation makes continuous mapping a complex task. Imaging spectroscopy has proven to be an useful technique for mapping of soil properties, but the applicability decreases rapidly when fields are partially covered with vegetation. In this paper we show that with only a few percent fractional maize cover the accuracy of a Partial Least Square Regression (PLSR) based SOC prediction model drops dramatically. However, this problem can be solved with the use of spectral unmixing techniques. First, the fractional maize cover is determined with linear spectral unmixing, taking the illumination and observation angles into account. In a next step the influence of maize is filtered out from the spectral signal by a new procedure termed Residual Spectral Unmixing (RSU). The residual soil spectra resulting from this procedure are used for mapping of SOC using PLSR, which could be done with accuracies comparable to studies performed on bare soil surfaces (Root Mean Standard Error of Calibration = 1.34 g/kg and Root Mean Standard Error of Prediction = 1.65 g/kg). With the presented RSU approach it is possible to filter out the influence of maize from the mixed spectra, and the residual soil spectra contain enough information for mapping of the SOC distribution within agricultural fields. This can improve the applicability of airborne imaging spectroscopy for soil studies in temperate climates, since the use of the RSU approach can extend the flight-window which is often constrained by the presence of vegetation.     

矿化蚀变区典型地物光谱特征分析与空间结构研究——以内蒙古突泉县—扎鲁特旗成矿带为例

基于野外实测光谱数据,分析岩石、土壤和植被的光谱反射特征,提出在光谱特征空间研究典型地物的几何结构与空间关系,分析总结中等植被覆盖的矿物蚀变区典型地物在光谱特征空间中的分布形态与空间分布关系,为遥感矿化蚀变信息提取提供科学依据。     

Integration of full-waveform LiDAR and hyperspectral data to enhance tea and areca classification

Hyperspectral image and full-waveform light detection and ranging (LiDAR) data provide useful spectral and geometric information for classifying land cover. Hyperspectral images contain a large number of bands, thus providing land-cover discrimination. Waveform LiDAR systems record the entire time-varying intensity of a return signal and supply detailed information on geometric distribution of land cover. This study developed an efficient multi-sensor data fusion approach that integrates hyperspectral data and full-waveform LiDAR information on the basis of minimum noise fraction and principal component analysis. Then, support vector machine was used to classify land cover in mountainous areas. Results showed that using multi-sensor fused data achieved better accuracy than using a hyperspectral image alone, with overall accuracy increasing from 83% to 91% using population error matrices, for the test site. The classification accuracies of forest and tea farms exhibited significant improvement when fused data were used. For example, classification results were more complete and compact in tea farms based on fused data. Fused data considered spectral and geometric land-cover information, and increased the discriminability of vegetation classes that provided similar spectral signatures.     

Use of field reflectance data for crop mapping using airborne hyperspectral image

Recent developments in hyperspectral remote sensing technologies enable acquisition of image with high spectral resolution, which is typical to the laboratory or in situ reflectance measurements. There has been an increasing interest in the utilization of in situ reference reflectance spectra for rapid and repeated mapping of various surface features. Here we examined the prospect of classifying airborne hyperspectral image using field reflectance spectra as the training data for crop mapping. Canopy level field reflectance measurements of some important agricultural crops, i.e. alfalfa, winter barley, winter rape, winter rye, and winter wheat collected during four consecutive growing seasons are used for the classification of a HyMAP image acquired for a separate location by (1) mixture tuned matched filtering (MTMF), (2) spectral feature fitting (SFF), and (3) spectral angle mapper (SAM) methods. In order to answer a general research question “what is the prospect of using independent reference reflectance spectra for image classification”, while focussing on the crop classification, the results indicate distinct aspects. On the one hand, field reflectance spectra of winter rape and alfalfa demonstrate excellent crop discrimination and spectral matching with the image across the growing seasons. On the other hand, significant spectral confusion detected among the winter barley, winter rye, and winter wheat rule out the possibility of existence of a meaningful spectral matching between field reflectance spectra and image. While supporting the current notion of “non-existence of characteristic reflectance spectral signatures for vegetation”, results indicate that there exist some crops whose spectral signatures are similar to characteristic spectral signatures with possibility of using them in image classification.     

Comparison of Different Mapping Techniques for Classifying Hyperspectral Data

Hyperion is a space borne sensor which provides powerful tool in discriminating land cover features including urban area and in preparation of urban maps. It gives hyperspectral images in 242 bands within 400?nm to 2,500?nm wavelength range with 10?nm band-width. The Hyperion image in raw form is badly affected with several atmospheric effects which cause haziness. In this study hyperspectral image is atmospherically corrected by using FLAASH model of ENVI. After atmospheric correction the urban area was mapped using the spectral endmember collected by the procedure which includes minimum noise fraction (MNF), pixel purity index (PPI) and n-dimensional visualization in ENVI software. The aim of this study is to map the urban area using several mapping techniques such as Spectral Angle Mapper (SAM), Mixture Tune Matched Filtering (MTMF) and Linear Spectral Unmixing. The urban land covers displayed noticeable differences from one another in the spectral responses in the Hyperion image. The overall accuracy of the SAM classified map was 89.41%, which indicated good potential of Hyperion image for Classification. Use of the other approaches, linear spectral unmixing and MTMF have improved the classification results.     

The use of single-date MODIS imagery for estimating large-scale urban impervious surface fraction with spectral mixture analysis and machine learning techniques

Urban impervious surface information is essential for urban and environmental applications at the regional/national scales. As a popular image processing technique, spectral mixture analysis (SMA) has rarely been applied to coarse-resolution imagery due to the difficulty of deriving endmember spectra using traditional endmember selection methods, particularly within heterogeneous urban environments. To address this problem, we derived endmember signatures through a least squares solution (LSS) technique with known abundances of sample pixels, and integrated these endmember signatures into SMA for mapping large-scale impervious surface fraction. In addition, with the same sample set, we carried out objective comparative analyses among SMA (i.e. fully constrained and unconstrained SMA) and machine learning (i.e. Cubist regression tree and Random Forests) techniques. Analysis of results suggests three major conclusions. First, with the extrapolated endmember spectra from stratified random training samples, the SMA approaches performed relatively well, as indicated by small MAE values. Second, Random Forests yields more reliable results than Cubist regression tree, and its accuracy is improved with increased sample sizes. Finally, comparative analyses suggest a tentative guide for selecting an optimal approach for large-scale fractional imperviousness estimation: unconstrained SMA might be a favorable option with a small number of samples, while Random Forests might be preferred if a large number of samples are available.     

Combining QuickBird,LiDAR, and GIS topography indices to identify a single native tree species in a complex landscape using an object-based classification approach

There are now a wide range of techniques that can be combined for image analysis. These include the use of object-based classifications rather than pixel-based classifiers, the use of LiDAR to determine vegetation height and vertical structure, as well terrain variables such as topographic wetness index and slope that can be calculated using GIS. This research investigates the benefits of combining these techniques to identify individual tree species. A QuickBird image and low point density LiDAR data for a coastal region in New Zealand was used to examine the possibility of mapping Pohutukawa trees which are regarded as an iconic tree in New Zealand. The study area included a mix of buildings and vegetation types. After image and LiDAR preparation, single tree objects were identified using a range of techniques including: a threshold of above ground height to eliminate ground based objects; Normalised Difference Vegetation Index and elevation difference between the first and last return of LiDAR data to distinguish vegetation from buildings; geometric information to separate clusters of trees from single trees, and treetop identification and region growing techniques to separate tree clusters into single tree crowns. Important feature variables were identified using Random Forest, and the Support Vector Machine provided the classification. The combined techniques using LiDAR and spectral data produced an overall accuracy of 85.4% (Kappa 80.6%). Classification using just the spectral data produced an overall accuracy of 75.8% (Kappa 67.8%). The research findings demonstrate how the combining of LiDAR and spectral data improves classification for Pohutukawa trees.     

机载激光雷达及高光谱的森林乔木物种多样性遥感监测

利用机载LiDAR和高光谱数据并结合37个地面调查样本数据,基于结构差异与光谱变异理论,通过相关分析法分别筛选了3个最优林冠结构参数和6个最优光谱指数,在单木尺度上利用自适应C均值模糊聚类算法,在神农架国家自然保护区开展森林乔木物种多样性监测,实现了森林乔木物种多样性的区域成图。研究结果表明,(1)基于结合形态学冠层控制的分水岭算法可以获得较高精度的单木分割结果(R~2=0.88,RMSE=13.17,P0.001);(2)基于LiDAR数据提取的9个结构参数中,95%百分位高度、冠层盖度和植被穿透率为最优结构参数,与Shannon-Wiener指数的相关性达到R~2=0.39—0.42(P0.01);(3)基于机载高光谱数据筛选的16个常用的植被指数中,CRI、OSAVI、Narrow band NDVI、SR、Vogelmann index1、PRI与Shannon-Wiener指数的相关性最高(R~2=0.37—0.45,P0.01);(4)在研究区,利用以30 m×30 m为窗口的自适应模糊C均值聚类算法可预测的最大森林乔木物种数为20,物种丰富度的预测精度为R~2=0.69,RMSE=3.11,Shannon-Wiener指数的预测精度为R~2=0.70,RMSE=0.32。该研究在亚热带森林开展乔木物种多样性监测,是在区域尺度上进行物种多样性成图的重要实践,可有效补充森林生物多样性本底数据的调查手段,有助于实现生物多样性的长期动态监测及科学分析森林物种多样性的现状和变化趋势。     

Incorporating band selection in the spatial selection of spectral endmembers

The impact of band selection on endmember selection is seldom explored in the analysis of hyperspectral imagery. This study incorporates the N-dimensional Spectral Solid Angle (NSSA) band selection tool into the Spectral-Spatial Endmember Extraction (SSEE) tool to determine a band set that can be used to better define endmembers classes used in spectral mixture analysis. The incorporation aims to define a band set that improves the spectral contrast between endmembers at each step of the spatial-spectral endmember search and ultimately captures key features for discriminating spectrally similar materials. The proposed method (NSSA-SSEE) was evaluated for lithological mapping using a hyperspectral image encompassing a range of spectrally similar mafic and ultramafic rock units. The band selected by NSSA-SSEE showed a good agreement with known features of scene components identified by experts. Results showed an improvement in the selection of detailed endmembers, endmembers that are similar and that can be significant for mapping. The incorporation of NSSA into SSEE was feasible because both methods are well suited for this process. NSSA is one of the few methods of band selection that is suitable for the analysis of a small number of endmembers and SSEE provides such endmember sets via spatial subsetting. The automated NSSA-SSEE approach can reduce the need for field-based information to guide the feature selection process.