Monitoring rubber plantation distribution on Hainan Island using Landsat OLI imagery

Rubber expansion has been widely explored at various scales due to its increasing environmental and socio-economic impact and the global demand for natural rubber products. However, accurate and up-to-date maps and large-scale monitoring of rubber plantation spatial distribution are not yet available. In this article, we developed a simple algorithm for rapidly and accurately mapping rubber plantations on Hainan Island, China, by combining survey samples and Landsat 8 Operational Land Imager data from 2014 in order to understand their spatial distribution. The results showed that rubber plantations are distinguishable from other land-cover types by band value changes, vegetation index changes, and phenological phase changes (defoliation and foliation) and can be accurately extracted from multi-temporal Landsat images using a decision tree method and Google Earth Engine. This method results in a high overall classification accuracy of 92.17% with a corresponding κ of 84.33%. Rubber plantations are concentrated in the northwest of Hainan Island and gradually decrease to the south and east. Most plantations are found in the 50–500-m elevation range, with few found outside this range. We believe that the proposed approach will have significant implications for mapping and monitoring rubber spatial distribution at regional scales.     

Mapping burn severity and burning efficiency in California using simulation models and Landsat imagery

Uncertainties in burning efficiency (BE) estimates can lead to large errors in fire emission quantification (from 23% to 46%). One of the main causes of these errors is the spatial variability of fuel consumption within burned areas. This paper studies whether burn severity (BS) maps can be used to improve BE assessment. A burn severity map of two large fires in California was obtained by inverting a simulation model constrained by post-fire observations from Landsat TM imagery. Model output values of BS were validated against field measurements, obtaining a high correlation (R² = 0.85) and low errors (Root Mean Square Error, RMSE = 0.14) throughout a wide range of BS levels. The BS map obtained was then used to adjust BE reference values per vegetation type found in the area before the fire. The adjusted burning efficiency (BEadj) was compared to the burned biomass, which was estimated by subtracting vegetation indices from pre- and post-fire images. Results showed a high correlation for conifers (R² = 0.75) and hardwoods (R² = 0.73), and a moderate correlation (R² ∼ 0.5) for shrubs and grasslands. In general, for all vegetation types BEadj performed better (R² = 0.4-0.75) than literature-based BE (R² < 0.0001). This study demonstrates: (i) the consistency of the simulation model inversion for BS estimation in temperate ecosystems, and (ii) the improvement of BE estimation when the spatial variability of the combustion was quantified in terms of BS.     

Surface soil moisture estimation over dense crop using Envisat ASAR and Landsat TM imagery: an approach

This study focuses on developing a new method of surface soil moisture estimation over wheat fields using Environmental Satellite Advanced Synthetic Aperture Radar (Envisat ASAR) and Landsat Thematic Mapper (TM) data. The Michigan Microwave Canopy Scattering (MIMICS) model was used to simulate wheat canopy backscattering coefficients from experiment plots at incidence angles of 23° (IS2) and 43.9° (IS7). Based on simulated data, the scattering characteristics of a wheat canopy were first investigated in order to derive an optimal configuration of polarization (HH) and incidence angle (IS2) for soil moisture estimation. Then a parametric model was developed to describe wheat canopy backscattering at the optimal configuration. In addition, direct backscattering and two-way transmissivity of wheat crowns were derived from the TM normalized difference vegetation index (NDVI); direct ground backscattering was derived from surface soil moisture and TM NDVI; and backscattering from double scattering was derived from total backscattering. A semi-empirical model for soil moisture estimation was derived from the parametric model. Coefficients in the semi-empirical model were obtained using a calibration approach based on measured soil moisture, ASAR, and TM data. A validation of the model was performed over the experimental area. In this study, the root mean square error (RMSE) for the estimated soil moisture was 0.041 m³ m^?3, and the correlation coefficient between the measured and estimated soil moisture was 0.84. The experimental results indicate that the semi-empirical model could improve soil moisture estimation compared to an empirical model.     

An automated approach for segmenting and classifying a large sample of multi-date Landsat imagery for pan-tropical forest monitoring

