Integrating LIDAR data and multispectral imagery for enhanced classification of rangeland vegetation: A meta analysis

This study compared the suitability of LIDAR (LIght Detection And Ranging) data, three-band multispectral data, and LIDAR data integrated with multispectral information, for classifying spatially complex vegetation in the Aspen Parkland of western Canada. Classifications were performed for both a) general vegetation classes limited to three major formations of deciduous forest, shrubland and grassland, and b) eight detailed vegetation classes including upland mixed prairie and fescue grasslands, closed and semi-open aspen forests, western snowberry and silverberry shrublands, and fresh and saline riparian (lowland) meadows. A Digital Elevation Model (DEM) and Surface Elevation Model (SEM) developed from LIDAR data incorporated both topographic and biological biases in community positioning across the landscape. Using multispectral data, the original digital image mosaic, its hybrid color composite, and an intensity-hue-saturation (IHS) image were each tested. Final vegetation classification was done through integration of information from both digital images and LIDAR data to evaluate the improvement in classification accuracy. Among the land cover schedules with three and eight classes of vegetation, classification from the multispectral imagery, specifically the hybrid color composite image, had the highest accuracy, peaking at 74.6% and 59.4%, respectively. In contrast, the LIDAR classification schedules led to an average classification accuracy of 64.8% and 52.3%, respectively, for the general and detailed vegetation data. Subsequent integration of the LIDAR and digital image classification schedules resulted in accuracy improvements of 16 to 20%, resulting in a superior final accuracy of 91% and 80.3%, respectively, for the three and eight classes of vegetation. A final land cover map including 8 classes of vegetation, fresh and saline water, as well as bare ground, was created for the study area with an overall accuracy of 83.9%, highlighting the benefit of integrating LIDAR and multispectral imagery for enhanced vegetation classification in heterogenous rangeland environments.     

A comparison of local variance,fractal dimension,and Moran's I as aids to multispectral image classification

The accuracy of traditional multispectral maximum‐likelihood image classification is limited by the multi‐modal statistical distributions of digital numbers from the complex, heterogenous mixture of land cover types in urban areas. This work examines the utility of local variance, fractal dimension and Moran's I index of spatial autocorrelation in segmenting multispectral satellite imagery with the goal of improving urban land cover classification accuracy. Tools available in the ERDAS Imagine^TM software package and the Image Characterization and Modeling System (ICAMS) were used to analyse Landsat ETM?+ imagery of Atlanta, Georgia. Images were created from the ETM?+ panchromatic band using the three texture indices. These texture images were added to the stack of multispectral bands and classified using a supervised, maximum likelihood technique. Although each texture band improved the classification accuracy over a multispectral only effort, the addition of fractal dimension measures is particularly effective at resolving land cover classes within urbanized areas, as compared to per‐pixel spectral classification techniques.     

The potential for integrating Sentinel 2 MSI with SPOT 5 HRG and Landsat 8 OLI imagery for monitoring semi-arid savannah woody cover

Multitemporal archived imagery enables the monitoring of savannah woody cover, for ecological purposes. Compatibility in multitemporal, multiple sensor image data would facilitate the monitoring. The decommissioning of SPOT 5 (Système Pour l’Observation de la Terre 5) left a void in multispectral imagery at the 10 m spatial resolution of its high-resolution geometric (HRG) sensor. The subsequent launch of Sentinel 2 presented an opportunity for data continuity to monitor the savannah woody cover, using equivalent 10 m resolution multispectral instrument (MSI) bands. This study examined the integration potential of Sentinel 2 MSI with the longer archive HRG and Landsat 8 (Land Satellite 8) Operational Land Imager (OLI) imagery, in assessing savannah woody cover. Images of three semi-arid savannah sites acquired on same season dates that excluded herbaceous vegetation from the spectral signature were used: November 2014 (HRG) and December 2015 (MSI, OLI). Using equivalent green (G), red (R), and near infrared (NIR) bands at 10 m (MSI, HRG) and 30 m (OLI) resolution, the woody cover was mapped through subpixel classification. The mapped woody cover was compared for statistical differences using χ² analysis at 10 m resolution (MSI, HRG) and at a degradation of the MSI and HRG images to the 30 m OLI pixel size. Conversion to top-of-atmosphere reflectance values facilitated inter-sensor correlation of G, R, and NIR reflectance for field sampling sites where woody cover was quantified. Inter-sensor regression functions in G, R, and NIR band MSI and HRG images were developed. The 10 m resolution classifications of woody cover were not statistically different. Due to spatial resolution similarity, SPOT 5 HRG multispectral imagery was established as suitable for integration with equivalent band MSI imagery in mapping the woody cover in a multitemporal analysis. For dense woody cover, Landsat 8 OLI imagery was more suitable for integration with MSI than HRG images due to higher radiometric sensitivity, which can permit monitoring physiology-related woody reflectance.     

