(Acest articol a fost publicat pentru prima dată pe r.iresmi.netși cu amabilitate a contribuit la R-bloggeri). (Puteți raporta problema legată de conținutul acestei pagini aici)
Doriți să vă distribuiți conținutul pe R-bloggeri? dați clic aici dacă aveți un blog, sau aici dacă nu aveți.
Ziua 5 din 30DayMapChallenge: „Pământ” (anterior).
Ce se întâmplă dacă batimetria Mării Mediterane este inversată? Am obține o insulă de vest și o insulă de est, un golf Sardinia-Corsica, laguna Baleare… Să vedem!
library(terra) library(sf) terraOptions(progress=0)
Un set de date de batimetrie globală (7.5 Go zipped NetCDF) este disponibil la GEBCO. O decupăm în jurul Mării Mediterane.
med <- c(xmin = -7, ymin = 30, xmax = 39, ymax = 45) |>
st_bbox() |>
st_as_sfc()
gebco <- rast("~/data/gebco/GEBCO_2025_sub_ice.nc") |>
crop(med)
Setăm toate valorile pozitive la 0 și toate valorile negative sunt inversate. Relieful este mai pronunțat cu o rezoluție ușoară mai mică, apoi calculăm umbrirea dealului.
mountain_med <- aggregate((gebco <= 0) * -1 * gebco, 10, mean) slope <- terrain(mountain_med, "slope", unit = "radians") aspect <- terrain(mountain_med, "aspect", unit = "radians") hill <- shade(slope, aspect, angle = 45, direction = 315, normalize = TRUE)
Hartă
plot(hill, col = grey(0:255/255),
main = "Mediterranean Sea upside down",
legend = FALSE, mar = c(2,2,2,4))
plot(mountain_med, col = topo.colors(30, alpha = 0.3),
plg = list(title = "Elevation (m)"), add = TRUE)
mtext(paste("data: GEBCOnhttps://r.iresmi.net", Sys.Date()),
side = 1, line = 3, cex = 0.5, adj = 1)
Figura 1: Marea Mediterană cu susul în jos