The TREES-3 project of the Joint Research Centre aims at assessing tropical forest cover changes for the periods 1990-2000 and 2000-2010 using a sample-based approach. This paper refers to the 1990-2000 assessment. Extracts of Landsat satellite imagery (20 km × 20 km) are analyzed for these reference dates for more than 4000 sample sites distributed systematically across the tropical belt. For the processing and analysis of such a large amount of satellite imagery a robust methodological approach for forest mapping and change detection has been developed. This approach comprises two automated steps of multi-date image segmentation and object-based land cover classification (based on a supervised spectral library), followed by an intense phase of visual control and expert refinement. Image segmentation is done at two spatial scales, introducing the concept of a minimum mapping unit via the automated selection of a site-specific scale parameter. The automated segmentation of land cover polygons and the pre-classification of land cover types mainly aim at avoiding manual delineation and at reducing the efforts of visual interpretation of land cover to a reasonable level, making the analysis of 4000 sample sites feasible. The importance of visual control and correction can be perceived when comparing to the initial automatic classification result: about 20% of the polygon labels were changed through expert knowledge by visual interpretation. The component of visual refinement of the mapping results had also a notable impact on forest area and change estimates: for a set of sample sites in Southeast Asia (~ 90% of all sites of SE-Asia) the rate of change in tree cover (deforestation) was assessed at 0.9% and 1.6% before and after visual control, respectively.     

Population estimation from Landsat imagery

The paper introduces a method of population estimation using the Landsat MSS data. The radiance in the four spectral bands, detected by the multi spectral scanner (MSS) depends upon the ground covering materials, albeit the land use of the area. A mathematical model is set up to express the relation between the reflected electromagnetic energy of sample areas and their population distribution. Landsat 1 and Landsat 3 data of the Kanto area (including Tokyo Metropolitan), acquired in 1972 and 1979, are used along with ground-based census data of 1970 and 1975 to monitor the population distribution and its temporal changes. The method provided a reliable assessment of the population density in residential zones, however land-use classification using MSS imagery previous to the modeling is expected to improve the results.     

Mapping forest vegetation using Landsat TM imagery and a canopy reflectance model

Estimates of mean tree size and cover for each forest stand from an invertible forest canopy reflectance model are part of a new forest vegetation mapping system. Image segmentation defines stands which are sorted into general growth forms using per-pixel image classifications. Ecological models based on terrain relations predict species associations for the conifer, hardwood, and brush growth forms. The combination of the model-based estimates of tree size and cover with species associations yields general-purpose vegetation maps useful for a variety of land management needs. Results of timber inventories in the Tahoe and Stanislaus National Forests indicate the vegetation maps form a useful basis for stratification. Patterns in timber volumes for the strata reveal that the cover estimates are more reliable than the tree size estimates. A map accuracy assessment of the Stanislaus National Forest shows high overall map accuracy and also illustrates the problems in estimating tree size.     

Predicting categorical forest variables using an improved k-Nearest Neighbour estimator and Landsat imagery

The k-Nearest Neighbour (k-NN) estimation and prediction technique is widely used to produce pixel-level predictions and areal estimates of continuous forest variables such as area and volume, often by sub-categories such as species. An advantage of k-NN is that the same parameters (e.g., k-value, distance metric, weight vector for the feature space variables) can be used for all variables, whether continuous or categorical. An obvious question is the degree to which accuracy can be improved if the k-NN estimation parameters are tailored for specific variable groups such as volumes by tree species or categorical variables. We investigated prediction of categorical forest attribute variables from satellite image spectral data using k-NN with optimisation of the weight vector for the ancillary variables obtained using a genetic algorithm. We tested several genetic algorithm fitness functions, all derived from well-known accuracy measures. For a Finnish test site, the categorical forest attribute variables were site fertility and tree species dominance, and for an Italian test site, the variables were forest type and conifer/broad-leaved dominance. The results for both test sites were validated using independent data sets. Our results indicate that use of the genetic algorithm to optimize the weight vector for prediction of a single forest attribute variable had a slight positive effect on the prediction accuracies for other variables. Errors can be further decreased if the optimisation is done by variable groups.     

