Assessment of Landsat 7 Scan Line Corrector-off data gap-filling methods for seagrass distribution mapping

Methods to predict and fill Landsat 7 Scan Line Corrector (SLC)-off data gaps are diverse and their usability is case specific. An appropriate gap-filling method that can be used for seagrass mapping applications has not been proposed previously. This study compared gap-filling methods for filling SLC-off data gaps with images acquired from different dates at similar mean sea-level tide heights, covering the Sungai Pulai estuary area inhabited by seagrass meadows in southern Peninsular Malaysia. To assess the geometric and radiometric fidelity of the recovered pixels, three potential gap-filling methods were examined: (a) geostatistical neighbourhood similar pixel interpolator (GNSPI); (b) weighted linear regression (WLR) algorithm integrated with the Laplacian prior regularization method; and (c) the local linear histogram matching method. These three methods were applied to simulated and original SLC-off images. Statistical measures for the recovered images showed that GNSPI can predict data gaps over the seagrass, non-seagrass/water body, and mudflat site classes with greater accuracy than the other two methods. For optimal performance of the GNSPI algorithm, cloud and shadow in the primary and auxiliary images had to be removed by cloud removal methods prior to filling data gaps. The gap-filled imagery assessed in this study produced reliable seagrass distribution maps and should help with the detection of spatiotemporal changes of seagrasses from multi-temporal Landsat imagery. The proposed gap-filling method can thus improve the usefulness of Landsat 7 ETM+ SLC-off images in seagrass applications.     

Evaluation of gap-filling methods for Landsat 7 ETM+ SLC-off image for LULC classification in a heterogeneous landscape of West Africa

ABSTRACT

The Landsat mission which has existed over five decades has remained at the forefront of providing consistent moderate spatial and temporal resolution optical images of the earth. The failure of the scan line corrector (SLC) on board the Landsat 7 Enhanced Thematic Mapper Plus (ETM+) in May 2003 has permanently resulted in data gaps on each Landsat 7 scene. Due to the obvious negative impacts on the image usability, a number of methods have been developed to fill the no-data areas in the image. This study assessed the performance of four Landsat 7 ETM+ SLC-off gap-filling methods in a highly heterogeneous landscape of West Africa for two different seasons (dry and rainy). The methods considered are: (1) Weighted Linear Regression (WLR) integrated with Laplacian Prior Regularization Method (LPRM), (2) Localised Linear Histogram Matching (LLHM), (3) Neighbourhood Similar Pixel Interpolator (NSPI) and (4) Geostatistical Neighbourhood Similar Pixel Interpolator (GNSPI). All the images used were Landsat 7 ETM+ SLC-off images, temporally close and from the same season for each set of time step. Visual comparison, mean, and standard deviations of the histograms of all bands of only the filled areas were used to assess the results. Additionally, overall accuracy (OA), kappa coefficient (κ), and balanced accuracy (BA) per class were used to evaluate a land use/cover (LULC) classification based on the gap-filled images. Visually, all the four methods were able to completely fill the gaps in the Landsat 7 ETM+ SLC-off image. They all look similar and spatially continuous with no anomalies or artefacts on them. The histograms from each band for only the filled areas for all the four methods also gave similar means and standard deviations in most cases. All the four gap-filling methods provided satisfactory results (OA >96% and κ> 0.937 in all methods for images in the dry season and OA >93% and κ> 0.877 for the image in the rainy season) in the land cover classification considering the complexity of the study area. But the GNSPI was superiority in all cases with the highest OA of 97.1% and κ of 0.947 in the dry season and OA of 94.6% and κ of 0.899 in the rainy season. This implies that the GNSPI is more robust in gap filling of Landsat 7 ETM+ SLC-off images than the other three methods in a heterogeneous landscape of West Africa regardless of the season. This study suggests that gap filling of Landsat 7 ETM+ SLC-off images will help to increase the number of Landsat images needed to build time-series data for a data-scarce region such as West Africa.     

A comparison of gap-filling approaches for Landsat-7 satellite data

