R/get-cnc-obs.R

In iozeroff/cncpointR: What the Package Does (One Line, Title Case)

library(units)

cnc_area <- st_read("data/inputs/CNC_boundary_NEA_Edit_4.9.2019.kml") # Reads in cnc area shapefile as sf object. 

buffer <- set_units(30, m)

# Filters inat for all observations within cnc_area. 
inat_sf <- inat_dat %>% st_as_sf(coords = c("decimalLongitude", "decimalLatitude"), crs = st_crs(cnc_area))
cnc_buffered <- cnc_area %>% 
  st_transform(crs = 3857) %>% st_buffer(dist = buffer)
inat_sf_met <- inat_sf %>% st_transform(crs = 3857)
test_mat <- st_within(inat_sf_met, cnc_buffered, sparse = FALSE)

# Filters lvl2 according to test matrix.  
inat_cnc <- inat_dat[test_mat[, 1] == 1,]
saveRDS(inat_cnc, "data/inat_cnc.rds")