在栅格中采样非相邻像元-R

Question

我想使用栅格像元值作为包含概率从栅格中采样非相邻像元。

下面是从栅格中采样像元的示例，这可能会导致选择相邻像元。

set.seed(1)
library(raster)
r <- raster(matrix(runif(2500), nrow = 50, ncol = 50))
plot(r)

cells = sample(1:ncell(r), 100, prob=r[], replace=FALSE)

cell.centers = xyFromCell(r, cells)

plot(r)
points(cell.centers, pch = 19)

我意识到，对于大样本量，避免选择相邻单元几乎是不可能的。实际上，我将要制作的样本数量总是比可用单元格的数量少得多。我只想避免包含相邻的单元格。我想我可以迭代采样并在选定单元格后将其删除..但这似乎很钝。任何帮助表示赞赏。