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

Question

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

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

底数的栅格

Answer 1

你应该总是提供一些例子data/code:

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)

Calculate logarithm function in a raster in R