Estimating impervious surfaces from medium spatial resolution imagery: a comparison between fuzzy classification and LSMA

Impervious surfaces are important environmental indicators and are related to many environmental issues, such as water quality, stream health and the urban heat island effect. Therefore, detailed impervious surface information is crucial for urban planning and environment management. To extract impervious surfaces from remote sensing imagery, many algorithms and techniques have been developed. However, there are still debates over the strengths and limitations of linear versus nonlinear algorithms in handling mixed pixels in the urban landscapes. In the meantime, although many previous studies have compared various techniques, few comparisons were made between linear and nonlinear techniques. The objective of this study is to compare the performance between nonlinear and linear methods for impervious surface extraction from medium spatial resolution imagery. A linear spectral mixture analysis (LSMA) and a fuzzy classifier were applied to three Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) images acquired on 5 April 2004, 16 June 2001 and 3 October 2000, which covered Marion County, Indiana, United States. An aerial photo of Marion County with a spatial resolution of 0.14 m was used for validation of estimation results. Six impervious surface maps were yielded, and an accuracy assessment was performed. The root mean square error (RMSE), the mean average error (MAE), and the coefficient of determination (R ²) were calculated to indicate the accuracy of impervious surface maps. The results show that the fuzzy classification outperformed LSMA in impervious surface estimation in all seasons. For the June image, LSMA yielded a result with an RMSE of 13.2%, while the fuzzy classifier yielded an RMSE of 12.4%. For the April image, LSMA yielded an accuracy of 21.1% and the fuzzy classifier yielded 17.0%. For the October image, LSMA yielded a result with an RMSE of 19.8%, but the fuzzy classifier yielded an RMSE of 17.5%. Moreover, a subset image of the commercial, high-density and low-density residential areas was selected in order to compare the effectiveness of the developed algorithms for estimating impervious surfaces of different land use types. The result shows that the fuzzy classification was more effective than LSMA in both high-density and low-density residential areas. These areas prevailed with mixed pixels in the medium resolution imagery, such as ASTER. The results from the tested commercial area had a very high RMSE value due to the prevalence of shade in the area. It is suggested that the fuzzy classifier based on the nonlinear assumption can handle mixed pixels more effectively than LSMA.     

Estimating impervious surfaces from medium spatial resolution imagery using the self-organizing map and multi-layer perceptron neural networks

The studies of impervious surfaces are important because they are related to many environmental problems, such as water quality, stream health, and the urban heat island effect. Previous studies have discussed that the self-organizing map (SOM) can provide a promising alternative to the multi-layer perceptron (MLP) neural networks for image classification at both per-pixel and sub-pixel level. However, the performances of SOM and MLP have not been compared in the estimation and mapping of urban impervious surfaces. In mid-latitude areas, plant phenology has a significant influence on remote sensing of the environment. When the neural networks approaches are applied, how satellite images acquired in different seasons impact impervious surface estimation of various urban surfaces (such as commercial, residential, and suburban/rural areas) remains to be answered. In this paper, an SOM and an MLP neural network were applied to three ASTER images acquired on April 5, 2004, June 16, 2001, and October 3, 2000, respectively, which covered Marion County, Indiana, United States. Six impervious surface maps were yielded, and an accuracy assessment was performed. The root mean square error (RMSE), the mean average error (MAE), and the coefficient of determination (R²) were calculated to indicate the accuracy of impervious surface maps. The results indicated that the SOM can generate a slightly better estimation of impervious surfaces than the MLP. Moreover, the results from three test areas showed that, in the residential areas, more accurate results were yielded by the SOM, which indicates that the SOM was more effective in coping with the mixed pixels than the MLP, because the residential area prevailed with mixed pixels. Results obtained from the commercial area possessed very high RMSE values due to the prevalence of shade, which indicates that both algorithms cannot handle the shade problem well. The lowest RMSE value was obtained from the rural area due to containing of less mixed pixels and shade. This research supports previous observations that the SOM can provide a promising alternative to the MLP neural network. This study also found that the impact of different map sizes on the impervious surface estimation is significant.     

Comprehensive urban area extraction from multispectral medium spatial resolution remote-sensing imagery based on a novel structural feature

