The utility of multi-sensor data for mapping eroded lands

The availability of remote sensing data with improved spatial, spectral and radiometric resolution is now available to fully exploit their potential for a specific application subject to the relative merits and the limitations of each sensor's data. Presented here is a case study where Landsat MSS and TM; and SPOT MLA data for part of the Bijapur district, southern India, which were acquired on the same day, have been evaluated for mapping eroded lands. The approach involves the geometric registration of all three data to a common map grid using tie points and third order polynomial transform; and resampling the MSS and TM data to a 20m by 20 m pixel dimension and radiometric normalization. Thematic maps showing eroded lands were generated on a micro-VAXbased DIPIX system using a maximum likelihood classifier. Accuracy estimates were made for the thematic maps following stratified unaligned random sampling technique, and subsequently, computing overall accuracy and Kappa coefficient. Spectral separability and classification accuracy was maximum from SPOT-MLA data followed by a combination of Landsat MSS band 1, SPOT-MLA band 2 and Landsat TM band 4; Landsat TM, a combination of Landsat MSS, TM and SPOT MLA; and Landsat MSS data.     

Cover. ERS-1 SAR interferogram of Mount Vesuvio (Italy)

A pixel block intensity modulation (PRIM) method has been developed to add spatial detail to Landsat Thematic Mapper (TM) thermal band TM6 images in regions with sufficient topography. The method uses 30 m resolution TM reflective spectral band images (TM1-5 and 7) to modulate the relevant TM6 image on the basis of its 120m resolution thermal pixel blocks. Topographic detail in each 120m resolution pixel block of the TM6 image is thus recovered, without altering the average thermal digital number level of the block, by the spatial information recorded in the reflective spectral bands at 30m resolution. Tests confirm that the PBIM can effectively integrate topographic spatial detail from reflective spectral bands with TM6 images while retaining the fidelity of the original thermal spectral information. PBIM is also applicable, as a general method, for data fusion of multispectral and panchromatic images with different spatial resolutions. Bearing in mind that, for space-borne remote sensing, the spatial resolution of the thermal band will continue to be lower than that of VISSWIR and panchromatic bands in the multispectral sensor systems of the next generation, such as Enhanced Thematic Mapper Plus, the PBIM method will remain a useful technique for enhancing thermal imagery data for some time.     

Assessing the potential of VEGETATION sensor data for mapping snow and ice cover: A Normalized Difference Snow and Ice Index

The VEGETATION (VGT) sensor in SPOT 4 has four spectral bands that are equivalent to Landsat Thematic Mapper (TM) bands (blue, red, near-infrared and mid-infrared spectral bands) and provides daily images of the global land surface at a 1-km spatial resolution. We propose a new index for identifying and mapping of snow/ice cover, namely the Normalized Difference Snow/Ice Index (NDSII), which uses reflectance values of red and mid-infrared spectral bands of Landsat TM and VGT. For Landsat TM data, NDSII is calculated as NDSII_TM=(TM3-TM5)/(TM3+TM5); for VGT data, NDSII is calculated as NDSII_VGT=(B2-MIR)/(B2+MIR). As a case study we used a Landsat TM image that covers the eastern part of the Qilian mountain range in the Qinghai-Xizang (Tibetan) plateau of China. NDSIITM gave similar estimates of the area and spatial distribution of snow/ice cover to the Normalized Difference Snow Index (NDSI=(TM2-TM5)/(TM2+TM5)) which has been proposed by Hall et al. The results indicated that the VGT sensor might have the potential for operational monitoring and mapping of snow/ice cover from regional to global scales, when using NDSII_VGT.     

The effects of spatial resolution on the classification of Thematic Mapper data

Actual and degraded LANDSAT-4 Thematic Mapper (TM) data were analysed to examine the effect of spatial resolution on the performance of a per pixel, maximum-likelihood classification algorithm. Analysis of variance (ANOVA) and a balanced, three-factor, eight-treatment, fixed-effects model were used to investigate the interactions between spatial resolution and two other TM refinements, spectral band configuration and data quantization. The goal was to evaluate quantitatively the effects of these attributes on classification accuracies obtained with all pixels (pure pixels plus mixed pixels) and on accuracies obtained with pure pixels alone.

