Rename project to reflect all landbouwstreken

source/targets/mas_steekproef_zl_lm_za_kp_po/.Rprofile

@@ -0,0 +1,2 @@
+source("renv/activate.R")
+library(sf)

source/targets/mas_steekproef_zl_lm_za_kp_po/_targets.R

@@ -0,0 +1,32 @@
+# Generated by targets: do not edit by hand
+library(targets)
+lapply(
+  X = list.files(
+    path = file.path(targets::tar_path_script_support(), "globals"),
+    pattern = "\\.R$",
+    full.names = TRUE,
+    all.files = TRUE,
+    recursive = TRUE
+  ),
+  FUN = function(path) {
+    eval(
+      expr = parse(file = path, keep.source = TRUE),
+      envir = targets::tar_option_get(name = "envir")
+    )
+  }
+)
+lapply(
+  X = list.files(
+    path = file.path(targets::tar_path_script_support(), "targets"),
+    pattern = "\\.R$",
+    full.names = TRUE,
+    all.files = TRUE,
+    recursive = TRUE
+  ),
+  FUN = function(path) {
+    eval(
+      expr = parse(file = path, keep.source = TRUE),
+      envir = targets::tar_option_get(name = "envir")
+    )
+  }
+)

source/targets/mas_steekproef_zl_lm_za_kp_po/_targets/.gitignore

