# NOT RUN {
## Plum Island Ecosystems
## Load observed land use maps
obs <- ObsLulcRasterStack(x=pie,
pattern="lu",
categories=c(1,2,3),
labels=c("forest","built","other"),
t=c(0,6,14))
summary(obs[[1]])
summary(obs[[1:2]])
# }
Run the code above in your browser using DataLab