One of the most important applications of remote sensing is urban area analysis in multispectral images. This article addresses a comprehensive method for urban area extraction in the mentioned images using a combination of classic spectral features and a new structural feature. The spectral features are used for eliminating the non-urban land covers including vegetation, water, shadows, and bright soil. Then the proposed structural feature based on the density of the spectral gradient is utilized for the final separation of urban points (UPs) from non-urban points (NUPs). The proposed method is insensitive to the spectral variation of urban areas caused by variant geographical conditions and coordinates. The proposed approach is applied to a wide variety of geographical areas from Landsat 7 Enhanced Thematic Mapper Plus (ETM+) multispectral data including arid and desert land, mountainous land, and plain area. Furthermore, a reference land-cover map from National Land Cover Data 2001 of USA is used as the ground reference for the accuracy assessment of the proposed method. Results analysis shows better performance of the proposed method in all study sites compared with the other spectral indices and the other methods. In addition, the proposed structural feature has shown promising results compared with the spectral indices when used in its own.     

Enhancing multispectral imagery spatial resolution using optimized adaptive PCA and high-pass modulation

ABSTRACT

The pan-sharpening scheme combines high-resolution panchromatic imagery (HRPI) data and low-resolution multispectral imagery (LRMI) data to get a single merged high-resolution multispectral image (HRMI). The pan-sharpened image has extensive information that will promote the efficiency of image analysis methods. Pan-sharpening technique is considered as a pixel-level fusion scheme utilized for enhancing LRMI using HRPI while keeping LRMI spectral information. In this article, an efficient optimized integrated adaptive principal component analysis (APCA) and high-pass modulation (HPM) pan-sharpening method is proposed to get excellent spatial resolution within fused image with minimal spectral distortion. The proposed method is adjusted with multi-objective optimizationto determine the optimal window size and σfor the Gaussian low-pass filter (GLPF) and gain factor utilized for adding the high-pass details extracted from the HRPI to the LRMI principlecomponent of maximum correlation. Optimization results show that if the spatial resolution ratio of HRPI to LRMI is 0.50, then a GLPF of 5 × 5 window size and σ = 1.640 yields HRMI with low spectral distortion and high spatial quality. If the HRPI/LRMI spatial resolution ratio is 0.25, then a GLPF of 7 × 7 window size and σ = 1.686 yields HRMI with low spectral distortion and high spatial quality. Simulation tests demonstrated that the proposed optimized APCA–HPM fusion scheme gives adjustment between spectral quality and spatial quality and has small computational and memory complexity.     

Using very high spatial resolution multispectral satellite sensor imagery to monitor refugee camps

Detailed geographic information is a key factor in decision making processes during refugee relief operations. The forthcoming commercial very high spatial resolution (VHSR) satellite sensors will be capable of acquiring multispectral (MS) images at spatial resolutions of 1m (panchromatic) and 4m (multispectral) of refugee camps and their environment. This work demonstrates how refugee camp environment, area and population can be analysed using a VHSR MS satellite sensor image from the Russian KVR-1000 sensor. This image, with a spatial resolution of 3.3m, was used to study Thailand's Site 2 refugee camps, which were established to accommodate Khmer refugees on the Thai-Kampuchean border. At the time of image acquisition, the total population of Site 2's five refugee camps was close to 143000. The VHSR MS image was found to be suitable for mapping the refugee camp environment and area. A statistically significant linear relationship between camp area and population was determined. Accordingly, the study suggests that VHSR MS images in general may be useful for refugee camp planning and management and points toward the utilization of forthcoming commercial VHSR MS satellite sensor images in humanitarian relief operations.     

Detection and classification of invasive saltcedar through high spatial resolution airborne hyperspectral imagery

