## load packages
library("testthat")
library("neotoma2")

context("Run Neotoma `test_sites` only when not on CRAN")

test_that("get_sites runs as expected.", {
  ## we don't want this to run on CRAN
  skip_on_cran()

  alexander_lake <- get_sites(sitename = "Alexander Lake")

  # Check that sites' names are only Alexander Lake

  expect_identical(alexander_lake[[1]]@sitename, "Alexander Lake")

})

test_that("get_sites runs as expected.", {
  ## we don't want this to run on CRAN
  skip_on_cran()

  alex <- get_sites(sitename = "Alex%")

  # Check that sites' names are only Alexander Lake

  for (i in seq_len(length(alex))) {
    alex_char <- substring(alex[[i]]@sitename, 1,4)
    alex_char <- tolower(alex_char)
    expect_identical(alex_char, "alex")
  }

})

test_that("get_sites runs as expected.", {
  ## we don't want this to run on CRAN
  skip_on_cran()

  site.1001 <- get_sites(1001)

  # Check that sites' id matches 1001

  expect_equal(site.1001[[1]]@siteid, 1001)

})

test_that("get_sites runs as expected.", {
  ## we don't want this to run on CRAN
  skip_on_cran()

  sites.ob <- get_sites(c(1001, 2001, 15, 24))
  sites.vec <- getids(sites.ob)
  sites.vec <- unique(sites.vec$siteid)

  # Check that sites' id matches 1001

  expect_setequal(sites.vec, c(1001, 2001, 15, 24))

})

test_that("All Czech sites work with different spatial bounds:", {
  skip_on_cran()
  cz_json <- '{"type": "Polygon",
        "coordinates": [[
            [12.40, 50.14],
            [14.10, 48.64],
            [16.95, 48.66],
            [18.91, 49.61],
            [15.24, 50.99],
            [12.40, 50.14]]]}'
  cz_WKT <- 'POLYGON ((12.4 50.14,
                       14.1 48.64,
                       16.95 48.66,
                       18.91 49.61,
                       15.24 50.99,
                       12.4 50.14))'
  cz_bbox <- c(12.4, 48.64, 18.91, 50.99)

  testthat::expect_true(all.equal(get_sites(loc=cz_json),
    get_sites(loc=cz_WKT)))

  # Now, we know that all sites in cz_sites[[1]] should be in cz_sites[[3]],
  # but the bounding box strategy means that the reverse is not true:
  cz_ids <- getids(get_sites(loc=cz_json[1]))
  testthat::expect_true(all(cz_ids$siteid %in%
    getids(get_sites(loc = cz_bbox, limit = 50))$siteid))
})

test_that("All Data + loc work", {
  skip_on_cran()
  europe_json <- '{"type": "Polygon",
        "coordinates": [
          [[-23.5546875, 70.8446726342528],
           [-25.3125, 39.36827914916014],
           [61.87499999999999, 39.36827914916014],
           [62.57812500000001, 74.01954331150228],
           [-23.5546875, 70.8446726342528]]
        ]}'
  data <- get_sites(loc = europe_json[1], all_data = TRUE)
  testthat::expect_gt(length(data), 50)

  # Now, we know that all sites in cz_sites[[1]] should be in cz_sites[[3]],
  # but the bounding box strategy means that the reverse is not true:
})