st_coordinates()
sf
inla()
ˆZ(s0)=k∑i=1(Z(si)⋅wi)
wi=d0i/k∑i=1d0i=1/k
nn <- gstat( formula = vble ~ 1, data = d2, locations = ~ x + y, nmax = 5, set = list(idp = 0) )