(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: „O călătorie” (anterior).
Treisprezece balene albastre (Balaenoptera musculus) au fost etichetate prin satelit în 2015 în Australia (Mӧller și colab. 2020). Le vom urmări vara în Great Southern Australian Coastal Upwelling System (GSACUS).
Datele sunt disponibile de la Australian Antarctic Data Center (Andrews-Goff și colab. 2020)direct ca fișier R RDS.
library(dplyr)
library(readr)
library(sf)
library(leaflet)
unzip("AAS_4101_pygmy_blue_whale_SSSM.zip",
files = "bw_3h_ssm.RDS")
bm <- read_rds("bw_3h_ssm.RDS")
# str(bm)
journeys <- bm$summary |>
st_as_sf(coords = c("lon", "lat"), crs = "EPSG:4326") |>
group_by(id) |>
arrange(date) |>
summarise(do_union = FALSE) %>%
st_cast("LINESTRING") |>
group_by(indiv = as.factor(parse_number(as.character(id)))) |>
summarise()
După unirea punctelor GPS ca linii și linii de la aceeași balenă cu mai multe linii, putem folosi {leaflet} pentru o hartă interactivă.
pal <- colorFactor("viridis", sort(unique(journeys$indiv)))
journeys |>
leaflet() |>
addTiles() |>
addPolylines(color = ~ pal(indiv),
popup = ~ indiv) |>
addLegend(pal = pal,
values = sort(unique(journeys$indiv)))
