R - 結合索引從矢量和柵格堆棧層的函數

Question

我想寫一個高效的腳本來校準數百個Landsat 8圖像。在校準步驟的某個點，我需要在柵格堆棧的每一層中應用一些係數。

這是一個示例堆棧：

fn <- system.file("external/test.grd", package="raster") 
s <- stack(fn, fn) 

而且這是樣本系數：

mult <- c(0.0003342, 0.0005534) 
add <- c(0.1, 0.2) 

我需要什麼是係數的各指標適用於堆棧層的相應指數，就像這個例子：

s[[1]] <- (s[[1]] * mult[1]) + add[1] 
s[[2]] <- (s[[2]] * mult[2]) + add[2] 

這是我可憐的嘗試，這顯然不工作：

cal.fun <- function(x) { 
x <- (x * mult) + add 
} 

s.cal <- calc(s, cal.fun, progress='text') 

有關如何做到這一點的任何想法？

非常感謝。

Answer 1

光柵是一個驚人構建良好的包，你可以簡單地做：

s2 <- s * mult + add 

對於快速的視覺確認，該矢量稱之爲「只是工程」，做這樣的事情：

library(gridExtra) 
library(rasterVis) 
grid.arrange(levelplot(s), levelplot(s2), nrow=2)