在一個點（光柵）距離內獲取細胞數

Question

我需要使用R獲得一點距離（例如10公里）內的細胞數量，但我沒有弄清楚如何處理它柵格數據。

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

如何獲得單點距離內的單元格數？

Answer 1

與緩衝選項和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)