From 2f328751a98b1e7a0535646f34dd320128d0b28a Mon Sep 17 00:00:00 2001 From: Ward Langeraert Date: Thu, 19 Dec 2024 09:16:04 +0100 Subject: [PATCH] cache --- .../markdown/test_densiteitsmodellering.Rmd | 21 ++++++++++++------- 1 file changed, 14 insertions(+), 7 deletions(-) diff --git a/source/markdown/test_densiteitsmodellering.Rmd b/source/markdown/test_densiteitsmodellering.Rmd index 77696ec..c79e3bb 100644 --- a/source/markdown/test_densiteitsmodellering.Rmd +++ b/source/markdown/test_densiteitsmodellering.Rmd @@ -2099,14 +2099,21 @@ prediction_grid <- st_make_grid( cellsize = units::set_units(300*300*pi, "m^2") ) prediction_grid_sf <- st_sf(geometry = prediction_grid) -joined_grid <- st_join(prediction_grid_sf, dsm1_pred_sf, join = st_nearest_feature) +joined_grid <- st_join(prediction_grid_sf, dsm1_pred_sf, + join = st_nearest_feature) -cropped_grid <- qgisprocess::qgis_run_algorithm( - "native:intersection", - INPUT = joined_grid, - OVERLAY = flanders_sf - ) %>% - st_as_sf() +if (file.exists(file.path(cache_dir, paste0("cropped_grid_dsm.Rds")))) { + cropped_grid <- readRDS(file.path(cache_dir, paste0("cropped_grid_dsm.Rds"))) +} else { + cropped_grid <- qgisprocess::qgis_run_algorithm( + "native:intersection", + INPUT = joined_grid, + OVERLAY = flanders_sf + ) %>% + st_as_sf() + + saveRDS(cropped_grid, file.path(cache_dir, paste0("cropped_grid_dsm.Rds"))) +} ```