st_coordinates()
sf
inla()
To predict values of a spatially continuous variable at unsampled locations s0: ˆZ(s0)