st_coordinates()
sf
inla()
poly2nb(<polygons>, queen = TRUE)
poly2nb(<polygons>, queen = FALSE)
dnearneigh(<points>, d1, d2)
knn2nb(<matrix of nearest neighbours>)