在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)