Assessment of tropical forest degradation by selective logging and fire using Landsat imagery

Many studies have assessed the process of forest degradation in the Brazilian Amazon using remote sensing approaches to estimate the extent and impact by selective logging and forest fires on tropical rain forest. However, only a few have estimated the combined impacts of those anthropogenic activities. We conducted a detailed analysis of selective logging and forest fire impacts on natural forests in the southern Brazilian Amazon state of Mato Grosso, one of the key logging centers in the country. To achieve this goal a 13-year series of annual Landsat images (1992-2004) was used to test different remote sensing techniques for measuring the extent of selective logging and forest fires, and to estimate their impact and interaction with other land use types occurring in the study region. Forest canopy regeneration following these disturbances was also assessed. Field measurements and visual observations were conducted to validate remote sensing techniques. Our results indicated that the Modified Soil Adjusted Vegetation Index aerosol free (MSAVI_af) is a reliable estimator of fractional coverage under both clear sky and under smoky conditions in this study region. During the period of analysis, selective logging was responsible for disturbing the largest proportion (31%) of natural forest in the study area, immediately followed by deforestation (29%). Altogether, forest disturbances by selective logging and forest fires affected approximately 40% of the study site area. Once disturbed by selective logging activities, forests became more susceptible to fire in the study site. However, our results showed that fires may also occur in undisturbed forests. This indicates that there are further factors that may increase forest fire susceptibility in the study area. Those factors need to be better understood. Although selective logging affected the largest amount of natural forest in the study period, 35% and 28% of the observed losses of forest canopy cover were due to forest fire and selective logging combined and to forest fire only, respectively. Moreover, forest areas degraded by selective logging and forest fire is an addition to outright deforestation estimates and has yet to be accounted for by land use and land cover change assessments in tropical regions. Assuming that this observed trend of land use and land cover conversion continues, we predict that there will be no undisturbed forests remaining by 2011 in this study site. Finally, we estimated that 70% of the total forest area disturbed by logging and fire had sufficiently recovered to become undetectable using satellite data in 2004.     

Comparison of automatic thresholding methods for snow-cover mapping using Landsat TM imagery

Knowledge of snow cover is essential to understanding the global water and energy cycle. Thresholding the normalized difference snow index (NDSI) image is a method frequently used to map snow cover from remotely sensed data. However, the threshold is dependent on the scenario and needs to be determined accordingly. In this study, nine automatic thresholding methods were tested on the NDSI. Comparisons of the automatic thresholding methods, optimal threshold, and support vector machine (SVM) classification show that Otsu's and Nie's methods appear to be the most robust among the nine automatic thresholding methods, achieving comparable accuracies with the latter two approaches. In addition, NDSI from the digital number (DN) can be an efficient substitution for NDSI obtained from atmospherically or topographically corrected data, with similar accuracy.     

Detection of forest harvest type using multiple dates of Landsat TM imagery

A simple and relatively accurate technique for classifying time-series Landsat Thematic Mapper (TM) imagery to detect levels of forest harvest is the topic of this research. The accuracy of multidate classification of the normalized difference vegetation index (NDVI) and the normalized difference moisture index (NDMI) were compared and the effect of the number of years (1–3, 3–4, 5–6 years) between image acquisition on forest change accuracy was evaluated. When Landsat image acquisitions were only 1–3 years apart, forest clearcuts were detected with producer's accuracy ranging from 79% to 96% using the RGB-NDMI classification method. Partial harvests were detected with lower producer's accuracy (55–80%) accuracy. The accuracy of both clearcut and partial harvests decreased as time between image acquisition increased. In all classification trials, the RGB-NDMI method produced significantly higher accuracies, compared to the RGB-NDVI. These results are interesting because the less common NDMI (using the reflected middle infrared band) outperformed the more popular NDVI. In northern Maine, industrial forest landowners have shifted from clearcutting to partial harvest systems in recent years. The RGB-NDMI change detection classification applied to Landsat TM imagery collected every 2–3 years appears to be a promising technique for monitoring forest harvesting and other disturbances that do not remove the entire overstory canopy.     

