空间连续性：编码Queen连续性为R中的虚拟变量

Question

我有空间时间序列数据，网格为基本单位。我现在需要为连续网格中的事件（由虚拟对象）的出现创建年度度量。

例如，在这里，根据女王连续性的规则，值为1/8 = 0,125。

0 0 1
0 x 0
0 0 0

我有ID变量，坐标以及数据中每个网格的行号和列号，但我不知道如何计算我的度量。

Answer 1

如果您正在处理栅格数据，则可以使用raster包计算焦点统计数据。

例如：

library(raster)

# Create fake binary data
r <- raster(matrix(rbinom(100, 1, 0.2), 10))

# Focal mean of the (binary) values of the 8-cell (fewer at grid's edges)
#  neighbourhood of each cell.
r_focal <- focal(r, matrix(c(1, 1, 1, 1, 0, 1, 1, 1, 1), 3), mean, 
                 pad=TRUE, na.rm=TRUE)

输入和输出如下所示：

plot(stack(r, r_focal), axes=FALSE, main='')

有关详细信息，请参阅?focal。