We evaluated the performance of airborne HyperSpecTIR (HST) images for detecting and classifying the invasive riparian vegetation saltcedar along the Muddy River in Clark County, Nevada. HyperSpecTIR image reflectance spectra (227 bands, 450–2450 nm) were acquired for the following four vegetation covers: invasive saltcedar, native honey mesquite, grassland patches and crops. We compared five feature reduction approaches: band selection based on Jeffreys–Matusita distance, principal component analysis (PCA), minimum noise fraction (MNF), segmented principal component transform (SPCT) and segmented minimum noise fraction (SMNF). In addition, maximum likelihood (ML) and two spectral angle mapper (SAM) classifiers were applied to all extracted bands or features. Classification accuracies were compared among all classification approaches. Although the overall accuracy of maximal likelihood classifiers generally surpassed that of SAM classifiers, the highest overall accuracy was achieved by a SMNF-SAM combination with adjusted angular thresholds for classes. We concluded that high spectral and spatial resolution imagery can be used to detect and classify invasive saltcedar in this arid area.     

Comparison of hyperspectral imagery with aerial photography and multispectral imagery for mapping broom snakeweed

Broom snakeweed (Gutierrezia sarothrae (Pursh) Britt. & Rusby) is one of the most widespread and abundant rangeland weeds in western North America. The objectives of this study were to evaluate airborne hyperspectral imagery and compare it with aerial colour-infrared (CIR) photography and multispectral digital imagery for mapping broom snakeweed infestations. Airborne hyperspectral imagery along with aerial CIR photographs and digital CIR images was acquired from a rangeland area in south Texas. The hyperspectral imagery was transformed using minimum noise fraction (MNF) and then classified using minimum distance, Mahalanobis distance, maximum likelihood, and spectral angle mapper (SAM) classifiers. The digitized aerial photographs and the digital images were respectively mosaicked as one photographic image and one digital image; these were then classified using the same classifiers. Accuracy assessment showed that the maximum likelihood classifier performed the best for the three types of images. The best overall accuracies for three-class classification maps (snakeweed, mixed woody and mixed herbaceous) were 91.0%, 92.5%, and 95.0%, respectively, for the CIR photographic image, the digital CIR image and the MNF-transformed hyperspectral image. Kappa analysis showed that there were no significant differences in maximum likelihood-based classifications among the three types of images. These results indicate that airborne hyperspectral imagery along with aerial photography and multispectral imagery can be used for monitoring and mapping broom snakeweed infestations on rangelands.     

Model development and comparison to predict softwood and hardwood per cent cover using high and medium spatial resolution imagery

The recent availability of high spatial resolution multispectral scanners provides an opportunity to adapt existing methods and test models to derive spatially explicit forest type and per cent cover information at the Landsat pixel level. A regression modelling methodology was applied for scaling‐up high resolution (IKONOS) to medium spatial resolution satellite imagery (Landsat) to predict softwood and hardwood forest type and density (per cent cover) in a northern Maine study area. Regression relationships (63 different models) were developed and compared. The model variables included vegetation indices and several date (season) combinations of Landsat Enhanced Thematic Mapper Plus (ETM+) imagery (August, September, October and May).

A model incorporating all variables from four dates of Landsat ETM+ imagery produced the highest coefficient of variation in predicting both softwood (0.655) and hardwood cover (0.66). The addition of vegetation indices with the six ETM+ reflected bands did not significantly improve or detract from the regression relationships for any of the multi‐date or single date models examined. A two‐date combination of October and May variables provided an acceptable (and arguably more cost‐effective) model as the adjusted R ² value was 0.645 for softwood and 0.649 for hardwood. A significant result was that all single‐date models produced inferior results with a sharp drop in adjusted R ², compared with the multi‐date seasonal models. This research has demonstrated that the regression models including multi‐date variables produce good results and can provide spatially explicit forest type and stand structure data that has been difficult or infeasible to obtain from medium spatial resolution imagery using traditional classification methods.     

Moment-based feature extraction from high spatial resolution hyperspectral images

In this article, we propose a method for extracting spatio-spectral features from high spatial resolution hyperspectral (HS) images. The method is based on extracting two-dimensional moments from neighbourhoods of pixels. Three different types of moments are considered: geometric, complex Zernike and Legendre. Moments of a given type are extracted from a few principal components (PC) of HS data, and are stacked on the original HS data to form a joint spatio-spectral feature space. These features are classified using a support vector machine (SVM) classifier. The influence of the moments orders and the size of the neighbourhood window on the quality of the extracted features are analysed. A few experiments are conducted on two widely used HS data sets, Pavia University and Salinas. The results demonstrate high capabilities of the proposed method in comparison with some state-of-the-art spatio-spectral HS classification methods.     

