scripts/landcover.R

In IVFL-BOKU/landsupport: Automates data processing for the Landsupport project

args = commandArgs(TRUE)
if (length(args) < 4) {
  stop('This scripts takes parameters: settingsFilePath regionId dateFrom dateTo')
}
names(args) = c('cfgFile', 'region', 'from', 'to')
t0 = Sys.time()
cat(paste0(c('Running landcover.R', args, as.character(Sys.time()), '\n'), collapse = '\t'))
source(args[1])

devtools::load_all(cubeRpath, quiet = TRUE)
library(sentinel2, quietly = TRUE)
library(dplyr, quietly = TRUE, warn.conflicts = FALSE)
library(mlr3, quietly = TRUE)
library(mlr3learners, quietly = TRUE)

data = suppressMessages(readr::read_csv(landcoverTrainDataPath))
period = paste0(substr(args['from'], 1, 4), 'y1')
modelFile = prepareLandcoverModel(modelsDir, tmpDir, landcoverModelName, data, landcoverTargetVar, landcoverFeaturesList, period, nCores, landcoverSkipExistingModel)
rm(data)

images = getCache(args['region'], args['from'], args['to'], args['cfgFile']) %>%
  dplyr::filter(.data$band == 'SCL') %>%
  dplyr::mutate(period = substr(.data$date, 1, 4)) %>%
  dplyr::group_by(utm, period) %>%
  dplyr::filter(dplyr::row_number() == 1) %>%
  dplyr::rename(tileFile = file) %>%
  dplyr::ungroup()
regionFile = getCachePath(cacheTmpl, args['region'], args['from'], args['to'], cloudCov, bands, 'geojson')
tiles = suppressMessages(mapTilesGrid(images, gridFile, regionFile)) %>%
  dplyr::select(.data$tile, .data$period) %>%
  dplyr::mutate(period = paste0(.data$period, 'y1'))
imgs = prepareLandcover(tiles, tilesDir, tmpDir, modelFile, landcoverModelName, nCores, 1000000, landcoverSkipExisting, landcoverGdalOpts)
logProcessingResults(imgs, t0)