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) 

을