我有一个包含120个多边形和一些相关数据的SpatialPolygonsDataFrame。现在,我想分别提取每个多边形内栅格上的值的平均值。我成功地绘制了单个多边形:
plot(SpatialPolygons(SPdataframe@polygons)[i])
但它不能以相同的方式提取值:
extract(raster, SpatialPolygons(SPdataframe@polygons)[i],fun="mean",na.rm=TRUE,method="simple")
有谁可以解释在这两种情况下使用相同索引的区别?使用索引选择SpatialPolygonsDataFrame的特定多边形的官方方法是什么?
非常感谢您的帮助!
答案 0 :(得分:2)
SpatialPolygonsDataFrame的单个多边形的正确索引是:SPdataframe[i,]
(Merci来自R-sig_geos用户RafaelWüest)