Mapping the land-cover distribution in arid and semiarid urban landscapes with Landsat Thematic Mapper imagery

Mapping the land-cover distribution in arid and semiarid urban landscapes using medium spatial resolution imagery is especially difficult due to the mixed-pixel problem in remotely sensed data and the confusion of spectral signatures among bare soils, sparse density shrub lands, and impervious surface areas (ISAs hereafter). This article explores a hybrid method consisting of linear spectral mixture analysis (LSMA), decision tree classifier, and cluster analysis for mapping land-cover distribution in two arid and semiarid urban landscapes: Urumqi, China, and Phoenix, USA. The Landsat Thematic Mapper (TM) imagery was unmixed into four endmember fraction images (i.e. high-albedo object, low-albedo object, green vegetation (GV), and soil) using the LSMA approach. New variables from these fraction images and TM spectral bands were used to map seven land-cover classes (i.e. forest, shrub, grass, crop, bare soil, ISA, and water) using the decision tree classifier. The cluster analysis was further used to modify the classification results. QuickBird imagery in Urumqi and aerial photographs in Phoenix were used to assess classification accuracy. Overall classification accuracies of 86.0% for Urumqi and 88.7% for Phoenix were obtained, much higher accuracies than those utilizing the traditional maximum likelihood classifier (MLC). This research demonstrates the necessity of using new variables from fraction images to distinguish between ISA and bare soils and between shrub and other vegetation types. It also indicates the different effects of spatial patterns of land-cover composition in arid and semiarid landscapes on urban land-cover classification.     

A weighted normalized difference water index for water extraction using Landsat imagery

Water extraction is one of challenging topics in studies on remote-sensing applications. Spectral profiles and experiments indicate that existing water indices often misclassified turbid water, small waterbodies, and some land features in a shadow area. In this study, a new water index called weighted normalized difference water index (WNDWI) was proposed to reduce those errors and improve the mapping accuracy of waterbodies by using Landsat imagery. To test the performance of the newly proposed water index, two test sites (Tampa Bay, FL, USA and Xiangshan Harbour, Zhejiang, China) were selected and the performances of three existing water indices including the normalized difference water index (NDWI), the modified NDWI (MNDWI), and the automated water extraction index (AWEI) were compared with that of the WNDWI. In addition, a default threshold 0 and automatically thresholding methods including Otsu threshold method and multiple thresholds identified by valley points in a histogram curve were tested to determine an optimal threshold that can be used to separate water and non-water features from grey images created by the four water indices. The experimental results indicate that the overall accuracies (OAs) created with WNDWI were all higher than those created with the three existing water indices: NDWI, MNDWI, and AWEI in both sites. Moreover, the results thresholded by 0 owned or shared the highest OAs with the results segmented by some of non-zero thresholds obtained from Otsu method and multiple thresholds method. Therefore, using an appropriate threshold, the proposed method could extract waterbodies from Landsat TM imagery with a high accuracy.     

Evaluation of photographic enhancements of Landsat imagery

In this study, color Landsat images were photographically enhanced and then evaluated for information content. Twenty-six color images were enhanced by changing gamma values on the images of bands 4 through 7 prior to generation of the color composites. This results in color-composited images with various color balances. These images were evaluated and compared to standard color-composite images. The evaluation of the photographic enhancements was conducted as a geologic and land-cover photo-interpretative analysis of scene no. 5470–19560 acquired on 1 August 1976 of the Denali Highway area, Alaska. The commonly used gamma values of 1.0 for all bands utilized during image reproduction does not result in the most interpretable image. A more informative image can be derived by using gamma values of 2.0 for band 4, 1.0 for band 5, and 2.0 for band 6 or 7. The remaining photographic enhancements that were attempted are generally found to enhance certain geologic or land-cover features at the expense of others.     

Classification of Amazonian primary rain forest vegetation using Landsat ETM+ satellite imagery

