siatka <- read_stars("simulation.tif")
paleta <-  hcl.colors(12, palette = "ag_Sunset")

# Kriging prosty - losowe
sk_random <-  gstat::krige(values ~ 1,
                     locations = random_sample,
                     newdata = siatka,
                     model = model_random,
                     beta = mean(random_sample$values))
plot(sk_random["var1.pred"], col = paleta)

# Kriging prosty - preferencyjne
sk_pref <-  gstat::krige(values ~ 1,
                     locations = pref_sample,
                     newdata = siatka,
                     model = fitted_pref,
                     beta = mean(pref_sample$values))
plot(sk_pref["var1.pred"], col = paleta)

# Kriging prosty - preferencyjne po rozgrupowaniu losowym
sk_pref_rand <-  gstat::krige(values ~ 1,
                         locations = decl_pref,
                         newdata = siatka,
                         model = fitted_pref_rand,
                         beta = mean(decl_pref$values))
plot(sk_pref_rand["var1.pred"], col = paleta)

# Kriging prosty - dla próbkowania preferencyjnego po rozgrupowaniu losowym na podstawie modelu próbkowania preferencyjnego
sk_model_pref <- gstat::krige(values~1,
                              locations = pref_sample,
                              newdata = siatka,
                              model = fitted_pref,
                              beta = mean(decl_pref$values))
plot(sk_model_pref["var1.pred"], col = paleta)