#Reclass SSP2 results for comparison with ESA
require(CLUvalid)
#define scenario
scenario <- "Scenario_SSP2"
outpath <- "C:/Bep/Projects/CLUMondo/SSP Validation/Data prep/Output/Maps/Scenario_SSP2/"
LandCover <- c("forest", "urb", "crop")
#create the reclass matrices
u <- c(0, 21, 0, 21, 22, 1, 22, 24, 2)
urb_rcl <- matrix(u, ncol = 3, byrow = TRUE)
c <- c(0, 5, 2, 5, 12, 1, 12, 24, 0)
crop_rcl <- matrix(c, ncol = 3, byrow = TRUE)
f <- c(0, 12, 0, 12, 13, 2, 13, 14, 1, 14, 24, 0)
forest_rcl <- matrix(f, ncol = 3, byrow = TRUE)
for(e in 1:length(LandCover)){
LC <- LandCover[e]
CLU_folders <- list_CLU_folders(inpath = "U:/CLUMondo/SSP Implementation/Scenario_", scenario = "SSP2")
for(i in 1:length(CLU_folders)){
region <- CLU_folders[i]
region_name <- basename(region)
region2000 <- raster(paste0(region, "cov_all.0.asc"))
region2015 <- raster(paste0(region, "cov_all.15.asc"))
rcl_LC <- get(paste0(LC, "_rcl"))
reclass2000 <- reclassify(region2000, rcl_LC, filename = paste0(outpath, scenario, "_", LC, "/Reclassified_", LC, "/2000/", region_name, "_", LC, "_rcl_2000.asc"), overwrite = TRUE)
reclass2015 <- reclassify(region2015, rcl_LC, filename = paste0(outpath, scenario, "_", LC, "/Reclassified_", LC, "/2015/", region_name, "_", LC, "_rcl_2015.asc"), overwrite = TRUE)
reclass2000[reclass2000 == 1 | reclass2000 == 2] <- 1
reclass2015[reclass2015 == 1 | reclass2015 == 2] <- 1
writeRaster(reclass2000, filename = paste0(outpath, scenario, "_", LC, "/Reclassified_", LC, "/2000/", region_name, "_", LC, "_bool_2000.asc"), overwrite = TRUE)
writeRaster(reclass2015, paste0(outpath, scenario, "_", LC, "/Reclassified_", LC, "/2015/", region_name, "_", LC, "_bool_2015.asc"), overwrite = TRUE)
}
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.