scholarly journals Harmonizing Surface Reflectance Between Landsat-7 ETM+, Landsat-8 OLI, and Sentinel-2 MSI

Author(s):  
Hongye Cao ◽  
Ling Han ◽  
Liangzhi Li

Abstract Remote sensing dynamic monitoring methods often benefit from a dense time series of observations. To enhance these time series, it is sometimes necessary to integrate data from multiple satellite systems. For more than 40 years, Landsat has provided the longest time record of space-based land surface observations, and the successful launch of the Landsat-8 Operational Land Imager (OLI) sensor in 2013 continues this tradition. However, the 16-day observation period of Landsat images has challenged the ability to measure subtle and transient changes like never before. The European Space Agency (ESA) launched the Sentinel-2A satellite in 2015. The satellite carries a Multispectral Instrument (MSI) sensor that provides a 10-20m spatial resolution data source providing an opportunity to complement the Landsat data record. The collection of Sentinel-2A MSI, Landsat-7 ETM+, and Landsat-8 OLI data provide multispectral global coverage from 10m to 30m with further reduced data revisit intervals. There are many differences between sensor data that need to be taken into account to use these data together reliably. The purpose of this study is to evaluate the potential of integrating surface reflectance data from Landsat-7, Landsat-8 and Sentinel-2 archived in the Google Earth Engine (GEE) cloud platform. To test and quantify the differences between these sensors, hundreds of thousands of surface reflectance data from sensor pairs were collected over China. In this study, some differences in the surface reflectance of the sensor pairs were identified, based upon which a cross-sensor conversion model was proposed, i.e., a suitable adjustment equation was fitted using an ordinary least squares (OLS) linear regression method to convert the Sentinel-2 reflectance values closer to the Landsat-7 or Landsat-8 values. The regression results show that the Sentinel MSI data are spectrally comparable to both types of Landsat image data, just as the Landsat sensors are comparable to each other. The root mean square error (RMSE) values between MSI and Landsat spectral values before coordinating the sensors ranged from 0.014 to 0.037, and the RMSE values between OLI and ETM + ranged from 0.019 to 0.039. After coordination, RMSE values between MSI and Landsat spectral values ranged from 0.011 to 0.026, and RMSD values between OLI and ETM + ranged from 0.013 to 0.034. The fitted adjustment equations were also compared to the HLS (Harmonized Landsat-8 Sentinel-2) global fitted equations (Sentinel-2 to Landsat-8) published by the National Aeronautics and Space Administration (NASA) and were found to be significantly different, increasing the likelihood that such adjustments would need to be fitted on a regional basis. This study believes that despite the differences in these datasets, it appears feasible to integrate these datasets by applying a linear regression correction between the bands.

2020 ◽  
Vol 9 (4) ◽  
pp. 257 ◽  
Author(s):  
Kiwon Lee ◽  
Kwangseob Kim ◽  
Sun-Gu Lee ◽  
Yongseung Kim

Surface reflectance data obtained by the absolute atmospheric correction of satellite images are useful for land use applications. For Landsat and Sentinel-2 images, many radiometric processing methods exist, and the images are supported by most types of commercial and open-source software. However, multispectral KOMPSAT-3A images with a resolution of 2.2 m are currently lacking tools or open-source resources for obtaining top-of-canopy (TOC) reflectance data. In this study, an atmospheric correction module for KOMPSAT-3A images was newly implemented into the optical calibration algorithm in the Orfeo Toolbox (OTB), with a sensor model and spectral response data for KOMPSAT-3A. Using this module, named OTB extension for KOMPSAT-3A, experiments on the normalized difference vegetation index (NDVI) were conducted based on TOC reflectance data with or without aerosol properties from AERONET. The NDVI results for these atmospherically corrected data were compared with those from the dark object subtraction (DOS) scheme, a relative atmospheric correction method. The NDVI results obtained using TOC reflectance with or without the AERONET data were considerably different from the results obtained from the DOS scheme and the Landsat-8 surface reflectance of the Google Earth Engine (GEE). It was found that the utilization of the aerosol parameter of the AERONET data affects the NDVI results for KOMPSAT-3A images. The TOC reflectance of high-resolution satellite imagery ensures further precise analysis and the detailed interpretation of urban forestry or complex vegetation features.


2020 ◽  
Vol 12 (19) ◽  
pp. 3121
Author(s):  
Roya Mourad ◽  
Hadi Jaafar ◽  
Martha Anderson ◽  
Feng Gao

