R/spatial.R

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)
}
moecampos/425-project documentation built on May 14, 2019, 2:03 a.m.