Convert selected layers of a Satellite object to a RasterStack
# S4 method for Satellite
stack(x, layer = names(x), ...)
an object of class 'Satellite'
character vector (bcde codes) or integer vector (index) of the layers to be stacked
additional arguments passed on to stack
path <- system.file("extdata", package = "satellite")
files <- list.files(path, pattern = glob2rx("LC08*.TIF"), full.names = TRUE)
sat <- satellite(files)
stck <- stack(sat, c("B001n", "B002n", "B003n"))
stck
Run the code above in your browser using DataLab