Deriving Land Surface Temperature Using the Earth Engine Landsat LST Toolbox

In the previous section, we explored an LST retrieval algorithm to give an example of the standard steps to get to LST from the satellite measurements in the thermal bands. In this section, we will use an Earth Engine module developed for this purpose to calculate LST.

// Link to the module that computes the Landsat LST.
var landsatLST = require('projects/gee-edu/book:Part A - Applications/A1 - Human Applications/A1.5 Heat Islands/modules/Landsat_LST.js');

// Select region of interest, date range, and Landsat satellite.
var geometry = regionInt.geometry();
var satellite = 'L8';
var dateStart = '2014-01-01';
var dateEnd = '2019-01-01';
var useNdvi = true;

// Get Landsat collection with additional necessary variables.
var landsatColl = landsatLST.collection(satellite, dateStart, dateEnd, geometry, useNdvi);

// Create composite, clip, filter to summer, mask, and convert to degree Celsius.
var landsatComp = landsatColl
  .select('LST')
  .filter(sumFilter)
  .median()
  .clip(regionInt)
  .updateMask(notWater)
  .subtract(273.15);

Map.addLayer(landsatComp, {
    min: 25,
    max: 38,
    palette: ['blue', 'white', 'red']
  },
  'LST_SMW');
LST SMW screenshot

As an aside, the Landsat Collection 2 products have recently incorporated LST bands, which can be processed similarly to the MODIS data, but with the bands’ own specific offsets and scaling factors.