A lack of spatially and thematically accurate vegetation maps complicates conservation and management planning, as well as ecological research, in tropical rain forests. Remote sensing has considerable potential to provide such maps, but classification accuracy within primary rain forests has generally been inadequate for practical applications. Here we test how accurately floristically defined forest types in lowland tropical rain forests in Peruvian Amazonia can be recognized using remote sensing data (Landsat ETM+ satellite image and STRM elevation model). Floristic data and a vegetation classification with four forest classes were available for eight line transects, each 8 km long, located in an area of ca 800 km². We compared two sampling unit sizes (line transect subunits of 200 and 500 m) and several image feature combinations to analyze their suitability for image classification. Mantel tests were used to quantify how well the patterns in elevation and in the digital numbers of the satellite image correlated with the floristic patterns observed in the field. Most Mantel correlations were positive and highly significant. Linear discriminant analysis was used first to build a function that discriminates between forest classes in the eight field-verified transects on the basis of remotely sensed data, and then to classify those parts of the line transects and the satellite image that had not been visited in the field. Classification accuracy was quantified by 8-fold crossvalidation. Two of the tierra firme (non-inundated) forest types were combined because they were too often misclassified. The remaining three forest types (inundated forest, terrace forest and Pebas formation/intermediate tierra firme forest) could be separated using the 500-m sampling units with an overall classification accuracy of 85% and a Kappa coefficient of 0.62. For the 200-m sampling units, the classification accuracy was clearly lower (71%, Kappa 0.35). The forest classification will be used as habitat data to study wildlife habitat use in the same area. Our results show that remotely sensed data and relatively simple classification methods can be used to produce reasonably accurate forest type classifications, even in structurally homogeneous primary rain forests.     

Object-based satellite image analysis applied for crop residue estimating using Landsat OLI imagery

ABSTRACT

The main objective of this study is to apply an object-based image analysis (OBIA) approach to satellite image processing and determining crop residue cover (CRC) and tillage intensity. To achieve this goal, we collected ground truth data using line-transect method from 35 plots of farmlands with an area of 528 ha. Accordingly, Landsat Operational Land Imager (OLI) satellite image together with global positioning system (GPS)-based survey data set were considered for applying the OBIA methods and deriving CRC. To process the data, object-based image processing steps including segmentation and classification were applied to develop intelligent objects and establish classification using spectral and spatial characteristics of CRC. We developed three categories of rule sets including mean indices, tillage indices, and grey-level co-occurrence matrix (GLCM) texture features using the OBIA algorithms and assign class method. Results were validated against of ground control data set and were collected by GPS in field survey. Results of this study indicated that the brightness, normalised difference tillage index, and GLCM texture feature mean performed out as effective techniques. Overall accuracy and kappa coefficient (κ) were computed to be about 0.91 and 0.86; 0.93 and 0.90; 0.60 and 0.35, respectively, for the above-mentioned indices. The foregoing discussion has attempted to demonstrate that the remotely sensed data can be effective approach and substitute for ground methods, especially in large areas.     

Object-based cloud and cloud shadow detection in Landsat imagery

A new method called Fmask (Function of mask) for cloud and cloud shadow detection in Landsat imagery is provided. Landsat Top of Atmosphere (TOA) reflectance and Brightness Temperature (BT) are used as inputs. Fmask first uses rules based on cloud physical properties to separate Potential Cloud Pixels (PCPs) and clear-sky pixels. Next, a normalized temperature probability, spectral variability probability, and brightness probability are combined to produce a probability mask for clouds over land and water separately. Then, the PCPs and the cloud probability mask are used together to derive the potential cloud layer. The darkening effect of the cloud shadows in the Near Infrared (NIR) Band is used to generate a potential shadow layer by applying the flood-fill transformation. Subsequently, 3D cloud objects are determined via segmentation of the potential cloud layer and assumption of a constant temperature lapse rate within each cloud object. The view angle of the satellite sensor and the illuminating angle are used to predict possible cloud shadow locations and select the one that has the maximum similarity with the potential cloud shadow mask. If the scene has snow, a snow mask is also produced. For a globally distributed set of reference data, the average Fmask overall cloud accuracy is as high as 96.4%. The goal is development of a cloud and cloud shadow detection algorithm suitable for routine usage with Landsat images.     