The purpose of this study is to assess the relative performance of four different gap-filling approaches across a range of land-surface conditions, including both homogeneous and heterogeneous areas as well as in scenes with abrupt changes in landscape elements. The techniques considered in this study include: (1) Kriging and co-Kriging; (2) geostatistical neighbourhood similar pixel interpolator (GNSPI); (3) a weighted linear regression (WLR) algorithm; and (4) the direct sampling (DS) method. To examine the impact of image availability and the influence of temporal distance on the selection of input training data (i.e. time separating the training data from the gap-filled target image), input images acquired within the same season (temporally close) as well as in different seasons (temporally far) to the target image were examined, as was the case of using information only within the target image itself. Root mean square error (RMSE), mean spectral angle (MSA), and coefficient of determination (R²) were used as the evaluation metrics to assess the prediction results. In addition, the overall accuracy (OA) and kappa coefficient (kappa) were used to assess a land-cover classification based on the gap-filled images. Results show that all of the gap-filling approaches provide satisfactory results for the homogeneous case, with R² > 0.93 for bands 1 and 2 in all cases and R² > 0.80 for bands 3 and 4 in most cases. For the heterogeneous example, GNSPI performs the best, with R² > 0.85 for all tested cases. WLR and GNSPI exhibit equivalent accuracy when a temporally close input image is used (i.e. WLR and GNSPI both have an R² equal to 0.89 for band 1). For the case of abrupt changes in scene elements or in the absence of ancillary data, the DS approach outperforms the other tested methods.     

Robust prediction of time-integrated NDVI

The time-integrated normalized difference vegetation index (iNDVI) provides key remote-sensing-derived information on the interactions between vegetation growth, climatic and soil conditions, and land use. Using a time-series of Landsat imagery obtained for Queensland, Australia, it has been demonstrated how robust geostatistics can be used to predict iNDVI. This approach is novel because it explicitly quantifies the uncertainty of prediction and uses Winsorizing, a data-censoring method, to minimize the distorting effects of outliers. Robust prediction of iNDVI, as opposed to non-robust prediction, was justifiable in 79% of the study area, highlighting the need for methods that deal with outliers in time-series analysis of remotely sensed imagery. There was a strong coarse-scale association between Queensland’s bioregions and iNDVI, and also between bioregion and the rain-induced difference in iNDVI through time (effects that were significant at p < 0.001 in both cases). At a finer spatial scale, prediction of iNDVI also appeared to be a promising way to distinguish long-term cropping land from adjacent long-term grazing land (effect significant at p < 0.001). The method is tied to a set of assumptions concerning image radiometry, cloud detection, variogram estimation, and variable additivity. The first two are fundamental remote-sensing issues that can be improved with additional labour; the last two can be improved statistically but would greatly increase the processing time per pixel. Robust geostatistical analysis of time-series has immediate relevance to gap-filling of SLC-off Landsat 7 Enhanced Thematic Mapper Plus (ETM+) imagery, and for generating novel covariates for digital soil mapping.     

Production of Landsat ETM+ reference imagery of burned areas within Southern African savannahs: comparison of methods and application to MODIS

Accurate production of regional burned area maps are necessary to reduce uncertainty in emission estimates from African savannah fires. Numerous methods have been developed that map burned and unburned surfaces. These methods are typically applied to coarse spatial resolution (1 km) data to produce regional estimates of the area burned, while higher spatial resolution (<30 m) data are used to assess their accuracy with little regard to the accuracy of the higher spatial resolution reference data. In this study we aimed to investigate whether Landsat Enhanced Thematic Mapper (ETM+)‐derived reference imagery can be more accurately produced using such spectrally informed methods. The efficacy of several spectral index methods to discriminate between burned and unburned surfaces over a series of spatial scales (ground, IKONOS, Landsat ETM+ and data from the MOderate Resolution Imaging Spectrometer, MODIS) were evaluated. The optimal Landsat ETM+ reference image of burned area was achieved using a charcoal fraction map derived by linear spectral unmixing (k = 1.00, a = 99.5%), where pixels were defined as burnt if the charcoal fraction per pixel exceeded 50%. Comparison of coincident Landsat ETM+ and IKONOS burned area maps of a neighbouring region in Mongu (Zambia) indicated that the charcoal fraction map method overestimated the area burned by 1.6%. This method was, however, unstable, with the optimal fixed threshold occurring at >65% at the MODIS scale, presumably because of the decrease in signal‐to‐noise ratio as compared to the Landsat scale. At the MODIS scale the Mid‐Infrared Bispectral Index (MIRBI) using a fixed threshold of >1.75 was determined to be the optimal regional burned area mapping index (slope = 0.99, r ² = 0.95, SE = 61.40, y = Landsat burned area, x = MODIS burned area). Application of MIRBI to the entire MODIS temporal series measured the burned area as 10 267 km² during the 2001 fire season. The char fraction map and the MIRBI methodologies, which both produced reasonable burned area maps within southern African savannah environments, should also be evaluated in woodland and forested environments.     

