inst/examples/compare_sentinel_modis/new_scripts/02_get_time_series/create_samples_3_labels.R

In albhasan/sits.starfm: Use the StarFM Fusion Model on Landsat 8 Images

#!/usr/bin/Rscript

# Recode the labels in the samples files from 5 to 3.

suppressMessages(library(dplyr))

base_dir <- "/home/alber/Documents/data/experiments/prodes_reproduction/papers/deforestation/data/validation"

# Read samples
recode_to_3_labels <- function(rds_file){
    rds_file %>%
        readRDS() %>%
        dplyr::filter(label != "Water") %>%
        dplyr::mutate(label = dplyr::recode(label,
                                            #"Forest"
                                            #"Deforestatio"
                                            #NonForest
                                            Pasture      = "NonForest",
                                            NatNonForest = "NonForest")) %>%
        ensurer::ensure_that(nrow(.) > 0, err_desc = "No rows!") %>%
        ensurer::ensure_that(length(unique(.$label)) == 3,
                             err_desc = "Label missmatch!") %>%
        return()
}

base_dir %>%
    file.path("samples_A_approx.rds") %>%
    recode_to_3_labels() %>%
    saveRDS(file.path(base_dir, "samples_A_approx_3l.rds"))

base_dir %>%
    file.path("samples_B_approx.rds") %>%
    recode_to_3_labels() %>%
    saveRDS(file.path(base_dir, "samples_B_approx_3l.rds"))

base_dir %>%
    file.path("samples_C_approx.rds") %>%
    recode_to_3_labels() %>%
    saveRDS(file.path(base_dir, "samples_C_approx_3l.rds"))

base_dir %>%
    file.path("samples_A_raw.rds") %>%
    recode_to_3_labels() %>%
    saveRDS(file.path(base_dir, "samples_A_raw_3l.rds"))

base_dir %>%
    file.path("samples_B_raw.rds") %>%
    recode_to_3_labels() %>%
    saveRDS(file.path(base_dir, "samples_B_raw_3l.rds"))

base_dir %>%
    file.path("samples_C_raw.rds") %>%
    recode_to_3_labels() %>%
    saveRDS(file.path(base_dir, "samples_C_raw_3l.rds"))