tests/testthat/test-tsf_process.R

In kongdd/rTIMESAT: Extract Remote Sensing Vegetation Phenology by TIMESAT Fortran library

context("test-tsf_process")


test_that("write_input works", {
    expect_silent({
      write_input(d$EVI/1e4 , file_y, nptperyear)
      write_input(d$SummaryQA, file_qc, nptperyear)
    })
})


test_that("write_setting works", {
    expect_silent({
        options <- update_setting(options)
        write_setting(options, file_set)
    })
})

test_that("TSF_main works", {
    nptperyear = 23
    options <- list(
        ylu         = c(0, 9999),     # Valid data range (lower upper)
        qc_1        = c(0, 0, 1),     # Quality range 1 and weight
        qc_2        = c(1, 1, 0.5),   # Quality range 2 and weight
        qc_3        = c(2, 3, 0.2),   # Quality range 3 and weight
        A           = 0.1,            # Amplitude cutoff value
        output_type = c(1, 1, 0),     # Output files (1/0 1/0 1/0), 1: seasonality data; 2: smoothed time-series; 3: original time-series
        seasonpar   = 1.0,            # Seasonality parameter (0-1)
        iters       = 2,              # No. of envelope iterations (3/2/1)
        FUN         = 2,              # Fitting method (1/2/3): (SG/AG/DL)
        half_win    = 7,              # half Window size for Sav-Gol.
        meth_pheno  = 1,              # (1: seasonal amplitude, 2: absolute value, 3: relative amplitude, 4: STL trend)
        trs         = c(0.5, 0.5)     # Season start / end values
    )
    
    data("MOD13A1")
    sitename <- "US-KS2"
    # sitename <- "CA-NS6"

    d <- subset(MOD13A1$dt, date >= as.Date("2004-01-01") & date <= as.Date("2010-12-31") & site == sitename)
    skip_on_os(os = c("mac", "linux", "solaris"))
    r <- TSF_main(
        y = d$EVI / 1e4, qc = d$SummaryQA, nptperyear,
        jobname = sitename, options, cache = FALSE
    )
    expect_true(mean(r$pheno$time_peak) > 40)
})