我有一个由值1、2和NA组成的栅格.我想使用带有3x3窗口的焦点函数来重新分配具有新值的单元格.附图显示了数据集和3x3窗口的示例.我希望新分配的值是窗口中不等于中心单元格的相邻单元格的计数，不包括NA个值.因此，基本上，如果中心单元格=1，我希望它计算相邻单元格的个数=2，反之亦然.在这张附加的照片中，值将是4，因为中心值=1，并且有四个相邻的单元格=2.
我正在使用下面的代码，它与基本功能(mean,min,max等)一起工作，但没有任何功能来做我正在寻找的，就我所能找到的.
focal.statistic <- focal(raster, window=window, function=function)个
我需要创建一个执行上述操作的新函数，然后可以使用"Function="代码调用该函数.我一直在考虑IF，THEN语句，但我正在努力编写这个函数……
这和很多年前的this post差不多，但我还是想不出来.
这里是用于对不等于中心的窗口的相邻单元进行计数的可能函数(细节如下).
f = function (x) sum(x != x[5L] & x != -999, na.rm = TRUE)
您提供的示例:
library(raster)
x <- matrix(c(
c(1, 1, NA, NA, 1, 2),
c(2, 2, 2, 1, 1, NA),
c(2, 2, 1, 1, 2, NA),
c(NA, NA, 2, 1, 2, NA),
c(NA, NA, 2, 2, 2, NA),
c(NA, NA, 1, 1, 2, NA)
),
nrow = 6)
r <- raster(x)
res <- as.matrix(focal(
x = r,
w = matrix(1, 3, 3),
fun = f,
pad = T,
padValue = -999
))
> x
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 1 2 2 NA NA NA
[2,] 1 2 2 NA NA NA
[3,] NA 2 1 2 2 1
[4,] NA 1 1 1 2 1
[5,] 1 1 2 2 2 2
[6,] 2 NA NA NA NA NA
> res
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 2 2 0 0 0 0
[2,] 3 3 1 0 0 0
[3,] 0 4 4 3 3 2
[4,] 0 2 4 6 3 4
[5,] 1 2 4 2 2 1
[6,] 2 0 0 0 0 0
解释:
padValue = -999呼叫focal.我使用PAD值是为了正确考虑矩阵的边缘.f()中，我们取和，一个逻辑向量，其元素表示各自的值是否不等于中心或不等于虚值.请注意，我们将窗口矩阵的中心称为x[5]，如果您采用另一个窗口大小，则必须对其进行调整(例如，对于5x5的窗口，调整x[13]).