时间戳

首页 > 解决方案 > 基于类的加权栅格聚合

r - 基于类的加权栅格聚合

问题描述

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

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

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

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

谢谢你的帮助？

标签： raggregateweighted

解决方案

您可以使用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)

推荐阅读