Modeling Land-Cover Types Using Multiple Endmember Spectral Mixture Analysis in a Desert City

Spectral mixture analysis is probably the most commonly used approach among sub-pixel analysis techniques. This method models pixel spectra as a linear combination of spectral signatures from two or more ground components. However, spectral mixture analysis does not account for the absence of one of the surface features or spectral variation within pure materials since it utilizes an invariable set of surface features. Multiple endmember spectral mixture analysis (MESMA), which addresses these issues by allowing endmembers to vary on a per pixel basis, was employed in this study to model Landsat ETM+ reflectance in the Phoenix metropolitan area. Image endmember spectra of vegetation, soils, and impervious surfaces were collected with the use of a fine resolution Quickbird image and the pixel purity index. This study employed 204 (=3x17x4) total four-endmember models for the urban subset and 96 (=6x6x2x4) total five-endmember models for the non-urban subset to identify fractions of soil, impervious surface, vegetation, and shade. The Pearson correlation between the fraction outputs from MESMA and reference data from Quickbird 60 cm resolution data for soil, impervious, and vegetation were 0.8030, 0.8632, and 0.8496 respectively. Results from this study suggest that the MESMA approach is effective in mapping urban land covers in desert cities at sub- pixel level.

An evaluation of Multiple Endmember Spectral Mixture Analysis applied to Landsat 8 OLI images for mapping land cover in southern Africa\u27s Savanna.

Mapping land cover in southern Africa’s savannas using traditional pixel based remote sensing techniques can be very challenging due to the heterogeneity of its vegetation structure and the spectral difficulty in separating similar land covers across various land uses. In order to overcome these complications, a Multiple Endmember Spectral Mixture Analysis (MESMA) provides a potential remote sensing approach to quantify spectral variation in the physical environment at a subpixel level. The MESMA approach was applied in the study area of the Mayuni Conservancy, in Namibia. Results show that 32.3% of the study area is covered by photosynthetic vegetation (PV), 32.0% by non-photosynthetic vegetation (NPV), 25.2% by bare soil (B) and 10.6% by shade. Post-classification validation shows that MESMA presented a moderate performance in estimating the proportions of land cover types in the study area. However the validation process is limited to the available resources and carries great subjectivity. It is concluded that future research on the matter should include a more consistent investigation on the endmember selection methodology and expand the study area inside of the same ecosystem

Temporal mixture analysis for estimating impervious surface area from multi-temporal MODIS NDVI data in Japan

As a proxy measure of the human ecological footprint, impervious surface area (ISA) has recently become a key concept in the field of urban remote sensing, with a focus on estimation of the ISA at a city-scale by using Landsat-style satellite images. However, ISA estimation is also in demand in disciplines such as the environmental assessment and policy making at a national scale. This paper proposes a new method for estimating the ISA fraction in Japan based on a temporal mixture analysis (TMA) technique. The required inputs for the proposed method are rearranged MODIS NDVI time-series datasets at the temporal stable zone (i.e., the first to the sixth largest NDVI values in a year). Three ISA distribution maps obtained from Landsat-5 TM data were used as reference maps to evaluate the performance of the proposed method. The results showed that the proposed TMA-based method achieved a large reduction in the effects of endmember variability compared with the previous methods (e.g., SMA and NSMA), and thus the new method has promising accuracy for estimating ISA in Japan. The overall root mean square error (RMSE) of the proposed method was 8.7%, with a coefficient of determination of 0.86, and there was no obvious underestimation or overestimation for the whole ISA range

Assessing the role of EO in biodiversity monitoring: options for integrating in-situ observations with EO within the context of the EBONE concept

