st_coordinates()
sf
inla()
The implementation just needs some simple raster algebra in {terra}.
{terra}
# using equal weights: r_ens <- mean(r_closest, r_idw, r_nn)