if (FALSE) {
# Continental outline only (default)
cr <- get_cr_outline()
plot(sf::st_geometry(cr), main = "Costa Rica (continental)")
# Full outline including islands
cr_full <- get_cr_outline(continental = FALSE)
plot(sf::st_geometry(cr_full), main = "Costa Rica (all territory)")
# Use as AOI — consistent with h3sdm
h7 <- get_h3_grid(cr, res = 7)
bio <- get_chelsa_historic(var = "bio1", aoi = cr)
# Persistent cache
cr <- get_cr_outline(path = "data/gadm")
}
Run the code above in your browser using DataLab