The European Biodiversity Observation Network (EBONE) is a European contribution on terrestrial monitoring to GEO BON, the Group on Earth Observations Biodiversity Observation Network. EBONE’s aims are to develop a system of biodiversity observation at regional, national and European levels by assessing existing approaches in terms of their validity and applicability starting in Europe, then expanding to regions in Africa. The objective of EBONE is to deliver: 1. A sound scientific basis for the production of statistical estimates of stock and change of key indicators; 2. The development of a system for estimating past changes and forecasting and testing policy options and management strategies for threatened ecosystems and species; 3. A proposal for a cost-effective biodiversity monitoring system. There is a consensus that Earth Observation (EO) has a role to play in monitoring biodiversity. With its capacity to observe detailed spatial patterns and variability across large areas at regular intervals, our instinct suggests that EO could deliver the type of spatial and temporal coverage that is beyond reach with in-situ efforts. Furthermore, when considering the emerging networks of in-situ observations, the prospect of enhancing the quality of the information whilst reducing cost through integration is compelling. This report gives a realistic assessment of the role of EO in biodiversity monitoring and the options for integrating in-situ observations with EO within the context of the EBONE concept (cfr. EBONE-ID1.4). The assessment is mainly based on a set of targeted pilot studies. Building on this assessment, the report then presents a series of recommendations on the best options for using EO in an effective, consistent and sustainable biodiversity monitoring scheme. The issues that we faced were many: 1. Integration can be interpreted in different ways. One possible interpretation is: the combined use of independent data sets to deliver a different but improved data set; another is: the use of one data set to complement another dataset. 2. The targeted improvement will vary with stakeholder group: some will seek for more efficiency, others for more reliable estimates (accuracy and/or precision); others for more detail in space and/or time or more of everything. 3. Integration requires a link between the datasets (EO and in-situ). The strength of the link between reflected electromagnetic radiation and the habitats and their biodiversity observed in-situ is function of many variables, for example: the spatial scale of the observations; timing of the observations; the adopted nomenclature for classification; the complexity of the landscape in terms of composition, spatial structure and the physical environment; the habitat and land cover types under consideration. 4. The type of the EO data available varies (function of e.g. budget, size and location of region, cloudiness, national and/or international investment in airborne campaigns or space technology) which determines its capability to deliver the required output. EO and in-situ could be combined in different ways, depending on the type of integration we wanted to achieve and the targeted improvement. We aimed for an improvement in accuracy (i.e. the reduction in error of our indicator estimate calculated for an environmental zone). Furthermore, EO would also provide the spatial patterns for correlated in-situ data. EBONE in its initial development, focused on three main indicators covering: (i) the extent and change of habitats of European interest in the context of a general habitat assessment; (ii) abundance and distribution of selected species (birds, butterflies and plants); and (iii) fragmentation of natural and semi-natural areas. For habitat extent, we decided that it did not matter how in-situ was integrated with EO as long as we could demonstrate that acceptable accuracies could be achieved and the precision could consistently be improved. The nomenclature used to map habitats in-situ was the General Habitat Classification. We considered the following options where the EO and in-situ play different roles: using in-situ samples to re-calibrate a habitat map independently derived from EO; improving the accuracy of in-situ sampled habitat statistics, by post-stratification with correlated EO data; and using in-situ samples to train the classification of EO data into habitat types where the EO data delivers full coverage or a larger number of samples. For some of the above cases we also considered the impact that the sampling strategy employed to deliver the samples would have on the accuracy and precision achieved. Restricted access to European wide species data prevented work on the indicator ‘abundance and distribution of species’. With respect to the indicator ‘fragmentation’, we investigated ways of delivering EO derived measures of habitat patterns that are meaningful to sampled in-situ observations

Uncertainty Assessment of Spectral Mixture Analysis in Remote Sensing Imagery

