st_coordinates()
sf
inla()
par(mfrow = c(1, 2), mar = c(2, 2, 2, 2)) plot(variog_cloud) plot(variog_bin)