st_coordinates()
sf
inla()
plot(map_geom, border = "lightgray") plot.nb(nblagsc, map_geom, add = TRUE)