local_R/test_NW_monitors.R

In MazamaScience/MazamaLocationUtils: Manage Spatial Metadata for Known Locations

library(PWFSLSmoke)
library(MazamaLocationUtils)

mazama_initialize()

jon <- table_initialize()

# ----- Add Washington

monitor <- monitor_loadLatest() %>% monitor_subset(stateCodes = "WA")
lons <- monitor$meta$longitude
lats <- monitor$meta$latitude

jon <- jon %>% table_addLocation(lons, lats, distanceThreshold = 500)

table_save(jon, "jon")

# ----- Add Oregon

monitor <- monitor_loadLatest() %>% monitor_subset(stateCodes = c("WA","OR"))
lons <- monitor$meta$longitude
lats <- monitor$meta$latitude

jon <- jon %>% table_addLocation(lons, lats, distanceThreshold = 500)

table_save(jon, "jon")

# ----- Add Pac NW random sample

monitor <- monitor_loadLatest() %>% monitor_subset(stateCodes = c("WA","OR","ID"))
indices <- sort(sample(seq_along(monitor$meta$longitude), 10))
lons <- monitor$meta$longitude[indices]
lats <- monitor$meta$latitude[indices]

jon <- jon %>% table_addLocation(lons, lats, distanceThreshold = 500)

table_save(jon, "jon")

# ----- Add Idaho

monitor <- monitor_loadLatest() %>% monitor_subset(stateCodes = c("WA","OR","ID"))
lons <- monitor$meta$longitude
lats <- monitor$meta$latitude

jon <- jon %>% table_addLocation(lons, lats, distanceThreshold = 500)

table_save(jon, "jon")

# ----- Leaflet

table_leaflet(jon, radius = 5)

# ----- Map

maps::map("state", c("idaho", "oregon", "washington"))
points(jon$longitude, jon$latitude)