Landsat-8 Operational Land Imager Data Research Articles

Subpixel Inundation Mapping Using Landsat-8 OLI and UAV Data for a Wetland Region on the Zoige Plateau, China

Wetland inundation is crucial to the survival and prosperity of fauna and flora communities in wetland ecosystems. Even small changes in surface inundation may result in a substantial impact on the wetland ecosystem characteristics and function. This study presented a novel method for wetland inundation mapping at a subpixel scale in a typical wetland region on the Zoige Plateau, northeast Tibetan Plateau, China, by combining use of an unmanned aerial vehicle (UAV) and Landsat-8 Operational Land Imager (OLI) data. A reference subpixel inundation percentage (SIP) map at a Landsat-8 OLI 30 m pixel scale was first generated using high resolution UAV data (0.16 m). The reference SIP map and Landsat-8 OLI imagery were then used to develop SIP estimation models using three different retrieval methods (Linear spectral unmixing (LSU), Artificial neural networks (ANN), and Regression tree (RT)). Based on observations from 2014, the estimation results indicated that the estimation model developed with RT method could provide the best fitting results for the mapping wetland SIP (R2 = 0.933, RMSE = 8.73%) compared to the other two methods. The proposed model with RT method was validated with observations from 2013, and the estimated SIP was highly correlated with the reference SIP, with an R2 of 0.986 and an RMSE of 9.84%. This study highlighted the value of high resolution UAV data and globally and freely available Landsat data in combination with the developed approach for monitoring finely gradual inundation change patterns in wetland ecosystems.

Mapping and Characterization of Agricultural Systems from Time Series of Normalized Difference Vegetation Index (NDVI) in the Northeast Area of Tadla, Morocco

Introduction: Agriculture is the focus of sustainable development and food security in Morocco, contributing 14% to the national economy and its socio-economic impact depends on agro-climatic conditions and the availability of water resources. The objective of this work is to demonstrate that agriculture in Morocco, and especially in the Tadla region, is distributed according to the availability of underground water. Method: This work was carried out in the transition zone between the phosphate plateau and the plain of Tadla, a site of more than 124,000 ha. To achieve the objectives set, we used Landsat-8 OLI (Operational Land Imager) data sharpened to 15 m. The Support Vector Machine (SVM) classifications were applied to the 10-period Normalized Difference Vegetation Index (NDVI) series. Results: It is shown that agriculture in Morocco is highly dependent on irrigation. In the classified map, it can be observed that the southern region, a zone with a semi-arid climate, which generally relies on irrigation, is very developed in agriculture, whereas in the north most of the soil is bared with cereals that mainly depend on rainfall. Conclusion: The possible water shortage in the region puts food availability in Morocco at risk of inadequacy. Food security depends on the availability of water, a natural resource at risk of extinction. Therefore, strategies to reduce the risk of agricultural drought must be at the forefront of Moroccan policies.

Evaluating and comparing performances of topographic correction methods based on multi-source DEMs and Landsat-8 OLI data

ABSTRACTTopographic correction is a crucial and challenging step in interpreting optical remote-sensing images of extremely complex terrain environments due to the lack of universally suitable correction algorithms and digital elevation models (DEMs) of adequate resolution and quality. The free availability of open source global DEMs provides an unprecedented opportunity to remove topographic effects associated with remote-sensing data in remote and rugged mountain terrains. This study evaluated the performances of seven topographic correction methods including C-correction (C), Cosine C-correction (CC), Minnaert correction (M), Sun–canopy–sensor (SCS) correction (S), SCS+C-correction (SC), Teillet regression correction (TR), and the Terrain illumination correction model (TI) based on multi-source DEMs (local topographic map, Shuttle Radar Topography Mission (SRTM) DEM filled-finished A/B and the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) global digital elevation model (GDEM) data sets) and Landsat-8 Operational Land Imager (OLI) data using visual and statistical evaluation strategies. Overall, these investigated topographic correction methods removed topographic effects associated with Landsat-8 OLI data to varying degrees. However, the performances of these methods generally depend on the use of different DEMs and evaluation strategies. Among these correction methods, the SCS+C-correction performed best and was less sensitive to the use of different DEMs. The performances of topographic corrections based on free and open-access DEMs were generally better than or comparable to those based on local topographic maps. In particular, the topographic correction performance was greatly improved using the SRTM filled-finished B (FFB) data set with a resampling scheme based on the average value within a 3 × 3 pixel window. Nevertheless, further quantitative investigation is needed to determine the relative importance of DEMs and evaluation strategies used to select topographic correction methods.

