r - sp::over() 用于多边形分析中的点

Question

我有一个名为“ind_adm”的 shapefile 和一个名为“pnts”的 SpatialPointsDataFrame。“pnts”包含随机生成的点，其中一些点与多边形重叠。见下图。在此处输入图像描述

现在，我想在多边形分析中做一个点，即我想找出哪些点位于代表印度边界的灰色多边形内。为此，我使用 sp 库中的 over() 函数。

pt.in.poly <- sp::over(ind_adm, pnts, fn = mean) #do the join

但是，我得到的输出是

    >pt.in.poly
    values
    0 6.019467

我实际上应该得到“在”多边形中的点的索引。

我哪里错了？

score 28 · Accepted Answer

28

为 over 找到了这种简洁直观的语法：

   pnts[ind_adm,]

从这个介绍文档

于 2016-12-02T04:57:30.747 回答

score 26 · Accepted Answer

您不应该提供函数。您正在聚合多边形几何上的点的属性值（即返回的数字是mean落在多边形内的点的属性值）。此外，对于您想做的事情，您有自己的错误方式x。y应该...

over( pnts , ind_adm , fn = NULL)

score 8 · Accepted Answer

您可以使用point.in.polyfomspatialEco包。它“与点和多边形要素类相交并将多边形属性添加到点”。

library(spatialEco)

new_shape <- point.in.poly(pnts, ind_adm)

score 1 · Accepted Answer

您还可以使用包中的st_intersection功能sf：

library(sf)

ind_adm <- st_as_sf(ind_adm)

（24047 是印度的 EPSG 代码）

pnts <- st_as_sf(pnts) %>% st_set_crs(., 24047)

kept_points <- st_intersection(ind_adm, pnts)

4 回答 4