Evaluation of LANDSAT and airborne multispectral data and aerial photographs for mapping forest features and phenomena in a part of the Godavari basin

The automatic interpretation of multispectral digital data obtained from LANDSAT as well as from an airborne multispectral scanner using an interactive computer system and visual interpretation of colour composites of LANDSAT imagery and aerial photographs of a dry deciduous forest tract were used for evaluating the discrimination capabilities of each technique and for comparative evaluation. While visual interpretation of LANDSAT imagery could give only general information, such as contiguity of vegetation cover, digital analysis of the same yielded more detailed information, such as teak-bearing and non-teak-bearing regions. The analysis of airborne multispectral data, in the present state of the art, for performing forest surveys and making maps is limited. Aerial photographs are very useful for mapping forest land features and stock, which can be done more reliably than could be done by ground surveys. Infrared photographs show better promise in mapping forest features. The integration of multitemporal data and the incorporation of digitized additional information into the data stream for the improvement of digital analysis are suggested. Acquisition of data including aerial photographs for general surveys during a period prior to leaf fall in a deciduous forest is also recommended.     

Maximum likelihood classification combined with spectral angle mapper algorithm for high resolution satellite imagery

The generation of precise land cover classification maps is an important application of high resolution satellite multispectral imagery. In this study, Spectral Angle Mapper algorithm (SAM) was used to extract the spectral characteristics from multispectral imagery. The spectral angle between neighbouring pixels was calculated. The distribution of spectral characteristics was derived from the average and variance of the calculated spectral angle in a 3×3 window of the image. The extracted spectral characteristics were combined with original multispectral imagery, and the data were classified by the maximum likelihood method. This approach was applied to Quickbird multispectral imagery. The extracted spectral characteristics highlighted boundaries between different types of land cover. The method proposed in this study exhibits an increase in overall classification accuracy relative to the original maximum likelihood method.     

Sub-pixel estimation of urban land cover components with linear mixture model analysis and Landsat Thematic Mapper imagery

We examine the utility of linear mixture modelling in the sub-pixel analysis of Landsat Enhanced Thematic Mapper (ETM) imagery to estimate the three key land cover components in an urban/suburban setting: impervious surface, managed/unmanaged lawn and tree cover. The relative effectiveness of two different endmember sets was also compared. The interior endmember set consisted of the median pixel value of the training pixels of each land cover and the exterior endmember set was the extreme pixel value. As a means of accuracy assessment, the resulting land cover estimates were compared with independent estimates obtained from the visual interpretation of digital orthophotography and classified IKONOS imagery. Impervious surface estimates from the Landsat ETM showed a high degree of similarity (RMS error (RMSE) within approximately ±10 to 15%) to that obtained using high spatial resolution digital orthophotography and IKONOS imagery. The partition of the vegetation component into tree vs grass cover was more problematic due to the greater spectral similarity between these land cover types with RMSE of approximately ±12 to 22%. The interior endmember set appeared to provide better differentiation between grass and urban tree cover than the exterior endmember set. The ability to separate the grass vs tree components in urban vegetation is of major importance to the study of the urban/suburban ecosystems as well as watershed assessment.     

Large area vegetation mapping in the Gisborne district,New Zealand,from Landsat TM

Debate over a forestry incentive scheme in the Gisborne district, New Zealand, highlighted the need for up to date information on the vegetation cover. Maps of vegetation at a scale of 1:100 000 were produced by automatically classifying Landsat Thematic Mapper (TM) imagery. The classified imagery was compared with existing vegetation information (20-years-old) from a GIS database to identify gross errors. Through field checking the discrepancies were identified as either real changes or errors in classification. Correction of errors increased the overall classification accuracy from 84 to 90 per cent.