An index and approach for water extraction using Landsat–OLI data

ABSTRACTThe extraction of water distribution is extremely useful in research and planning activities, including those associated with water resources, environments, disasters, local climates, and other factors. Remote-sensing images with moderate resolution have been the main data source due to the vast distribution of water and the high cost, access difficulty, and massive size of high-resolution images. Although some water indices and methods for water extraction have been proposed, there is still a lack of these resources to easily, accurately, efficiently, and automatically extract water. This paper focused on some improvements that mainly used the most traditional but also the newest Operational Land Imager (OLI) images in Landsat 8. This study first analysed the variation features of previous water indices. Secondly, taking the city of Beijing and its surrounding area as the experimental site, a spectral curve analysis was performed and a new water index was proposed. This index was compared to three typical indices. Thirdly, a new approach was proposed to accurately and easily extract water. It included four major steps: background partitioning, thresholding and preliminary segmentation, noise removal by patch size, and local region growth. Next, the stricter and more effective stratified random sampling method was used to test the accuracy. Then, we tested the generality of the proposed water index and extraction method using nine typical test sites from around the world and tried to simplify the workflow. Finally, this paper discusses threshold optimization issues, such as automatic selection and reduction of the number of thresholds. The results show that the normalized water index (NDWI), modified normalized water index (MNDWI), and normalized difference built-up index (NDBI) may fail in some situations due to the complex spectrum of the impervious surface class. Some shadow pixels were impossible to remove using only spectral analysis because both the digital number (DN) trends and values were similar to those of water. The proposed water index was easy and simple, but it corresponded better to water bodies. Additionally, it was more accurate and universal and showed greater potential for extracting water. This method relatively accurately and completely extracted various water bodies from plain city, plain country, and natural mountainous regions in many typical climate zones, eliminating interference caused by dark impervious surfaces, plants, sand, suspended sediments, snow, ice, bedrock, reservoir drawdown areas, shadows from mountains and buildings, mixed pixels, etc. The mean kappa coefficients were 0.988, 0.982, and 0.984 in plain city, plain country, and natural mountainous regions, respectively. This paper suggests that thresholds can be automatically determined by comparing the accuracy changes of different thresholds according to preselected sample and test points. Furthermore, the combined use of the maximum class square error method (also known as the Ostu algorithm) and the adaptive thresholding method exhibits great potential for automatic determination of thresholds in regions without many noises with higher water index values. In addition, water bodies could also be accurately extracted by setting these thresholds to fixed values based on the results at more test sites.

Timber production assessment of a plantation forest: An integrated framework with field-based inventory, multi-source remote sensing data and forest management history

Timber production is the purpose for managing plantation forests, and its spatial and quantitative information is critical for advising management strategies. Previous studies have focused on growing stock volume (GSV), which represents the current potential of timber production, yet few studies have investigated historical process-harvested timber. This resulted in a gap in a synthetical ecosystem service assessment of timber production. In this paper, we established a Management Process–based Timber production (MPT) framework to integrate the current GSV and the harvested timber derived from historical logging regimes, trying to synthetically assess timber production for a historical period. In the MPT framework, age-class and current GSV determine the times of historical thinning and the corresponding harvested timber, by using a “space-for-time” substitution. The total timber production can be estimated by the historical harvested timber in each thinning and the current GSV. To test this MPT framework, an empirical study on a larch plantation (LP) with area of 43,946ha was conducted in North China for a period from 1962 to 2010. Field-based inventory data was integrated with ALOS PALSAR (Advanced Land-Observing Satellite Phased Array L-band Synthetic Aperture Radar) and Landsat-8 OLI (Operational Land Imager) data for estimating the age-class and current GSV of LP. The random forest model with PALSAR backscatter intensity channels and OLI bands as input predictive variables yielded an accuracy of 67.9% with a Kappa coefficient of 0.59 for age-class classification. The regression model using PALSAR data produced a root mean square error (RMSE) of 36.5m3ha−1. The total timber production of LP was estimated to be 7.27×106m3, with 4.87×106m3 in current GSV and 2.40×106m3 in harvested timber through historical thinning. The historical process-harvested timber accounts to 33.0% of the total timber production, which component has been neglected in the assessments for current status of plantation forests. Synthetically considering the RMSE for predictive GSV and misclassification of age-class, the error in timber production were supposed to range from −55.2 to 56.3m3ha−1. The MPT framework can be used to assess timber production of other tree species at a larger spatial scale, providing crucial information for a better understanding of forest ecosystem service.

