围绕坐标创建圆形区域

Question

I need to create circular areas around a coordinate to identify points inside this circular area. Any ideas how can I do this the easiest way in R?

Answer 1

You can do this with the sf package. Will need to bring in your data, this differs if you have a spatail data file, or are getting your data from tables. Then buffer your co-ordinate to get your circular polygon, then intersect that against your points ot get only those that overlap.

library(sf)
library(tidyverse)

##if from table data###
CO_ORD <- st_as_sf(x =  "co-ord-data-file.csv", coords=c("Long", "Lat"))

POINTS <- st_as_sf(x = "points-data-file.csv", coords=c("Long", "Lat"))


###if from spatail data file####

CO_ORD <- st_read(dsn="co-ord-spatial-data.gdb", layer="co-ord-layer")

POINTS <- st_read(dsn="poins-spatial-data.gdb", layer="points-layer")



CO_ORD_BUFF <- st_buffer(x = CO_ORD, dist=500)

OVERLAPPING_POINTS <- st_intersection(POINTS, CO_ORD_BUFF)