如何解决在R中的缓冲区内提取点ID
我有一个具有纬度/经度坐标和ID的文件。我试图在每个点周围创建一个1km的缓冲区，然后获取每个缓冲区内所有点的ID。我不能做后者。这是我到目前为止所拥有的：
library(maps)
library(mapdata)
library(sp)
library(raster)
map('world2Hires',"Canada")
Geocoded_Data <- read.csv("filepath/Geocoded_Data.csv",header = TRUE)
# Lay the points over the map:
points(Geocoded_Data$long,Geocoded_Data$lat,col = "red",cex = .1)
# Remove NA values
Geocoded_Data <- Geocoded_Data[complete.cases(data),]
# Make a spatial data frame of locations
coords<-data.frame(x = Geocoded_Data$long,y= Geocoded_Data$lat)
crs<-"+proj=longlat +datum=wgs84 +no_defs +ellps=wgs84 +towgs84=0,0"
neighborhood.spdf <- SpatialPointsDataFrame(coords=coords,data=Geocoded_Data,proj4string=CRS(crs))
# Create 1km buffer and plot **strong text**it
buffer_1km <- buffer(neighborhood.spdf,width=1000)
plot(buffer_1km,add=TRUE,col='light blue',border='light blue')
# How do I Now extract the IDs of the points inside each buffer?
版权声明：本文内容由互联网用户自发贡献，该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容， 请发送邮件至 dio@foxmail.com 举报，一经查实，本站将立刻删除。