Leaf area index (LAI) is an essential indicator of crop development and growth. For many agricultural applications, satellite-based LAI estimates at the farm-level often require near-daily imagery at medium to high spatial resolution. The combination of data from different ongoing satellite missions, Sentinel 2 (ESA) and Landsat 8 (NASA), provides this opportunity. In this study, we evaluated the leaf area index generated from three methods, namely, existing vegetation index (VI) relationships applied to Harmonized Landsat-8 and Sentinel-2 (HLS) surface reflectance produced by NASA, the SNAP biophysical model, and the THEIA L2A surface reflectance products from Sentinel-2. The intercomparison was conducted over the agricultural scheme in Bekaa (Lebanon) using a large set of in-field LAIs and other biophysical measurements collected in a wide variety of canopy structures during the 2018 and 2019 growing seasons. The major studied crops include herbs (e.g., cannabis: Cannabis sativa, mint: Mentha, and others), potato (Solanum tuberosum), and vegetables (e.g., bean: Phaseolus vulgaris, cabbage: Brassica oleracea, carrot: Daucus carota subsp. sativus, and others). Additionally, crop-specific height and above-ground biomass relationships with LAIs were investigated. Results show that of the empirical VI relationships tested, the EVI2-based HLS models statistically performed the best, specifically, the LAI models originally developed for wheat (RMSE:1.27), maize (RMSE:1.34), and row crops (RMSE:1.38). LAI derived through European Space Agency’s (ESA) Sentinel Application Platform (SNAP) biophysical processor underestimated LAI and provided less accurate estimates (RMSE of 1.72). Additionally, the S2 SeLI LAI algorithm (from SNAP biophysical processor) produced an acceptable accuracy level compared to HLS-EVI2 models (RMSE of 1.38) but with significant underestimation at high LAI values. Our findings show that the LAI-VI relationship, in general, is crop-specific with both linear and non-linear regression forms. Among the examined indices, EVI2 outperformed other vegetation indices when all crops were combined, and therefore it can be identified as an index that is best suited for a unified algorithm for crops in semi-arid irrigated regions with heterogeneous landscapes. Furthermore, our analysis shows that the observed height-LAI relationship is crop-specific and essentially linear with an R2 value of 0.82 for potato, 0.79 for wheat, and 0.50 for both cannabis and tobacco. The ability of the linear regression to estimate the fresh and dry above-ground biomass of potato from both observed height and LAI was reasonable, yielding R2: ~0.60.


2020 ◽  
Vol 12 (12) ◽  
pp. 2065 ◽  
Author(s):  
Feng Xu ◽  
Zhaofu Li ◽  
Shuyu Zhang ◽  
Naitao Huang ◽  
Zongyao Quan ◽  
...  

Winter wheat is one of the major cereal crops in China. The spatial distribution of winter wheat planting areas is closely related to food security; however, mapping winter wheat with time-series finer spatial resolution satellite images across large areas is challenging. This paper explores the potential of combining temporally aggregated Landsat-8 OLI and Sentinel-2 MSI data available via the Google Earth Engine (GEE) platform for mapping winter wheat in Shandong Province, China. First, six phenological median composites of Landsat-8 OLI and Sentinel-2 MSI reflectance measures were generated by a temporal aggregation technique according to the winter wheat phenological calendar, which covered seedling, tillering, over-wintering, reviving, jointing-heading and maturing phases, respectively. Then, Random Forest (RF) classifier was used to classify multi-temporal composites but also mono-temporal winter wheat development phases and mono-sensor data. The results showed that winter wheat could be classified with an overall accuracy of 93.4% and F1 measure (the harmonic mean of producer’s and user’s accuracy) of 0.97 with temporally aggregated Landsat-8 and Sentinel-2 data were combined. As our results also revealed, it was always good to classify multi-temporal images compared to mono-temporal imagery (the overall accuracy dropped from 93.4% to as low as 76.4%). It was also good to classify Landsat-8 OLI and Sentinel-2 MSI imagery combined instead of classifying them individually. The analysis showed among the mono-temporal winter wheat development phases that the maturing phase’s and reviving phase’s data were more important than the data for other mono-temporal winter wheat development phases. In sum, this study confirmed the importance of using temporally aggregated Landsat-8 OLI and Sentinel-2 MSI data combined and identified key winter wheat development phases for accurate winter wheat classification. These results can be useful to benefit on freely available optical satellite data (Landsat-8 OLI and Sentinel-2 MSI) and prioritize key winter wheat development phases for accurate mapping winter wheat planting areas across China and elsewhere.


2020 ◽  
Vol 12 (16) ◽  
pp. 2597
Author(s):  
Cibele Teixeira Pinto ◽  
Xin Jing ◽  
Larry Leigh

