From 3d8c7be442feadc68a6779aa875bd761d44f453d Mon Sep 17 00:00:00 2001
From: "WUR\\simon083" <wolfram.simon@wur.nl>
Date: Thu, 12 Aug 2021 13:35:02 +0200
Subject: [PATCH] externalized the functions. Fixed the forage.
---
.../SPAM2010_ADM0_SOILagri_AEZ_SpatialEx.R | 54 +++----------------
1 file changed, 7 insertions(+), 47 deletions(-)
diff --git a/SPAM_Data_Transformation/SPAM2010_ADM0_SOILagri_AEZ_SpatialEx.R b/SPAM_Data_Transformation/SPAM2010_ADM0_SOILagri_AEZ_SpatialEx.R
index bc6eb7f..d26c47a 100644
--- a/SPAM_Data_Transformation/SPAM2010_ADM0_SOILagri_AEZ_SpatialEx.R
+++ b/SPAM_Data_Transformation/SPAM2010_ADM0_SOILagri_AEZ_SpatialEx.R
@@ -11,49 +11,8 @@ easypackages::packages(c("tidyverse", 'sf', 'raster','sp', 'spatialEco',
"rlist", "largeList", "snow", "foreign",
"reshape2", "margrittr", "stringi"))
-
-# Functions ---------------------------------------------------------------
-# Functions loading dbf SPAM files
-FUNload_list_dbf_SPAM <- function(path, type) {
- pattern <- switch (type,
- "harv" = "_H_",
- "phys" = "_A_",
- "prod" = "_P_")
- path <- list.files(path = path,
- full.names = TRUE,
- pattern = pattern)
- l_dbf <- lapply(path, foreign::read.dbf)
- names(l_dbf) <- basename(path)
- return(l_dbf)
-}
-
-# Functions loading dbf veggie and nuts Monfreda files
-FUNload_list_monfreda <- function(path, type) {
- pattern <- switch (type,
- "harv" = "_HarvestedAreaHectares",
- "phys" = "_HarvestedAreaHectares",
- "prod" = "_Production")
- path <- list.files(path = path,
- full.names = TRUE,
- pattern = pattern)
- l_csv <- lapply(path, read.csv)
- names(l_csv) <- basename(path)
- return(l_csv)
-}
-
-# # Functions loading dbf forage Monfreda files
-# FUNload_list_forage <- function(path, type) {
-# pattern <- switch (type,
-# "harv" = "_HarvestedAreaHectares",
-# "phys" = "_HarvestedAreaHectares",
-# "prod" = "_Production")
-# path <- list.files(path = path,
-# full.names = TRUE,
-# pattern = pattern)
-# l_csv <- lapply(path, read.csv)
-# names(l_csv) <- basename(path)
-# return(l_csv)
-# }
+# Sourcing functions from external R script
+source("Functions_cifos/function_collection.R")
# Load data --------------------------------------------------------------
SPAMdbf_harv = FUNload_list_dbf_SPAM(path="C:/Users/simon083/OneDrive - Wageningen University & Research/PhD_WJS/Academic/RQ1/Data_analysis/Spatial_mod/Modified_scripts/R/CifoS_Crop/Input_data/DBF_ADM0_AEZ_Soilagri_exctraction/SPAM2010", type = "harv")
@@ -93,8 +52,8 @@ SPAM2010_ZONES_exdbf = df_rec %>%
'H.dbf' = 'H'))) %>%
dplyr::select(-Rec_type)
-write_csv(SPAM2010_ZONES_exdbf, here("Input_data", "SPAM2010_ADM0_AEZ_Soilagri_exdbf_global.csv"))
-# SPAM2010_ZONES_exdbf =read_csv(here("Input_data", "SPAM2010_adm1_exdbf_global.csv"))
+write_csv(SPAM2010_ZONES_exdbf, here::here("Input_data", "SPAM2010_ADM0_AEZ_Soilagri_exdbf_global.csv"))
+SPAM2010_ZONES_exdbf =read_csv(here::here("Input_data", "SPAM2010_ADM0_AEZ_Soilagri_exdbf_global.csv"))
# Monfreda data ----------------------------------------------------------------
# Load data --------------------------------------------------------------
@@ -128,8 +87,8 @@ monf_dat = monf_rec %>%
dplyr::select(-c(Rec_type, stats)) %>%
mutate(tech_type = 'A') #add a tech type all together ("A") -> Monfreda has no distinction btw tech types like SPAM does.
-write_csv(monf_dat, here("Input_data", "Monfreda2008_ADM0_AEZ_Soilagri_excsv_global.csv"))
-# SPAM2010_ZONES_exdbf =read_csv(here("Input_data", "SPAM2010_adm1_exdbf_global.csv"))
+write_csv(monf_dat, here::here("Input_data", "Monfreda2008_ADM0_AEZ_Soilagri_excsv_global.csv"))
+monf_dat = read_csv(here::here("Input_data", "Monfreda2008_ADM0_AEZ_Soilagri_excsv_global.csv"))
# vector containing all different category names
# nuts
@@ -219,6 +178,7 @@ Forage_nuts_veg_spam = SPAM2010_global_ZONES %>%
bind_rows(., forageCrop_dat, forageGrass_dat)
write_csv(Forage_nuts_veg_spam, here("Input_data", "SPAM2010_global_nuts_veg_adm0_aez_soilagri_forage.csv"))
+Forage_nuts_veg_spam = read_csv(here::here("Input_data", "SPAM2010_global_nuts_veg_adm0_aez_soilagri_forage.csv"))
# Subsetting -------------------------------------------------------------
EU_iso3 = c("AUT","BEL", "BGR", "HRV", "CZE", "DNK", "EST","FIN", "FRA", "DEU", "GRC", "HUN", "IRL", "ITA","LVA",
--
GitLab