R/spatial.R

In moecampos/425-project: Venice Pirates (STAT 425 Project)

library(geosphere)

within_radius <- function(lon, lat, coordinates, radius = 1000) {
  point <- cbind(lon, lat)
  distances <- geosphere::distm(point, coordinates, 
                                fun = geosphere::distVincentyEllipsoid)
  which(distances <= radius)
}