在R中的栅格中计算对数函数

Question

如何在R中的栅格中计算对数函数？

类似下面的代码片段     t = log(3-raster1/raster2-raster1,rasterBase) 其中rasterBase是用作对数底数的栅格

Answer 1

您应该始终提供一些示例数据/代码：

library(raster)
r <- raster(ncol=10, nrow=10)
r1 <- setValues(r, 1000)
r2 <- setValues(r,  1:100)

解决方案：

x <- overlay(r1, r2, fun=function(x, y) log(x, y))
x

替代方法：

s <- stack(r1, r2)
y <- calc(s, fun=function(x) log(x[1], x[2]))

出于完整性考虑，请注意，对于固定基数，例如10，你可以做

x <- log(r1, 10)