Вы всегда должны предоставлять пример данных / код:
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)