Comparison of Landsat 8 OLI and Landsat 7 ETM+ for estimating grassland LAI using model inversion and spectral indices: case study of Mpumalanga,South Africa

The leaf area index (LAI) is the key biophysical indicator used to assess the condition of rangeland. In this study, we investigated the implications of narrow spectral response, high radiometric resolution (12 bits), and higher signal-to-noise ratio of the Landsat 8 Operational Land Imager (OLI) sensor for the estimation of LAI. The Landsat 8 LAI estimates were compared to that of its predecessors, namely Landsat 7 Enhanced Thematic Mapper Plus (ETM+) (8 bits). Furthermore, we compared the radiative transfer model (RTM) and spectral indices approaches for estimating LAI on rangeland systems in South Africa. The RTM was inverted using artificial neural network (ANN) and lookup table (LUT) algorithms. The accuracy of the models was higher for Landsat 8 OLI, where ANN (root mean squared error, RMSE = 0. 13; R² = 0. 89), LUT (RMSE = 0. 25; R² = 0. 50), compared to Landsat 7 ETM+, where ANN (RMSE = 0. 35; R² = 0. 60), LUT (RMSE = 0. 38; R² = 0. 50). Compared to an empirical approach, the RTM provided higher accuracy. In conclusion, Landsat 8 OLI provides an improvement for the estimation of LAI over Landsat 7 ETM+. This is useful for rangeland monitoring.     

Multi-temporal MODIS-Landsat data fusion for relative radiometric normalization, gap filling, and prediction of Landsat data

A semi-physical fusion approach that uses the MODIS BRDF/Albedo land surface characterization product and Landsat ETM+ data to predict ETM+ reflectance on the same, an antecedent, or subsequent date is presented. The method may be used for ETM+ cloud/cloud shadow and SLC-off gap filling and for relative radiometric normalization. It is demonstrated over three study sites, one in Africa and two in the U.S. (Oregon and Idaho) that were selected to encompass a range of land cover land use types and temporal variations in solar illumination, land cover, land use, and phenology. Specifically, the 30 m ETM+ spectral reflectance is predicted for a desired date as the product of observed ETM+ reflectance and the ratio of the 500 m surface reflectance modeled using the MODIS BRDF spectral model parameters and the sun-sensor geometry on the predicted and observed Landsat dates. The difference between the predicted and observed ETM+ reflectance (prediction residual) is compared with the difference between the ETM+ reflectance observed on the two dates (temporal residual) and with respect to the MODIS BRDF model parameter quality. For all three scenes, and all but the shortest wavelength band, the mean prediction residual is smaller than the mean temporal residual, by up to a factor of three. The accuracy is typically higher at ETM+ pixel locations where the MODIS BRDF model parameters are derived using the best quality inversions. The method is most accurate for the ETM+ near-infrared (NIR) band; mean NIR prediction residuals are 9%, 12% and 14% of the mean NIR scene reflectance of the African, Oregon and Idaho sites respectively. The developed fusion approach may be applied to any high spatial resolution satellite data, does not require any tuning parameters and so may be automated, is applied on a per-pixel basis and is unaffected by the presence of missing or contaminated neighboring Landsat pixels, accommodates for temporal variations due to surface changes (e.g., phenological, land cover/land use variations) observable at the 500 m MODIS BRDF/Albedo product resolution, and allows for future improvements through BRDF model refinement and error assessment.     

A multi-sensor approach for improving biodiversity estimation in the Hyrcanian mountain forest,Iran

ABSTRACT

