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)
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.