如何过滤SpatialPolygonDataFrame的部分内容?

时间:2016-02-16 16:29:20

标签: r rgdal sp

我想如果有一个简单的解决方案来过滤R中的SpatialPolygonDataFrame。想象一下,我有很多多边形,但我只想选择其中一些来使用leaflet来绘制它们

我的数据来自data.gouv.fr

这是我的gist,用于将数据加载到R中。

如果我只想在leaflet中绘制一个多边形,我会过滤SpatialPolygonDataFrame的@data部分,获取ID

> ign_shape_iris@data %>% 
+   filter(DCOMIRIS == "606120301")
 DEPCOM NOM_COM IRIS  DCOMIRIS               NOM_IRIS TYP_IRIS  id 
1  60612  Senlis 0301 606120301 Vald'Aunette-Gateliere        H 790

然后我只绘制我想要的多边形:

leaflet(ign_shape_iris@polygons[[790]]) %>% 
 addTiles() %>%
 addPolygons()

我确信有更好的解决方案。

1 个答案:

答案 0 :(得分:1)

尝试将其视为通常的data.frame,并考虑子集而不是过滤器:

ign_shape_iris[ign_shape_iris$DCOMIRIS == "606120301",]

subset(ign_shape_iris, DCOMIRIS == "606120301")