draft/demo-forestgap.R

In fdschneider/caspr: Cellular Automata for Spatial Pressure in R

# 
# A file for testing the use of the forest gap model
# 

initial <- init_landscape(states=c("+","0"), 
                          cover=c(.5,.5))

parms <- forestgap$parms # Import defaults
parms$delta <- .15

output <- ca(initial, parms, model=forestgap, t_max=1000)

library(ggplot2)
output_fmt <- fortify(output)
ggplot(output_fmt) + 
  geom_raster(aes(x,y,fill=state)) + 
  facet_wrap( ~ time ) + 
  scale_fill_manual(values=c('#111111','#FFFFFF'))