如何查找多個柵格的像素最大值？

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。）