## ----setup--------------------------------------------------------------------
library(potentiomap)
data("synthetic_wells")
data("synthetic_candidate_sites")
p <- ps_make_points(synthetic_wells[1:12, ], "x", "y", "gw_elevation",
                    "well_id", "EPSG:26916")

## ----network------------------------------------------------------------------
thin <- ps_network_thinning(p, retain = 0.75, repeats = 2,
                            method = "IDW", grid_res = 400, seed = 3)
thin$summary
candidates <- terra::vect(subset(synthetic_candidate_sites, !excluded),
                          geom = c("x", "y"), crs = "EPSG:26916")
design <- ps_candidate_network(p, candidates, n_select = 2)
design$selected_sequence