Spectral mixture analysis (SMA), a scheme of sub-pixel-based classifications, is one of the widely used models to map fractional land use and land cover information in remote sensing imagery. It assumes that: 1) a mixed pixel is composed by several pure land cover classes (endmembers) linearly or nonlinearly, and 2) the spectral signature of each endmember is a constant within the entire spatial extent of analysis. SMA has been commonly applied to impervious surface area extraction, vegetation fraction estimation, and land use and land cover change (LULC) mapping. Limitations of SMA, however, still exist. First, the existence of between- and within-class variability prevents the selection of accurate endmembers, which results in poor accuracy of fractional land cover estimates. Weighted spectral mixture analysis (WSMA) and transformed spectral mixture analysis (TSMA) are alternate means to address the within- and between- class variability. These methods, however, have not been analyzed systematically and comprehensively. The effectiveness of each WSMA and TSMA scheme is still unknown, in particular within different urban areas. Second, multiple endmember SMA (MESMA) is a better alternative to address spectral mixture model uncertainties. It, nonetheless, is time consuming and inefficient. Further, incorrect endmember selections may still limit model performance as the best-fit endmember model might not be the optimal model due to the existence of spectral variability. Therefore, this study aims 1) to explore endmember uncertainties by examining WSMA and TSMA modeling comprehensively, and 2) to develop an improved MESMA model in order to address the uncertainties of spectral mixture models. Results of the WSMA examination illustrated that some weighting schemes did reduce endmember uncertainties since they could improve the fractional estimates significantly. The results also indicated that spectral class variance played a key role in addressing the endmember uncertainties, as the better performing weighting schemes were constructed with spectral class variance. In addition, the results of TSMA examination demonstrated that some TSMAs, such as normalized spectral mixture analysis (NSMA), could effectively solve the endmember uncertainties because of their stable performance in different study areas. Results of Class-based MEMSA (C-MESMA) indicated that it could address spectral mixture model uncertainties by reducing a lot of the calculation burden and effectively improving accuracy. Assessment demonstrated that C-MEMSA significantly improving accuracy. Major contributions of this study can be summarized as follow. First, the effectiveness of addressing endmember uncertainties have been fully discussed by examining: 1) the effectiveness of ten weighted spectral mixture models in urban environments; and 2) the effectiveness of 26 transformed spectral mixture models in three locations. Constructive guidance regarding handling endmember uncertainties using WSMA and TSMA have been provided. Second, the uncertainties of spectral mixture model were reduced by developing an improved MESMA model, named C-MESMA. C-MESMA could restrict the distribution of endmembers and reduce the calculation burden of traditional MESMA, increasing SMA accuracy significantly

Urban Image Classification: Per-Pixel Classifiers, Sub-Pixel Analysis, Object-Based Image Analysis, and Geospatial Methods

Remote sensing methods used to generate base maps to analyze the urban environment rely predominantly on digital sensor data from space-borne platforms. This is due in part from new sources of high spatial resolution data covering the globe, a variety of multispectral and multitemporal sources, sophisticated statistical and geospatial methods, and compatibility with GIS data sources and methods. The goal of this chapter is to review the four groups of classification methods for digital sensor data from space-borne platforms; per-pixel, sub-pixel, object-based (spatial-based), and geospatial methods. Per-pixel methods are widely used methods that classify pixels into distinct categories based solely on the spectral and ancillary information within that pixel. They are used for simple calculations of environmental indices (e.g., NDVI) to sophisticated expert systems to assign urban land covers. Researchers recognize however, that even with the smallest pixel size the spectral information within a pixel is really a combination of multiple urban surfaces. Sub-pixel classification methods therefore aim to statistically quantify the mixture of surfaces to improve overall classification accuracy. While within pixel variations exist, there is also significant evidence that groups of nearby pixels have similar spectral information and therefore belong to the same classification category. Object-oriented methods have emerged that group pixels prior to classification based on spectral similarity and spatial proximity. Classification accuracy using object-based methods show significant success and promise for numerous urban 3 applications. Like the object-oriented methods that recognize the importance of spatial proximity, geospatial methods for urban mapping also utilize neighboring pixels in the classification process. The primary difference though is that geostatistical methods (e.g., spatial autocorrelation methods) are utilized during both the pre- and post-classification steps. Within this chapter, each of the four approaches is described in terms of scale and accuracy classifying urban land use and urban land cover; and for its range of urban applications. We demonstrate the overview of four main classification groups in Figure 1 while Table 1 details the approaches with respect to classification requirements and procedures (e.g., reflectance conversion, steps before training sample selection, training samples, spatial approaches commonly used, classifiers, primary inputs for classification, output structures, number of output layers, and accuracy assessment). The chapter concludes with a brief summary of the methods reviewed and the challenges that remain in developing new classification methods for improving the efficiency and accuracy of mapping urban areas

