(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 7 din 30DayMapChallenge: «Accesibilitate» (anterior).
Ei bine, să fim răzvrătiți și să căutăm în schimb inaccesibilitatea; mai exact polul de inaccesibilitate al Frantei (Hexagonul): locatia cea mai indepartata de limita. A nu se confunda cu centroidul.
library(sf) library(dplyr) library(ggplot2) library(glue) library(purrr) library(polylabelr)
Date
Vom folosi din nou unitățile administrative franceze (obțineți datele din această postare).
# France boundary
fr <- read_sf("~/data/adminexpress/adminexpress_cog_simpl_000_2022.gpkg",
layer = "region") |>
filter(insee_reg > "06",
insee_reg != "94") |>
st_transform("EPSG:2154") |>
st_union()
# French communes to get the point name
com <- read_sf("~/data/adminexpress/adminexpress_cog_simpl_000_2022.gpkg",
layer = "commune") |>
filter(insee_reg > "06",
insee_reg != "94") |>
st_transform("EPSG:2154")
Calculați POI
Obțineți polul de inaccesibilitate al Franței cu {polylabelr} și se intersectează cu comuna strat pentru a găsi cel mai apropiat oraș.
fr_poi <- poi(fr) |>
pluck(1) |>
as_tibble() |>
st_as_sf(coords = c("x", "y"), crs = "EPSG:2154") |>
st_join(com)
fr_poi_circle <- fr_poi |>
mutate(geometry = st_buffer(geometry, dist))
fr_centroid <- fr |>
st_centroid()
Se pare că este în Saint-Palais în Cher departament.
Hartă
fr_poi |>
ggplot() +
geom_sf(data = fr) +
geom_sf(data = fr_poi_circle, linewidth = 1, linetype = 3) +
geom_sf(data = fr_centroid, color = "darkgrey") +
geom_sf() +
geom_sf_text(aes(label = nom), vjust = -.5) +
labs(title = "Pole of inaccessibility",
subtitle = "France",
x = "", y = "",
caption = glue("https://r.iresmi.net/ - {Sys.Date()}
data from IGN Adminexpress 2022")) +
theme_minimal() +
theme(plot.caption = element_text(size = 6,
color = "darkgrey"))
Figura 1: Stâlpul de inaccesibilitate al Franței (negru). Punctul gri este centroidul