Radiometric cross Calibration of Gaofen-1 WFV Cameras Using Landsat-8 OLI Images: A Simple Image-Based Method

WFV (Wide Field of View) cameras on-board Gaofen-1 satellite (gaofen means high resolution) provide unparalleled global observations with both high spatial and high temporal resolutions. However, the accuracy of the radiometric calibration remains unknown. Using an improved cross calibration method, the WFV cameras were re-calibrated with well-calibrated Landsat-8 OLI (Operational Land Imager) data as reference. An objective method was proposed to guarantee the homogeneity and sufficient dynamic coverage for calibration sites and reference sensors. The USGS spectral library was used to match the most appropriate hyperspectral data, based on which the spectral band differences between WFV and OLI were adjusted. The TOA (top-of-atmosphere) reflectance of the cross-calibrated WFV agreed very well with that of OLI, with the mean differences between the two sensors less than 5% for most of the reflectance ranges of the four spectral bands, after accounting for the spectral band difference between the two sensors. Given the calibration error of 3% for Landsat-8 OLI TOA reflectance, the uncertainty of the newly-calibrated WFV should be within 8%. The newly generated calibration coefficients established confidence when using Gaofen-1 WFV observations for their further quantitative applications, and the proposed simple cross calibration method here could be easily extended to other operational or planned satellite missions.

Monitoring Plastic-Mulched Farmland by Landsat-8 OLI Imagery Using Spectral and Textural Features

In recent decades, plastic-mulched farmland has expanded rapidly in China as well as in the rest of the world because it results in marked increases of crop production. However, plastic-mulched farmland significantly influences the environment and has so far been inadequately investigated. Accurately monitoring and mapping plastic-mulched farmland is crucial for agricultural production, environmental protection, resource management, and so on. Monitoring plastic-mulched farmland using moderate-resolution remote sensing data is technically challenging because of spatial mixing and spectral confusion with other ground objects. This paper proposed a new scheme that combines spectral and textural features for monitoring the plastic-mulched farmland and evaluates the performance of a Support Vector Machine (SVM) classifier with different kernel functions using Landsat-8 Operational Land Imager (OLI) imagery. The textural features were extracted from multi-bands OLI data using a Grey Level Co-occurrence Matrix (GLCM) algorithm. Then, six combined feature sets were developed for classification. The results indicated that Landsat-8 OLI data are well suitable for monitoring plastic-mulched farmland; the SVM classifier with a linear kernel function is superior both to other kernel functions and to two other widely used supervised classifiers: Maximum Likelihood Classifier (MLC) and Minimum Distance Classifier (MDC). For the SVM classifier with a linear kernel function, the highest overall accuracy was derived from combined spectral and textural features in the 90° direction (94.14%, kappa 0.92), followed by the combined spectral and textural features in the 45° (93.84%, kappa 0.92), 135° (93.73%, kappa 0.92), 0° (93.71%, kappa 0.92) directions, and the spectral features alone (93.57%, kappa 0.91). Spectral features make a more significant contribution to monitoring the plastic-mulched farmland; adding textural features from medium resolution imagery provide only limited improvement in accuracy.

Estimating wetland vegetation abundance from Landsat-8 operational land imager imagery: a comparison between linear spectral mixture analysis and multinomial logit modeling methods

Mapping vegetation abundance by using remote sensing data is an efficient means for detecting changes of an eco-environment. With Landsat-8 operational land imager (OLI) imagery acquired on July 31, 2013, both linear spectral mixture analysis (LSMA) and multinomial logit model (MNLM) methods were applied to estimate and assess the vegetation abundance in the Wild Duck Lake Wetland in Beijing, China. To improve mapping vegetation abundance and increase the number of endmembers in spectral mixture analysis, normalized difference vegetation index was extracted from OLI imagery along with the seven reflective bands of OLI data for estimating the vegetation abundance. Five endmembers were selected, which include terrestrial plants, aquatic plants, bare soil, high albedo, and low albedo. The vegetation abundance mapping results from Landsat OLI data were finally evaluated by utilizing a WorldView-2 multispectral imagery. Similar spatial patterns of vegetation abundance produced by both fully constrained LSMA algorithm and MNLM methods were observed: higher vegetation abundance levels were distributed in agricultural and riparian areas while lower levels in urban/built-up areas. The experimental results also indicate that the MNLM model outperformed the LSMA algorithm with smaller root mean square error (0.0152 versus 0.0252) and higher coefficient of determination (0.7856 versus 0.7214) as the MNLM model could handle the nonlinear reflection phenomenon better than the LSMA with mixed pixels.