A comparison of results from these separate analyses supported previous explanations of the effects of increasing spatial resolution. First, the difficulty in classifying mixed pixels was demonstrated by an average 21 per cent decrease in percentage accuracy from the pure-pixel case to the pure-plus-mixed-pixel case for the eight ANOVA treatments. In the pure-pixel case, an increase in spatial resolution from 80 to 30 m caused an average 6·1 per cent decrease in percentage accuracy when the other factors were held constant. This decrease was attributed to increased within-class spectral variability at the TM resolution. Finally, in the pure-plus-mixed-pixel case, increasing the spatial resolution did not significantly affect accuracy. This insignificance was attributed to a reduced proportion of mixed pixels at the TM resolution which counteracted the detrimental effects of increased spectral variability. These results point to a need for the development of new approaches to classification which take full advantage of the TM spatial resolution.     

An operational radiometric calibration procedure for the Landsat sensors based on pseudo-invariant target sites

The Statewide Landcover and Trees Study (SLATS) use both Landsat-7 ETM+ and Landsat-5 TM imagery to monitor short-term woody vegetation changes throughout Queensland, Australia. In order to analyse more subtle long-term vegetation change, time-based trends resulting from artefacts introduced by the sensor system must be removed. In this study, a reflectance-based vicarious calibration approach using high-reflectance, pseudo-invariant targets in western Queensland was developed. This calibration procedure was used to test the existing calibration models for ETM+ and TM, and develop a consistent operational calibration procedure which provides calibration information for the MSS sensors. Ground based data, sensor spectral response functions and atmospheric variables were used as input to MODTRAN radiative transfer code to estimate top-of-atmosphere radiance. The estimated gains for Landsat-7 ETM+ (1999-2003), -5 TM (1987-2004), -5 MSS (1984-1993) and -2 MSS (1979-1982) are presented. Results confirm the stability and accuracy of the ETM+ calibration, and the suitability of this data as a radiometric standard for cross-calibration with TM. Vicarious data support the use of the existing TM calibration model for the red and two shortwave-infrared bands. However, alternative models for blue, green and near-infrared bands are presented. The models proposed differ most noticeably at dates prior to 1995, with differences in estimated gains of up to 9.7%, 10.8% and 6.9% for the blue, green and near-infrared bands respectively. Vicarious gains for Landsat-2 MSS and Landsat-5 MSS are presented and are compared with those applied by the on-board calibration system. Updated calibration coefficients to scale MSS data to the SLATS vicarious measurements are given. The removal of time based calibration trends in the SLATS data archive will enable the measurement of vegetation changes over the 26 year period covered by Landsat -2, -5 and -7.     

An evaluation of Landsat TM spectral data and SAR-derived textural information for lithological discrimination in the Red Sea Hills,Sudan

The effectiveness of spectral and textural information in the identification of surface rock types in an arid region, the Red Sea Hills of Sudan, is evaluated using spectral information from the six Landsat TM optical bands and textural features derived from Shuttle Imaging Radar-C (SIR-C) C-band HH polarization data. An initial classification is derived from Landsat TM data alone using three classification algorithms, Gaussian maximum likelihood, a multi-layer feed-forward neural network and a Kohonen self-organizing feature map (SOM), to generate lithological maps, with classification accuracy being measured using a confusion matrix approach. The feed-forward neural net produced the highest overall classification accuracy of 57 per cent and was, therefore, selected for the second experiment, in which texture measures from SIR-C C-band HH-polarized synthetic aperture radar (SAR) data are added to selected TM spectral features. Four methods of measuring texture are employed, based on the Fourier power spectrum, grey level co-occurrence matrix (GLCM), multi-fractal measures, and the multiplicative autoregressive random field (MAR) model. The use of textural information together with a subset of the TM spectral features leads to an increase in classification accuracy to almost 70 per cent. Both the MAR model and the GLCM matrix approach perform better than Fourier and multi-fractal based methods of texture characterization.     

Image segmentation for humid tropical forest classification in Landsat TM data

Humid tropical forest types have low spectral separability in Landsat TM data due to highly textured reflectance patterns at the 30m spatial resolution. Two methods of reducing local spectral variation, low-pass spatial filtering and image segmentation, were examined for supervised classification of 10 forest types in TM data of Peruvian Amazonia. The number of forest classes identified at over 90% accuracy increased from one in raw imagery to three in filtered imagery, and six in segmented imagery. The ability to derive less generalised tropical forest classes may allow greater use of classified imagery in ecology and conservation planning.     

