## Necessary pre-specified arguments required for the function to work.
## image resoluation + number of spatial predictors
im.res <- c(5, 5)
J <- prod(im.res)
## create predictor indices w/ coordinates
row.id <-rep(1, im.res[2])
for (i in 2:im.res[1]) {
row.id <- c(row.id, rep(i, im.res[2]))
}
coords <- data.frame(index = 1:J,
row.id = row.id,
col.id = rep(c(1:im.res[2]), im.res[1]) )
neighbors_by_dist(x = 2, y = 2, im.res = im.res, coords = coords, r = 2)
Run the code above in your browser using DataLab