── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
✖ tidyr::extract() masks raster::extract()
✖ dplyr::filter() masks stats::filter()
✖ dplyr::lag() masks stats::lag()
✖ dplyr::select() masks raster::select()
ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errors
Load the bounding box polygon and create a bounding box.
#Loading bounding box for the area of interestfil <- here::here("data", "region", "BoundingBox.shp")extent_polygon <- sf::read_sf(fil)bbox <- sf::st_bbox(extent_polygon)wkt_geometry <- extent_polygon$geometry %>%st_as_text()
Make a map of our region so we know we have the right area.
world <-ne_countries(scale ="medium", returnclass ="sf")ggplot(data = world) +geom_sf() +geom_sf(data = extent_polygon, color ="red", fill=NA)
Deleting layer `pts_absence' using driver `CSV'
Writing layer `pts_absence' to data source
`/Users/eli.holmes/Documents/GitHub/ohw23_proj_marinesdms/data/raw-bio/pts_absence.csv' using driver `CSV'
options: GEOMETRY=AS_XY
Updating existing layer pts_absence
Writing 1000 features with 0 fields and geometry type Point.