@@ -0,0 +1,5 @@
+*
+!.gitignore
+!meta
+meta/*
+!meta/meta

source/targets/mas_steekproef_zl_lm_za_kp_po/_targets_packages.R

@@ -0,0 +1,27 @@
+# Generated by targets::tar_renv(): do not edit by hand
+library(GVI)
+library(arrow)
+library(bs4Dash)
+library(clustermq)
+library(dplyr)
+library(exactextractr)
+library(future)
+library(gt)
+library(httr)
+library(markdown)
+library(nngeo)
+library(osmextract)
+library(pingr)
+library(qs)
+library(readr)
+library(rstudioapi)
+library(sf)
+library(shiny)
+library(shinyWidgets)
+library(shinybusy)
+library(terra)
+library(tibble)
+library(tidyr)
+library(visNetwork)
+library(SamplingBigData)
+library(duckdb)

source/targets/mas_steekproef_zl_lm_za_kp_po/_targets_r/globals/mas-globals.R

@@ -0,0 +1,31 @@
+options(tidyverse.quiet = TRUE)
+tar_option_set(packages = c("dplyr",
+                            "readr",
+                            "tidyr",
+                            "tibble",
+                            "nngeo",
+                            "osmextract",
+                            "terra",
+                            "exactextractr",
+                            "arrow",
+                            "httr",
+                            "sf",
+                            "GVI"),
+               format = "qs",
+               memory = "transient",
+               garbage_collection = TRUE,
+               workspace_on_error = TRUE)
+
+target_dir <- rprojroot::find_root_file(
+    "source", "targets", "mas_steekproef_zn_lm",
+    criterion = rprojroot::is_git_root)
+mbag_dir <- rprojroot::find_root_file(
+    criterion = rprojroot::is_git_root)
+
+source(file.path(mbag_dir, "source", "R", "landuse_maps.R"))
+source(file.path(mbag_dir, "source", "R", "geocomputations.R"))
+source(file.path(mbag_dir, "source", "R", "draw_sample.R"))
+source(file.path(mbag_dir, "source", "R", "steekproefkader.R"))
+source(file.path(mbag_dir, "source", "R", "berekening_hulpvariabelen.R"))
+source(file.path(mbag_dir, "source", "R", "steekproeftrekking_nabehandeling.R"))
+source(file.path(mbag_dir, "source", "R", "wfs_wcs.R"))

source/targets/mas_steekproef_zl_lm_za_kp_po/_targets_r/targets/hulpvariabelen.R

@@ -0,0 +1,44 @@
+list(
+  tar_target(
+    name = plus_sb,
+    command = add_bo_to_frame(
+      punten_df = punten_selectie_zichtbaarheid,
+      path_bo = bo_file
+      ),
+    pattern = map(punten_selectie_zichtbaarheid),
+    iteration = "list"
+  ),
+  tar_target(
+    name = plus_openheid_landschap,
+    command = add_openheid_landschap_to_frame(
+      path = openheid_landschap_file,
+      punten_sf = plus_sb,
+      gebied = perimeters_data,
+      cutlevels = c(1.25, 1.35, 1.51),
+      class_labels = c("GL", "HGL", "HOL", "OL")
+      ),
+    pattern = map(perimeters_data, plus_sb),
+    iteration = "list"
+  ),
+  tar_target(
+    name = sbp_akkervogels,
+    command = read_sbp_akkervogels(
+      path = sbp_akkervogels_file,
+      gebied = perimeters_data,
+      path_extra_soorten = sbp_overige_file,
+      extra_soorten = c("hamster", "bruine kiekendief", "zomertortel",
+                        "grauwe kiekendief")
+      ),
+    pattern = map(perimeters_data),
+    iteration = "list"
+  ),
+  tar_target(
+    name = steekproefkader_finaal,
+    command = add_stratum_sbp(
+      punten_sf = plus_openheid_landschap,
+      sbp = do.call(rbind.data.frame, sbp_akkervogels)
+      ),
+    pattern = map(plus_openheid_landschap),
+    iteration = "list"
+  )
+)

source/targets/mas_steekproef_zl_lm_za_kp_po/_targets_r/targets/nabehandeling.R

@@ -0,0 +1,22 @@
+list(
+  tar_target(
+    steekproef_thinned,
+    thin_sample(
+      sample = steekproef,
+      thin_dist = 600
+      ),
+    pattern = map(steekproef),
+    iteration = "list"
+  ),
+  tar_target(
+    steekproef_final,
+    replace_by_existing(
+      sample = steekproef_thinned,
+      existing_points = existing_data,
+      overlap_prop = 0.5,
+      sbp_file = do.call(rbind.data.frame, sbp_akkervogels)
+      ),
+    pattern = map(steekproef_thinned),
+    iteration = "list"
+  )
+)

source/targets/mas_steekproef_zl_lm_za_kp_po/_targets_r/targets/raw-data.R

@@ -0,0 +1,51 @@
+list(
+  tarchetypes::tar_file(
+    perimeters_file,
+    path_to_perimeters("zl_lm_za_kp_po_perimeters.gpkg")
+    ),
+  tar_target(
+    perimeters_data,
+    st_read(perimeters_file) %>%
+      arrange(Naam)
+    ),
+  tarchetypes::tar_file(
+    osm_belgium,
+    path_to_osm_download()
+    ),
+  tarchetypes::tar_file(
+    vito_lum_2019_file,
+    path_to_lum(jaar = 2019)
+    ),
+  tarchetypes::tar_file(
+    legend_lum_file,
+    path_to_legend_lum()
+    ),
+  tar_target(
+    legend_lum,
+    read_legend_lum(file = legend_lum_file)
+  ),
+  tarchetypes::tar_file(
+    openheid_landschap_file,
+    path_to_openheid_landschap()
+    ),
+  tarchetypes::tar_file(
+    bo_file,
+    path_to_bo(jaar = 2022)
+    ),
+  tarchetypes::tar_file(
+    sbp_akkervogels_file,
+    path_to_sbp_akkervogels(file = "akkervogelgebieden2022.shp")
+  ),
+  tarchetypes::tar_file(
+    sbp_overige_file,
+    path_to_sbp_akkervogels(file = "sbp_overige_soorten.shp")
+  ),
+  tarchetypes::tar_file(
+    existing_file,
+    path_to_existing(file = "steekproef_piloot_avimap.geojson")
+  ),
+  tar_target(
+    existing_data,
+    st_read(existing_file)
+  )
+)

source/targets/mas_steekproef_zl_lm_za_kp_po/_targets_r/targets/steekproefkader.R

@@ -0,0 +1,114 @@
+list(
+  tar_target(
+    name = selectie_openheid_klasses,
+    command = selectie_openheid(
+      gebied = perimeters_data,
+      ol_strata = c("OL", "HOL")
+      ),
+    pattern = map(perimeters_data),
+    iteration = "list"
+  ),
+  tar_target(
+    name = exclusie_osm_landgebruiken,
+    command = exclusie_landgebruik_osm(
+      gebied = selectie_openheid_klasses,
+      osmdata = osm_belgium,
+      landuse = c("residential", "military", "industrial", "cemetery",
+                  "railway", "commercial", "farmyard"),
+      leisure = c("park"),
+      buffer_poly = 0,
+      layer_poly = list(aeroway = c("aerodrome")),
+      buffer_line = 100,
+      layer_line = list(highway = c("motorway", "motorway_link")),
+      update_osm_layer = FALSE
+      ),
+    pattern = map(selectie_openheid_klasses),
+    iteration = "list"
+  ),
+  tar_target(
+    name = paden,
+    command = extract_osm_paden(
+      gebied = selectie_openheid_klasses,
+      exclusie = exclusie_osm_landgebruiken,
+      osmdata = osm_belgium,
+      paths_include = c("track", "tertiary", "tertiary_link", "unclassified"),
+      cutting_exclude = NULL,
+      historic_exclude = NULL,
+      waterway = NULL,
+      update_osm_layer = FALSE
+      ),
+    pattern = map(selectie_openheid_klasses, exclusie_osm_landgebruiken),
+    iteration = "list"
+  ),
+  tar_target(
+    name = punten,
+    command = paden_naar_punten(
+      data_paden = paden,
+      gebieden = perimeters_data,
+      interpoint_distance = 50,
+      border_distance = 300
+      ),
+    pattern = map(perimeters_data, paden),
+    iteration = "list"
+  ) ,
+  tar_target(
+    name = telcirkels_landgebruik,
+    command = punten_lum_buffer(
+      punten_sf = punten,
+      radius = 300,
+      file = vito_lum_2019_file,
+      legend = legend_lum
+      ),
+    pattern = map(punten),
+    iteration = "list"
+  ),
+  tar_target(
+    name = telcirkels_selectie_landgebruik,
+      command = punten_selectie_landgebruik(
+        lum_extract_result = telcirkels_landgebruik,
+        legend_rast = legend_lum,
+        max_prop_overige = 0.5,
+        min_prop_akker = 0.3,
+        min_prop_akker_grasland = 0.4
+        ),
+    pattern = map(telcirkels_landgebruik)
+  ),
+  tar_target(
+    name = selectie_landgebruik,
+    command = selectie_landgebruik_vito(
+      punten_sf = punten,
+      selectie_df = telcirkels_selectie_landgebruik
+      ),
+    pattern = map(punten, telcirkels_selectie_landgebruik),
+    iteration = "list"
+  ),
+  tarchetypes::tar_group_size(
+    name = selectie_landgebruik_per_size,
+    command = do.call(rbind.data.frame, selectie_landgebruik),
+    size = 200
+  ),
+  tar_target(
+    name = punten_zichtbaarheid,
+    command = add_visibility_to_frame(
+      punten_sf = selectie_landgebruik_per_size,
+      resolution = 5,
+      spacing = 10
+    ),
+    pattern = map(selectie_landgebruik_per_size),
+    iteration = "list"
+  ),
+  tarchetypes::tar_group_by(
+    name = punten_zichtbaarheid_per_regio,
+    command = do.call(rbind.data.frame, punten_zichtbaarheid),
+    Naam
+  ),
+  tar_target(
+    name = punten_selectie_zichtbaarheid,
+    command = filter_zichtbaarheid(
+      punten_sf = punten_zichtbaarheid_per_regio,
+      min_cvvi = 0.1
+    ),
+    pattern = map(punten_zichtbaarheid_per_regio),
+    iteration = "list"
+  )
+)

source/targets/mas_steekproef_zl_lm_za_kp_po/_targets_r/targets/steekproeftrekking.R

@@ -0,0 +1,47 @@
+list(
+  tar_target(
+    name = target_size,
+    command = calc_target_samplesize(gebied = selectie_openheid_klasses),
+    pattern = map(selectie_openheid_klasses)
+  ),
+  tar_target(
+    allocatie_df,
+    allocatie(
+      steekproefkader = steekproefkader_finaal,
+      min_samplesize = 200,
+      target_samplesize = target_size,
+      popsize_minimum = 200,
+      allocatie_binnen_sbp = 0.5,
+      allocatie_leemstreek = 0,
+      ol_strata = c("OL", "HOL")
+      ),
+    pattern = map(steekproefkader_finaal, target_size)
+  ),
+  tarchetypes::tar_group_by(
+    allocatie_per_stratum,
+    allocatie_df,
+    Naam,
+    is_sbp,
+    openheid_klasse
+  ),
+  tarchetypes::tar_group_by(
+    steekproefkader_per_stratum,
+    do.call(rbind.data.frame, steekproefkader_finaal) %>%
+      semi_join(allocatie_df %>%
+                  select(Naam, is_sbp, openheid_klasse),
+                by = c("Naam", "is_sbp", "openheid_klasse")),
+    Naam,
+    is_sbp,
+    openheid_klasse
+  ),
+  tar_target(
+    steekproef,
+    draw_sample(
+      sampling_frame = steekproefkader_per_stratum,
+      sample_size_multiplication = 1,
+      balance = c("X", "Y", "area_prop_sb")
+      ),
+    pattern = map(steekproefkader_per_stratum, allocatie_per_stratum),
+    iteration = "list"
+  )
+)

source/targets/mas_steekproef_zl_lm_za_kp_po/_targets_r/targets/version-control.R

@@ -0,0 +1,10 @@
+list(
+  tar_target(
+    output_objecten,
+    output_finaal(
+      files = list(steekproefkader = steekproefkader_finaal,
+                   steekproef = steekproef_final_600),
+      write_out = FALSE
+    )
+  )
+)