An automated approach for reconstructing recent forest disturbance history using dense Landsat time series stacks

A highly automated algorithm called vegetation change tracker (VCT) has been developed for reconstructing recent forest disturbance history using Landsat time series stacks (LTSS). This algorithm is based on the spectral-temporal properties of land cover and forest change processes, and requires little or no fine tuning for most forests with closed or near close canopy cover. It was found very efficient, taking 2-3 h on average to analyze an LTSS consisting of 12 or more Landsat images using an average desktop PC. This LTSS-VCT approach has been used to examine disturbance patterns with a biennial temporal interval from 1984 to 2006 for many locations across the conterminous U.S. Accuracy assessment over 6 validation sites revealed that overall accuracies of around 80% were achieved for disturbances mapped at individual year level. Average user's and producer's accuracies of the disturbance classes were around 70% and 60% in 5 of the 6 sites, respectively, suggesting that although forest disturbances were typically rare as compared with no-change classes, on average the VCT detected more than half of those disturbances with relatively low levels of false alarms. Field assessment revealed that VCT was able to detect most stand clearing disturbance events, including harvest, fire, and urban development, while some non-stand clearing events such as thinning and selective logging were also mapped in western U.S. The applicability of the LTSS-VCT approach depends on the availability of a temporally adequate supply of Landsat imagery. To ensure that forest disturbance records can be developed continuously in the future, it is necessary to plan and develop observational capabilities today that will allow continuous acquisition of frequent Landsat or Landsat-like observations.     

A preliminary assessment of Landsat TM imagery for mapping vegetation and sediment distribution in the Wash estuary

Abstract

Analysis of Landsat Thematic Mapper (TM) image of 14 May 1984 has shown that such data can be used to survey vegetation and sediment distributions in the intertidal zone of the Wash estuary at a spatial detail comparable with current methods practised by the.Nature Conservancy Council. Multispectral classification of this TM image showed good separation of salt-marsh vegetation communities which had recently been surveyed by the Nature Conservancy Council and for which reliable training data could be taken. The sensitivity of classification performance, using both parametric and non-parametric algorithms, to apparently minor differences in phenology at training site locations demonstrates two requirements for improved salt-marsh classification. They are the need for strict definition of training data and that TM wave bands 2, 3, 4 and 5 provide suitable spectral vectors for classifying intertidal environments.     

Assessment of deforestation in the Lower Amazon floodplain using historical Landsat MSS/TM imagery

The floodplain forests bordering the Amazon River have outstanding ecological, economic, and social importance for the region. However, the original distribution of these forests is not well known, since they have suffered severe degradation since the 16th century. The previously published vegetation map of the Amazon River floodplain (Hess et al., 2003), based on data acquired in 1996, shows enormous difference in vegetation cover classes between the regions upstream and downstream of the city of Manaus. The upper floodplain is mostly covered by forests, while the lower floodplain is predominantly occupied by grasses and shrubs.This study assesses deforestation in the Lower Amazon floodplain over a ~ 30 year period by producing and comparing a historical vegetation map based on MSS/Landsat images acquired in the late 1970s with a recent vegetation map produced from TM/Landsat images obtained in 2008. The maps were generated through the following steps: 1) normalization and mosaicking of images for each decade; 2) application of a linear mixing model transformation to produce vegetation, soil and shade fraction-images; and 3) object-oriented image analysis and classification. For both maps, the following classes were mapped: floodplain forest, non-forest floodplain vegetation, bare soil and open water. The two maps were combined using object-level Boolean operations to identify time transitions among the mapped classes, resulting in a map of the land cover change occurred over ~ 30 years. Ground information collected at 168 ground points was used to build confusion matrices and calculate Kappa indices of agreement. A survey strategy combining field observations and interviews allowed the collection of information about both recent and historical land cover for validation purposes. Kappa values (0.77, 0.75 and 0.75) indicated the good quality of the maps, and the error estimates were used to adjust the estimated deforested area to a value of 3457 km² ± 1062 km² (95% CI) of floodplain deforestation over the ~ 30 years.