The digital vegetation map was intersected with land use suitability data to provide a two-way table that provided land managers with quantitative information suitable for making regional planning decisions. Although the 90 per cent accuracy is high enough to permit the calculation of vegetation areas and to achieve an adequate representation of regional vegetation patterns, it is not high enough to permit the digital vegetation map to be used as a vegetation database where point queries are important.     

遥感影像超分辨率制图研究进展

尺度问题是土地覆盖分类中的一个核心问题,向下尺度转换又是其中的难点。混合像元分解可以得到亚像元尺度的类别组分百分比,但无法求得亚像元的具体位置。遥感影像超分辨率制图是由粗空间分辨率的影像得到高空间分辨率分类结果图的技术,可用于地表分类向下尺度转换,近年来该技术已成为遥感影像分类和尺度转换领域的研究热点。对超分辨率制图研究进展做了详细论述,从超分辨率制图的发展和研究现状、主要方法、精度评价等几方面进行了详细阐述,并分析了当前超分辨率制图算法存在的主要问题,以及可能的研究重点和发展空间。     

Approaches to fractional land cover and continuous field mapping: A comparative assessment over the BOREAS study region

Subpixel land cover mapping involves the estimation of surface properties using sensors whose spatial sampling is coarse enough to produce mixtures of the properties within each pixel. This study evaluates five algorithms for mapping subpixel land cover fractions and continuous fields of vegetation properties within the BOREAS study area. The algorithms include a conventional “hard”, per-pixel classifier, a neural network, a clustering/look-up-table approach, multivariate regression, and linear least squares inversion. A land cover map prepared using a Landsat TM mosaic was adopted as the source of fine scale calibration and validation data. Coarse scale mixtures of five basic land cover classes and continuous vegetation fields, both corresponding to the field of view of SPOT-VEGETATION imagery (1.15-km pixel size), were synthesised from the TM mosaic using a modelled point spread function. Two measures of land cover distribution were used, fractions of fine scale land cover categories and continuous fields of vegetation structural characteristics. The subpixel algorithms were applied using both proximate (<100 km) and distant (>400 km) separation between training and validation regions. “Hard” classification performed poorly in estimating proportions or continuous fields. The neural network, look-up-table and multivariate regression algorithms produced good matches of spatial patterns and regional land cover composition for the proximate treatment. However, all three methods exhibited substantial biases with the distant treatment due to the characteristics of the training data. Linear least squares inversion offers a relatively unbiased but less precise alternative for subpixel proportion and fraction mapping as it avoids calibration to the a priori distribution of land cover in the training data. In general, a combination of multivariate regression for proximate training data and linear least squares inversion for distant training data resulted in woody fraction estimates within 20% of the Landsat TM classification-based estimates.     

Integrating airborne multispectral imagery and airborne LiDAR data for enhanced lithological mapping in vegetated terrain

Practical and financial constraints associated with traditional field-based lithological mapping are often responsible for the generation of maps with insufficient detail and inaccurately located contacts. In arid areas with well exposed rocks and soils, high-resolution multi- and hyperspectral imagery is a valuable mapping aid as lithological units can be readily discriminated and mapped by automatically matching image pixel spectra to a set of reference spectra. However, the use of spectral imagery in all but the most barren terrain is problematic because just small amounts of vegetation cover can obscure or mask the spectra of underlying geological substrates. The use of ancillary information may help to improve lithological discrimination, especially where geobotanical relationships are absent or where distinct lithologies exhibit inherent spectral similarity. This study assesses the efficacy of airborne multispectral imagery for detailed lithological mapping in a vegetated section of the Troodos ophiolite (Cyprus), and investigates whether the mapping performance can be enhanced through the integration of LiDAR-derived topographic data. In each case, a number of algorithms involving different combinations of input variables and classification routine were employed to maximise the mapping performance. Despite the potential problems posed by vegetation cover, geobotanical associations aided the generation of a lithological map - with a satisfactory overall accuracy of 65.5% and Kappa of 0.54 - using only spectral information. Moreover, owing to the correlation between topography and lithology in the study area, the integration of LiDAR-derived topographic variables led to significant improvements of up to 22.5% in the overall mapping accuracy compared to spectral-only approaches. The improvements were found to be considerably greater for algorithms involving classification with an artificial neural network (the Kohonen Self-Organizing Map) than the parametric Maximum Likelihood Classifier. The results of this study demonstrate the enhanced capability of data integration for detailed lithological mapping in areas where spectral discrimination is complicated by the presence of vegetation or inherent spectral similarities.     

