R提取不同空间多边形属性的栅格和直方图

Question

我想基于不同的属性组合为多边形创建栅格值的直方图。以下可重复数据：

library(raster)
library(sp)   
poly <- raster(nrow=10, ncol=10)
poly[] <- runif(ncell(poly)) * 10
poly <- rasterToPolygons(poly, fun=function(x){x > 9})
r <- raster(nrow=100, ncol=100)
r[] <- runif(ncell(r))
poly@data$place<-sample(letters[1:3], length(poly), TRUE)
poly@data$rank<-sample.int(3, length(poly), replace = TRUE)
plot(r)
plot(poly, add=TRUE, lwd=4) 
v <- raster::extract(r, poly, df=TRUE)

我可以用v绘制ggplot中所有ID（即多边形）的直方图

ggplot(v, aes(layer)) + geom_histogram(aes(y = stat(count / sum(count))), binwidth = 0.25)

但是，我想基于rank属性（即1,2,3）和基于place属性的另一组三个直方图（即a，b，c）创建一组三个直方图。也许在ggplot中使用facet，但我不确定如何将v中的ID链接到poly中的属性。

Answer 1

你的例子：

library(raster)
#Loading required package: sp
pr <- raster(nrow=10, ncol=10)
set.seed(1)
values(pr) <- runif(ncell(pr)) * 10
poly <- rasterToPolygons(pr, fun=function(x){x > 9})
poly$place <- sample(letters[1:3], length(poly), TRUE)
poly$rank <- sample.int(3, length(poly), replace = TRUE)

r <- raster(nrow=100, ncol=100)
values(r) <- runif(ncell(r))
v <- raster::extract(r, poly, df=TRUE)

为多边形分配显式ID，仅保留感兴趣的变量，并从SpatialPolygonsDataFrame中提取data.frame。

poly$ID <- 1:length(poly)
poly$layer <- NULL
d <- data.frame(poly)

去

vd <- merge(d, v, by="ID")

选择一个子集并制作直方图

x <- vd[vd$place == "a",]
hist(x$layer)