From f17f2846ae5e23f1c177f5b4b6844f8a6ba13738 Mon Sep 17 00:00:00 2001
From: nsaby <nicolas.saby@inrae.fr>
Date: Wed, 14 Feb 2024 15:31:36 +0000
Subject: [PATCH] tidy things
---
1- DataPreparationJRC.R | 35 +++++++++++++++++++
...ableJRCTests.R => 2-SMLFilltableJRCTests.R | 0
README.md | 3 +-
3 files changed, 37 insertions(+), 1 deletion(-)
create mode 100644 1- DataPreparationJRC.R
rename 1-SMLFilltableJRCTests.R => 2-SMLFilltableJRCTests.R (100%)
diff --git a/1- DataPreparationJRC.R b/1- DataPreparationJRC.R
new file mode 100644
index 0000000..8aee231
--- /dev/null
+++ b/1- DataPreparationJRC.R
@@ -0,0 +1,35 @@
+library(sf)
+library(terra)
+library(ggplot2)
+library(dplyr)
+library(tidyr)
+library(clhs)
+library(viridis)
+library(tmap)
+
+# 1 Data preparation -----
+
+setwd("~/serena/data/sml/europe/")
+stacka <- rast(c("FR_OC.tif", "FR_pH.tif" , "FR_TXT.tif","FR_N.tif","FR_P.tif",
+ "FR_CEC.tif","FR_BD010.tif") )
+
+stackb <- rast(c("FR_CLC.tif", "FR_NUTS2.tif","FR_NUTS1.tif"))
+stack = c(stacka,stackb)
+
+change = FALSE # set true to resample rasters
+
+if(change) {
+ # Change resolution
+ stac500a <- terra::aggregate(stacka, 5, fun=mean)
+ stac500b <- terra::aggregate(stackb, 5, fun=modal)
+ stack500 = c(stac500a,stac500b)
+
+ terra::writeRaster(stack500,filename = "~/serena/data/sml/europe/stack500.tiff",
+ overwrite = TRUE)
+
+} else {
+ stack500 = rast("~/serena/data/sml/europe/stack500.tiff")
+
+}
+
+
diff --git a/1-SMLFilltableJRCTests.R b/2-SMLFilltableJRCTests.R
similarity index 100%
rename from 1-SMLFilltableJRCTests.R
rename to 2-SMLFilltableJRCTests.R
diff --git a/README.md b/README.md
index 81fb802..bbbce66 100644
--- a/README.md
+++ b/README.md
@@ -4,7 +4,8 @@ this project to share the statistical work on sampling exercise
# code
-use smltest.R to run the optimisation algorithm
+use 2-SMLFilltableJRCTests.R to run different tests required.
+use 1-DataPreparationJRC.R to prepare the rasters from JRC sharepoint
# Questions
--
GitLab