Evaluation of MK-4 multispectral satellite photography in land cover classification of eastern Ontario

Russian MK-4 multispectral satellite photography has been investigated for potential in land cover classification. Thematic maps were generated using maximum likelihood, neural network and context classifiers. Classifications of the raw spectral data, of spectral transforms, and of combined spectral/textural data were evaluated. Low point-based class accuracies resulted for land cover types exhibiting high spatial variability at the given pixel spacing of 7.5m, while more spatially homogeneous cover types were well classified. Several issues arose which need to be addressed for effective future use of high-resolution satellite sensors in regional land cover mapping. They include the need for further research in techniques for classification and accuracy assessment which are sensitive to the spatial variance of such high resolution imagery, and optimization of class attribute definitions.     

A scalable approach to mapping annual land cover at 250 m using MODIS time series data: A case study in the Dry Chaco ecoregion of South America

Land use and land cover (LULC) maps from remote sensing are vital for monitoring, understanding and predicting the effects of complex human-nature interactions that span local, regional and global scales. We present a method to map annual LULC at a regional spatial scale with source data and processing techniques that permit scaling to broader spatial and temporal scales, while maintaining a consistent classification scheme and accuracy. Using the Dry Chaco ecoregion in Argentina, Bolivia and Paraguay as a test site, we derived a suite of predictor variables from 2001 to 2007 from the MODIS 250 m vegetation index product (MOD13Q1). These variables included: annual statistics of red, near infrared, and enhanced vegetation index (EVI), phenological metrics derived from EVI time series data, and slope and elevation. For reference data, we visually interpreted percent cover of eight classes at locations with high-resolution QuickBird imagery in Google Earth. An adjustable majority cover threshold was used to assign samples to a dominant class. When compared to field data, we found this imagery to have georeferencing error < 5% the length of a MODIS pixel, while most class interpretation error was related to confusion between agriculture and herbaceous vegetation. We used the Random Forests classifier to identify the best sets of predictor variables and percent cover thresholds for discriminating our LULC classes. The best variable set included all predictor variables and a cover threshold of 80%. This optimal Random Forests was used to map LULC for each year between 2001 and 2007, followed by a per-pixel, 3-year temporal filter to remove disallowed LULC transitions. Our sequence of maps had an overall accuracy of 79.3%, producer accuracy from 51.4% (plantation) to 95.8% (woody vegetation), and user accuracy from 58.9% (herbaceous vegetation) to 100.0% (water). We attributed map class confusion to limited spectral information, sub-pixel spectral mixing, georeferencing error and human error in interpreting reference samples. We used our maps to assess woody vegetation change in the Dry Chaco from 2002 to 2006, which was characterized by rapid deforestation related to soybean and planted pasture expansion. This method can be easily applied to other regions or continents to produce spatially and temporally consistent information on annual LULC.     

A method for mapping Australian woody vegetation cover by linking continental-scale field data and long-term Landsat time series