Modelling tree biodiversity in mountainous forests using remote-sensing data is challenging because forest composition and structure change along elevation. Topographic variations also affect vegetation’s spectral and backscattering behaviour. We demonstrate the potential of multi-source integration to tackle this challenge in a mountainous part of the Hyrcanian forest in Iran. This forest is a remnant of a deciduous broadleaved forest with heterogeneous structure affected by natural and anthropogenic factors. The multi-source approach (i.e. Landsat Enhanced Thematic Mapper Plus (ETM +), Advanced Land Observing Satellite/ Phased Array type L-band Synthetic Aperture Radar (ALOS/PALSAR), and topographic variables) allows us to propose a biodiversity estimation model using partial least square regression (PLSR) calibrated and validated with limited field data. The effective number of species was calculated based on field measurements of the biodiversity in the study area. In order to model species diversity in more homogeneous extrinsic environmental conditions, we divided data into two groups with relatively uniform slope values. In each slope group, we modelled the correlation between observed biodiversity and satellite-derived data. For that, we followed three scenarios: (A) multispectral Landsat ETM + alone, (B) ALOS/PALSAR alone, and (C) inclusion of both sensors. In each scenario, elevation and slope data were also considered as predictors. We observed that in all scenarios, coefficient of determination (R²) in gentler slopes was higher than that in areas with steeper slopes (average difference in R²: ?R² = 0.21). The highest correlation was achieved by inclusion of synthetic aperture radar (SAR) and ETM + (R² = 0.87). The results clearly confirm that the multi-source remote-sensing approach can provide a practical estimate of biodiversity across the Hyrcanian forest and potentially in other deciduous broadleaved forests in complex terrain.     

Recovering missing pixels for Landsat ETM + SLC-off imagery using HJ-1A /1B as auxiliary data

To improve the usability of Enhanced Thematic Mapper Plus (ETM+) scan line corrector (SLC)-off data, this article proposes using HJ-1A/1B imagery as auxiliary (i.e. reference) data to recover the SLC-off ETM+ data. The least-median-of-squares (LMedS) method is newly proposed to recover missing pixels of Landsat 7 by removing the variant or abnormal digital number values. In particular, for the visible and near-infrared bands, using HJ-1A/1B for recovery has three clear advantages: the same spatial resolution, similar spectral resolution, and approximate temporal resolution. The experiments show that all of the reference-recovery methods are better than the non-reference-recovery method. The results of using of auxiliary data in reference-recovery methods, from best to worst, are Landsat 8, HJ-1A/1B, and Landsat 7. However, for recovering missing pixels, HJ-1A/1B is superior to the ETM+ auxiliary data due to the shorter time interval in Landsat 7 (a few hours). Hence, HJ-1A/1B should be considered a useful auxiliary data to recover ETM+ SLC-off imagery data.     

A simple and effective method for filling gaps in Landsat ETM+ SLC-off images

The scan-line corrector (SLC) of the Landsat 7 Enhanced Thematic Mapper Plus (ETM+) sensor failed in 2003, resulting in about 22% of the pixels per scene not being scanned. The SLC failure has seriously limited the scientific applications of ETM+ data. While there have been a number of methods developed to fill in the data gaps, each method has shortcomings, especially for heterogeneous landscapes. Based on the assumption that the same-class neighboring pixels around the un-scanned pixels have similar spectral characteristics, and that these neighboring and un-scanned pixels exhibit similar patterns of spectral differences between dates, we developed a simple and effective method to interpolate the values of the pixels within the gaps. We refer to this method as the Neighborhood Similar Pixel Interpolator (NSPI). Simulated and actual SLC-off ETM+ images were used to assess the performance of the NSPI. Results indicate that NSPI can restore the value of un-scanned pixels very accurately, and that it works especially well in heterogeneous regions. In addition, it can work well even if there is a relatively long time interval or significant spectral changes between the input and target image. The filled images appear reasonably spatially continuous without obvious striping patterns. Supervised classification using the maximum likelihood algorithm was done on both gap-filled simulated SLC-off data and the original “gap free” data set, and it was found that classification results, including accuracies, were very comparable. This indicates that gap-filled products generated by NSPI will have relevance to the user community for various land cover applications. In addition, the simple principle and high computational efficiency of NSPI will enable processing large volumes of SLC-off ETM+ data.     

