demTools: Compute terrain characteristics from digital elevation models

In satellite: Handling and Manipulating Remote Sensing Data

Description

Compute terrain characteristics from digital elevation models (DEM) using raster::terrain or raster::hillShade.

Usage

## S4 method for signature 'Satellite'
demTools(x, method = "hillShade", bcde = "DEM")

## S4 method for signature 'RasterLayer'
demTools(x, sunElev, sunAzim, method = "hillShade")

Arguments

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 method = "hillShade", the elevation angle of the sun in degrees. See parameter angle in hillShade.
sunAzim	If method = "hillShade", the sun azimuth angle in degree. See parameter direction in hillShade.

Value

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.

See Also

raster::terrain, raster::hillShade.

Examples

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)

satellite documentation built on Oct. 12, 2021, 5:07 p.m.