If the sensor does not provide linear conversion coefficients for reflectance
computation, the reflectance is calculated using the solar irradiance
following the functions taken from USGS' Landsat 7 manual, chapter 11.3.2,
which is available online at
http://landsathandbook.gsfc.nasa.gov/data_prod/prog_sect11_3.html.
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).
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 <- convRad2Ref(sat)
# If you use a raster layer, supply required meta information
bcde <- "B002n"
convRad2Ref(x = getSatDataLayer(sat, bcde),
mult = getSatRADM(sat, bcde),
add = getSatRADA(sat, bcde))