Comparison of optical sensors discrimination ability using spectral libraries

In remote sensing, the ability to discriminate different land covers or material types is directly linked with the spectral resolution and sampling provided by the optical sensor. Previous studies have shown that spectral resolution is a critical issue, especially in complex environments. In spite of the increasing availability of hyperspectral data, multispectral optical sensors onboard various satellites are acquiring every day a massive amount of data with a relatively poor spectral resolution (i.e. usually about four to seven spectral bands). These remotely sensed data are intensively used for Earth observation regardless of their limited spectral resolution. In this article, we studied seven of these optical sensors: Pleiades, QuickBird, Système Pour l'Observation de la Terre 5 (SPOT5), IKONOS, Landsat Thematic Mapper (TM), FORMOSAT, and Medium Resolution Imaging Spectrometer (MERIS). This study focuses on the ability of each sensor to discriminate different materials according to its spectral resolution. We used four different spectral libraries which contain around 2500 spectra of materials and land covers with a fine spectral resolution. These spectra were convolved with the relative spectral responses (RSRs) of each sensor to create spectra at the sensors' resolutions. Then, these reduced spectra were compared using separability indices (divergence, transformed divergence (TD), Bhattacharyya, and Jeffreys-Matusita) and machine learning tools. In the experiments, we highlighted that the configuration of spectral bands could lead to important differences in classification accuracy according to the context of application (e.g. urban area).     

Multi-scale GEOBIA with very high spatial resolution digital aerial imagery: scale,texture and image objects

This study used geographic object-based image analysis (GEOBIA) with very high spatial resolution (VHR) aerial imagery (0.3 m spatial resolution) to classify vegetation, channel and bare mud classes in a salt marsh. Three classification issues were investigated in the context of segmentation scale: (1) a comparison of single- and multi-scale GEOBIA using spectral bands, (2) the relative benefit of incorporating texture derived from the grey-level co-occurrence matrix (GLCM) in classifying the salt marsh features in single- and multi-scale GEOBIA and (3) the effect of quantization level of GLCM texture in the context of multi-scale GEOBIA. The single-scale GEOBIA experiments indicated that the optimal segmentation was both class and scale dependent. Therefore, the single-scale approach produced an only moderately accurate classification for all marsh classes. A multi-scale approach, however, facilitated the use of multiple scales that allowed the delineation of individual classes with increased between-class and reduced within-class spectral variation. With only spectral bands used, the multi-scale approach outperformed the single-scale GEOBIA with an overall accuracy of 82% vs. 76% (Kappa of 0.71 vs. 0.62). The study demonstrates the potential importance of ancillary data, GLCM texture, to compensate for limited between-class spectral discrimination. For example, gains in classification accuracies ranged from 3% to 12% when the GLCM mean texture was included in the multi-scale GEOBIA. The multi-scale classification overall accuracy varied with quantization level of the GLCM texture matrix. A quantization level of 2 reduced misclassifications of channel and bare mud and generated a statistically higher classification than higher quantization levels. Overall, the multi-scale GEOBIA produced the highest classification accuracy. The multi-scale GEOBIA is expected to be a useful methodology for creating a seamless spatial database of marsh landscape features to be used for further geographic information system (GIS) analyses.     

Application of thematic mapper data to corn and soybean development stage estimation

A model, utilizing direct relationship between remotely sensed spectral data and the development stage of both corn and soybeans has been proposed and published previously (Badhwar and Henderson, 1981; and Henderson and Badhwar, 1984). This model was developed using data acquired by instruments mounted on trucks over field plots of corn and soybeans as well as satellite data from Landsat. In all cases, the data was analyzed in the spectral bands equivalent to the four bands of Landsat multispectral scanner (MSS). In this study the same model has been applied to corn and soybeans using Landsat-4 Thematic Mapper (TM) data combined with simulated TM data to provide a multitemporal data set in TM band intervals. All data (five total acquisitions) were acquired over a test site in Webster County, Iowa from June to October 1982. The use of TM data for determining development state is as accurate as with Landsat MSS and field plot data in MSS bands. The maximum deviation of 0.6 development stage for corn and 0.8 development stage for soybeans is well within the uncertainty with which a field can be estimated with procedures used by observers on the ground in 1982.     

