获取一个点（栅格）距离内的单元格数

Question

我需要使用R获取一个点的距离（例如，10 km）内的单元格数，但我还没有弄清楚如何处理栅格数据。

library(raster)
r <- raster(ncols=10, nrows=10)
cellFromXY(r, c(2,2)) # get the cell number of one point

如何在一个点的距离内获取细胞数？

Answer 1

使用extract和buffer选项，cellnumbers = TRUE

r <- raster(ncols=100, nrows=100)
  r[]<-runif(ncell(r))

xy <- cbind(-50, seq(-80, 80, by=20))
extract(r, xy[1:3,], buffer=1000000,cellnumbers=TRUE)