Extracting forest canopy structure from spatial information of high resolution optical imagery: tree crown size versus leaf area index

Leaves are the primary interface where energy, water and carbon exchanges occur between the forest ecosystems and the atmosphere. Leaf area index (LAI) is a measure of the amount of leaf area in a stand, and the tree crown size characterizes how leaves are clumped in the canopy. Both LAI and tree crown size are of essential ecological and management value. There is a lot of interest in extracting both canopy structural parameters from remote sensing. The LAI is generally estimated with spectral information from remotely sensed images at relatively coarse spatial resolution. There has been much less success in estimating tree crown size with remote sensing. The recent availability of abundant high spatial resolution imagery from space offers new potential for extracting LAI and tree crown size, particularly in the spatial domain. This study found that the spatial information in Ikonos imagery is highly valuable in estimating both tree crown size and LAI. When the conifer‐ and hardwood‐dominated stands are pooled, tree crown sizes of conifer stands relate best to the ratio of image variance at 2×2 m spatial resolution to that at 3×3 m spatial resolution, while LAI relates best to image variance at 4×4 m spatial resolution. When the conifer‐ and hardwood‐dominated stands are separated, image spatial information estimates tree crown size much better for conifer‐dominated stands than for the hardwood‐dominated stands, while the relationship between image spatial information and LAI is strengthened after the two types of stands are combined. Tree crown size is more sensitive to image spatial resolution than LAI. Image variance is more useful in estimating LAI than normalized difference vegetation index (NDVI) and simple ratio vegetation index (SRVI). Combining both spatial and spectral information provides some improvement in estimating LAI compared with using spatial information alone. Therefore, future efforts to estimate canopy structure with high resolution imagery should also use image spatial information.     

A comparison of three methods for automatic tree crown detection and delineation from high spatial resolution imagery

This article compares the performance of three algorithms representative of published methods for tree crown detection and delineation from high spatial resolution imagery, and demonstrates a standardized accuracy assessment framework. The algorithms – watershed segmentation, region growing and valley-following – were tested on softwood and hardwood sites using Emerge natural colour vertical aerial imagery with 60 cm ground sampled distance and QuickBird panchromatic imagery with an 11? look angle. The evaluation considered both plot-level and individual tree crown detection and delineation results. The study shows that while all three methods reasonably delineate crowns in the softwood stand on the Emerge image, region growing provided the highest accuracies, with producer's and user's accuracy for tree detection reaching 70% and root mean square error for crown diameter estimation of 15%. Crown detection accuracies were lower on the QuickBird image. No algorithm proved accurate for the hardwood stand on either image set (both producer's and user's accuracies < 30%).     

Integrating visible,near-infrared and short-wave infrared hyperspectral and multispectral thermal imagery for geological mapping at Cuprite,Nevada: a rule-based system

Previous research has shown that integrating hyperspectral visible and near-infrared (VNIR) / short-wave infrared (SWIR) with multispectral thermal infrared (TIR) data can lead to improved mineral and rock identification. However, inconsistent results were found regarding the relative accuracies of different classification methods for dealing with the integrated data set. In this study, a rule-based system was developed for integration of VNIR/SWIR hyperspectral data with TIR multispectral data and evaluated using a case study of Cuprite, Nevada. Previous geological mapping, supplemented by field work and sample spectral measurements, was used to develop a generalized knowledge base for analysis of both spectral reflectance and spectral emissivity. The characteristic absorption features, albedo and the location of the spectral emissivity minimum were used to construct the decision rules. A continuum removal algorithm was used to identify absorption features from VNIR/SWIR hyperspectral data only; spectral angle mapper (SAM) and spectral feature fitting (SFF) algorithms were used to estimate the most likely rock type. The rule-based system was found to achieve a notably higher performance than the SAM, SFF, minimum distance and maximum likelihood classification methods on their own.     

Semi-automatic detection and counting of oil palm trees from high spatial resolution airborne imagery