An experimental comparison between KELM and CART for crop classification using Landsat-8 OLI data

The operational land imager (OLI) is the latest instrument in the Landsat series of satellite imagery, which officially began normal operations on 30 May 2013. The OLI includes two bands that are not on the thematic mapper series of sensors aboard Landsat-5 and 7; a cirrus band and a coastal/aerosol band. This paper compares the classification and regression tree and the kernel-based extreme learning machine (KELM) for mapping crops in Hokkaido, Japan, using OLI data, except the cirrus band and the pan band. The OLI data acquired on 8 July 2013 was used for crop classification of beans, beets, grassland, maize, potatoes and winter wheat. The KELM algorithm performed better in this study and achieved overall accuracies of 90.1%. According to the Jeffries–Matusita (J–M) distances, the short wavelength infrared band provides the greater contribution (the highest value was observed for band 6 in OLI data).

Extent and Area of Swidden in Montane Mainland Southeast Asia: Estimation by Multi-Step Thresholds with Landsat-8 OLI Data

Information on the distribution, area and extent of swidden agriculture landscape is necessary for implementing the program of Reducing Emissions from Deforestation and Forest Degradation (REDD), biodiversity conservation and local livelihood improvement. To our knowledge, explicit spatial maps and accurate area data on swidden agriculture remain surprisingly lacking. However, this traditional farming practice has been transforming into other profit-driven land use, like tree plantations and permanent cash agriculture. Swidden agriculture is characterized by a rotational and dynamic nature of agroforestry, with land cover changing from natural forests, newly-cleared swiddens to different-aged fallows. The Operational Land Imager (OLI) onboard the Landsat-8 satellite has visible, near-infrared and shortwave infrared bands, which are sensitive to the changes in vegetation cover, land surface moisture content and soil exposure, and therefore, four vegetation indices (VIs) were calculated, including the Normalized Difference Vegetation Index (NDVI), the Normalized Difference Moisture Index (NDMI), the Normalized Burn Ratio (NBR) and the Soil Adjusted Vegetation Index (SAVI). In this study, we developed a multi-step threshold approach that uses a combination of thresholds of four VIs and local elevation range (LER) and applied it to detect and map newly-opened swiddens and different-aged fallows using OLI imagery acquired between 2013 and 2015. The resultant Landsat-derived swidden agriculture maps have high accuracy with an overall accuracy of 86.9% and a Kappa coefficient of 0.864. The results of this study indicated that the Landsat-based multi-step threshold algorithms could potentially be applied to monitor the long-term change pattern of swidden agriculture in montane mainland Southeast Asia since the late 1980s and also in other tropical regions, like insular Southeast Asia, South Asia, Latin America and Central Africa, where swidden agriculture is still common.

An Effective Method for Snow-Cover Mapping of Dense Coniferous Forests in the Upper Heihe River Basin Using Landsat Operational Land Imager Data

The Normalized Difference Snow Index (NDSI) is an effective index for snow-cover mapping at large scales, but in forested regions the identification accuracy for snow using the NDSI is low because of forest cover effects. In this study, typical evergreen coniferous forest zones on Qilian Mountain in the Upper Heihe River Basin (UHRB) were chosen as example regions. By analyzing the spectral signature of snow-covered and snow-free evergreen coniferous forests with Landsat Operational Land Imager (OLI) data, a novel spectral band ratio using near-infrared (NIR) and shortwave infrared (SWIR) bands, defined as (ρnir − ρswir)/(ρnir + ρswir), is proposed. Our research shows that this band ratio, named the normalized difference forest snow index (NDFSI), can be used to effectively distinguish snow-covered evergreen coniferous forests from snow-free evergreen coniferous forests in UHRB.

Radiometric cross-calibration of Gaofen-1 WFV cameras using Landsat-8 OLI images: A solution for large view angle associated problems