Land-cover classification in a moist tropical region of Brazil with Landsat Thematic Mapper imagery

This research aims to improve land-cover classification accuracy in a moist tropical region in Brazil by examining the use of different remote-sensing-derived variables and classification algorithms. Different scenarios based on Landsat Thematic Mapper (TM) spectral data and derived vegetation indices and textural images and different classification algorithms, maximum likelihood classification (MLC), artificial neural network (ANN), classification tree analysis (CTA) and object-based classification (OBC), were explored. The results indicate that a combination of vegetation indices as extra bands into Landsat TM multi-spectral bands did not improve the overall classification performance, but the combination of textural images was valuable for improving vegetation classification accuracy. In particular, the combination of both vegetation indices and textural images into TM multi-spectral bands improved the overall classification accuracy (OCA) by 5.6% and the overall kappa coefficient (OKC) by 6.25%. Comparison of the different classification algorithms indicated that CTA and ANN have poor classification performance in this research, but OBC improved primary forest and pasture classification accuracies. This research indicates that use of textural images or use of OBC are especially valuable for improving the vegetation classes such as upland and liana forest classes that have complex stand structures and large patch sizes.     

Effects of spatial and spectral resolution in estimating ecosystem α-diversity by satellite imagery

Remote sensing represents a powerful tool to derive quantitative and qualitative information about ecosystem biodiversity. In particular, since plant species richness is a fundamental indicator of biodiversity at the community and regional scales, attempts were made to predict species richness (spatial heterogeneity) by means of spectral heterogeneity. The possibility of using spectral variance of satellite images for predicting species richness is known as Spectral Variation Hypothesis. However, when using remotely sensed data, researchers are limited to specific scales of investigation. This paper aims to investigate the effects of scale (both as spatial and spectral resolution) when searching for a relation between spectral and spatial (related to plant species richness) heterogeneity, by using satellite data with different spatial and spectral resolution. Species composition was sampled within square plots of 100 m² nested in macroplots of 10,000 m². Spectral heterogeneity of each macroplot was calculated using satellite images with different spatial and spectral resolution: a Quickbird multispectral image (4 bands, spatial resolution of 3 m), an Aster multispectral image (first 9 bands used, spatial resolution of 15 m for bands 1 to 3 and 30 m for bands 4 to 9), an ortho-Landsat ETM+ multispectral image (bands 1 to 5 and band 7 used; spatial resolution, 30 m), a resampled 60 m Landsat ETM+ image.Quickbird image heterogeneity showed a statistically highly significant correlation with species richness (r = 0.69) while coarse resolution images showed contrasting results (r = 0.43, r = 0.67, and r = 0.69 considering the Aster, Landsat ETM+, and the resampled 60 m Landsat ETM+ images respectively). It should be stressed that spectral variability is scene and sensor dependent. Considering coarser spatial resolution images, in such a case even using SWIR Aster bands (i.e. the additional spectral information with respect to Quickbird image) such an image showed a very low power in catching spectral and thus spatial variability with respect to Landsat ETM+ imagery. Obviously coarser resolution data tend to have mixed pixel problems and hence less sensitive to spatial complexity. Thus, one might argue that using a finer pixel dimension should inevitably result in a higher level of detail. On the other hand, the spectral response from different land-cover features (and thus different species) in images with higher spectral resolution should exhibit higher complexity.Spectral Variation Hypothesis could be a basis for improving sampling designs and strategies for species inventory fieldwork. However, researchers must be aware on scale effects when measuring spectral (and thus spatial) heterogeneity and relating it to field data, hence considering the concept of scale not only related to a spatial framework but even to a spectral one.     

Land-cover classification in a moist tropical region of Brazil with Landsat TM imagery

