satellite (version 0.1.0)

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

## S3 method for class 'Satellite':
convSC2Rad(x, szen_correction = "TRUE")
## S3 method for class 'RasterStack':
convSC2Rad(x, mult, add, szen)
## S3 method for class 'RasterLayer':
convSC2Rad(x, mult, add, szen)

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.

mult

Multiplicative coefficient for value transformation (i.e. slope).

add

Additive coefficient for value transformation (i.e. offset).

szen

Cosine of solar zenith angle.

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).

Details

The conversion functions are taken from USGS' Landsat 8 manual which is available online at http://landsat.usgs.gov/Landsat8_Using_Product.php.

See Also

calcAtmosCorr for conversions of scaled counts to physical units including a scene-based atmospheric correction.

Examples

path <- system.file("extdata", package = "satellite")
files <- list.files(path, pattern = glob2rx("LC8*.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))