st_coordinates()
sf
inla()
ggplot() + geom_sf(data = map_geom) + geom_sf(data = centroids) + theme_bw()