Creating a dataset with rangeland, arable land, arable pasture, natural...

Creating a dataset with rangeland, arable land, arable pasture, natural pasture per zone (adm0, aez, soilagri). draft, not finished.

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