将栅格值提取到R中的多边形时包括多边形ID

Question

我遵循了如何从多边形数据中提取栅格值，然后合并到空间数据框中？ （这很有帮助）创建平均栅格值到多边形的矩阵（然后是数据框）。 现在的问题是我想知道哪个多边形。 我的SpatialPolygonsDataFrame在p$Block_ID具有ID值。 有没有办法把它带到extract()代码中？

或者， extract()函数是否按照输入顺序报告输出（这会有意义）？ 即p$Block_ID的顺序将保留在输出中？ 我浏览了文档，但不清楚一种方法还是另一种方法。 如果是这样，将ID列添加到extract()输出就很容易了。

这是我的通用代码供参考。 NOTE note具有可重现性，因为我认为此时并不需要。 其中r是栅格，多边形中的p

extract(r, p, small = TRUE, fun = mean, na.rm = TRUE, df = TRUE, nl = 1)

思考？

Answer 1

值将按顺序返回，就像在R中所期望的那样，并按手册中的说明（ ?extract ）：*返回值的顺序对应于对象y的顺序*

这样你就可以做

e <- extract(r, p)
ee <- data.frame(ID=p$Block_ID, e)

您应该提供了一个可重现的示例，并且很容易制作一个示例，例如，将其基于?extract的示例。

Answer 2

我无法得到R. Hijmans的答案为我工作。 我发现这可行。

e = extract(r, p)
e$ID = as.factor(e$ID)
levels(e$ID) = levels(p$Block_ID)