小编Kie*_*fer的帖子

我有一个空间点数据框和一个空间多边形数据框.例如,我的多边形将是曼哈顿中每个块的多边形.分数是人,分散在各处,有时落在街道中间,不是多边形的一部分.

我知道如何检查多边形中是否包含一个点,但是如何将点分配给它们最接近的多边形？

## Make some example  data
set.seed(1)
library(raster)
library(rgdal)
library(rgeos)
p <- shapefile(system.file("external/lux.shp", package="raster"))
p2 <- as(1.5*extent(p), "SpatialPolygons")
proj4string(p2) <- proj4string(p)
pts <- spsample(p2-p, n=10, type="random")

## Plot to visualize
plot(pts, pch=16, cex=.5,col="red")
plot(p, col=colorRampPalette(blues9)(12), add=TRUE)