The following example shows how to read a Landsat-8 data set downloaded from the USGS Earth Explorer and how to convert the bands’ scaled counts to radiance, reflectance and/or brightness temperature.
We will use the test data set supplied with this package. In order to read the data set, just supply the filename (or list of filenames) to the
## Loading required package: raster
## Loading required package: sp
<- system.file("extdata", package = "satellite") path <- list.files(path, pattern = glob2rx("LC08*.TIF"), full.names = TRUE) files <- satellite(files)sat
The meta data will automatically be loaded since (i) Landsat-8 is a known platform to the
satellite package and (ii) the name of the meta data file is in accordance with the USGS conventions.
In order to convert scaled counts of a Satellite object to radiance, reflectance and/or temperature, just supply it to the
<- convSC2Rad(sat) sat <- convSC2Ref(sat) sat <- convRad2BT(sat)sat