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)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.