st_coordinates()
sf
inla()
Basic conversion from the neighbours list to a neighbourhood matrix:
nb2mat(nb1, style = "B") |> dim()
## [1] 49 49