library(sf)
# get the projection of the epmGrid object
proj <- summary(tamiasEPM)$crs
# define some points
pts <- rbind.data.frame(
c(-120.5, 38.82),
c(-84.02, 42.75),
c(-117.95, 55.53))
colnames(pts) <- c('x', 'y')
ptsSF <- st_as_sf(pts, coords = 1:2, crs = "epsg:4326")
pts <- st_coordinates(st_transform(ptsSF, crs = proj))
# extract with table of coordinates
extractFromEpmGrid(tamiasEPM, pts)
# extract with spatial points object
extractFromEpmGrid(tamiasEPM, ptsSF)
# extract with spatial polygon
hull <- st_convex_hull(st_union(ptsSF))
extractFromEpmGrid(tamiasEPM, hull)
# returns each cell's contents
extractFromEpmGrid(tamiasEPM, hull, collapse=FALSE)
# collapses results to unique set of species
extractFromEpmGrid(tamiasEPM, hull, collapse=TRUE)
Run the code above in your browser using DataLab