Estimating east Mediterranean forest parameters using Landsat ETM

The conservation of Jordan's Mediterranean forest requires the use of remote sensing. Among the most important parameters needed are the crown-cover percentage (C) and above-ground biomass (A). This study aims to: (1) identify the best predictor(s) of C using Landsat Enhanced Thematic Mapper (ETM) bands and the derived transformed normalized difference vegetation index (TNDVI); (2) determine if C is a good predictor of A, volume (V), Shannon diversity index (S) and basal area (B); and (3) generate maps of all these parameters. A Landsat ETM image, aerial photographs and ground surveys are used to model C using multiple regression. C is then modelled to A, V, S and B using linear regression. The relationship between C and Landsat ETM bands (1 and 7) plus the TNDVI is significantly high (coefficient of determination R ² = 0.8) and is used to produce the C map. The generated C map is used to predict A (R ² = 0.56), V (R ² = 0.58), S (R ² = 0.50) and B (R ² = 0.43). Cross validation for the predicted C map (cross-validation error = 5.3%) and for the predicted forest-parameter maps (cross-validation error = 13.7%–19.9%) shows acceptable error levels. Results indicate that Jordan's east Mediterranean forest parameters can be mapped and monitored for biomass accumulation and carbon dioxide (CO₂) flux using Landsat ETM images.     

Defining the Landsat TM/ETM+ and CHRIS/PROBA spectral regions in which turbidity can be retrieved in inland waterbodies using field spectroscopy

This article describes the results obtained by an existing campaign in which in situ spectroradiometric measurements using a GER1500 field spectroradiometer, Secchi disk depth, and turbidity measurements (using a portable turbidity meter) were acquired at Asprokremmos Reservoir in Paphos District, Cyprus. Field spectroradiometric and water quality data span 18 sampling campaigns during the period May 2010–October 2010. By applying several regression analyses between ‘In-Band’ mean reflectance values against turbidity values for all spectral bands corresponding to Landsat TM/ETM+ (Bands 1 to 4) and CHRIS/PROBA (Bands A1 to A62), the highest correlation was found for Landsat TM/ETM+ Band 3 (R² = 0.85) and for CHRIS/PROBA Bands A30 to A32 (R² = 0.90).     

Monitoring boreal forest leaf area index across a Siberian burn chronosequence: a MODIS validation study

Landscapes containing differing amounts of ecological disturbance provide an excellent opportunity to validate and better understand the emerging Moderate Resolution Imaging Spectrometer (MODIS) vegetation products. Four sites, including 1‐year post‐fire coniferous, 13‐year post‐fire deciduous, 24‐year post‐fire deciduous, and >100 year old post‐fire coniferous forests, were selected to serve as a post‐fire chronosequence in the central Siberian region of Krasnoyarsk (57.3°N, 91.6°E) with which to study the MODIS leaf area index (LAI) and vegetation index (VI) products. The collection 4 MODIS LAI product correctly represented the summer site phenologies, but significantly underestimated the LAI value of the >100 year old coniferous forest during the November to April time period. Landsat 7‐derived enhanced vegetation index (EVI) performed better than normalized difference vegetation index (NDVI) to separate the deciduous and conifer forests, and both indices contained significant correlation with field‐derived LAI values at coniferous forest sites (r ² = 0.61 and r ² = 0.69, respectively). The reduced simple ratio (RSR) markedly improved LAI prediction from satellite measurements (r ² = 0.89) relative to NDVI and EVI. LAI estimates derived from ETM+ images were scaled up to evaluate the 1 km resolution MODIS LAI product; from this analysis MODIS LAI overestimated values in the low LAI deciduous forests (where LAI<5) and underestimated values in the high LAI conifer forests (where LAI>6). Our results indicate that further research on the MODIS LAI product is warranted to better understand and improve remote LAI quantification in disturbed forest landscapes over the course of the year.     

Settlement location and population density estimation in rugged terrain using information derived from Landsat ETM and SRTM data

