From 71ad9a4b13f5a2de7a2e5382bdee4cb43506c4b8 Mon Sep 17 00:00:00 2001
From: "WUR\\simon083" <wolfram.simon@wur.nl>
Date: Wed, 11 Aug 2021 16:35:59 +0200
Subject: [PATCH] Creating a dataset with rangeland, arable land, arable
 pasture, natural pasture per zone (adm0, aez, soilagri). draft, not finished.

---
 Zones/Zones_AEZ.R | 12 +++++++-----
 1 file changed, 7 insertions(+), 5 deletions(-)

diff --git a/Zones/Zones_AEZ.R b/Zones/Zones_AEZ.R
index 64c0b2f..c83987e 100644
--- a/Zones/Zones_AEZ.R
+++ b/Zones/Zones_AEZ.R
@@ -148,7 +148,8 @@ dat_soil = dat_moist %>%
                            DOMS =="WR"~ "other",
                            DOMS =="ND"~ "other",
                           DOMS =="GG"~ "other",
-                          DOMS =="RK"~ "other"),
+                          DOMS =="RK"~ "other",
+                          DOMS =="ST"~ "other"),
          Soil = case_when(DOMS =="ORG"~ "organic",TRUE ~ Soil)) %>% 
   dplyr::select(AEZ,DOMS, ADM1_CODE, ADM0_CODE, Adm0_Adm1, Moist, Soil,geometry)
 
@@ -163,12 +164,13 @@ dat_agrisoil = dat_soil %>%
                           DOMS =="WR"~ "other",
                           DOMS =="ND"~ "other",
                           DOMS =="GG"~ "other",
-                          DOMS =="RK"~ "other"),
+                          DOMS =="RK"~ "other",
+                          DOMS =="ST"~ "other"),
          Soil_agri = case_when(DOMS =="ORG"~ "organic",TRUE ~ Soil_agri)) %>% 
-  dplyr::select(AEZ,DOMS, ADM1_CODE, ADM0_CODE, Adm0_Adm1, Moist, Soil, Soil_agri, geometry)
-
-
+  dplyr::select(AEZ,DOMS, ADM1_CODE, ADM0_CODE, Adm0_Adm1, Moist, Soil, Soil_agri, geometry) %>% 
+  mutate_at(vars(Soil_agri, Soil, DOMS), ~replace_na(., "other"))
 
+soil = dat_agrisoil %>% st_set_geometry(NULL)
 # 1         HAC - Low activity clay soils
 # 2          WR - Water 
 # 3         POD - Spodic soils
-- 
GitLab