计算栅格单元和邻域之间的差异

Question

我正在尝试计算栅格像元与 3 x 3 邻域中的平均值之间的差异。

有了这个光栅

r = raster(ncol=10, nrow=10, xmn=0, xmx=100, ymn=0, ymx=100)
set.seed(123)
values(r) = round(runif(ncell(r),1,100))

我可以使用类似于

的东西

r2 <- focal(r2, w=matrix(1,nrow=3,ncol=3), fun=)

但我认为我需要一个自定义函数作为 fun 的参数来计算 3 x 3 邻域的平均值，然后从中心单元格中减去该值。

我知道我可以使用两个不同的栅格图层来做到这一点，但我怀疑有一个更好的方法，如上所述。

任何建议将不胜感激。

Answer 1

您可以分两步完成：

#First, get the mean in a 3x3 neighbourhood
r2 <- focal(r, w=matrix(1,nrow=3,ncol=3), fun=mean)

#Then subtract the focal mean from the original cell values using simple raster arithmetic:
r3 <- r - r2

然后，您可以轻松地将上述步骤包装在一个更方便的函数中，干杯。