Sedimentary evidence of the Late Holocene tsunami in the Shetland Islands (UK) at Loch Flugarth, northern Mainland

Tsunami deposits around the North Sea basin are needed to assess the long-term hazard of tsunamis. Here, we present sedimentary evidence of the youngest tsunami on the Shetland Islands from Loch Flugarth, a coastal lake on northern Mainland. Three gravity cores show organic-rich background sedimentation with many sub-centimetre-scale sand layers, reflecting recurring storm overwash and a sediment source limited to the active beach and uppermost subtidal zone. A basal 13-cm-thick sand layer, dated to 426–787 cal. a CE based on 14C, 137Cs and Bayesian age–depth modelling, was found in all cores. High-resolution grain-size analysis identified four normally graded or massive sublayers with inversely graded traction carpets at the base of two sublayers. A thin organic-rich ‘mud’ drape and a ‘mud’ cap cover the two uppermost sublayers, which also contain small rip-up clasts. Grain-size distributions show a difference between the basal sand layer and the coarser and better sorted storm layers above. Multivariate statistical analysis of X-ray fluorescence core scanning data also distinguishes both sand units: Zr, Fe and Ti dominate the thick basal sand, while the thin storm layers are high in K and Si. Enriched Zr and Ti in the basal sand layer, in combination with increased magnetic susceptibility, may be related to higher heavy mineral content reflecting an additional marine sediment source below the storm-wave base that is activated by a tsunami. Based on reinterpretation of chronological data from two different published sites and the chronostratigraphy of the present study, the tsunami seems to date to c. 1400 cal. a BP. Although the source of the tsunami remains unclear, the lack of evidence for this event outside of the Shetland Islands suggests that it had a local source and was smaller than the older Storegga tsunami (8.15 cal. ka BP), which affected most of the North Sea basin.</p

A Comparison of the Classification of Vegetation Characteristics by Spectral Mixture Analysis and Standard Classifiers on Remotely Sensed Imagery within the Siberia Region

As an alternative to the traditional method of inferring vegetation cover characteristics from satellite data by classifying each pixel into a specific land cover type based on predefined classification schemes, the Spectral Mixture Analysis (SMA) method is applied to images of the Siberia region. A linear mixture model was applied to determine proportional estimates of land cover for, (a) agriculture and floodplain soils, (b) broadleaf, and (c) conifer classes, in pixels of 30 m resolution Landsat data. In order to evaluate the areal estimates, results were compared with ground truth data, as well as those estimates derived from more sophisticated method of image classification, providing improved estimates of endmember values and subpixel areal estimates of vegetation cover classes than the traditional approach of using predefined classification schemes with discrete numbers of cover types. This technique enables the estimation of proportional land cover type in a single pixel and could potentially serve as a tool for deriving improved estimates of vegetation parameters that are necessary for modeling carbon processes

Spatial and Temporal Dust Source Variability in Northern China Identified Using Advanced Remote Sensing Analysis

The aim of this research is to provide a detailed characterization of spatial patterns and temporal trends in the regional and local dust source areas within the desert of the Alashan Prefecture (Inner Mongolia, China). This problem was approached through multi-scale remote sensing analysis of vegetation changes. The primary requirements for this regional analysis are high spatial and spectral resolution data, accurate spectral calibration and good temporal resolution with a suitable temporal baseline. Landsat analysis and field validation along with the low spatial resolution classifications from MODIS and AVHRR are combined to provide a reliable characterization of the different potential dust-producing sources. The representation of intra-annual and inter-annual Normalized Difference Vegetation Index (NDVI) trend to assess land cover discrimination for mapping potential dust source using MODIS and AVHRR at larger scale is enhanced by Landsat Spectral Mixing Analysis (SMA). The combined methodology is to determine the extent to which Landsat can distinguish important soils types in order to better understand how soil reflectance behaves at seasonal and inter-annual timescales. As a final result mapping soil surface properties using SMA is representative of responses of different land and soil cover previously identified by NDVI trend. The results could be used in dust emission models even if they are not reflecting aggregate formation, soil stability or particle coatings showing to be critical for accurately represent dust source over different regional and local emitting areas