This research aims to improve land-cover classification accuracy in a moist tropical region in Brazil by examining the use of different remote sensing-derived variables and classification algorithms. Different scenarios based on Landsat Thematic Mapper (TM) spectral data and derived vegetation indices and textural images, and different classification algorithms - maximum likelihood classification (MLC), artificial neural network (ANN), classification tree analysis (CTA), and object-based classification (OBC), were explored. The results indicated that a combination of vegetation indices as extra bands into Landsat TM multispectral bands did not improve the overall classification performance, but the combination of textural images was valuable for improving vegetation classification accuracy. In particular, the combination of both vegetation indices and textural images into TM multispectral bands improved overall classification accuracy by 5.6% and kappa coefficient by 6.25%. Comparison of the different classification algorithms indicated that CTA and ANN have poor classification performance in this research, but OBC improved primary forest and pasture classification accuracies. This research indicates that use of textural images or use of OBC are especially valuable for improving the vegetation classes such as upland and liana forest classes having complex stand structures and having relatively large patch sizes.     

Sensor‐independent analysis method for hyperspectral data based on the pattern decomposition method

This paper describes a modified pattern decomposition method with a supplementary pattern. The proposed approach can be regarded either as a type of spectral mixing analysis or as a kind of multivariate analysis; the later explanation is more suitable considering the presence of the additional supplementary patterns. The sensor‐independent method developed herein uses the same normalized spectral patterns for any sensor: fixed multi‐band (1260 bands) spectra serve as the universal standard spectral patterns. The resulting pattern decomposition coefficients showed sensor independence. That is, regardless of sensor, the three coefficients had nearly the same values for the same samples. The estimation errors for pattern decomposition coefficients depended on the sensor used. The estimation errors for Landsat/MSS and ALOS/AVNIR‐2 were larger than those of Landsat/TM (ETM+), Terra/MODIS and ADEOS‐II/GLI. The latter three sensors had negligibly small errors.     

面向亚像元级积雪信息提取的MODIS影像波段融合方法研究

MODIS影像因其共享性和时间序列的完整性而成为大区域积雪监测研究广泛使用的数据源,进行MODIS影像波段间融合,能够为积雪研究提供较高分辨率的影像数据源。为了充分利用MODIS影像250 m分辨率波段的空间和光谱信息,提取亚像元级的积雪面积,使用两种具有高光谱保真度的影像融合方法：基于SFIM变换和基于小波变换的融合方法,采取不同的波段组合策略,对MODIS影像bands 1~2和bands 3~7进行融合,并以Landsat TM影像的积雪分类图作为“真值”,对融合后影像进行混合像元分解得到的积雪丰度图的精度进行评价。结果表明：利用基于SFIM变换和小波变换方法融合后影像提取的积雪分类图精度较高,数量精度为75%,比未融合影像积雪分类图的精度提高了6%,表明MODIS影像波段融合是一种提取高精度积雪信息的有效方法。     

Mapping of land degradation from space: a comparative study of Landsat ETM+ and ASTER data

The purpose of this study is to compare the role of spectral and spatial resolutions in mapping land degradation from space‐borne imagery using Landsat ETM+ and ASTER data as examples. Land degradation in the form of salinization and waterlogging in Tongyu County, western Jilin Province of northeast China was mapped from an ETM+ image of 22 June 2002 and an ASTER image recorded on 24 June 2001 using supervised classification, together with several other land covers. It was found that the mapping accuracy was achieved at 56.8% and higher for moderately degraded (e.g. salinized) farmland, and over 80% for severely degraded land (e.g. barren) from both ASTER and ETM+ data. The spatial resolution of the ASTER data exerts only a negligible effect on the mapping accuracy. The 30 m ETM+ outperforms the ASTER image of both 15 m and 30 m resolution in consistently generating a higher overall accuracy as well as a higher user's accuracy for barren land. The inferiority of ASTER data is attributed to the highly repetitive spectral content of its six shortwave infrared bands. It is concluded that the spectral resolution of an image is not as important as the information content of individual bands in accurately mapping land covers automatically.     

Spatial discrimination of salt- and sodium-affected soil surfaces

Salinization-alkalinization is a time- and space-dynamic soil degradation process in semiarid regions. This study implements a synergistic approach to map salt- and sodium affected surfaces, combining digital image classification with field observation of soil degradation features and laboratory determinations. Salinity-alkalinity classes were established using the electrical conductivity (EC) and pH values. A neighbourhood operator, with spatial and spectral user-defined constraints determined the spectral objects constituting the training set. Six combined Landsat TM bands (1,2,4,5,6,7) provided the highest separability between salt- and sodium-affected soil classes. Although the overall accuracy was slightly low (64 per cent), accuracies of 100 per cent were obtained for some classes. Main causes of spectral confusions, masking different salinity-alkalinity degrees were the type and abundance of salt-tolerant vegetation cover, the topsoil textures, and the mixture of topsoil properties under field conditions.     

