## rgeedim v0.4.0 -- using geedim 2.0.0 w/ earthengine-api 1.7.10
## terra 1.8.97
## Using Application Default Credentials (ADC)
ee <- earthengine()
p <- ee$FeatureCollection("WCMC/WDPA/current/polygons")$filter(
ee$Filter$And(
ee$Filter$eq("NAME", "Yellowstone"),
ee$Filter$eq("DESIG_ENG", "National Park")
)
)
y <- gd_region(p) |>
gd_region_to_vect()
b <- gd_bbox(p) |>
gd_region() |>
gd_region_to_vect()
r <- gd_image_from_id("USGS/SRTMGL1_003") |>
gd_download(
"dem.tif",
region = b,
overwrite = TRUE,
resampling = "bilinear",
scale = 100,
crs = "EPSG:4326",
bands = list("elevation")
) |>
terra::rast()
plot(r)
plot(as.lines(y), add = TRUE, col = "red", lwd = 3)
