# NOT RUN {
## Plum Island Ecosystems
obs <- ObsLulcRasterStack(x=pie,
pattern="lu",
categories=c(1,2,3),
labels=c("forest","built","other"),
t=c(0,6,14))
## Sibuyan Island
obs <- ObsLulcRasterStack(x=sibuyan$maps,
pattern="lu",
categories=c(1,2,3,4,5),
labels=c("forest","coconut","grass","rice","other"),
t=c(0,14))
# }
Run the code above in your browser using DataLab