rnaturalearth
geodata
osmdata
We are limiting the extend:
srtm_cropped = terra::crop(srtm, zion) plot(srtm_cropped, col = grDevices::terrain.colors(50)) plot(zion$geom, add = TRUE)