path <- system.file("extdata", package = "satellite")
files <- list.files(path, pattern = glob2rx("LC08*.TIF"), full.names = TRUE)
sat <- satellite(files)
brck <- brick(sat, c("B001n", "B002n", "B003n"))
brck
Run the code above in your browser using DataLab