A Software Tool for Atmospheric Correction and Surface Temperature Estimation of Landsat Infrared Thermal Data

Land surface temperature (LST) is an important variable involved in the Earth's surface energy and water budgets and a key component in many aspects of environmental research. The Landsat program, jointly carried out by NASA and the USGS, has been recording thermal infrared data for the past 40 years. Nevertheless, LST data products for Landsat remain unavailable. The atmospheric correction (AC) method commonly used for mono-window Landsat thermal data requires detailed information concerning the vertical structure (temperature, pressure) and the composition (water vapor, ozone) of the atmosphere. For a given coordinate, this information is generally obtained through either radio-sounding or atmospheric model simulations and is passed to the radiative transfer model (RTM) to estimate the local atmospheric correction parameters. Although this approach yields accurate LST data, results are relevant only near this given coordinate. To meet the scientific community's demand for high-resolution LST maps, we developed a new software tool dedicated to processing Landsat thermal data. The proposed tool improves on the commonly-used AC algorithm by incorporating spatial variations occurring in the Earth's atmosphere composition. The ERA-Interim dataset (ECMWFmeteorological organization) was used to retrieve vertical atmospheric conditions, which are available at a global scale with a resolution of 0.125 degrees and a temporal resolution of 6 h. A temporal and spatial linear interpolation of meteorological variables was performed to match the acquisition dates and coordinates of the Landsat images. The atmospheric correction parameters were then estimated on the basis of this reconstructed atmospheric grid using the commercial RTMsoftware MODTRAN. The needed surface emissivity was derived from the common vegetation index NDVI, obtained from the red and near-infrared (NIR) bands of the same Landsat image. This permitted an estimation of LST for the entire image without degradation in resolution. The software tool, named LANDARTs, which stands for Landsat automatic retrieval of surface temperatures, is fully automatic and coded in the programming language Python. In the present paper, LANDARTs was used for the local and spatial validation of surface temperature obtained from a Landsat dataset covering two climatically contrasting zones: southwestern France and central Tunisia. Long-term datasets of in situ surface temperature measurements for both locations were compared to corresponding Landsat LST data. This temporal comparison yielded RMSE values ranging from 1.84 • C–2.55 • C. Landsat surface temperature data obtained with LANDARTs were then spatially compared using the ASTER data products of kinetic surface temperatures (AST08) for both geographical zones. This comparison yielded a satisfactory RMSE of about 2.55 • C. Finally, a sensitivity analysis for the effect of spatial validation on the LST correction process showed a variability of up to 2 • C for an entire Landsat image, confirming that the proposed spatial approach improved the accuracy of Landsat LST estimations.

[1]  Albert Olioso,et al.  Estimating the difference between brightness and surface temperatures for a vegetal canopy , 1995 .

[2]  S. Idso.  A set of equations for full spectrum and 8- to 14-μm and 10.5- to 12.5-μm thermal radiation from cloudless skies , 1981 .

[3]  Miquel Ninyerola,et al.  Revision of the Single-Channel Algorithm for Land Surface Temperature Retrieval From Landsat Thermal-Infrared Data , 2009, IEEE Transactions on Geoscience and Remote Sensing.

[4]  Albert Olioso,et al.  Impact of surface emissivity and atmospheric conditions on surface temperatures estimated from top of canopy brightness temperatures derived from Landsat 7 data , 2013, 2013 IEEE International Geoscience and Remote Sensing Symposium - IGARSS.

[5]  Jeff Dozier,et al.  A generalized split-window algorithm for retrieving land-surface temperature from space , 1996, IEEE Trans. Geosci. Remote. Sens..

[6]  Zhao-Liang Li,et al.  Angular effect of MODIS emissivity products and its application to the split-window algorithm , 2011 .

[7]  Michaela Buenemann,et al.  ATMOSPHERIC CORRECTION OF LANDSAT THERMAL INFRARED DATA : A CALCULATOR BASED ON NORTH AMERICAN REGIONAL REANALY SIS ( NARR ) DATA , 2011 .

[8]  Z. Wan,et al.  Quality assessment and validation of the MODIS global land surface temperature , 2004 .

[9]  Vincent Rivalland,et al.  Crops’ water use efficiencies in temperate climate: Comparison of stand, ecosystem and agronomical approaches , 2013 .

[10]  T. J. Majumdar,et al.  Study of land surface temperature and spectral emissivity using multi-sensor satellite data , 2010 .

[11]  Ji Zhou,et al.  Intercomparison of methods for estimating land surface temperature from a Landsat-5 TM image in an arid region with low water vapour in the atmosphere , 2012 .

[12]  H. Roquet,et al.  Use of satellite observations for operational oceanography: recent achievements and future prospects , 2015 .

[13]  K. Moffett,et al.  Remote Sens , 2015 .

[14]  José A. Sobrino,et al.  Land surface temperature retrieval from LANDSAT TM 5 , 2004 .

