(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 12 din 30DayMapChallenge: « 2125 » (anterior).
Vidaller et al. (2021) indică faptul că ghețarii din Pirinei se află într-un dezechilibru clar cu clima regională și probabil că vor dispărea în următoarele câteva decenii
.
library(dplyr) library(ggplot2) library(glue) library(sf) library(elevatr) library(rnaturalearth) library(rnaturalearthhires) library(osmdata) library(terra) library(ggspatial) library(ggnewscale) library(ggrepel) sf_use_s2(FALSE)
Date
pyrenees <- getbb("Pyrénées", format_out = "sf_polygon")$polygon
if (!file.exists("elevation_pyrenees.tif")) {
get_elev_raster(pyrenees, z = 9) |>
writeRaster("elevation_pyrenees.tif")
}
elevation_pyrenees <- rast("elevation_pyrenees.tif")
elevation_pyrenees(elevation_pyrenees <=0) <- NA
# map background
world <- ne_countries(scale = 10) |>
st_intersection(elevation_pyrenees |>
st_bbox() |>
st_as_sfc())
slope <- terrain(elevation_pyrenees, "slope", unit = "radians") aspect <- terrain(elevation_pyrenees, "aspect", unit = "radians") hillshade <- shade(slope, aspect, angle = 45, direction = 315, normalize = TRUE)
Hartă
# we create a void layer to train the scale and have a legend
glaciers <- st_polygon() |>
st_sfc() |>
st_sf() |>
st_set_crs("EPSG:4326") |>
mutate(glaciers = NA_character_)
ggplot() +
geom_sf(data = world, fill = "#c4c4c4", color = "#c4c4c4") +
layer_spatial(hillshade) +
scale_fill_distiller(palette = "Greys", guide = "none",
na.value = "transparent") +
geom_sf(data = world, fill = NA, color = "#333333") +
new_scale_fill() +
geom_sf(data = glaciers, aes(fill = glaciers)) +
scale_fill_brewer(palette="RdYlBu", na.value="blue") +
geom_text_repel(data = world, aes(label = sovereignt, geometry = geometry),
stat = "sf_coordinates", color = "#333333",
bg.colour = "#ffffffaa", bg.r = .2, alpha = 0.4) +
coord_sf() +
labs(title = "Pyrenean glaciers",
subtitle = "2125",
x = "", y = "",
caption = glue("https://r.iresmi.net/ - {Sys.Date()}
Elevation data from //
https://registry.opendata.aws/terrain-tiles
Natural Earth")) +
theme_minimal() +
theme(plot.caption = element_text(size = 6,
color = "darkgrey"),
legend.position = "bottom",
legend.text = element_blank())
Figura 1: Harta ghețarilor din Pirinei în 2025
O hartă tristă cu siguranță…
Referințe
Vidaller, I., J. Revuelto, E. Izagirre, F. Rojas-Heredia, E. Alonso-González, S. Gascoin, P. René, et al. 2021. „Spre un Lanțul muntos fără gheață: Deces de Ghetarii Pirinei In timpul 2011–2020.” Scrisori de cercetare geofizică 48 (18): e2021GL094339. https://doi.org/10.1029/2021GL094339.
