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

【问题标题】：Calculate logarithm function in a raster in R在R中的栅格中计算对数函数
【发布时间】：2018-09-03 21:31:17
【问题描述】：

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

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

【问题讨论】：

【解决方案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)

【讨论】：