[15]  P. Béziat,et al.  Carbon balance of a three crop succession over two cropland sites in South West France , 2009 .

[16]  Vicente Caselles,et al.  Validation of Landsat-7/ETM+ Thermal-Band Calibration and Atmospheric Correction With Ground-Based Measurements , 2010, IEEE Transactions on Geoscience and Remote Sensing.

[17]  Zhao-Liang Li,et al.  Quantitative Remote Sensing in Thermal Infrared: Theory and Applications , 2014 .

[18]  M. Abrams.  The Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER): Data products for the high spatial resolution imager on NASA's Terra platform , 2000 .

[19]  M. Lawrence.  The relationship between relative humidity and the dewpoint temperature in moist air - A simple conversion and applications , 2005 .

[20]  Gérard Dedieu,et al.  A multi-temporal method for cloud detection, applied to FORMOSAT-2, VENµS, LANDSAT and SENTINEL-2 images , 2010 .

[21]  I. F. Trigo,et al.  Kalman filter physical retrieval of surface emissivity and temperature from SEVIRI infrared channels: a validation and intercomparison study , 2015 .

[22]  Alan R. Gillespie,et al.  Temperature / Emissivity Separation Algorithm Theoretical Basis Document , Version 2 . 4 , 1999 .

[23]  James A. Gardner,et al.  MODTRAN5: a reformulated atmospheric band model with auxiliary species and practical multiple scattering options , 2004, SPIE Asia-Pacific Remote Sensing.

[24]  S. Hook,et al.  The ASTER spectral library version 2.0 , 2009 .

[25]  John R. Schott,et al.  Validation of a web-based atmospheric correction tool for single thermal band instruments , 2005, SPIE Optics + Photonics.

[26]  J. Sobrino,et al.  A generalized single‐channel method for retrieving land surface temperature from remote sensing data , 2003 .

[27]  Martha C. Anderson,et al.  Use of Landsat thermal imagery in monitoring evapotranspiration and managing water resources , 2012 .

[28]  K. Wittich,et al.  Some simple relationships between land-surface emissivity, greenness and the plant cover fraction for use in satellite remote sensing , 1997 .

[29]  B. Markham,et al.  Summary of Current Radiometric Calibration Coefficients for Landsat MSS, TM, ETM+, and EO-1 ALI Sensors , 2009 .

[30]  Antonio J. Plaza,et al.  Land Surface Emissivity Retrieval From Different VNIR and TIR Sensors , 2008, IEEE Transactions on Geoscience and Remote Sensing.

[31]  Martha C. Anderson,et al.  Advances in thermal infrared remote sensing for land surface modeling , 2009 .

[32]  R. Eskridge,et al.  Improved Magnus form approximation of saturation vapor pressure , 1996 .

[33]  Julia A. Barsi,et al.  An Atmospheric Correction Parameter Calculator for a single thermal band earth-sensing instrument , 2003, IGARSS 2003. 2003 IEEE International Geoscience and Remote Sensing Symposium. Proceedings (IEEE Cat. No.03CH37477).

[34]  Zhilin Zhu,et al.  Study of emissivity scaling and relativity of homogeneity of surface temperature , 2004 .

[35]  Shuichi Rokugawa,et al.  A temperature and emissivity separation algorithm for Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) images , 1998, IEEE Trans. Geosci. Remote. Sens..

[36]  John R. Schott,et al.  Landsat-8 Thermal Infrared Sensor (TIRS) Vicarious Radiometric Calibration , 2014, Remote. Sens..

[37]  Martha C. Anderson,et al.  Evaluation of Drought Indices Based on Thermal Remote Sensing of Evapotranspiration over the Continental United States , 2011 .

[38]  Albert Olioso,et al.  Quantifying uncertainties in land surface temperature due to atmospheric correction: Application to Landsat-7 data over a Mediterranean agricultural region , 2014, 2014 IEEE Geoscience and Remote Sensing Symposium.

[39]  Ji Zhou,et al.  Validation and Performance Evaluations of Methods for Estimating Land Surface Temperatures from ASTER Data in the Middle Reach of the Heihe River Basin, Northwest China , 2015, Remote. Sens..

[40]  A. Karnieli,et al.  A mono-window algorithm for retrieving land surface temperature from Landsat TM data and its application to the Israel-Egypt border region , 2001 .

[41]  J. Thepaut,et al.  The ERA‐Interim reanalysis: configuration and performance of the data assimilation system , 2011 .

[42]  Ning Wang,et al.  Land surface emissivity retrieval from satellite data , 2013 .

[43]  I. D. Feis,et al.  Kalman filter physical retrieval of surface emissivity and temperature from geostationary infrared radiances , 2013 .

[44]  Z. Wan,et al.  Evaluation of Six Methods for Extracting Relative Emissivity Spectra from Thermal Infrared Images , 1999 .

[45]  Paul E. Lewis,et al.  MODTRAN 5: a reformulated atmospheric band model with auxiliary species and practical multiple scattering options: update , 2005 .