st_coordinates()
sf
inla()
ggplot() + geom_spatraster(data = grid) + coord_sf(datum = st_crs(grid)) + facet_wrap(~lyr) + scale_fill_viridis_c(na.value = "transparent", direction = -1) + theme_minimal()