Skip to content
Snippets Groups Projects
Commit c14cb441 authored by lopotel's avatar lopotel
Browse files

modifications du coe pour le transformer en CLC 

(rmqs204-25, sasn rmqs)
parent 625bd389
No related branches found
No related tags found
1 merge request!1Resolve "adapter les scripts dirsol pour le calcul avec CLC"
Hide whitespace changes
Inline Side-by-side
dirsol/Simul_DIRSOL_avec_RMQS_2024_2025_CLC.R
+ 25
− 142
View file @ c14cb441
# nouveau CLC
# Ce script permet de faire tourner les 4 scénarios aux deux résolutions attendues, en prenant en compte
# les points du RMQS prospéctés en 2024-2025 et en utilisant les données d'occupation du sol Corinne Land Cover
# les points du RMQS prospéctés en 2024-2025 et en utilisant les données d'occupation du sol CLC
# initialisations --------
......@@ -43,8 +43,8 @@ sample3 <- sample3 %>%
sample3 <- filter(sample3, !FR_WRB %in% excludeWRB )
# saveRDS( sample3, "~/smlsamplingeu_louis/data/octobre_avec_rmqs2024_2025_sample3_lulucf.rds")
# sample3 <- readRDS("~/smlsamplingeu_louis/data/octobre_avec_rmqs2024_2025_sample3_lulucf.rds")
# saveRDS( sample3, "~/smlsamplingeu_louis/data/octobre_avec_rmqs2024_2025_sample3_CLC.rds")
# sample3 <- readRDS("~/smlsamplingeu_louis/data/octobre_avec_rmqs2024_2025_sample3_CLC.rds")
# en forçant la prise en compte de tous les points RMQS
......@@ -52,9 +52,11 @@ sample3 <- filter(sample3, !FR_WRB %in% excludeWRB )
# Suppose we have a regular grid we want to preserve in total
echant = "RMQS2024-25"
## (S1_CLC_JRC) ------
scenario = "S1_CLC_JRC"
resolution = "100m"
ListVarScenario <- c("FR_CLC",
"FR_regions",
"FR_WRB",
......@@ -63,7 +65,7 @@ ListVarScenario <- c("FR_CLC",
## 100m-------
domainsS1 <- 1e6 * stackb[["FR_CLC"]] +
domains <- 1e6 * stackb[["FR_CLC"]] +
1e4 * stackb[["code_supra"]] +
1e2 *stackb[["WRBLEV1"]] +
stackb[["climate_metzeger"]]
......@@ -74,26 +76,13 @@ sampleRand <- sample1[rnd.lhs$index_samples,]
source("dirsol/ComputeSampleOptim.R")
# création du masque pour avoir les pixels pris en compte
mask_values <- terra::ifel( domainsS1 %in% as.numeric(tab1$domains[tab1$n>10]),
1, 0
)
titreCartetot = paste(titreCarte,"\n ",scenario , "100m")
r = plotResuRMQS( mask_values,
sample7,
sampleRand ,
sample9 ,
titreCartetot ,
stack500)
tmap_save(r , file = paste0("dirsol/figures/",scenario,"100.jpeg") )
## 500m ----
domainsS1500 <- 1e6 * stack500[["FR_CLC"]] +
resolution = "500m"
domains <- 1e6 * stack500[["FR_CLC"]] +
1e4 * stack500[["code_supra"]] +
1e2 *stack500[["WRBLEV1"]] +
stack500[["climate_metzeger"]]
......@@ -104,24 +93,6 @@ sampleRand <- sample1500[rnd.lhs500$index_samples,]
source("dirsol/ComputeSampleOptim.R")
# création du masque pour avoir les pixels pris en compte
mask_values <- terra::ifel( domainsS1500 %in% as.numeric(tab1$domains[tab1$n>10]),
1, 0
)
titreCartetot = paste(titreCarte,"\n ",scenario , "500m")
r = plotResuRMQS( mask_values,
sample7,
sampleRand ,
sample9 ,
titreCartetot ,
stack500)
tmap_save(r , file = paste0("dirsol/figures/",scenario,"500.jpeg") )
# (S2_CLC_JRC) ---------
scenario = "S2_CLC_JRC"
ListVarScenario <- c("FR_CLC",
......@@ -129,8 +100,9 @@ ListVarScenario <- c("FR_CLC",
"FR_WRB")
## 100m------
resolution = "100m"
domainsS2 <- 1e4 * stackb[["FR_CLC"]] +
domains <- 1e4 * stackb[["FR_CLC"]] +
1e2 * stackb[["code_supra"]] +
stackb[["WRBLEV1"]]
......@@ -141,29 +113,12 @@ sampleRand <- sample1[rnd.lhs$index_samples,]
source("dirsol/ComputeSampleOptim.R")
# création du masque pour avoir les pixels pris en compte
mask_values <- terra::ifel( domainsS2 %in% as.numeric(tab1$domains[tab1$n>10]),
1, 0
)
titreCartetot = paste(titreCarte,"\n ",scenario , "100m")
r = plotResuRMQS( mask_values,
sample7,
sampleRand ,
sample9 ,
titreCartetot ,
stack500)
tmap_save(r , file = paste0("dirsol/figures/",scenario,"100.jpeg") )
##500m----
resolution = "500m"
domainsS2500 <- 1e4 * stack500[["FR_CLC"]] +
domains <- 1e4 * stack500[["FR_CLC"]] +
1e2 * stack500[["code_supra"]] +
stack500[["WRBLEV1"]]
......@@ -173,22 +128,6 @@ sampleRand <- sample1500[rnd.lhs500$index_samples,]
source("dirsol/ComputeSampleOptim.R")
# création du masque pour avoir les pixels pris en compte
mask_values <- terra::ifel( domainsS2500 %in% as.numeric(tab1$domains[tab1$n>10]),
1, 0
)
titreCartetot = paste(titreCarte,"\n ",scenario , "500m")
r = plotResuRMQS( mask_values,
sample7,
sampleRand ,
sample9 ,
titreCartetot ,
stack500)
tmap_save(r , file = paste0("dirsol/figures/",scenario,"500.jpeg") )
# (S1_CLC_IGCS) -----
......@@ -196,8 +135,9 @@ scenario = "S1_CLC_IGCS"
ListVarScenario <- c("FR_CLC", "FR_regions", "FR_IGCS", "FR_metzger")
### 100m-------
resolution = "100m"
domainsS3 <- 1e6 * stackb[["FR_CLC"]] + 1e4 * stackb[["code_supra"]] +
domains <- 1e6 * stackb[["FR_CLC"]] + 1e4 * stackb[["code_supra"]] +
1e2 *stackb[["ger"]] + stackb[["climate_metzeger"]]
sampleRand <- sample1[rnd.lhs$index_samples,]
......@@ -205,25 +145,10 @@ sampleRand <- sample1[rnd.lhs$index_samples,]
source("dirsol/ComputeSampleOptim.R")
# création du masque pour avoir les pixels pris en compte
mask_values <- terra::ifel( domainsS3 %in% as.numeric(tab1$domains[tab1$n>10]),
1, 0
)
titreCartetot = paste(titreCarte,"\n ",scenario , "100m")
r = plotResuRMQS( mask_values,
sample7,
sampleRand ,
sample9 ,
titreCartetot ,
stack500)
tmap_save(r , file = paste0("dirsol/figures/",scenario,"100.jpeg") )
## 500m-----
resolution = "500m"
domainsS3500 <- 1e6 * stack500[["FR_CLC"]] + 1e4 * stack500[["code_supra"]] +
domains <- 1e6 * stack500[["FR_CLC"]] + 1e4 * stack500[["code_supra"]] +
1e2 *stack500[["ger"]] + stack500[["climate_metzeger"]]
sampleRand <- sample1500[rnd.lhs500$index_samples,]
......@@ -231,30 +156,16 @@ sampleRand <- sample1500[rnd.lhs500$index_samples,]
source("dirsol/ComputeSampleOptim.R")
# création du masque pour avoir les pixels pris en compte
mask_values <- terra::ifel( domainsS3500 %in% as.numeric(tab1$domains[tab1$n>10]),
1, 0
)
titreCartetot = paste(titreCarte,"\n ",scenario , "500m")
r = plotResuRMQS( mask_values,
sample7,
sampleRand ,
sample9 ,
titreCartetot ,
stack500)
tmap_save(r , file = paste0("dirsol/figures/",scenario,"500.jpeg") )
# (S1_CLC_IGCS) ----------
# (S2_CLC_IGCS) ----------
scenario = "S2_CLC_IGCS"
ListVarScenario <- c("FR_CLC", "FR_regions", "FR_IGCS")
## 100m-------
resolution = "100m"
domainsS4 <- 1e4 * stackb[["FR_CLC"]] +
domains <- 1e4 * stackb[["FR_CLC"]] +
1e2 * stackb[["code_supra"]] +
stackb[["ger"]]
......@@ -267,25 +178,12 @@ sampleRand <- sample1[rnd.lhs$index_samples,]
source("dirsol/ComputeSampleOptim.R")
# création du masque pour avoir les pixels pris en compte
mask_values <- terra::ifel( domainsS3500 %in% as.numeric(tab1$domains[tab1$n>10]),
1, 0
)
titreCartetot = paste(titreCarte,"\n ",scenario , "100m")
r = plotResuRMQS( mask_values,
sample7,
sampleRand ,
sample9 ,
titreCartetot ,
stack500)
tmap_save(r , file = paste0("dirsol/figures/",scenario,"100.jpeg") )
## 500m------------
resolution = "500m"
domainsS4500 <- 1e4 * stack500[["FR_CLC"]] + 1e2 * stack500[["code_supra"]] +
domains <- 1e4 * stack500[["FR_CLC"]] +
1e2 * stack500[["code_supra"]] +
stack500[["ger"]]
sampleRand <- sample1500[rnd.lhs500$index_samples,]
......@@ -293,18 +191,3 @@ sampleRand <- sample1500[rnd.lhs500$index_samples,]
source("dirsol/ComputeSampleOptim.R")
# création du masque pour avoir les pixels pris en compte
mask_values <- terra::ifel( domainsS3500 %in% as.numeric(tab1$domains[tab1$n>10]),
1, 0
)
titreCartetot = paste(titreCarte,"\n ",scenario , "500m")
r = plotResuRMQS( mask_values,
sample7,
sampleRand ,
sample9 ,
titreCartetot ,
stack500)
tmap_save(r , file = paste0("dirsol/figures/",scenario,"500.jpeg") )
dirsol/Simul_DIRSOL_sans_RMQS_CLC.R
+ 32
− 22
View file @ c14cb441
......@@ -2,7 +2,7 @@
# Script pour faire tourner les simulations SANS RMQS, avec les données CLC----
# initialisations --------
source("dirsol/0_InitialisationsCLC.R")
source("dirsol/0_InitialisationsLULCF.R")
......@@ -92,9 +92,10 @@ mask_values <- terra::ifel( domainsS1500 %in% tab1$domains[tab1$n>10], 1, 0 )
plot(mask_values)
titreCartetot = paste(titreCarte,"\n",scenario ,resolution)
plotResu(testA1,
nom = "S1500_sansRmqs_lulucf_test_function",
nom = "dirsol/figures/S1_CLC_JRC500m",
crs = crs(stack500),
sample2 = sample2,
stack500 = stack500,
......@@ -105,6 +106,7 @@ plotResu(testA1,
## (S1_CLC_JRC) 100m ----------
# A JRC scenario with default layers (A)
# calcul d'de la grille representant les domaines
......@@ -188,7 +190,7 @@ mask_values <- terra::ifel( domainsS1%in% tab1100$domains[tab1100$n>10], 1, 0 )
plotResu(testA1100,
nom = "testA1100_lulucf_test_function_1millions",
nom = "dirsol/figures/S1_CLC_JRC100m",
crs = crs(stack500),
sample2 = sample2,
stack500 = stack500,
......@@ -196,10 +198,11 @@ plotResu(testA1100,
mask_values = mask_values
)
##(S2_CLC_JRC) 500metres --------------
##(S2_CLC_JRC) 500m --------------
# calcul d'de la grille representant les domaines
domainsS2500 <- 1e4 * stack500[["FR_CLC"]] + 1e2 * stack500[["code_supra"]] +
domainsS2500 <- 1e4 * stack500[["FR_CLC"]] +
1e2 * stack500[["code_supra"]] +
stack500[["WRBLEV1"]]
sample2 <- sample1500[rnd.lhs500$index_samples,]
......@@ -271,7 +274,7 @@ mask_values <- terra::ifel( domainsS2500 %in% tab2500$domains[tab2500$n>10], 1,
plotResu(testA2500,
nom = "testA2500_lulucf_test_function",
nom = "dirsol/figures/S2_CLC_JRC500m",
crs = crs(stack500),
sample2 = sample2,
stack500 = stack500,
......@@ -279,10 +282,11 @@ plotResu(testA2500,
mask_values = mask_values
)
## (S2_CLC_JRC) 100metres --------------
## (S2_CLC_JRC) 100m --------------
# calcul d'de la grille representant les domaines
domainsS2 <- 1e4 * stackb[["FR_CLC"]] + 1e2 * stackb[["code_supra"]] +
domainsS2 <- 1e4 * stackb[["FR_CLC"]] +
1e2 * stackb[["code_supra"]] +
stackb[["WRBLEV1"]]
# selection of the pixels using clhs
......@@ -356,7 +360,7 @@ saveRDS(testA2100, "/home/lpotel/smlsamplingeu_louis/export_data/S2_CLC_JRC_100m
mask_values <- terra::ifel( domainsS2 %in% tab2100$domains[tab2100$n>10], 1, 0 )
plotResu(testA2100,
nom = "testA2100_lulucf_test_function",
nom = "dirsol/figures/S2_CLC_JRC100m",
crs = crs(stack500),
sample2 = sample2,
stack500 = stack500,
......@@ -367,8 +371,10 @@ plotResu(testA2100,
## (S1_CLC_IGCS) 500m --------------
# calcul d'de la grille representant les domaines
domainsS3500 <- 1e6 * stack500[["FR_CLC"]] + 1e4 * stack500[["code_supra"]] +
1e2 *stack500[["ger"]] + stack500[["climate_metzeger"]]
domainsS3500 <- 1e6 * stack500[["FR_CLC"]] +
1e4 * stack500[["code_supra"]] +
1e2 *stack500[["ger"]] +
stack500[["climate_metzeger"]]
# selection of the pixels using clhs
sample2 <- sample1500[rnd.lhs500$index_samples,]
......@@ -439,7 +445,7 @@ mask_values <- terra::ifel( domainsS3500 %in% tab3500$domains[tab3500$n>10], 1,
plotResu(testA3500,
nom = "test3500_lulucf_test_function",
nom = "dirsol/figures/S1_CLC_IGCS500m",
crs = crs(stack500),
sample2 = sample2,
stack500 = stack500,
......@@ -449,8 +455,10 @@ plotResu(testA3500,
## (S1_CLC_IGCS) 100m --------------
# calcul d'de la grille representant les domaines
domainsS3 <- 1e6 * stackb[["FR_CLC"]] + 1e4 * stackb[["code_supra"]] +
1e2 *stackb[["ger"]] + stackb[["climate_metzeger"]]
domainsS3 <- 1e6 * stackb[["FR_CLC"]] +
1e4 * stackb[["code_supra"]] +
1e2 *stackb[["ger"]] +
stackb[["climate_metzeger"]]
# selection of the pixels using clhs
sample2 <- sample1[rnd.lhs$index_samples,]
......@@ -515,7 +523,7 @@ testA3100 = optimisesampleSML(
cv = cv
)
saveRDS(testA3100, "/home/lpotel/smlsamplingeu_louis/export_data/S1_CLC_IGCS_500m.rds" )
saveRDS(testA3100, "/home/lpotel/smlsamplingeu_louis/export_data/S1_CLC_IGCS_100m.rds" )
......@@ -524,7 +532,7 @@ mask_values <- terra::ifel( domainsS3%in% tab3100$domains[tab3100$n>10], 1, 0 )
plot(mask_values)
plotResu(testA3100,
nom = "test3100_lulucf_test_function",
nom = "dirsol/figures/S1_CLC_IGCS100m",
crs = crs(stack500),
sample2 = sample2,
stack500 = stack500,
......@@ -532,10 +540,11 @@ plotResu(testA3100,
mask_values = mask_values
)
## (S2_CLC_IGCS) 500 --------------
## (S2_CLC_IGCS) 500m --------------
# calcul d'de la grille representant les domaines
domainsS4500 <- 1e4 * stack500[["FR_CLC"]] + 1e2 * stack500[["code_supra"]] +
domainsS4500 <- 1e4 * stack500[["FR_CLC"]] +
1e2 * stack500[["code_supra"]] +
stack500[["ger"]]
......@@ -609,17 +618,18 @@ mask_values <- terra::ifel( domainsS4500 %in% tab4500$domains[tab4500$n>10], 1,
plot(mask_values)
plotResu(testA4500,
nom = "test4500_lulucf_test_function",
nom = "dirsol/figures/S2_CLC_IGCS500m",
crs = crs(stack500),
sample2 = sample2,
stack500 = stack500,
sample3 = sample3,
mask_values = mask_values
)
## (S2_CLC_IGCS) 100 --------------
## (S2_CLC_IGCS) 100m --------------
# calcul d'de la grille representant les domaines
domainsS4 <- 1e4 * stackb[["FR_CLC"]] + 1e2 * stackb[["code_supra"]] +
domainsS4 <- 1e4 * stackb[["FR_CLC"]] +
1e2 * stackb[["code_supra"]] +
stackb[["ger"]]
......@@ -699,7 +709,7 @@ mask_values <- terra::ifel( domainsS4 %in% tab4100$domains[tab4100$n>10], 1, 0 )
plot(mask_values)
plotResu(testA4100,
nom = "test4100_lulucf_test_function",
nom = "dirsol/figures/S2_CLC_IGCS100m",
crs = crs(stack500),
sample2 = sample2,
stack500 = stack500,
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment