Convert a band's scaled counts or radiance values to reflectance
Convert a band's scaled counts to reflectance using a simple linear conversion without any kind of atmospheric correction etc.
## S4 method for signature 'Satellite' convRad2Ref(x, szen_correction = "TRUE") ## S4 method for signature 'RasterStack' convRad2Ref(x, mult, add, szen) ## S4 method for signature 'RasterLayer' convRad2Ref(x, mult, add, szen)
x | 
 An object of class Satellite, raster::RasterStack or raster::RasterLayer providing radiance values.  | 
szen_correction | 
 Logical; if   | 
mult | 
 Multiplicative coefficient for value transformation (i.e. slope).  | 
add | 
 Additive coefficient for value transformation (i.e. offset)  | 
szen | 
 Cosine of solar zenith angle.  | 
The conversion functions are taken from USGS' Landsat 8 Data Users Handbook which is available online at https://www.usgs.gov/land-resources/nli/landsat/landsat-8-data-users-handbook.
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 https://landsat.gsfc.nasa.gov/wp-content/uploads/2016/08/Landsat7_Handbook.pdf.
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).
calcAtmosCorr for conversions of scaled counts 
to physical units including a scene-based atmospheric correction.
path <- system.file("extdata", package = "satellite")
files <- list.files(path, pattern = glob2rx("LC08*.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))Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.