As spdep::poly2nb: queen neighbours share at least one boundary
vertex within snap in both coordinates, rook neighbours at least
two. Polygons are two-column coordinate matrices; a repeated closing
vertex is ignored.
Usage
PolygonContiguity(polygons, queen = TRUE, snap = sqrt(.Machine$double.eps))