Ajout de la préparation des domaines climat, ger et wrb

1-DataPreparationJRC.R

@@ -32,13 +32,16 @@ if(change) {
   
 }
 
+
 # 2 prepare soil region data --------
 
-library(sf)
+## soil region
 soil <- st_read("~/serena/data/sml/europe/eusr5000.shp") %>%
   st_transform(crs(stacka))
 
 soilR <- rasterize(soil, stacka, field = "SOILNR")
+tst = as.data.frame(soilR)
+
 
 stacka <- rast(c("FR_OC.tif") ) 
 
@@ -53,3 +56,52 @@ plot(stack)
 
 terra::writeRaster(soilRextr,filename = "~/serena/data/sml/europe/soilRegion.tiff",
                    overwrite = TRUE)
+
+
+## WRB 1
+### 100
+wrb1 <- rast("/home/lpotel/smlsamplingeu/data/shp_pour_raster/soil_map_WRB_final.tif")
+
+wrb2 <- terra::project(wrb1,
+                       stacka,
+                       method= "near")
+
+# wbr3 <- crop(wrb2 , stacka[[1]] ,mask=TRUE)
+wrb3 <- mask(wrb2 , stacka[[1]])
+
+test = c(stacka,wrb3)
+
+terra::writeRaster(wrb3,filename = "~/serena/data/sml/europe/wrb3.tiff",
+                   overwrite = TRUE)
+
+rm(wrb1,wrb2)
+
+### 500
+
+wrb2 <- terra::project(wrb1,
+                       stacka500,
+                       method= "near")
+
+
+test = c(stacka,wrb2)
+
+## GER IGCS
+### 100 m
+soil <- st_read("/home/lpotel/smlsamplingeu/data/shp_pour_raster/carte_sol_France.shp") %>%
+  st_transform(crs(stacka))
+
+soilR <- rasterize(soil, stacka, field = "ger")
+
+test = c(stacka,soilR)
+
+
+
+freq(soilR)
+terra::writeRaster(soilR,filename = "~/serena/data/sml/europe/",
+                   overwrite = TRUE)
+
+# 3 Climate data
+## Metzger
+
+
+