Using MERIS fused images for land-cover mapping and vegetation status assessment in heterogeneous landscapes

In this paper we evaluate the potential of ENVISAT–Medium Resolution Imaging Spectrometer (MERIS) fused images for land-cover mapping and vegetation status assessment in heterogeneous landscapes. A series of MERIS fused images (15 spectral bands; 25 m pixel size) is created using the linear mixing model and a Landsat Thematic Mapper (TM) image acquired over the Netherlands. First, the fused images are classified to produce a map of the eight main land-cover types of the Netherlands. Subsequently, the maps are validated using the Dutch land-cover/land-use database as a reference. Then, the fused image with the highest overall classification accuracy is selected as the best fused image. Finally, the best fused image is used to compute three vegetation indices: the normalized difference vegetation index (NDVI) and two indices specifically designed to monitor vegetation status using MERIS data: the MERIS terrestrial chlorophyll index (MTCI) and the MERIS global vegetation index (MGVI).

Results indicate that the selected data fusion approach is able to downscale MERIS data to a Landsat-like spatial resolution. The spectral information in the fused images originates fully from MERIS and is not influenced by the TM data. Classification results for the TM and for the best fused image are similar and, when comparing spectrally similar images (i.e. TM with no short-wave infrared bands), the results of the fused image outperform those of TM. With respect to the vegetation indices, a good correlation was found between the NDVI computed from TM and from the best fused image (in spite of the spectral differences between these two sensors). In addition, results show the potential of using MERIS vegetation indices computed from fused images to monitor individual fields. This is not possible using the original MERIS full resolution image. Therefore, we conclude that MERIS–TM fused images are very useful to map heterogeneous landscapes.     

A wavelet transform method to merge Landsat TM and SPOT panchromatic data

To take advantage of the high spectral resolution of Landsat TM images and the high spatial resolution of SPOT panchromatic images (SPOT PAN), we present a wavelet transform method to merge the two data types. In a pyramidal fashion, each TM reflective band or SPOT PAN image was decomposed into an orthogonal wavelet representation at a given coarser resolution, which consisted of a low frequency approximation image and a set of high frequency, spatially-oriented detail images. Band-by-band, the merged images were derived by performing an inverse wavelet transform using the approximation image from each TM band and detail images from SPOT PAN. The spectral and spatial features of the merged results of the wavelet methods were compared quantitatively with those of intensity-hue-saturation (IHS), principal component analysis (PCA), and the Brovey transform. It was found that multisensor data merging is a trade-off between the spectral information from a low spatial-high spectral resolution sensor and the spatial structure from a high spatial-low spectral resolution sensor. With the wavelet merging method, it is easy to control this trade-off. Experiments showed that the simultaneous best spectral and spatial quality can only be achieved with wavelet transform methods, compared with the three other approaches examined.     

Using Landsat TM data to estimate carbon release from burned biomass in an Alaskan spruce forest complex

Fire disturbance in boreal forests can release carbon to the atmosphere stored in both the aboveground vegetation and the organic soil layer. Estimating pyrogenic emissions of carbon released during biomass burning in these forests is useful for understanding and estimating global carbon budgets. In this work, we have developed a method to estimate carbon efflux for the burned black spruce in an Alaskan forest by combining information derived from Landsat Thematic Mapper (TM) data and field measurements. We have used the spatial and spectral information of TM data to identify and measure two important factors: pre-burn black spruce density and burn severity. Field measurements provided estimates of aboveground and ground layer carbon per unit area for the pre-burn Landsat spectral classes, and percentage of carbon consumed for the post-burn Landsat spectral classes. Carbon release estimates for the burned black spruce were computed using field data and the co-occurrence of the pre-burn and post-burn spectral classes. The estimated carbon released was 39.9tha-1, which is 57% greater than an estimate computed using AVHRR data and estimates of pre-burn biomass and carbon fractions consumed that were not site specific or spatially varying. We conclude that the spectral bands and spatial resolution of Landsat TM data provide the potential for improved estimates of pyrogenic carbon efflux relative to the coarser spectral and spatial resolution of other multispectral sensors.