This function calculates the TOA (Top Of Atmosphere) reflectance considering only the image metadata.
calcTOAr(image.DN, sat = "auto", aoi, incidence.rel, MTL)raw image in digital numbers
"L7" for Landsat 7, "L8" for Landsat 8 or "auto" to guess from filenames
area of interest to crop images, if waterOptions("autoAoi") == TRUE will look for any object called aoi on .GlobalEnv
solar incidence angle, considering the relief
Landsat Metadata File
R. G. Allen, M. Tasumi, and R. Trezza, "Satellite-based energy balance for mapping evapotranspiration with internalized calibration (METRIC) - Model" Journal of Irrigation and Drainage Engineering, vol. 133, p. 380, 2007
LPSO. (2004). Landsat 7 science data users handbook, Landsat Project Science Office, NASA Goddard Space Flight Center, Greenbelt, Md., (http://landsathandbook.gsfc.nasa.gov/) (Feb. 5, 2007)
Other remote sensing support functions: calcRadiance,
calcSR, cfmask,
loadImageSR, loadImage