Running examples in 'ClimaRep-Ex.R' failed
The error most likely occurred in:

> ### Name: mh_rep
> ### Title: Multivariate climate representativeness analysis
> ### Aliases: mh_rep
> 
> ### ** Examples
> 
> library(terra)
terra 1.8.86
> library(sf)
Linking to GEOS 3.13.0, GDAL 3.10.1, PROJ 9.5.1; sf_use_s2() is TRUE
> set.seed(2458)
> n_cells <- 100 * 100
> r_clim_present <- terra::rast(ncols = 100, nrows = 100, nlyrs = 7)
> values(r_clim_present) <- c(
+    (terra::rowFromCell(r_clim_present, 1:n_cells) * 0.2 + rnorm(n_cells, 0, 3)),
+    (terra::rowFromCell(r_clim_present, 1:n_cells) * 0.9 + rnorm(n_cells, 0, 0.2)),
+    (terra::colFromCell(r_clim_present, 1:n_cells) * 0.15 + rnorm(n_cells, 0, 2.5)),
+    (terra::colFromCell(r_clim_present, 1:n_cells) +
+      (terra::rowFromCell(r_clim_present, 1:n_cells)) * 0.1 + rnorm(n_cells, 0, 4)),
+    (terra::colFromCell(r_clim_present, 1:n_cells) /
+      (terra::rowFromCell(r_clim_present, 1:n_cells)) * 0.1 + rnorm(n_cells, 0, 4)),
+    (terra::colFromCell(r_clim_present, 1:n_cells) *
+      (terra::rowFromCell(r_clim_present, 1:n_cells) + 0.1 + rnorm(n_cells, 0, 4))),
+    (terra::colFromCell(r_clim_present, 1:n_cells) *
+      (terra::colFromCell(r_clim_present, 1:n_cells) + 0.1 + rnorm(n_cells, 0, 4)))
+ )
> names(r_clim_present) <- c("varA", "varB", "varC", "varD", "varE", "varF", "varG")
> terra::crs(r_clim_present) <- "EPSG:4326"
Warning: PROJ: proj_create_from_database: Cannot find proj.db (GDAL error 1)
Warning: [crs<-] Cannot set raster SRS: empty srs
> 
> vif_result <- ClimaRep::vif_filter(r_clim_present, th = 5)
Starting iterative VIF filtering process
 - Removing variable with highest VIF ( 15.55 ): varD
 - Removing variable with highest VIF ( 7.07 ): varB
VIF filtering process completed
> print(vif_result$summary)
$kept_layers
[1] "varA" "varC" "varE" "varF" "varG"

$excluded_layers
[1] "varD" "varB"

$original_correlation_matrix
        varA    varB   varC   varD    varE    varF   varG
varA  1.0000  0.8870 0.0087 0.0938 -0.0745  0.5816 0.0028
varB  0.8870  1.0000 0.0043 0.0997 -0.0769  0.6513 0.0003
varC  0.0087  0.0043 1.0000 0.8523  0.0517  0.5676 0.8341
varD  0.0938  0.0997 0.8523 1.0000  0.0576  0.7077 0.9529
varE -0.0745 -0.0769 0.0517 0.0576  1.0000 -0.0256 0.0653
varF  0.5816  0.6513 0.5676 0.7077 -0.0256  1.0000 0.6305
varG  0.0028  0.0003 0.8341 0.9529  0.0653  0.6305 1.0000

$final_vif_values
  varA   varC   varE   varF   varG 
2.2832 3.3473 1.0121 3.8325 4.3545 

> r_clim_present_filtered <- vif_result$filtered_raster
> hex_grid <- sf::st_sf(
+    sf::st_make_grid(
+      sf::st_as_sf(
+        terra::as.polygons(
+          terra::ext(r_clim_present_filtered))),
+      square = FALSE)
+      )
> sf::st_crs(hex_grid) <- "EPSG:4326"
> polygons <- hex_grid[sample(nrow(hex_grid), 2), ]
> polygons$name <- c("Pol_A", "Pol_B")
> study_area_polygon <- sf::st_as_sf(terra::as.polygons(terra::ext(r_clim_present_filtered)))
> sf::st_crs(study_area_polygon) <- "EPSG:4326"
> terra::plot(r_clim_present_filtered[[1]])
> terra::plot(polygons, add = TRUE, color = "transparent", lwd = 3)
> terra::plot(study_area_polygon, add = TRUE, col = "transparent", lwd = 3, border = "red")
> 
> ClimaRep::mh_rep(
+    polygon = polygons,
+    col_name = "name",
+    climate_variables = r_clim_present_filtered,
+    study_area = study_area_polygon,
+    th = 0.95,
+    dir_output = file.path(tempdir(), "ClimaRep"),
+    save_raw = TRUE
+    )
Establishing output file structure
Validating and adjusting Coordinate Reference Systems (CRS)
Adjusting CRS of climate_variables to match the polygon's system (EPSG:4326).
Warning: PROJ: proj_create: no database context specified (GDAL error 1)
Warning: PROJ: proj_create_from_database: Cannot find proj.db (GDAL error 1)
Warning: PROJ: proj_create: no database context specified (GDAL error 1)
Error: [project] cannot get output boundaries for the target crs
Execution halted