[GIS] Land Surface Temperature

landsatlandsat 8remote sensing

Estimation of Land Surface Temperature

Contact measurements of emissivity values of the surfaces using thermal radiometers and spectral assessment of samples supplementedLST estimation procedures in this study.Since the temperature values obtained as atsatellite temperature Tʙ,referenced to black bodies,the land cover classes were assigned emissivity values derived through field measurements.Thus the emissivity corrected land surface temperatures can be calculated as [8]:

St = Tʙ ∕1+ [λ* Tʙ∕ρ]Ιnɛ
λ = Wavelength of emitted radiance,ρ =h * c / σ (1.438 * 10^-2 m K)
h = Planck’s constant (6.626 * 10^-34 Js)
σ = Boltzmann constant (1.38 * 10^-23 J/K, c = velocity of light (2.998 * 10^8 m/s)

Source: Estimation of Spatial Variability of Land Surface Temperature using Landsat 8 Imagery, P Rajendran, K Mani, 2015 [PDF]

Can anyone tell which algorithm is used for Land Surface Temperature retrieval?

Best Answer

Conversion to At-Satellite Brightness Temperature

http://landsat.usgs.gov/Landsat8_Using_Product.php

TIRS band data can be converted from spectral radiance to brightness temperature using the thermal constants provided in the metadata file:

T = K2 ln( K1 +1) Lλ where:

T = At-satellite brightness temperature (K)

Lλ = TOA spectral radiance (Watts/( m2 * srad * μm))

K1 = Band-specific thermal conversion constant from the metadata

(K1_CONSTANT_BAND_x, where x is the thermal band number)

K2 + Band-specific thermal conversion constant from the metadata (K2_CONSTANT_BAND_x, where x is the thermal band number)

Standard Landsat 8 data products provided by the USGS EROS Center consist of quantized and calibrated scaled Digital Numbers (DN) representing multispectral image data acquired by both the Operational Land Imager (OLI) and Thermal Infrared Sensor (TIRS).

This is the tutorial. This website provided a full step-by-step tutorial.

https://blogs.esri.com/esri/arcgis/2014/01/06/deriving-temperature-from-landsat-8-thermal-bands-tirs/

You would then need to convert to land temperature using https://www.youtube.com/watch?v=uDQo2a5e7dM

Related Solutions

[GIS] how i calculate land surface temperature using landsat 8 in Envi 5.1

Step 1: Convert from digital numbers (DN) to radiance
This is done by applying the multiplier and addition numbers as found in the metadata (.MTL) file. For the thermal bands (B10 and B11), the values are usually, but you should check the file:
Add: 0.1
Multiply by: 0.0003342 (3.3420E-04)

In ENVI you can apply this correction using 'band math':
float(b10)*0.0003342+0.1
This gives you the radiance value.

Step 2: Convert from radiance to kelvin
The formula needed here is
K2 / ln(K1/TOA_r + 1)
Again, the important values can be found in the metadata file. Usually, the K1 and K2 values are as follows:
K1_CONSTANT_BAND_10 = 774.89
K1_CONSTANT_BAND_11 = 480.89
K2_CONSTANT_BAND_10 = 1321.08
K2_CONSTANT_BAND_11 = 1201.14
In ENVI band math the formula becomes:
1321.08 / alog(774.89/B1+1)
Where alog is the ENVI band math version of the natural log.

This could be combined into one step - for example, band 11 becomes:
1201.14 / alog(480.89/(float(b11)*0.0003342+0.1)+1)

R Programming – Using Conditional Statements for Raster in R for Land Surface Temperature Analysis

I did not notice in my initial answer that NDVI was a vector created from the values of a raster. That is why you are receiving that warning message. It will not evaluate the condition for more than one element. Since you need to evaluate an entire raster and use other rasters in the calculations you can use the method shown here: https://gis.stackexchange.com/a/167514/70980

library(raster)
#Create sample rasters
NDVI <- raster(nrow=5, ncol=5)
NDVI <- setValues(NDVI, runif(ncell(NDVI)))
Pv <- setValues(NDVI, runif(25 ,0.6,0.9))
Palfa.B4 <- setValues(NDVI, runif(25 ,0.2,0.4))
C10 <- setValues(NDVI, runif(25 ,1,2))

#Stack rasters 
s <- stack(NDVI, Pv, Palfa.B4, C10)
#Sample values for constants
es10 <- 0.85
ev10 <- 0.54

#if-else reclassify statement
rc <- function(x1,x2,x3,x4) {
  #x1 - NDVI
  #x2 - Pv
  #x3 - Palfa.B4
  #x4 - C10
  ifelse(x1 <= 0.2, (-0.047 * x3 + 0.973), 
         ifelse((x1 > 0.2) & (x1 <= 0.5), (ev10 * x2 + es10 * (1 - x2) + x4), 
                ev10 + x4))
}

r.class <- overlay(s, fun=rc)

Best Answer

Related Solutions

[GIS] how i calculate land surface temperature using landsat 8 in Envi 5.1

R Programming – Using Conditional Statements for Raster in R for Land Surface Temperature Analysis

Related Question