convSC2Rad: Convert a band's scaled counts to radiance

Description

Convert a band's scaled counts to radiance using a simple linear conversion without any kind of atmospheric correction etc.

Usage

# S4 method for Satellite
convSC2Rad(x, szen_correction = "TRUE", subset = FALSE)
# S4 method for RasterStack
convSC2Rad(x, mult, add, szen)
# S4 method for RasterLayer
convSC2Rad(x, mult, add, szen)

Value

If x is a Satellite object, a Satellite object with added converted layers;

if x is a raster::Raster* object, a raster::Raster* object with converted layer(s).

Arguments

x: An object of class Satellite, raster::RasterStack or raster::RasterLayer providing scaled counts (DNs).
szen_correction: Logical; if TRUE, sun zenith correction is being applied.
subset: Logical; if TRUE, all layers but the cropped ones are being dropped; if FALSE (default), cropped layers are appended to the Satellite object.
mult: Multiplicative coefficient for value transformation (i.e. slope).
add: Additive coefficient for value transformation (i.e. offset).
szen: Cosine of solar zenith angle.

Details

The conversion functions are taken from USGS' Landsat 8 Data Users Handbook which is available online at https://www.usgs.gov/landsat-missions/landsat-8-data-users-handbook.

Examples

Run this code

path <- system.file("extdata", package = "satellite")
files <- list.files(path, pattern = glob2rx("LC08*.TIF"), full.names = TRUE)
sat <- satellite(files)  
sat <- convSC2Rad(sat)

# If you use a raster layer, supply required meta information
bcde <- "B002n"
convSC2Rad(x = getSatDataLayer(sat, bcde),
           mult = getSatRADM(sat, bcde),
           add = getSatRADA(sat, bcde))