lstetm: Calculate Brightness Temperature as proxy for Surface...

In DGampe/TriangleMethod: Calculation of NDVI, LST and finally provide an estimation for actual evapotranspiration patterns using remote sensing data (default LANDSAT TM and ETM+ scenes).

Calculate Brightness Temperature as proxy for Surface Temperature (Ts) from Landsat ETM+ scenes.

Description

Brightness Temperature as proxy for Surface Temperature (Ts) will be calculated in degC or K from Landsat ETM+ thermal bands (as ETM+ has upper and lower thermal bands, the average radiation will be used.) here. Default parameters to calibrate the temperature at the sensor etc. are implemented for the ETM+ sensor. If TM scenes should be processed, see function lsttm. If LST should be calculated from other remote sensing data adjust parameters for calibration. NOTE: the parameters here match most of the scenes, however to be sure, check the corresponding calibration file for lmax and lmin!

Usage

  lstetm(TEMPL,TEMPH, ...)

Arguments

TEMPL	Lower thermal band (B61) for Landsat ETM+ or other sensors image, as raster data set - in memory, or as files (which will be rastered).
TEMPH	Upper thermal band (B62) for Landsat ETM+, as raster data set - in memory, or as files (which will be rastered).
write	Specify, if the resulting Ts grid should be written to file.Default = F (T/F)
outpath	User defined outpath to write the final Ts file.
outname	User defined name for the resulting file.(Default name / format: 'ts.tif'
unit	Define, if Ts should be calculated in degC or K.
LminLow,...	Calibration parameter files for the sensor. Landsat ETM+ is the default parameterization, if images originate from a different sensor, these parameters should be adjusted.

Details

Calibration for ETM+ Low gain thermal band (band B61)

• LminLow: 0.0*10^6 (can be found in calibration parameter file for Landsat scene)

• LmaxLow: 17.04*10^6

  # Calibration for ETM+ High gain thermal band (band B62)

• LminHigh: 3.2*10^6

• LmaxHigh: 12.65*10^6 # values here for Lmin+Lmax for ETM+ B5 Low Gain data

• QcalMin:Minimum digital number, default 1

• QcalMax:Maximum digital number, default 255

• Calculation of radiation:

  RadLow= ((LmaxLow-LminLow)/(QcalMax-QcalMin)) * (TEMPL - QcalMin)+LminLow

  RadHigh= ((LmaxHigh-LminHigh)/(QcalMax-QcalMin)) * (TEMPH - QcalMin)+LminHigh

  RadAvg = (RadLow+RadHigh)/2 # Average of the Low Gain and High Gain thermal bands

  see lsttm for additional information.

RasterLayer, .tif file

DGampe/TriangleMethod documentation built on March 18, 2022, 6:43 a.m.