Landsat Level-1 products are delivered as quantized and calibrated scaled Digital Numbers (DN). The Level-1 DN data can be rescaled to Top-of-Atmosphere (TOA) reflectance applying radiometric rescaling coefficients. Currently, the Level-1 product is the standard data product of the Landsat sensors. The more recent Level-2 data products contain surface reflectance values, i.e., reflectance as it would be measured at ground level in the absence of atmospheric effects; in the near future, these products are anticipated to become the standard products of the Landsat sensors. The purpose of this paper is to present a radiometric performance evaluation of Level-1 and Level-2 data products for the Landsat-7 Enhanced Thematic Mapper Plus (ETM+) and Landsat-8 Operational Land Imager (OLI) sensors. TOA reflectance and derived surface reflectance values from both data products were evaluated and compared to in situ measurements from eight test sites located in Turkey, Brazil, Chile, the United States, France, and Namibia. The results indicate an agreement between the ETM+ and OLI Level-1 TOA reflectance data and the in situ measurements of 3.9% to 6.5% and 3.9% to 6.0%, respectively, across all spectral bands. Agreement between the in situ measurements and both Level-2 surface reflectance data products were consistently decreased in the shorter wavelength bands, and slightly better in the longer wavelength bands. The mean percent absolute error for Level-2 surface reflectance data ranged from 3.3% to 10% for both Landsat sensors. The significant decay in agreement with the data collected in situ in the short wavelength spectral bands with Level-2 data suggests issues with retrieval of aerosol concentration at some sites. In contrast, the results indicate a reasonably accurate estimate of water vapor in the mid-infrared spectrum. Lastly, despite the less reliable performance of Level-2 data product in the visible spectral region (compared with Level-1 data) in both sensors, the Landsat-8 OLI Level-2 showed an improvement of surface reflectance product over all spectral bands in common with the Landsat-7 ETM+ Level-2 data.


2019 ◽  
Vol 11 (13) ◽  
pp. 1547 ◽  
Author(s):  
Najib Djamai ◽  
Detang Zhong ◽  
Richard Fernandes ◽  
Fuqun Zhou

Time series of vegetation biophysical variables (leaf area index (LAI), fraction canopy cover (FCOVER), fraction of absorbed photosynthetically active radiation (FAPAR), canopy chlorophyll content (CCC), and canopy water content (CWC)) were estimated from interpolated Sentinel-2 (S2-LIKE) surface reflectance images, for an agricultural region located in central Canada, using the Simplified Level 2 Product Prototype Processor (SL2P). S2-LIKE surface reflectance data were generated by blending clear-sky Sentinel-2 Multispectral Imager (S2-MSI) images with daily BRDF-adjusted Moderate Resolution Imaging Spectrometer images using the Prediction Smooth Reflectance Fusion Model (PSFRM), and validated using thirteen independent S2-MSI images (RMSE ≤ 6%). The uncertainty of S2-LIKE surface reflectance data increases with the time delay between the prediction date and the closest S2-MSI image used for training PSFRM. Vegetation biophysical variables from S2-LIKE products are validated qualitatively and quantitatively by comparison to the corresponding vegetation biophysical variables from S2-MSI products (RMSE = 0.55 for LAI, ~10% for FCOVER and FAPAR, and 0.13 g/m2 for CCC and 0.16 kg/m2 for CWC). Uncertainties of vegetation biophysical variables derived from S2-LIKE products are almost linearly related to the uncertainty of the input reflectance data. When compared to the in situ measurements collected during the Soil Moisture Active Passive Validation Experiment 2016 field campaign, uncertainties of LAI (0.83) and FCOVER (13.73%) estimates from S2-LIKE products were slightly larger than uncertainties of LAI (0.57) and FCOVER (11.80%) estimates from S2-MSI products. However, equal uncertainties (0.32 kg/m2) were obtained for CWC estimates using SL2P with either S2-LIKE or S2-MSI input data.


2020 ◽  
Vol 12 (18) ◽  
pp. 3109 ◽  
Author(s):  
Manjunatha Venkatappa ◽  
Sutee Anantsuksomsri ◽  
Jose Alan Castillo ◽  
Benjamin Smith ◽  
Nophea Sasaki

Although vegetation phenology thresholds have been developed for a wide range of mapping applications, their use for assessing the distribution of natural bamboo and the related carbon stocks is still limited, especially in Southeast Asia. Here, we used Google Earth Engine (GEE) to collect time-series of Landsat 8 Operational Land Imager (OLI) and Sentinel-2 images and employed a phenology-based threshold classification method (PBTC) to map the natural bamboo distribution and estimate carbon stocks in Siem Reap Province, Cambodia. We processed 337 collections of Landsat 8 OLI for phenological assessment and generated 121 phenological profiles of the average vegetation index for three vegetation land cover categories from 2015 to 2018. After determining the minimum and maximum threshold values for bamboo during the leaf-shedding phenology stage, the PBTC method was applied to produce a seasonal composite enhanced vegetation index (EVI) for Landsat collections and assess the bamboo distributions in 2015 and 2018. Bamboo distributions in 2019 were then mapped by applying the EVI phenological threshold values for 10 m resolution Sentinel-2 satellite imagery by accessing 442 tiles. The overall Landsat 8 OLI bamboo maps for 2015 and 2018 had user’s accuracies (UAs) of 86.6% and 87.9% and producer’s accuracies (PAs) of 95.7% and 97.8%, respectively, and a UA of 86.5% and PA of 91.7% were obtained from Sentinel-2 imagery for 2019. Accordingly, carbon stocks of natural bamboo by district in Siem Reap at the province level were estimated. Emission reductions from the protection of natural bamboo can be used to offset 6% of the carbon emissions from tourists who visit this tourism-destination province. It is concluded that a combination of GEE and PBTC and the increasing availability of remote sensing data make it possible to map the natural distribution of bamboo and carbon stocks.


Sign in / Sign up

Export Citation Format

Share Document