st_coordinates()
sf
inla()
closest observation method
ˆZ(s0)=Z(sclosest)
inverse distance weighting method (IDW)
ˆZ(s0)=∑ni=1(Z(si)⋅wi)