st_coordinates()
sf
inla()
r_closest <- vect(d) |> voronoi() |> rasterize(grid, field = "vble") |> mask(map)