tests/testthat/test_wtss.R

In wtss: An R Client for Web Time-Series Service

context("WTSS service")
test_that("Connection to a WTSS service", {
    URL <- "http://www.esensing.dpi.inpe.br/wtss/"
    wtss1 <-  wtss::WTSS(URL) 
    expect_true(wtss1$url == URL)
    coverages <- wtss1$coverages
    expect_true("MOD13Q1" %in% coverages)
})
test_that("Bad connection to a WTSS service", {
    URL <- "http://www.dpi.inpe.br/wtss/"
    expect_message(wtss1 <- wtss::WTSS(URL), 
                   "WTSS server not responding - please check URL")
    expect_true(purrr::is_null(wtss1))
})

test_that("List coverages", {
    URL   <- "http://www.esensing.dpi.inpe.br/wtss/"
    wtss2 <-  wtss::WTSS(URL) 
    output <- capture.output(wtss::list_coverages(wtss2))
    expect_true(as.logical(grep(wtss2$url, output[2])))
    expect_true(as.logical(grep(wtss2$coverages[1], output[3])))
})

test_that("Describe coverage", {
    URL   <- "http://www.esensing.dpi.inpe.br/wtss/"
    wtss3 <-  wtss::WTSS(URL) 
    output <- capture.output(wtss::describe_coverage(wtss3, wtss3$coverages[1]))
    expect_true(as.logical(grep("satellite", output[5])))
    expect_true(as.logical(grep("minimum_values", output[15])))
    expect_true(as.logical(grep("maximum_values", output[20])))
    expect_true(as.logical(grep("Timeline", output[29])))
})

test_that("Describe coverage", {
    URL   <- "http://www.esensing.dpi.inpe.br/wtss/"
    wtss4 <-  wtss::WTSS(URL) 
    output <- capture.output(wtss::describe_coverage(wtss4, wtss4$coverages[1]))
    expect_true(nrow(wtss4$description) == 1)
    output <- capture.output(wtss::describe_coverage(wtss4, wtss4$coverages[1]))
    expect_true(nrow(wtss4$description) == 1)
    output <- capture.output(wtss::describe_coverage(wtss4, wtss4$coverages[2]))
    expect_true(nrow(wtss4$description) == 2)
    output <- capture.output(wtss::describe_coverage(wtss4, wtss4$coverages[2]))
    expect_true(nrow(wtss4$description) == 2)
})

test_that("Time Series", {
    wtss5 <- wtss::WTSS("http://www.esensing.dpi.inpe.br/wtss/")
    ts    <- wtss::time_series(wtss5, "MOD13Q1", c("ndvi","evi"), 
                     longitude = -45.00, latitude  = -12.00,
                     start_date = "2000-02-18", end_date = "2016-12-18")
    expect_true(nrow(ts$time_series[[1]]) == 388)
    expect_true(ncol(ts$time_series[[1]]) == 3)
})

test_that("Time Series 2", {
    wtss6 <- wtss::WTSS("http://www.esensing.dpi.inpe.br/wtss/")
    ts    <- wtss::time_series(wtss6, "MOD13Q1", 
                               longitude = -45.00, latitude  = -12.00)
    expect_true(ncol(ts$time_series[[1]]) == 7)
    expect_true(all(c("Index", "mir", "blue", "nir", "red", "evi", "ndvi") 
                    %in% names(ts$time_series[[1]])))
    expect_true(nrow(ts$time_series[[1]]) >= 452)
})

test_that("Time Series - errors", {
    wtss7 <- wtss::WTSS("http://www.esensing.dpi.inpe.br/wtss/")
    expect_error(ts <- wtss::time_series(wtss7, "MOD13Q1", 
                                         longitude = 45.00, latitude  = -12.00))
    expect_error(ts <- wtss::time_series(wtss7, "MOD13Q1", 
                                         longitude = -45.00, latitude  = 12.00))
    expect_error(ts <- wtss::time_series(wtss7, "MOD13Q1", c("swir", "ndvi"),
                                        longitude = -45.00, latitude  = -12.00))
    expect_error(ts <- wtss::time_series(wtss7, "MOD13Q1",
                                         longitude = -45.00, latitude  = -12.00,
                                         start_date = "1999-01-01"))    
    expect_error(ts <- wtss::time_series(wtss7, "MOD13Q1",
                                         longitude = -45.00, latitude  = -12.00,
                                         end_date = "2045-01-01")) 
    
    expect_error(ts <- wtss::time_series(wtss7, "MOD13Q1",
                                         longitude = -45.00, latitude  = -12.00,
                                         start_date = "2010-01-01", 
                                         end_date   = "2005-01-01"))

})

test_that("Time Series - conversion to ts and zoo", {
    wtss8 <- wtss::WTSS("http://www.esensing.dpi.inpe.br/wtss/")
    ts    <- wtss::time_series(wtss8, "MOD13Q1", 
                               longitude = -45.00, latitude  = -12.00,
                             start_date = "2002-01-01", end_date = "2002-12-31")
    ts_zoo <- wtss::wtss_to_zoo(ts)
    expect_true(nrow(ts_zoo) == nrow(ts$time_series[[1]]))
    expect_true(ncol(ts_zoo) == (ncol(ts$time_series[[1]]) - 1))
    expect_true(all(as.vector(ts_zoo[,1]) == 
                        dplyr::pull(ts$time_series[[1]][,2])))
    
    ts_ts <- wtss::wtss_to_ts(ts, band = "ndvi")
    
    ts_start <- c(as.numeric(lubridate::year(ts$start_date)), 
                  as.numeric(lubridate::week(ts$start_date)))
    ts_end   <- c(as.numeric(lubridate::year(ts$end_date)), 
                  as.numeric(lubridate::week(ts$end_date)))
    
    expect_true(all(stats::start(ts_ts) == ts_start))
    expect_true(length(ts_ts) == 53)
    expect_true(ts_ts[1] == as.numeric(ts_zoo[1,"ndvi"]))
})