基于类的加权栅格聚合

Question

假设我有一个以特定分辨率表示土地利用类别的栅格。我必须将此栅格与 R 聚合为较粗的分辨率和模态值方法，以便在较粗的栅格中具有最主要的像元值。这很容易通过

实现

m <- aggregate(r, fact = 3, fun = modal, na.rm = TRUE)

但是，我想对不同的土地利用类别进行加权——例如森林类（代码 1）的权重为 4，水类（代码 2）的权重为 2，街道类的权重为 1。

是否有一个函数可以遍历栅格单元并为每个单元应用权重？

感谢您的帮助？

Answer 1

您可以使用reclassify 来应用权重，但是然后呢？之后你还想计算模态值吗？

我认为你想要的是你自己提供的聚合函数。也许是这样的

library(raster)
f <- function(x, ...) {
    y <- c(
       rep(x[x==1], 4),
       rep(x[x==2], 2),
       x[x==3]
    )
    modal(y, ...)
}

r <- raster(res=5)
values(r) <- sample(c(1:3,2,3,3), ncell(r), replace=TRUE)

a <- aggregate(r, fact=10, fun=f)