Merge pull request #4 from leymanan/open_data

Open data
inbo · Jun 13, 2023 · 4508fe6 · 4508fe6
2 parents 3b7e5f8 + 9f52f8c
commit 4508fe6
Show file tree

Hide file tree

Showing 26 changed files with 5,663 additions and 710 deletions.
diff --git a/.gitignore b/.gitignore
@@ -1,6 +1,7 @@
 .Rproj.user
 .Rhistory
-.RData
+*.RData
+*.Rdata
 .Ruserdata
 Data/
 Output/

diff --git a/...R_VBI/meanNGV_per_bms_levend_dood_PNV.Rmd → ...R_VBI/meanNGV_per_bms_levend_dood_PNV.Rmd b/...R_VBI/meanNGV_per_bms_levend_dood_PNV.Rmd → ...R_VBI/meanNGV_per_bms_levend_dood_PNV.Rmd
@@ -70,8 +70,8 @@ Folder op gdrive:
 # deze worden ook in "Setup.R" gedefinieerd (aanroepen van "Paths.R")
 # hier ev. extra "paths" toevoegen indien nodig 
 
-path_to_request <- paste0(path_to_datarequests, "INBO_biomassa_vglBR_VBI/")
-path_to_gdrive <- paste0(path_to_datarequests_gdrive, "INBO_biomassa_vglBR_VBI/")
+path_to_request <- paste0(path_to_datarequests, "INBO_artikel_vglCarbon_BR_VBI/")
+path_to_gdrive <- paste0(path_to_datarequests_gdrive, "INBO_artikel_vglCarbon_BR_VBI/")
 
 # rechtstreeks naar data-folder van artikel
 path_to_article <- "G:/.shortcut-targets-by-id/1GY3uYeL5JjG14h_BAkzsQwse31el-yhC/artikel_FrontiersForestsGlobalChange/data/"