It is useful to have a disaggregated population database at uniform grid units in disaster situations. This study presents a method for settlement location probability and population density estimations at a 90 m resolution for northern Iraq using the Shuttle Radar Topographic Mission (SRTM) digital terrain model and Landsat Enhanced Thematic Mapper satellite imagery. A spatial model each for calculating the probability of settlement location and for estimating population density is described. A randomly selected subset of field data (equivalent to 50%) is first analysed for statistical links between settlement location probability and population density; and various biophysical features which are extracted from Landsat or SRTM data. The model is calibrated using this subset. Settlement location probability is attributed to the distance from roads and water bodies and land cover. Population density can be estimated based upon land cover and topographic features. The Landsat data are processed using a segmentation and subsequent feature–based classification approach making this method robust to seasonal variations in imagery and therefore applicable to a time series of images regardless of acquisition date. The second half of the field data is used to validate the model. Results show a reasonable estimate of population numbers (r = 0.205, p<0.001) for both rural and urban settlements. Although there is a strong overall correlation between the results of this and the LandScan model (r = 0.464, p<0.001), this method performs better than the 1 km resolution LandScan grid for settlements with fewer than 1000 people, but is less accurate for estimating population numbers in urban areas (LandScan rural r = 0.181, p<0.001; LandScan urban r = 0.303, p<0.001). The correlation between true urban population numbers is superior to that of LandScan however when the 90 m grid values are summed using a filter which corresponds to the LandScan spatial resolution (r = 0.318, p<0.001).     

Multitemporal assessment of the geomorphologic evolution of the Restinga of Marambaia,Rio de Janeiro,Brazil

The Restinga of Marambaia is an emerged sand bar located between the Sepetiba Bay and the South Atlantic Ocean, on the south‐east coast of Brazil. The objective of this study was to observe the geomorphologic evolution of the coastal zone of the Restinga of Marambaia using multitemporal satellite images acquired by multisensors from 1975 to 2004. The images were digitally segmented by a region growth algorithm and submitted to an unsupervised classification procedure (ISOSEG) followed by a raster edit based on visual interpretation. The image time‐series showed a general trend of decrease in the total sand bar area with values varying from 80.61 km² in 1975 to 78.15 km² in 2004. The total area calculation based on the 1975 and 1978 Landsat MSS data was shown to be super‐estimated in relation to the Landsat TM, Landsat ETM+, and CBERS‐2 CCD data. These differences can also be associated to the relatively poorer spatial resolution of the MSS data, nominally 79 m, against the 20 m of the CCD data and 30 m of the TM and ETM+ data. For the estimates of the width in the central portion of the sand bar the variation was from 158 m (1975) to 100 m (2004). The formation of a spit in the northern region of the study area was visually observed. The area of the spit was estimated, with values varying from 0.82 km² (1975) to 0.55 km² (2004).     

Distorted spatial warping to compress large,high-resolution rasters for remote-sensing classification

Decisions made early in the data preparation phases of remote-sensing classification projects set fundamental limits on the value to society of the final products. The often-used approach of degrading/down-sampling high-resolution (e.g. 1 m pixel size) imagery to match lower-resolution data (e.g. Landsat 30 m) through averaging or majority-rule solves the problem of aligning pixels across bands of differing resolution, but does so by forgoing all ability to detect features smaller than 30 m in addition to potentially discarding up to 99% of the information content of the high-resolution data. The alternative of up-sampling coarser-resolution data into smaller-sized synthetic pixels creates its own set of problems, including potentially enormous file sizes, likely absence of meaningful variation over small spatial scales (which may generate matrix singularities fatal to the maximum likelihood classifier), and no assurance of meaningful improvement in classification accuracy despite guaranteed increases in computational time and resource requirements. We propose a new ‘warped space compression technique’ as a variation of vector quantization that analyses local variability in the finest-resolution data available to define acceptable pixel-based neighbourhood (N × N) sizes over which data can be averaged while minimizing overall information loss. Alternative neighbourhoods are aligned so that nine smaller ones nest within each progressively larger one as 3 × 3 squares, resulting in local data compression options of 3 × 3 (ninefold), 9 × 9 (81-fold), 27 × 27 (729-fold), and 81 × 81 (6561-fold). Our transformation process to ‘warped space’ created spatially distorted images with jagged east edges and little visually discernible relationship to the original data. We achieved compressions of 48- to 138-fold in disc storage and 292- to 785-fold in actual numbers of non-null pixels through our choice of cut-off values for accepting 3 × 3, 9 × 9, 27 × 27, or 81 × 81 neighbourhoods of tolerable variability, while otherwise retaining full (1 m) resolution data in regions three cells wide by three cells high. Medium-resolution data (e.g. Landsat 30 m) can be translated into the warped space defined by high-resolution data and composited with it for conducting remote-sensing classifications. When applied to a 71-band, 55-class remote-sensing classification of a 25,500 km² region centred on the Willamette Valley of Oregon, USA, classification accuracy increased from 64.4% in normal space to 71.3% in warped space. Unsupervised classification in warped space identified several additional categories that could be appended to the 55 existing ground-truth classes, leading to further increases in accuracy. Warped-space compression may be particularly beneficial for ecological studies where it could maintain high resolution in features of interest such as riparian buffers without creating exorbitantly large data files.     

