# Wczytanie danych
base_sample <- sf::read_sf(dsn = "base_sample.gpkg")
random_sample <-  sf::read_sf(dsn = "random_sample.gpkg")
pref_sample <-  sf::read_sf(dsn = "pref_sample.gpkg")
decl_pref <- declusteringr::declstr_random(pref_sample, 100000, 2)

set.seed(75757561)

# Semiwariogram - losowe
vario_random <-  gstat::variogram(values~ 1, locations = random_sample)
model_random <-  gstat::vgm(psill = 10.3, model = "Sph", range = 295000)

# Semiwariogram - preferencyjne
vario_pref <-  gstat::variogram(values~ 1, locations = pref_sample)
model_pref <-  gstat::vgm(psill = 8, model = "Sph", range = 280000)
fitted_pref <-  gstat::fit.variogram(vario_pref, model_pref)

# Semiwariogram - preferencyjne po rozgrupowaniu losowym
vario_pref_rand <-  gstat::variogram(values~ 1, locations = decl_pref)
model_pref_rand <-  gstat::vgm(nugget = 0.25, psill = 12, model = "Sph", range = 340000)
fitted_pref_rand <-  gstat::fit.variogram(vario_pref_rand, model_pref_rand)