diff --git a/source/markdown/exploratie_duinen_weidestreek.Rmd b/source/markdown/exploratie_duinen_weidestreek.Rmd index 02836b8..9fe264a 100644 --- a/source/markdown/exploratie_duinen_weidestreek.Rmd +++ b/source/markdown/exploratie_duinen_weidestreek.Rmd @@ -122,7 +122,10 @@ selectie_openheid_klasses <- lapply(kleine_landbouwstreken, function(streek) { ### Duinen ```{r} -mapview(kleine_landbouwstreken$Duinen, color = "red", alpha.regions = 0, legend = FALSE) + +mapview(kleine_landbouwstreken$Duinen, + color = "red", + alpha.regions = 0, + legend = FALSE) + mapview(openheid_sf$Duinen, zcol = "openheid", layer = "openheid") ``` @@ -150,7 +153,10 @@ openheid_sf$Duinen %>% ### Weidestreek ```{r} -mapview(kleine_landbouwstreken$Weidestreek, color = "red", alpha.regions = 0, legend = FALSE) + +mapview(kleine_landbouwstreken$Weidestreek, + color = "red", + alpha.regions = 0, + legend = FALSE) + mapview(openheid_sf$Weidestreek, zcol = "openheid", layer = "openheid") ``` @@ -783,7 +789,8 @@ sbp_akkervogels_duinen <- read_sbp_akkervogels( sbp_akkervogels_weidestreek <- read_sbp_others( path_to_sbp_akkervogels(file = "sbp_overige_soorten.shp"), - soorten = c("hamster", "bruine kiekendief", "zomertortel", "grauwe kiekendief"), + soorten = c("hamster", "bruine kiekendief", + "zomertortel", "grauwe kiekendief"), gebied = kleine_landbouwstreken$Weidestreek ) %>% summarise() @@ -1244,11 +1251,13 @@ balansvergelijking_duinen <- steekproef_thinned$Duinen %>% filter(batch == "eerste set") %>% select_at(names(steekproefkader_finaal$Duinen)) %>% mutate(fill_var = "steekproef") %>% - bind_rows(steekproefkader_finaal$Duinen %>% mutate(fill_var = "steekproefkader")) %>% + bind_rows(steekproefkader_finaal$Duinen %>% + mutate(fill_var = "steekproefkader")) %>% mutate(sbp_akkervogels = ifelse(is_sbp, "binnen SBP", "buiten SBP")) %>% cbind(st_coordinates(.)) %>% st_drop_geometry() %>% - mutate(fill_var = factor(fill_var, levels = c("steekproefkader", "steekproef"))) + mutate(fill_var = factor(fill_var, + levels = c("steekproefkader", "steekproef"))) count_df_duinen <- balansvergelijking_duinen %>% group_by(fill_var, openheid_klasse, sbp_akkervogels) %>% @@ -1354,11 +1363,13 @@ balansvergelijking_weidestreek <- steekproef_thinned$Weidestreek %>% filter(batch == "eerste set") %>% select_at(names(steekproefkader_finaal$Weidestreek)) %>% mutate(fill_var = "steekproef") %>% - bind_rows(steekproefkader_finaal$Weidestreek %>% mutate(fill_var = "steekproefkader")) %>% + bind_rows(steekproefkader_finaal$Weidestreek %>% + mutate(fill_var = "steekproefkader")) %>% mutate(sbp_akkervogels = ifelse(is_sbp, "binnen SBP", "buiten SBP")) %>% cbind(st_coordinates(.)) %>% st_drop_geometry() %>% - mutate(fill_var = factor(fill_var, levels = c("steekproefkader", "steekproef"))) + mutate(fill_var = factor(fill_var, + levels = c("steekproefkader", "steekproef"))) count_df_weidestreek <- balansvergelijking_weidestreek %>% group_by(fill_var, openheid_klasse, sbp_akkervogels) %>% @@ -1447,7 +1458,7 @@ sizes <- 40:100 ``` De vervolgvraag stelt zich of we een grotere finale steekproef kunnen bekomen door geen reservepunten te voorzien. -We analiseren groottes van `r min(sizes)` tot `r max(sizes)`. +We analyseren groottes van `r min(sizes)` tot `r max(sizes)`. ```{r} calc_chisq <- function(observed, expected) { @@ -1673,11 +1684,13 @@ steekproef_extended_thinned %>% balansvergelijking_weidestreek2 <- steekproef_extended_thinned %>% select_at(names(steekproefkader_finaal$Weidestreek)) %>% mutate(fill_var = "steekproef") %>% - bind_rows(steekproefkader_finaal$Weidestreek %>% mutate(fill_var = "steekproefkader")) %>% + bind_rows(steekproefkader_finaal$Weidestreek %>% + mutate(fill_var = "steekproefkader")) %>% mutate(sbp_akkervogels = ifelse(is_sbp, "binnen SBP", "buiten SBP")) %>% cbind(st_coordinates(.)) %>% st_drop_geometry() %>% - mutate(fill_var = factor(fill_var, levels = c("steekproefkader", "steekproef"))) + mutate(fill_var = factor(fill_var, + levels = c("steekproefkader", "steekproef"))) count_df_weidestreek2 <- balansvergelijking_weidestreek2 %>% group_by(fill_var, openheid_klasse, sbp_akkervogels) %>% diff --git a/source/markdown/steekproefontwerp_kempen.Rmd b/source/markdown/steekproefontwerp_kempen.Rmd index cccac9d..a787534 100644 --- a/source/markdown/steekproefontwerp_kempen.Rmd +++ b/source/markdown/steekproefontwerp_kempen.Rmd @@ -1213,11 +1213,12 @@ intersect %>% ```{r} # visualisatie -mapview(oude_cirkels_buffer %>% filter(definitief_punt %in% intersect$definitief_punt), - col.regions = "red", layer = "bestaande punten" +mapview(oude_cirkels_buffer %>% + filter(definitief_punt %in% intersect$definitief_punt), + col.regions = "red", layer = "bestaande punten" ) + mapview(cirkels_kempen %>% filter(pointid %in% intersect$pointid), - layer = "nieuwe punten" + layer = "nieuwe punten" ) ``` @@ -1349,7 +1350,9 @@ kempen_avimap_df <- read_delim( "naam" = Naam, "pointid_from_targets" = pointid ) %>% - mutate(pointid = ifelse(is.na(vervanging_id), pointid_from_targets, vervanging_id)) + mutate(pointid = ifelse(is.na(vervanging_id), + pointid_from_targets, + vervanging_id)) kempen_avimap_sf <- st_as_sf(kempen_avimap_df, wkt = "WKT", diff --git a/source/markdown/steekproefontwerp_polders.Rmd b/source/markdown/steekproefontwerp_polders.Rmd index 2d3cced..4e5d2d8 100644 --- a/source/markdown/steekproefontwerp_polders.Rmd +++ b/source/markdown/steekproefontwerp_polders.Rmd @@ -1211,11 +1211,12 @@ intersect %>% ```{r} # visualisatie -mapview(oude_cirkels_buffer %>% filter(definitief_punt %in% intersect$definitief_punt), - col.regions = "red", layer = "bestaande punten" +mapview(oude_cirkels_buffer %>% + filter(definitief_punt %in% intersect$definitief_punt), + col.regions = "red", layer = "bestaande punten" ) + mapview(cirkels_polders %>% filter(pointid %in% intersect$pointid), - layer = "nieuwe punten" + layer = "nieuwe punten" ) ``` @@ -1344,7 +1345,9 @@ polders_avimap_df <- read_csv(file.path( "naam" = Naam, "pointid_from_targets" = pointid ) %>% - mutate(pointid = ifelse(is.na(vervanging_id), pointid_from_targets, vervanging_id)) + mutate(pointid = ifelse(is.na(vervanging_id), + pointid_from_targets, + vervanging_id)) polders_avimap_sf <- st_as_sf(polders_avimap_df, diff --git a/source/markdown/steekproefontwerp_zandstreek.Rmd b/source/markdown/steekproefontwerp_zandstreek.Rmd index 5371ed4..fe7ac74 100644 --- a/source/markdown/steekproefontwerp_zandstreek.Rmd +++ b/source/markdown/steekproefontwerp_zandstreek.Rmd @@ -1213,11 +1213,12 @@ intersect %>% ```{r} # visualisatie -mapview(oude_cirkels_buffer %>% filter(definitief_punt %in% intersect$definitief_punt), - col.regions = "red", layer = "bestaande punten" +mapview(oude_cirkels_buffer %>% + filter(definitief_punt %in% intersect$definitief_punt), + col.regions = "red", layer = "bestaande punten" ) + mapview(cirkels_zandstreek %>% filter(pointid %in% intersect$pointid), - layer = "nieuwe punten" + layer = "nieuwe punten" ) ``` @@ -1346,7 +1347,9 @@ zandstreek_avimap_df <- read_csv(file.path( "naam" = Naam, "pointid_from_targets" = pointid ) %>% - mutate(pointid = ifelse(is.na(vervanging_id), pointid_from_targets, vervanging_id)) + mutate(pointid = ifelse(is.na(vervanging_id), + pointid_from_targets, + vervanging_id)) zandstreek_avimap_sf <- st_as_sf(zandstreek_avimap_df, wkt = "WKT",