r - 如何找到多个栅格的像素最大值？

Question

我有数百个具有相同分辨率和范围的栅格。这是一个时间序列，每个栅格代表一个时间点。

我知道如何在栅格中找到绝对最大值。

但是如何在整个时间序列的栅格中找到每个单元格中的最大值？

如果 a1,a2,......a1000 是光栅，我想创建一个光栅 x，其中每个像素是 a1-a1000 所有对应像素的最大值。

Answer 1

如果您首先将栅格放入堆栈中，则可以简单地将 min() 或 max() 应用于堆栈以获取摘要 RasterLayer 你在追求

## Example rasters and stack
r1 <- raster(matrix(1:4,ncol=4))
r2 <- -2*r1
r3 <- 2*r1
rr <- list(r1,r2,r3)
s <- stack(rr)

## Extract the pixel-wise min and max values
min(s)
max(s)

(要应用一些其他，为堆栈中的每个像素返回一个标量的更复杂的函数，您可能需要使用 calc()，如所演示的(对于例子) here .)