按坐标子集空间多边形 df

Posted 2023-02-23

【中文标题】按坐标子集空间多边形 df

【英文标题】：subset spatial polygon df by coordinates

【发布时间】：2017-10-26 18:13:48

【问题描述】：

我正在使用 spdf 构建带有传单和闪亮的交互式地图。我想根据当前视图中地图边界的坐标对我的数据进行反应性子集化（以便绘制相应变化的图形）。

可以通过数据属性对 spdf 进行子集化，例如

spdf2 <- spdf1[spdf1$Indicator == "pop" & spdf1$Country == "US",]

但是是否可以按多边形属性进行子集化 - 特别是 labpt 槽（多边形的中心）的经纬度坐标？

【参考方案1】：

由于我只对数据槽感兴趣，我的解决方案是使用 apply 提取 labpt 坐标，将生成的矩阵绑定到 sp 数据部分。例如：

my_df2 <- cbind(my_spdf@data, as.data.frame(t(sapply(slot(my_spdf, 'polygons'), 
                                                    function(i) slot(i, 'labpt')))))

然后使用传统的子集/过滤器方法。