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