inst/template/02_3_pre_processing_cropland/select_esacci.R

In michielvandijk/mapspamc: An R package to create crop distribution maps

#'========================================================================================================================================
#' Project:  mapspamc
#' Subject:  Code to select ESACCI cropland map
#' Author:   Michiel van Dijk
#' Contact:  michiel.vandijk@wur.nl
#'========================================================================================================================================

# SOURCE PARAMETERS ----------------------------------------------------------------------
source(here::here("01_model_setup/01_model_setup.r"))


# LOAD DATA ------------------------------------------------------------------------------
load_data(c("adm_map", "grid"), param)


# PROCESS --------------------------------------------------------------------------------
temp_path <- file.path(param$model_path, glue("processed_data/maps/cropland/{param$res}"))
dir.create(temp_path, showWarnings = FALSE, recursive = TRUE)

# Warp and mask
# If needed change the year of the cropland map
input <- file.path(param$db_path, glue("esacci/esacci_{param$year}.tif"))
output <- align_raster(input, grid, adm_map, method = "bilinear")
names(output) <- "esacci"
plot(output)
writeRaster(output, file.path(temp_path, glue("cropland_esacci_{param$res}_{param$year}_{param$iso3c}.tif")),
            overwrite = TRUE)


# CLEAN UP -------------------------------------------------------------------------------
rm(input, output, grid, adm_map, temp_path)