Neighbours list based on distance bounds

ggplot() + 
  geom_sf(data = map_geom) +
  geom_sf(data = centroids) +
  theme_bw()