There is a significant need to provide nationwide consistent information for land managers and scientists to assist with property planning, vegetation monitoring applications, risk assessment, and conservation activities at an appropriate spatial scale. We created maps of woody vegetation cover of Australia using a consistent method applied across the continent, and made them accessible. We classified pixels as woody or not woody, quantified their foliage projective cover, and classed them as forest or other wooded lands based on their cover density. The maps provide, for the first time, cover density estimates of Australian forests and other wooded lands with the spatial detail required for local-scale studies. The maps were created by linking field data, collected by a network of collaborators across the continent, to a time series of Landsat-5 TM and Landsat-7 ETM+ images for the period 2000–2010. The fractions of green vegetation cover, non-green vegetation cover, and bare ground were calculated for each pixel using a previously developed spectral unmixing approach. Time series statistics, for the green vegetation cover, were used to classify each pixel as either woody or not using a random forest classifier. An estimate of woody foliage projective cover was made by calibration with field measurements, and woody pixels classified as forest where the foliage cover was at least 0.1. Validation of the foliage projective cover with field measurements gave a coefficient of determination, R²,of 0.918 and root mean square error of 0.070. The user’s and producer’s accuracies for areas mapped as forest were high at 92.2% and 95.9%, respectively. The user’s and producers’s accuracies were lower for other wooded lands at 75.7% and 61.3%, respectively. Further research into methods to better separate areas with sparse woody vegetation from those without woody vegetation is needed. The maps provide information that will assist in gaining a better understanding of our natural environment. Applications range from the continental-scale activity of estimating national carbon stocks, to the local scale activities of assessing habitat suitability and property planning.     

Mapping sub-pixel proportional land cover with AVHRR imagery

A problem with NOAA AVHRR imagery is that the intrinsic scale of spatial variation in land cover in the U.K. is usually finer than the scale of sampling imposed by the image pixels. The result is that most NOAA AVHRR pixels contain a mixture of land cover types (sub-pixel mixing). Three techniques for mapping the sub-pixel proportions of land cover classes in the New Forest, U.K. were compared: (i) artificial neural networks (ANN); (ii) mixture modelling; and (iii) fuzzy c -means classification. NOAA AVHRR imagery and SPOT HRV imagery, both for 28 June 1994, were obtained. The SPOT HRV images were classified using the maximum likelihood method, and used to derive the 'known' sub-pixel proportions of each land cover class for each NOAA AVHRR pixel. These data were then used to evaluate the predictions made (using the three techniques and the NOAA AVHRR imagery) in terms of the amount of information provided, the accuracy with which that information is provided, and the ease of implementation. The ANN was the most accurate technique, but its successful implementation depended on accurate co-registration and the availability of a training data set. Supervised fuzzy c -means classification was slightly more accurate than mixture modelling.     

Land use studies in drylands: an evaluation of object‐oriented classification of very high resolution panchromatic imagery

Object‐oriented classification approaches offer an alternative to per‐pixel methods for assessment of land use and land cover. Combining object‐oriented approaches with very high resolution imagery may provide enhanced possibilities for applications requiring land use and land cover data. The aim of this study is to evaluate the application of object‐oriented classification of panchromatic very high resolution data in African drylands, where sizes and shapes of fields are varied, and intercropping practised, which might lead to difficulties in image segmentation. The results show that region‐based segmentation is sensitive to the proportion of spectral and shape information and the best results were gained when the segmentation was based on predominately spectral information. The accuracy (Kappa value of 0.6) for the object‐oriented classification was significantly higher than that for per‐pixel classification. However, both the segmentation and the classification were time‐consuming based on a trial and error process.     

Sub-pixel land cover mapping for per-field classification

A method was developed to transform a soft land cover classification into hard land cover classes at the sub-pixel scale for subsequent per-field classification. First, image pixels were segmented using vector boundaries. Second, the pixel segments (ranked by area) were labelled with a land cover class (ranked by class typicality). Third, a hard per-field classification was generated by examining each polygon (representing a land cover parcel, or field) in its entirety (by grouping the fragments of the polygon contained within different image pixels) and assigning to it the modal land cover class. The accuracy of this technique was considerably higher than that of both a corresponding hard per-pixel classification and a perfield classification based on hard per-pixel classified imagery.     

ICA mixture model algorithm for unsupervised classification of remote sensing imagery

Conventional remote sensing classification algorithms assume that the data in each class can be modelled using a multivariate Gaussian distribution. As this assumption is often not valid in practice, conventional algorithms do not perform well. In this paper, we present an independent component analysis (ICA)‐based approach for unsupervised classification of multi/hyperspectral imagery. ICA used for a mixture model estimates the data density in each class and models class distributions with non‐Gaussian (sub‐ and super‐Gaussian) probability density functions, resulting in the ICA mixture model (ICAMM) algorithm. Independent components and the mixing matrix for each class are found using an extended information‐maximization algorithm, and the class membership probabilities for each pixel are computed. The pixel is allocated to the class having maximum class membership probability to produce a classification. We apply the ICAMM algorithm for unsupervised classification of images obtained from both multispectral and hyperspectral sensors. Four feature extraction techniques are considered as a preprocessing step to reduce the dimensionality of the hyperspectral data. The results demonstrate that the ICAMM algorithm significantly outperforms the conventional K‐means algorithm for land cover classification produced from both multi‐ and hyperspectral remote sensing images.     

