存储大量栅格的栅格名称和表示

Question

对于非常多的栅格，我需要在栅格上循环，并获得栅格的名称和栅格的平均单元值。

这段代码工作正常，但只存储循环中最后一个元素的结果，当然我需要为每个i (=for each raster)存储它。

 for(i in ras){
      dfff$A <- cellStats(i,"mean")
      dfff$B <- i 
    }

我需要的是一个代码来存储每个栅格的结果，包括名称和均值(这样就可以很容易地看到哪个均值对应于哪个栅格)。这个代码会导致致命的错误。

   for(i in ras){
      dfff[i]$A <- cellStats(i,"mean")
      dfff[i]$B <- i 
    }

有谁知道怎么解决这个问题吗？

Answer 1

如果栅格具有相同的原点和范围，则只需执行以下操作

library(raster)
s <- stack(ras)
cellStats(s)

否则，您可以在循环中执行cellStats

d <- rep(NA, length(ras))
nms <- rep(NA, length(ras))
for(i in 1:length(ras)){
     d[i] <- cellStats(ras[[i]], "mean")
     nms[i] <- names(ras[[i]])
}

我不知道你说的栅格的“名字”是什么意思。在您的代码中，您将整个对象分配给"B“。但无论哪种方式，您都不需要在循环中这样做，因为您已经有了一个栅格列表(我假设这就是ras )。所以你也可以这样做

nms <- sapply(ras, names)

或者你指的是文件名？

nms <- sapply(ras, filename)

现在你可以做

df <- data.frame(name=nms, mean=d)

第二个循环的问题是i不是一个数字索引，它是一个光栅对象。所以你不能做diff[i]