# NOT RUN {
## Load raster landscape
tiny <- raster::raster(system.file("extdata/tiny.asc", package = "grainscape"))
## Create a resistance surface from a raster using an is-becomes reclassification
tinyCost <- raster::reclassify(tiny, rcl = cbind(c(1, 2, 3, 4), c(1, 5, 10, 12)))
## Produce a patch-based MPG where patches are resistance features=1
tinyPatchMPG <- MPG(cost = tinyCost, patch = tinyCost == 1)
if (interactive()) {
library(ggplot2)
## Plot the patches in a minimum planar graph
ggplot() +
geom_tile(ggGS(tinyPatchMPG, "patchId"), mapping = aes(x = x, y = y, fill = value)) +
theme_grainscape()
}
# }
Run the code above in your browser using DataLab