Sub-pixel mapping of urban land cover using multiple endmember spectral mixture analysis: Manaus, Brazil

The spatial and spectral variability of urban environments present fundamental challenges to deriving accurate remote sensing products for urban areas. Multiple endmember spectral mixture analysis (MESMA) is a technique that potentially addresses both challenges. MESMA models spectra as the linear sum of spectrally pure endmembers that vary on a per-pixel basis. Spatial variability is addressed by mapping sub-pixel components of land cover as a combination of endmembers. Spectral variability is addressed by allowing the number and type of endmembers to vary from pixel to pixel. This paper presents an application of MESMA to map the physical components of urban land cover for the city of Manaus, Brazil, using Landsat Enhanced Thematic Mapper (ETM+) imagery.We present a methodology to build a regionally specific spectral library of urban materials based on generalized categories of urban land-cover components: vegetation, impervious surfaces, soil, and water. Using this library, we applied MESMA to generate a total of 1137 two-, three-, and four-endmember models for each pixel; the model with the lowest root-mean-squared (RMS) error and lowest complexity was selected on a per-pixel basis. Almost 97% of the pixels within the image were modeled within the 2.5% RMS error constraint. The modeled fractions were used to generate continuous maps of the per-pixel abundance of each generalized land-cover component. We provide an example to demonstrate that land-cover components have the potential to characterize trajectories of physical landscape change as urban neighborhoods develop through time. Accuracy of land-cover fractions was assessed using high-resolution, geocoded images mosaicked from digital aerial videography. Modeled vegetation and impervious fractions corresponded well with the reference fractions. Modeled soil fractions did not correspond as closely with the reference fractions, in part due to limitations of the reference data. This work demonstrates the potential of moderate-resolution, multispectral imagery to map and monitor the evolution of the physical urban environment.     

Separating grassland and shrub vegetation by multidate pixel‐adaptive spectral mixture analysis

Monitoring and assessment of land degradation and the processes driving it require effective change indicators at appropriate scales and spatial extent. In this context, the decomposition of Mediterranean rangeland vegetation into woody and herbaceous fractions is of great significance. This study demonstrates that a stratification of vegetation into woody and herbaceous components is possible with two satellite images of moderate spatial and spectral resolution. We used a pixel‐adaptive spectral mixture analysis to derive subpixel‐level vegetation abundances from satellite imagery representing two specific phenological stages of Mediterranean rangeland vegetation. The transferability of endmember models is often a problem of multidate spectral mixture analysis because of uneven spectral dimensionality within and among datasets. In our approach, the dimensionality of the mixture model was determined automatically, based on error calculations. This method enables the transfer of the mixture model to multiple scenes and allows for quantitative comparison of vegetation abundances. The results show that the woody vegetation fraction corresponds well with field data (R ² = 0.76–0.91) and vegetation cover mapped from a very high resolution satellite image. The herbaceous vegetation fraction displays a good correlation compared to field mapped cover but still implies a moderate level of uncertainty (R ² = 0.52–0.76). The approach pursued in this research may be valuable for the characterization of rangeland plant communities and for the derivation of vegetation‐related indicators useful for the monitoring and assessment of degradation.     

基于压缩数据维的城市建筑用地遥感信息提取

通过压缩数据维的方式，研究城市建筑用地信息准确提取的原理和方法。通过对城市土地利用类型的分析，选取了归一化差异建筑指数、修正归一化差异水体指数和土壤调节植被指数来代表城市建成区的3种最主要地类——建筑用地、水体和植被。通过将ETM 影像原有的7个波段压缩为由它们衍生的这3个采用比值运算构成的指数波段，大大压缩了数据维数、减少了数据的相关度并降低了不同地类的光谱混淆性。因此采用简单的最大似然分类和掩膜处理技术，就可以将城市建筑用地信息提取出来，其精度可达91．2％。