Plantation inventory and management require a range of fine-scale remote-sensing data. Remote-sensing images with high spatial and spectral resolution are an efficient source of such information. This article presents an approach to the extraction and counting of oil palm trees from high spatial resolution airborne imagery data. Counting oil palm trees is a crucial problem in specific agricultural areas, especially in Malaysia. The proposed scheme comprises six major parts: (1) discrimination of oil palms from non-oil palms using spectral analysis, (2) texture analysis, (3) edge enhancement, (4) segmentation process, (5) morphological analysis and (6) blob analysis. The average accuracy obtained was 95%, which indicates that high spatial resolution airborne imagery data with an appropriate assessment technique have the potential to provide us with vital information for oil palm plantation management. Information on the number of oil palm trees is crucial to the ability of plantation management to assess the value of the plantation and to monitor its production.     

Comparison of spectral and spatial windows for local anomaly detection in hyperspectral imagery

A new way of implementing two local anomaly detectors in a hyperspectral image is presented in this study. Generally, most local anomaly detector implementations are carried out on the spatial windows of images, because the local area of the image scene is more suitable for a single statistical model than for global data. These detectors are applied by using linear projections. However, these detectors are quite improper if the hyperspectral dataset is adopted as the nonlinear manifolds in spectral space. As multivariate data, the hyperspectral image datasets can be considered to be low-dimensional manifolds embedded in the high-dimensional spectral space. In real environments, the nonlinear spectral mixture occurs more frequently, and these manifolds could be nonlinear. In this case, traditional local anomaly detectors are based on linear projections and cannot distinguish weak anomalies from background data. In this article, local linear manifold learning concepts have been adopted, and anomaly detection algorithms have used spectral space windows with respect to the linear projection. Output performance is determined by comparison between the proposed detectors and the classic spatial local detectors accompanied by the hyperspectral remote-sensing images. The result demonstrates that the effectiveness of the proposed algorithms is promising to improve detection of weak anomalies and to decrease false alarms.     

Object-based land cover classification of shaded areas in high spatial resolution imagery of urban areas: A comparison study

A significant proportion of high spatial resolution imagery in urban areas can be affected by shadows. Considerable research has been conducted to investigate shadow detection and removal in remotely sensed imagery. Few studies, however, have evaluated how applications of these shadow detection and restoration methods can help eliminate the shadow problem in land cover classification of high spatial resolution images in urban settings. This paper presents a comparison study of three methods for land cover classification of shaded areas from high spatial resolution imagery in an urban environment. Method 1 combines spectral information in shaded areas with spatial information for shadow classification. Method 2 applies a shadow restoration technique, the linear-correlation correction method to create a “shadow-free” image before the classification. Method 3 uses multisource data fusion to aid in classification of shadows. The results indicated that Method 3 achieved the best accuracy, with overall accuracy of 88%. It provides a significantly better means for shadow classification than the other two methods. The overall accuracy for Method 1 was 81.5%, slightly but not significantly higher than the 80.5% from Method 2. All of the three methods applied an object-based classification procedure, which was critical as it provides an effective way to address the problems of radiometric difference and spatial misregistration associated with multisource data fusion (Method 3), and to incorporate thematic spatial information (Method 1).     

Geostatistical and local cluster analysis of high resolution hyperspectral imagery for detection of anomalies

This paper describes a new methodology to detect small anomalies in high resolution hyperspectral imagery, which involves successively: (1) a multivariate statistical analysis (principal component analysis, PCA) of all spectral bands; (2) a geostatistical filtering of noise and regional background in the first principal components using factorial kriging; and finally (3) the computation of a local indicator of spatial autocorrelation to detect local clusters of high or low reflectance values and anomalies. The approach is illustrated using 1 m resolution data collected in and near northeastern Yellowstone National Park. Ground validation data for tarps and for disturbed soils on mine tailings demonstrate the ability of the filtering procedure to reduce the proportion of false alarms (i.e., pixels wrongly classified as target), and its robustness under low signal to noise ratios. In almost all scenarios, the proposed approach outperforms traditional anomaly detectors (i.e., RX detector which computes the Mahalanobis distance between the vector of spectral values and the vector of global means), and fewer false alarms are obtained when using a novel statistic S₂ (average absolute deviation of p-values from 0.5 through all spectral bands) to summarize information across bands. Image degradation through addition of noise or reduction of spectral resolution tends to blur the detection of anomalies, increasing false alarms, in particular for the identification of the least pure pixels. Results from a mine tailings site demonstrate the approach performs reasonably well for highly complex landscape with multiple targets of various sizes and shapes. By leveraging both spectral and spatial information, the technique requires little or no input from the user, and hence can be readily automated.     

