tests/testthat/test-season.R

In phenofit: Extract Remote Sensing Vegetation Phenology

context("season")

source('helper_MOD13A1.R')
# source('tests/testthat/helper_MOD13A1.R')
lambda   <- init_lambda(INPUT$y) # lambda for whittaker

# `season` produce a pdf file here
param = listk(
    INPUT,
    rFUN = smooth_wWHIT, wFUN = wBisquare, iters = 2,
    lambda,
    rtrough_max = 0.6, ypeak_min = ypeak_min,
    MaxPeaksPerYear =2.5, MaxTroughsPerYear = 3.5
)


test_that("`season` with smooth_wWHIT and calendarYear", {
    param$calendarYear <- TRUE
    expect_silent(brks <- do.call(season, param))
})

test_that("`season` with smooth_wWHIT", {
    param$calendarYear <- FALSE
    expect_silent(brks <- do.call(season, param))
})

test_that("`season` with wHANTS", {
    param$rFUN <- smooth_wHANTS
    expect_silent(brks <- do.call(season, param))
})

test_that("`season` with smooth_wSG", {
    param$rFUN <- smooth_wSG
    expect_silent(brks <- do.call(season, param))
})