## sample data
path <- system.file("extdata", package = "satellite")
files <- list.files(path, pattern = glob2rx("LC08*.TIF"), full.names = TRUE)
sat <- satellite(files)
sat[[2:5]]
subset(sat, cid = "SC")
Run the code above in your browser using DataLab