Detecting Tamarisk species (Tamarix spp.) in riparian habitats of Southern California using high spatial resolution hyperspectral imagery

Exotic plant invasion is a major environmental and ecological concern and is a particular issue for Mediterranean-type ecosystems. Early detection of invasive plants is crucial for effective weed management. Several studies have explored hyperspectral imagery for mapping invasive plants with promising results. However, only a few extensive or comparative studies about image processing techniques for invasive plant detection have been reported, and even fewer studies have involved very high spatial and spectral resolution imagery. The primary goal of this study was to investigate the utility of very high spatial (0.5 m) and spectral (4 nm) resolution imagery and several classification approaches for detecting tamarisk (Tamarix spp.) infestations, the most problematic invasive plant species in the riparian habitats of southern California.Hierarchical clustering was a particularly effective and efficient statistical method for identifying wavebands and spectral transforms having the greatest discriminatory power. Products resulting from the classification of airborne hyperspectral image data varied by scene, input data type, classifier, and minimum patch size. Overall accuracy of image classification accuracy of products co-varied with commission error rates, such that products having strong agreement with reference data also had a high number of false detections. Integrating the findings from qualitative map analysis, areal proportion statistics, and object-based accuracy assessment indicates that the parallelepiped classifier with several narrow wavebands selected through hierarchical clustering yielded the most accurate and reliable tamarisk classification products.     

Integrating visible, near-infrared and short-wave infrared hyperspectral and multispectral thermal imagery for geological mapping at Cuprite, Nevada

This study investigated the potential value of integrating hyperspectral visible, near-infrared, and short-wave infrared imagery with multispectral thermal data for geological mapping. Two coregistered aerial data sets of Cuprite, Nevada were used: Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) hyperspectral data, and MODIS/ASTER Airborne Simulator (MASTER) multispectral thermal data. Four classification methods were each applied to AVIRIS, MASTER, and a combined set. Confusion matrices were used to assess the classification accuracy. The assessment showed, in terms of kappa coefficient, that most classification methods applied to the combined data achieved a marked improvement compared to the results using either AVIRIS or MASTER thermal infrared (TIR) data alone. Spectral angle mapper (SAM) showed the best overall classification performance. Minimum distance classification had the second best accuracy, followed by spectral feature fitting (SFF) and maximum likelihood classification. The results of the study showed that SFF applied to the combination of AVIRIS with MASTER TIR data are especially valuable for identification of silicified alteration and quartzite, both of which exhibit distinctive features in the TIR region. SAM showed some advantages over SFF in dealing with multispectral TIR data, obtaining higher accuracy in discriminating low albedo volcanic rocks and limestone which do not have unique, distinguishing features in the TIR region.     

Extracting crop radiometric responses from simulated low and high spatial resolution satellite data using a linear mixing model

The images from optical sensors with a broad path width (e.g. NOAA-AVHRR) are used for monitoring vegetation on a regional scale. The European agricultural land uses, which are generally heterogeneous, can be coarsely distinguished by these radiometers. Such sensors, however, do not allow the discrimination of seasonal radiometric changes of a given crop. Some future Earth observation platforms will carry two types of instruments on board. The first instrument will have moderate spatial resolution but a broad path width to allow almost daily observations of the emerged areas. The second will have high spatial resolution and a narrow path width to give the opportunity of making land use thematic maps from the few images recorded per year. The combination of these two types of data allows the medium resolution signal to be unmixed in order to restore the radiometric evolution of a particular crop or of a group of crops.

From the application of a linear mixing model to coarse spatial resolution data, this article presents an unmixing method based on the techniques of multiple regression. This approach has been applied to a simulated coarse resolution dataset to calculate the spectral response of the mixture components.

The most promising results of this first study encourage us to assess the method with real images (e.g. NOAA-AVHRR). Additionally, the results can be seen as an argument in favour of the complementary use of these two types of optical instrument.