@@ -834,13 +834,13 @@ write.xlsx2(x = Result_NGV_per_PNV_per_species,
 
 # idem - aandelen
 write.xlsx(result_G_species_pnv
-           ,paste0(path_to_request, "aandeel_NGV_per_pnv.xlsx")
+           ,paste0(path_to_request, "NGV_per_pnv_aandeel.xlsx")
            ,sheetName = "G")
-write.xlsx(result_N_species_pnv,paste0(path_to_request, "aandeel_NGV_per_pnv.xlsx")
+write.xlsx(result_N_species_pnv,paste0(path_to_request, "NGV_per_pnv_aandeel.xlsx")
            ,sheetName = "N", append = T)
-write.xlsx(result_V_species_pnv,paste0(path_to_request, "aandeel_NGV_per_pnv.xlsx")
+write.xlsx(result_V_species_pnv,paste0(path_to_request, "NGV_per_pnv_aandeel.xlsx")
            ,sheetName = "V", append = T)
-write.xlsx(result_Vdead_species_pnv,paste0(path_to_request, "aandeel_NGV_per_pnv.xlsx")
+write.xlsx(result_Vdead_species_pnv,paste0(path_to_request, "NGV_per_pnv_aandeel.xlsx")
            ,sheetName = "V_dead", append = T)
 
 ```
@@ -856,13 +856,13 @@ write.xlsx2(x = Result_NGV_per_PNV_per_species,
 
 # idem - aandelen
 write.xlsx(result_G_species_pnv
-           ,paste0(path_to_gdrive, "aandeel_NGV_per_pnv.xlsx")
+           ,paste0(path_to_gdrive, "NGV_per_pnv_aandeel.xlsx")
            ,sheetName = "G")
-write.xlsx(result_N_species_pnv,paste0(path_to_gdrive, "aandeel_NGV_per_pnv.xlsx")
+write.xlsx(result_N_species_pnv,paste0(path_to_gdrive, "NGV_per_pnv_aandeel.xlsx")
            ,sheetName = "N", append = T)
-write.xlsx(result_V_species_pnv,paste0(path_to_gdrive, "aandeel_NGV_per_pnv.xlsx")
+write.xlsx(result_V_species_pnv,paste0(path_to_gdrive, "NGV_per_pnv_aandeel.xlsx")
            ,sheetName = "V", append = T)
-write.xlsx(result_Vdead_species_pnv,paste0(path_to_gdrive, "aandeel_NGV_per_pnv.xlsx")
+write.xlsx(result_Vdead_species_pnv,paste0(path_to_gdrive, "NGV_per_pnv_aandeel.xlsx")
            ,sheetName = "V_dead", append = T)
 
 ```

diff --git a/...BI/plot_level_C_levend_dood_PNV_G_bms.Rmd → ...BI/plot_level_C_levend_dood_PNV_G_bms.Rmd b/...BI/plot_level_C_levend_dood_PNV_G_bms.Rmd → ...BI/plot_level_C_levend_dood_PNV_G_bms.Rmd
@@ -37,14 +37,27 @@ source(here::here("scripts/Setup.R"))
 ```
 
 
+```{r Path}
+# deze worden ook in "Setup.R" gedefinieerd (aanroepen van "Paths.R")
+# hier ev. extra "paths" toevoegen indien nodig 
+
+path_to_request <- paste0(path_to_datarequests, "INBO_artikel_vglCarbon_BR_VBI/")
+path_to_gdrive <- paste0(path_to_datarequests_gdrive, "INBO_artikel_vglCarbon_BR_VBI/")
+
+# rechtstreeks naar data-folder van artikel
+path_to_article <- "G:/.shortcut-targets-by-id/1GY3uYeL5JjG14h_BAkzsQwse31el-yhC/artikel_FrontiersForestsGlobalChange/data/"
+# !! oppassen om niks te overschrijven
+```
+
+
 # Vraagstelling
 
 Kadert in de opmaak van een artikel in themanummer Frontiers Forests Global Change.
 
 Gevraagd (Luc DK):
 
 - op plotniveau (cirkel of kernvlakte gebruikt voor de NARA gegevens)   
-- C voorraden in levende en dode biomassa !! bovengronds
+- C voorraden in levende en dode biomassa !! enkel bovengronds (ABG)
 - voor de twee tijdstippen
 - met PNV type 
 - globaal grondvlak
@@ -56,17 +69,15 @@ Folder op gdrive:
 - [artikel](https://drive.google.com/drive/folders/1GY3uYeL5JjG14h_BAkzsQwse31el-yhC)   
 - [artikel - data](https://drive.google.com/drive/folders/1J9ZOC5K6WqjxuV298dd-iNfcKzLvAWjp?usp=share_link)   
 
-```{r Path}
-# deze worden ook in "Setup.R" gedefinieerd (aanroepen van "Paths.R")
-# hier ev. extra "paths" toevoegen indien nodig 
+# Referentie C-berekening levende biomassa
 
-path_to_request <- paste0(path_to_datarequests, "INBO_biomassa_vglBR_VBI/")
-path_to_gdrive <- paste0(path_to_datarequests_gdrive, "INBO_biomassa_vglBR_VBI/")
+< S. Bauwens: 
+Note technique version 20/11/2017. Estimation du volume et biomasse de différents compartiments de l’arbre. 
+Accompagnement scientifique de l’IPRFW. ULiège - Wallonie environnement SPW.
+
+In artikel: 
+Perin, J., Bauwens, S., Pitchugin, M., Lejeune, P., Hébert, J. (2019). National forestry accounting plan of Belgium. Gembloux: University of Liège.
 
-# rechtstreeks naar data-folder van artikel
-path_to_article <- "G:/.shortcut-targets-by-id/1GY3uYeL5JjG14h_BAkzsQwse31el-yhC/artikel_FrontiersForestsGlobalChange/data/"
-# !! oppassen van nisk te overschrijven
-```
 
 
 # Inladen analysesets NARA
@@ -128,7 +139,7 @@ analyseSet_all <- analyseSet_alive %>%
 
 # Logs CP's Kersselaerspleyn, periode 1
 
-Zitten niet in fieldmap, maar er zijn wel gegevens uit het monitorinngrapport, zie 
+Zitten niet in fieldmap, maar er zijn wel gegevens uit het monitoringrapport, zie 
 "CP_Kerss_logs_2000.xlsx".
 
 ```{r}
@@ -316,14 +327,14 @@ vgl <- oude_data %>%
 # Export
 
 ```{r export_to_csv}
-write_excel_csv2(analyseSet_finaal,paste0(path_to_request, "carbon_plot_level.csv"))
-write_excel_csv2(analyseSet_finaal_wide,paste0(path_to_request, "carbon_plot_level_wide.csv"))
+write_excel_csv2(analyseSet_finaal,paste0(path_to_request, "carbon_ABG_plot_level.csv"))
+write_excel_csv2(analyseSet_finaal_wide,paste0(path_to_request, "carbon_ABG_plot_level_wide.csv"))
 
 ```
 
 ```{r export_to_csv_gdrive}
-write_excel_csv2(analyseSet_finaal,paste0(path_to_gdrive, "carbon_plot_level.csv"))
-write_excel_csv2(analyseSet_finaal_wide,paste0(path_to_gdrive, "carbon_plot_level_wide.csv"))
+write_excel_csv2(analyseSet_finaal,paste0(path_to_gdrive, "carbon_ABG_plot_level.csv"))
+write_excel_csv2(analyseSet_finaal_wide,paste0(path_to_gdrive, "carbon_ABG_plot_level_wide.csv"))
 
 ```
 

diff --git a/...NBO_biomassa_vglBR_VBI/volume_per_PNV.Rmd → ..._level_volume_en_total_volume_per_PNV.Rmd b/...NBO_biomassa_vglBR_VBI/volume_per_PNV.Rmd → ..._level_volume_en_total_volume_per_PNV.Rmd
@@ -54,12 +54,12 @@ Folder op gdrive:
 # deze worden ook in "Setup.R" gedefinieerd (aanroepen van "Paths.R")
 # hier ev. extra "paths" toevoegen indien nodig 
 
-path_to_request <- paste0(path_to_datarequests, "INBO_biomassa_vglBR_VBI/")
-path_to_gdrive <- paste0(path_to_datarequests_gdrive, "INBO_biomassa_vglBR_VBI/")
+path_to_request <- paste0(path_to_datarequests, "INBO_artikel_vglCarbon_BR_VBI/")
+path_to_gdrive <- paste0(path_to_datarequests_gdrive, "INBO_artikel_vglCarbon_BR_VBI/")
 
 # rechtstreeks naar data-folder van artikel
 path_to_article <- "G:/.shortcut-targets-by-id/1GY3uYeL5JjG14h_BAkzsQwse31el-yhC/artikel_FrontiersForestsGlobalChange/data/"
-# !! oppassen van nisk te overschrijven
+# !! oppassen om niks te overschrijven
 ```
 
 ```{r load_info_lkp}
@@ -91,7 +91,8 @@ n_distinct(analyseSet$plot_id) == nrow(analyseSet)/2
 
 plots <- unique(analyseSet$plot_id)
 ```
-
+Deze analyseset bevat geen "volume inclusief kroon", enkel "volume bole".
+Daarom `dendro_by_plot` importeren.
 
 ```{r load_dendro_by_plot}
 # data_dendro <- load_data_dendrometry(path_to_fieldmap, plottype = "Circular plot", extra_variables = TRUE)
@@ -118,6 +119,19 @@ analyseSet <- dendro_by_plot %>%
 
 ```
 
+# Export plot-level
+
+```{r eval = FALSE}
+write_excel_csv2(analyseSet,paste0(path_to_request, "volume_plot_level.csv"))
+
+write_excel_csv2(analyseSet,paste0(path_to_gdrive, "volume_plot_level.csv"))
+
+# oppassen om niet te overschrijven
+write_excel_csv2(analyseSet,paste0(path_to_article, "BR_volume_plot_level.csv"))
+
+
+```
+
 
 
 # Statistics per PNV