我有数百个具有相同分辨率和范围的栅格。这是一个时间序列，每个栅格代表一个时间点。
我知道如何在栅格中找到绝对最大值。
但是如何在整个栅格时间序列中找到每个单元格中的最大值？
如果 a1,a2,......a1000 是栅格，我想创建一个栅格 x，其中每个像素是 a1-a1000 的所有对应像素的最大值。
【问题讨论】：
?pmax。试试pmax(a1[],a2[],...)，其中a1、a2 等是您的raster 对象。
如果您首先将栅格放入堆栈中，则可以简单地将 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)
（要应用一些other，更复杂的函数，它为堆栈中的每个像素返回一个标量，您可能需要使用calc()，如（例如）here 所示。）
【讨论】：
raster list 而不是raster stack？就我而言，我的栅格没有完全相同的范围，我想避免raster::resample，因为处理时间显着增加