Abstract Four wide-field-of-view (WFV) instruments are onboard the Gaofen-1 (or GF-1) satellite, providing a combined swath of ~ 800 km. However, observations with large view angles pose new challenges for radiometric cross-calibration with a simple image-based method when using Landsat-8 Operational Land Imager (OLI) data as a reference. Based on radiative transfer modeling, a novel radiometric cross-calibration method has been proposed in this study to solve large view angle-associated problems. The Moderate Resolution Imaging Spectroradiometer (MODIS) aerosol products were used to simulate the top of atmosphere (TOA) signal of the reference and target instruments, and the unequal bidirectional effects were corrected using MODIS bi-directional reflectance distribution function (BRDF) products. Extensive validations with both satellite data and in situ measurements revealed an uncertainty of ~ 8% for the newly produced cross-calibration coefficients when they were used to calibrate TOA reflectance for both close-nadir and off-nadir instruments. The improvements are discernable when compared with the official provided coefficients and that were derived using the image-based cross-calibration method. This study demonstrated not only the usefulness of Landsat-8 OLI data in sensor radiometric calibration but also the impressive accuracy of the MODIS BRDF and aerosol products in radiative transfer simulations. The proposed method can be used in the future to monitor and correct potential radiometric degradations of the GF-1 WFV instruments, and it also can be easily extended to other similar satellite missions to conduct radiometric cross-calibrations.

Removal of Thin Clouds in Landsat-8 OLI Data with Independent Component Analysis

An approach to remove clouds in Landsat-8 operational land imager (OLI) data was developed with independent component analysis (ICA). Within cloud-covered areas, histograms were derived to quantify changes of the reflectance values before and after the use of the algorithm. Referred to a cloud-free image, changes of histogram curves validated the algorithm. Scatterplots were generated and linear regression performed for the reflectance values of each band before and after the algorithm, and compared to those of the reference image. Band-by-band, results in cloud removal were acceptable. The algorithm had little effect on pixels in cloud-free areas after the analyses of histograms, scatterplots, and linear regression equations. Finally, the algorithm was applied to various land use and land cover types and cloud conditions, and to a full Landsat-8 scene yielding satisfactory results efficiently.

An investigation of spatial variation of suspended sediment concentration induced by a bay bridge based on Landsat TM and OLI data

This study investigates the spatial variation of suspended sediment concentration (SSC) induced by Hangzhou Bay Bridge (HBB) in the coastal waters of the East China Sea. Based on Landsat Thematic Mapper (TM), Operational Land Imager (OLI) data and in-situ measurements, spectral characteristics of waters in Hangzhou Bay and SSC retrieved from near infrared single band have been analyzed. The results revealed significant difference of SSC in the water on two sides of HBB. SSC increases downstream of the bridge under conditions of low turbidity (SSC<300mgl−1) upstream water, while SSC decreases when high turbidity water (SSC>400mgl−1) upstream. This study shows that the interaction of bridge piers and currents has important influences on SSC distribution by inducing hydrodynamic factors and by changing suspended sediment transport. Remotely sensed data with high spatial resolution as TM and OLI can be applied to investigate SSC variations induced by a bridge in the bay area.

A new index for mapping built-up and bare land areas from Landsat-8 OLI data

Remote sensing is a useful technology for monitoring the spatial distribution and expansion of built-up and bare land areas. One effective approach, known as the normalized difference built-up index (NDBI), has been promoted for identifying built-up areas based on Landsat Thematic Mapper (TM)/Enhanced Thematic Mapper (ETM+) data. The successful launch of the Landsat-8 satellite has made possible the continued acquisition of high-quality data that meet requirements for observing land use and land cover (LULC), and whether or not NDBI approach can be used with Landsat-8 data needs further verification. In this study, we researched the reflectance spectral characteristics of different land cover types in different bands of Landsat-8 operational land imager (OLI) data and found that the trend of some built-up areas in the OLI data from the near-infrared band to the shortwave infrared band is different from that in the TM/ETM+ data. This different trend made the conventional NDBI approach unsuitable for Landsat-8 OLI data. We propose a new index, called the build-up and bare land areas index, for transforming Landsat-8 OLI data to map built-up and bare land areas automatically. This new index was used to detect the built-up and bare land areas in Zhengzhou (Henan, China). The accuracy assessment indicates that our index has much higher accuracy (90.8%) than the conventional NDBI approach does (57.4%).