landsat thermal imaging

July 1, 2014
By

(This article was first published on Dan Kelley Blog/R, and kindly contributed to R-bloggers)

Introduction

The tirs1 band of landsat-8 imagery (see [1]) is a near-thermal band, from which a quantity called ``brightness temperature’’ may be inferred.

As explained in [2], the brightness temperature in Kelvin is computed with where and are stored in the image metadata, and , where again and are stored in the metadata. In this last expression, is the image value, which when read with the oce package must be multiplied by a scale factor for the formulae to work.

Methods

Note in the code that the SST is computed in Celcius, and inserted into the data object before plotting.

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
library(oce)
l <- read.landsat("LC80080292014065LGN00", band="tirs1")
tirs1 <- l[["tirs1"]]
ML <- l@metadata$header$radiance_mult_band_10
AL <- l@metadata$header$radiance_add_band_10
K1 <- l@metadata$header$k1_constant_band_10
K2 <- l@metadata$header$k2_constant_band_10
d <- tirs1 * (2^16 - 1)            # convert from range 0 to 1
Llambda <- ML * d + AL
dd <- K2 / log(K1 / Llambda + 1)
SST <- dd - 273.15                 # convert Kelvin to Celcius
l@data$SST <- SST
plot(l, band="SST", col=oceColorsJet)
mtext(l[["time"]])

Results

The figure below shows the results. This image is from winter, so the land (Nova Scotia) is colder than the ocean.

The high land is especially cold, as those familiar with the topography can tell at a glance. Homogenous blue patches indicate frozen lakes.

Generally, the ocean temperature hovers around freezing (i.e. -2C). The anomalously cold patches in the Bay of Fundy are probably a signature of ice. Some of the orange colour suggests temperatures a much as 1C below freezing. Whether this indicates a difficulty in the algorithm or the interpretation is a matter for Landsat experts.

The beautiful swirls south of Nova Scotia are typical of these coastal waters.

center

Resources

  1. Landsat bands are explained at http://landsat.gsfc.nasa.gov/?page_id=5377.

  2. The conversion to brightness temperature is described at http://landsat.usgs.gov/Landsat8_Using_Product.php

  3. Source code: 2014-07-01-landsat.R.

To leave a comment for the author, please follow the link and comment on his blog: Dan Kelley Blog/R.

R-bloggers.com offers daily e-mail updates about R news and tutorials on topics such as: visualization (ggplot2, Boxplots, maps, animation), programming (RStudio, Sweave, LaTeX, SQL, Eclipse, git, hadoop, Web Scraping) statistics (regression, PCA, time series, trading) and more...



If you got this far, why not subscribe for updates from the site? Choose your flavor: e-mail, twitter, RSS, or facebook...

Comments are closed.