# Road building (stored in the landscape object and called "roads") acts on the landscape
# each year.
if (FALSE) {
road_effect <- disturbance(disturbance_layers = "roads", effect_time = 1)
ls <- landscape(population = egk_pop, suitability = egk_hab, "roads" = egk_road)
pd <- population_dynamics(change = growth(egk_mat))
sim <- simulation(landscape = ls,
population_dynamics = pd,
habitat_dynamics = list(road_effect),
timesteps = 20)
plot(sim, object = "suitability", type = "raster", timesteps = 1:9)
}
Run the code above in your browser using DataLab