Statistical model development and estimation of suspended particulate matter concentrations with Landsat 8 OLI images of Dongting Lake,China

Suspended particulate matter (SPM) is a dominant water constituent of case-II waters, and SPM concentration (C_SPM) is a key parameter describing water quality. This study, using Landsat 8 Operational Land Imager (OLI) images, aimed to develop the C_SPM retrieval models and further to estimate the C_SPM values of Dongting Lake. One Landsat 8 OLI image and 53 C_SPM measurements were employed to calibrate Landsat 8-based C_SPM retrieval models. The C_SPM values derived from coincident Landsat 8 OLI and Moderate Resolution Imaging Spectroradiometer (MODIS) images were compared to validate calibrated Landsat 8-based C_SPM models. After the best stable Landsat 8-based C_SPM retrieval model was further validated using an independent Landsat 8 OLI image and its coincident C_SPM measurements, it was applied to four Landsat 8 OLI images to retrieve the C_SPM values in the South and East Dongting Lake. Model calibration results showed that two exponential models of the red band explained 61% (estimated standard error (SE) = 7.96 mg l^–1) and 67% (SE = 6.79 mg l^–1) of the variation of C_SPM; two exponential models of the red:panchromatic band ratio obtained 81% (SE = 5.48 mg l^–1) and 77% (SE = 4.96 mg l^–1) fitting accuracy; and four exponential and quadratic models of the infrared band explained 72–83% of the variation of C_SPM (SE = 5.18–5.52 mg l^–1). By comparing the MODIS- and Landsat 8-based C_SPM values, an exponential model of the Landsat 8 OLI red band (C_SPM = 1.1034 × exp(23.61 × R)) obtained the best consistent C_SPM estimations with the MODIS-based model (r = 0.98, p < 0.01), and its further validation result using an independent Landsat 8 OLI image showed a significantly strong correlation between the measured and estimated C_SPM values at a significance level of 0.05 (r = 0.91, p < 0.05). The C_SPM spatiotemporal distribution derived from four Landsat 8 images revealed a clear spatial distribution pattern of C_SPM in the South and East Dongting Lake, which was caused by natural and anthropogenic factors together. This study confirmed the potential of Landsat 8 OLI images in retrieving C_SPM and provided a foundation for retrieving the spatial distribution of C_SPM accurately from this new data source in Dongting Lake.     

Detecting an invasive shrub in a deciduous forest understory using late‐fall Landsat sensor imagery

Landsat TM and ETM+ imagery was used to distinguish areas of high vs. low cover of Amur honeysuckle (Lonicera maackii), taking advantage of the late leaf retention of this invasive shrub. L. maackii cover was measured in eight stands and compared to 15 Landsat 5 TM and Landsat 7 ETM+ images from spring and autumn dates from 1999 to 2006. Jeffries–Matusita (JM) distance calculations showed potential separability between high vs. low/zero cover classes of L. maackii on some late fall images. The Soil Adjusted Atmospheric Resistant Vegetation Index (SARVI2) revealed higher levels of green biomass in high L. maackii cover plots than low/zero cover plots for November images only. These findings justify further investigation of the effectiveness of late fall images to map the historical spread of L. maackii and other forest understory invasives with similar phenology.     

