## ----setup--------------------------------------------------------------------
library(potentiomap)
data("synthetic_wells")
p <- ps_make_points(synthetic_wells[1:16, ], "x", "y", "gw_elevation",
                    "well_id", "EPSG:26916")
fit <- suppressWarnings(ps_interpolate(p, methods = "OK", grid_res = 350,
                                      support = TRUE, return = "result"))

## ----uncertainty--------------------------------------------------------------
fit$support$summary
u <- ps_surface_uncertainty(fit, approach = "kriging_variance")
u$method_manifest
band <- ps_contour_uncertainty(u, 168, method = "gaussian_pointwise",
                               accept_gaussian = TRUE)
band$level_manifest