Compute terrain characteristics from digital elevation models
Compute terrain characteristics from digital elevation models (DEM) using
raster::terrain or raster::hillShade.
## S4 method for signature 'Satellite' demTools(x, method = "hillShade", bcde = "DEM") ## S4 method for signature 'RasterLayer' demTools(x, sunElev, sunAzim, method = "hillShade")
x |
A DEM provided as an object of class Satellite or RasterLayer. |
method |
Currently "slope", "aspect" and "hillshade" are implemented. |
bcde |
The name of the DEM layer in the Satellite object. |
sunElev |
If |
sunAzim |
If |
If x is a Satellite object, a Satellite object with added layer containing calculated
terrain information; if x is a raster::RasterLayer object, a
raster::RasterLayer object with calculated terrain information.
raster::terrain, raster::hillShade.
path <- system.file("extdata", package = "satellite")
files <- list.files(path, pattern = glob2rx("LC08*.TIF"), full.names = TRUE)
sat <- satellite(files)
## dem
files_dem <- list.files(path, pattern = "DEM", full.names = TRUE)
DEM <- raster(files_dem)
sat <- addSatDataLayer(sat, data = DEM, info = NULL, bcde = "DEM", in_bcde="DEM")
sat <- demTools(sat)Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.