Semi‐automated methods for mapping wetlands using Landsat ETM+ and SRTM data

The overarching goal of this study was to develop a comprehensive methodology for mapping natural and human‐made wetlands using fine resolution Landsat enhanced thematic mapper plus (ETM+), space shuttle radar topographic mission digital elevation model (SRTM DEM) data and secondary data. First, automated methods were investigated in order to rapidly delineate wetlands; this involved using: (a) algorithms on SRTM DEM data, (b) thresholds of SRTM‐derived slopes, (c) thresholds of ETM+ spectral indices and wavebands and (d) automated classification techniques using ETM+ data. These algorithms and thresholds using SRTM DEM data either over‐estimated or under‐estimated stream densities (S _d) and stream frequencies (S _f), often generating spurious (non‐existent) streams and/or, at many times, providing glaring inconsistencies in the precise physical location of the streams. The best of the ETM+‐derived indices and wavebands either had low overall mapping accuracies and/or high levels of errors of omissions and/or errors of commissions.

Second, given the failure of automated approaches, semi‐automated approaches were investigated; this involved the: (a) enhancement of images through ratios to highlight wetlands from non‐wetlands, (b) display of enhanced images in red, green, blue (RGB) false colour composites (FCCs) to highlight wetland boundaries, (c) digitizing the enhanced and displayed images to delineate wetlands from non‐wetlands and (d) classification of the delineated wetland areas into various wetland classes. The best FCC RGB displays of ETM+ bands for separating wetlands from other land units were: (a) ETM+4/ETM+7, ETM+4/ETM+3, ETM+4/ETM+2, (b) ETM+4, ETM+3, ETM+5 and (c) ETM+3, ETM+2, ETM+1. In addition, the SRTM slope threshold of less than 1% was very useful in delineating higher‐order wetland boundaries. The wetlands were delineated using the semi‐automated methods with an accuracy of 96% as determined using field‐plot data.

The methodology was evaluated for the Ruhuna river basin in Sri Lanka, which has a diverse landscape ranging from sea shore to hilly areas, low to very steep slopes (0° to 50°), arid to semi‐arid zones and rain fed to irrigated lands. Twenty‐four per cent (145 733 ha) of the total basin area was wetlands as a result of a high proportion of human‐made irrigated areas, mainly under rice cropping. The wetland classes consisted of irrigated areas, lagoons, mangroves, natural vegetation, permanent marshes, salt pans, lagoons, seasonal wetlands and water bodies. The overall accuracies of wetland classes varied between 87% and 94% (K _hat = 0.83 to 0.92) with errors of omission less than 13% and errors of commission less than 1%.     

The effects of topographic correction and gap filling in imagery on the detection of tropical forest disturbances using a Landsat time series in Myanmar

In this study, we evaluated the effects of topographic correction and gap filling of Landsat Enhanced Thematic Mapper Plus (ETM+) images on the accuracy of forest change detection through a trajectory-based approach. Four types of Landsat time series stacks (LTSS) were generated. These stacks resulted from combinations of topographically corrected and uncorrected imagery combined with gap-filled and unfilled stacks. These combinations of stacks were then used as input into a trajectory-based change detection. The results of change detection from trajectory-based analysis using these LTSS were compared in order to assess the effects of both topographic correction and gap-filling procedures on the ability to detect forest disturbances. The results showed that overall accuracies of change detection were improved after gap filling (10.5% and 7.5%), but were only slightly improved after topographic correction (3.6% and 0.6%). Although the gap-filling process introduced some uncertainty that might have caused false change detection, the number of pixels whose detection of disturbance was enhanced after gap filling exceeded those detecting false change. The results also showed that the topographic correction did not contribute much to improve the change detection in this study area. However, topographic correction has a potential to increase the accuracy of change detection in areas of more rugged terrain and steep slopes. This is because a direct relationship between the slope of the topography with topographic correction and an enhanced detection of disturbance in pixels from year to year was observed in this study. For robust change detection, we recommend that a gap-filling process should be included in the trajectory-based analysis procedures such as the one used in this study where a single image per year is used to characterize change. We also recommend that in areas of rugged terrain, a topographic correction in the image pre-processing should be implemented.