st_coordinates()
sf
inla()
In the model formula, the SPDE model and the index set will be used in defining the random effect:
f(s, model = spde)