In dbca-wa/etlTurtleNesting: Data ETL of Turtle Nesting Census Data from ODK Central to WAStD

knitr::opts_chunk$set(echo = TRUE)
library(wastdr)
library(etlTurtleNesting)
library(drake)
library(tmap)
library(spData)
# library(spDataLarge)

Data

drake::loadd("wastd_data")
x <- wastd_data %>%
  wastdr::filter_wastd_turtledata(area_name = "Port Hedland")
site_cb <- x$sites %>% dplyr::filter(site_name == "Port Hedland Cemetery Beach")

From Tracks

nests_pth <- x$tracks %>%
  dplyr::filter(site_name == "Port Hedland Cemetery Beach") %>%
  dplyr::filter(nest_type == "successful-crawl") %>%
  wastdr::filter_realspecies()

nests_2021 <- nests_pth %>% wastdr::filter_wastd_season(2021)
nests_2020 <- nests_pth %>% wastdr::filter_wastd_season(2020)
nests_2019 <- nests_pth %>% wastdr::filter_wastd_season(2019)

nests_2021 %>% geojsonio::geojson_write(
  lat = "latitude", lon = "longitude", crs = 4326,
  file = here::here("inst/data/nests_2021.geojson")
)
nests_2020 %>% geojsonio::geojson_write(
  lat = "latitude", lon = "longitude", crs = 4326,
  file = here::here("inst/data/nests_2020.geojson")
)
nests_2019 %>% geojsonio::geojson_write(
  lat = "latitude", lon = "longitude", crs = 4326,
  file = here::here("inst/data/nests_2019.geojson")
)

nests_2021 %>% map_tracks(sites = site_cb)
nests_2020 %>% map_tracks(sites = site_cb)
nests_2019 %>% map_tracks(sites = site_cb)

# nests_2021_sf <- nests_2021 %>%
#   sf::st_as_sf(coords=c("longitude", "latitude"), crs = sf::st_crs(4326))
#
# tmap_mode("view")
# tm_basemap("Esri.WorldImagery") +
#   tm_shape(nests_2021_sf) +
#   tm_bubbles(size=0.5, col="red")

From tagging

tags_pth <- x$animals %>%
  dplyr::filter(site_name == "Port Hedland Cemetery Beach") %>%
  wastdr::filter_realspecies() %>%
  dplyr::filter(nesting_event %in% c("nest-unsure-of-eggs", "nest-with-eggs"))

tags_2021 <- tags_pth %>% wastdr::filter_wastd_season(2021)
tags_2020 <- tags_pth %>% wastdr::filter_wastd_season(2020)
tags_2019 <- tags_pth %>% wastdr::filter_wastd_season(2019)

tags_2021 %>% geojsonio::geojson_write(
  lat = "latitude", lon = "longitude", crs = 4326,
  file = here::here("inst/data/tags_2021.geojson")
)
tags_2020 %>% geojsonio::geojson_write(
  lat = "latitude", lon = "longitude", crs = 4326,
  file = here::here("inst/data/tags_2020.geojson")
)
tags_2019 %>% geojsonio::geojson_write(
  lat = "latitude", lon = "longitude", crs = 4326,
  file = here::here("inst/data/tags_2019.geojson")
)

tags_2021 %>% map_mwi(sites = site_cb)
tags_2020 %>% map_mwi(sites = site_cb)
tags_2019 %>% map_mwi(sites = site_cb)

dbca-wa/etlTurtleNesting documentation built on Nov. 18, 2022, 8:03 a.m.