Nothing
context("WTSS service")
test_that("List coverages", {
vcr::use_cassette("list_coverages", {
wtss2 <- "http://www.esensing.dpi.inpe.br/wtss/"
output <- capture.output(Rwtss::list_coverages(wtss2))
if (!purrr::is_null(output)) {
expect_true(as.logical(grep("MOD13Q1", output[1])))
}
})
})
test_that("Describe coverage", {
vcr::use_cassette("describe_coverage", {
wtss3 <- "http://www.esensing.dpi.inpe.br/wtss/"
output <- capture.output(Rwtss::describe_coverage(wtss3, "MOD13Q1"))
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("Time Series", {
vcr::use_cassette("time_series", {
wtss5 <- "http://www.esensing.dpi.inpe.br/wtss/"
ts <- Rwtss::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", {
vcr::use_cassette("time_series_2", {
wtss6 <- "http://www.esensing.dpi.inpe.br/wtss/"
ts <- Rwtss::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("Plot",{
vcr::use_cassette("plot", {
wtssg <- "http://www.esensing.dpi.inpe.br/wtss/"
tsg <- Rwtss::time_series(wtssg, "MOD13Q1",
longitude = -45.00,
latitude = -12.00)
g <- plot(tsg)
expect_true(unname(summary(g)[1,2]) == "gg")
})
})
test_that("Time Series - errors", {
vcr::use_cassette("time_series_errors", {
wtss7 <- "http://www.esensing.dpi.inpe.br/wtss/"
expect_message(ts <- Rwtss::time_series(wtss7, "MOD13Q1",
longitude = 45.00,
latitude = -12.00))
expect_true(purrr::is_null(ts))
expect_message(ts <- Rwtss::time_series(wtss7, "MOD13Q1",
longitude = -45.00,
latitude = 12.00))
expect_true(purrr::is_null(ts))
expect_message(ts <- Rwtss::time_series(wtss7, "MOD13Q1",
c("swir", "ndvi"),
longitude = -45.00,
latitude = -12.00))
expect_true(purrr::is_null(ts))
expect_message(ts <- Rwtss::time_series(wtss7, "MOD13Q1",
longitude = -45.00,
latitude = -12.00,
start_date = "1999-01-01"))
expect_true(purrr::is_null(ts))
expect_message(ts <- Rwtss::time_series(wtss7, "MOD13Q1",
longitude = -45.00,
latitude = -12.00,
end_date = "2045-01-01"))
expect_true(purrr::is_null(ts))
expect_message(ts <- Rwtss::time_series(wtss7, "MOD13Q1",
longitude = -45.00,
latitude = -12.00,
start_date = "2010-01-01",
end_date = "2005-01-01"))
expect_true(purrr::is_null(ts))
expect_message(ts <- Rwtss::time_series(wtss7, "MOD13Q12",
longitude = -45.00,
latitude = -12.00))
expect_true(!purrr::is_null(ts))
})
})
test_that("Time Series - conversion to ts and zoo", {
vcr::use_cassette("time_series_conversion_to_ts_and_zoo", {
wtss8 <- "http://www.esensing.dpi.inpe.br/wtss/"
ts <- Rwtss::time_series(wtss8, "MOD13Q1",
longitude = -45.00, latitude = -12.00)
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)))
ts_zoo <- Rwtss::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])))
ts1 <- Rwtss::time_series(wtss8, "MOD13Q1",
longitude = -45.00, latitude = -12.50)
ts2 <- dplyr::bind_rows(ts,ts1)
expect_message(Rwtss::wtss_to_zoo(ts2))
ts_week <- Rwtss::wtss_to_ts(ts, band = "ndvi", period = "week")
ts_weeks <- Rwtss::wtss_to_ts(ts, band = "ndvi", period = "weeks")
ts_week_52 <- Rwtss::wtss_to_ts(ts, band = "ndvi", period = 52)
expect_true(all(stats::start(ts_week) == ts_start))
expect_true(ts_week[1] == as.numeric(ts_zoo[1, "ndvi"]))
expect_true(all(ts_week == ts_weeks))
expect_true(all(ts_week == ts_week_52))
ts_day <- Rwtss::wtss_to_ts(ts, band = "ndvi", period = "day")
ts_days <- Rwtss::wtss_to_ts(ts, band = "ndvi", period = "days")
ts_day_1 <- Rwtss::wtss_to_ts(ts, band = "ndvi", period = 365)
expect_true(all(stats::start(ts_day) == ts_start))
expect_true(ts_day[1] == as.numeric(ts_zoo[1, "ndvi"]))
expect_true(all(ts_day == ts_days))
expect_true(all(ts_day == ts_day_1))
ts_month <- Rwtss::wtss_to_ts(ts, band = "ndvi", period = "month")
ts_months <- Rwtss::wtss_to_ts(ts, band = "ndvi", period = "months")
ts_month_12 <- Rwtss::wtss_to_ts(ts, band = "ndvi", period = 12)
expect_true(all(ts_month == ts_months))
expect_true(all(ts_month == ts_month_12))
expect_true(all(stats::start(ts_month) == ts_start))
expect_true(ts_month[1] == as.numeric(ts_zoo[1, "ndvi"]))
expect_message(Rwtss::wtss_to_ts(ts, band = "ndvi", period = "year"))
expect_message(Rwtss::wtss_to_ts(ts, band = c("ndvi", "evi")))
expect_message(Rwtss::wtss_to_ts(ts))
expect_message(ts2_ts <- Rwtss::wtss_to_ts(ts2, band = "ndvi"))
})
})
test_that("Bad connection to a WTSS service", {
wtss10 <- "http://www.dpi.inpe.br2/wtss/"
output <- suppressMessages(Rwtss::list_coverages(wtss10))
expect_true(purrr::is_null(output))
output <- suppressMessages(Rwtss::describe_coverage(wtss10, "MOD"))
expect_true(purrr::is_null(output))
ts <- suppressMessages(Rwtss::time_series(wtss10, "MOD13Q1",
longitude = 45.00,
latitude = -12.00))
expect_true(purrr::is_null(ts))
})
test_that("Guess satellite",{
sat_sensor <- Rwtss:::.wtss_guess_satellite(0.002)
expect_equal(unname(sat_sensor[1]), "TERRA")
expect_equal(unname(sat_sensor[2]), "MODIS")
sat_sensor <- Rwtss:::.wtss_guess_satellite(10)
expect_equal(unname(sat_sensor[1]), "SENTINEL-2")
expect_equal(unname(sat_sensor[2]), "MSI")
sat_sensor <- Rwtss:::.wtss_guess_satellite(70)
expect_equal(unname(sat_sensor[1]), "CBERS")
expect_equal(unname(sat_sensor[2]), "AWFI")
sat_sensor <- Rwtss:::.wtss_guess_satellite(30)
expect_equal(unname(sat_sensor[1]), "LANDSAT")
expect_equal(unname(sat_sensor[2]), "OLI")
sat_sensor <- Rwtss:::.wtss_guess_satellite(5)
expect_equal(unname(sat_sensor[1]), "UNKNOWN")
expect_equal(unname(sat_sensor[2]), "UNKNOWN")
})
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.