st_coordinates()
sf
inla()
Combining interpolation method 1, 2, …, j, …, M:
ˆZ(s0)=M∑j=1(ˆZj(s0)⋅wj) With wj the weight for each method; ∑Mj=1wj=1