test_that("esa global landcover works", {
  shp <- read_sf(
    system.file("extdata", "sierra_de_neiba_478140.gpkg",
      package = "mapme.biodiversity"
    )
  )
  esalandcover <- list.files(system.file("res", "esalandcover",
    package = "mapme.biodiversity"
  ), pattern = ".tif$", full.names = TRUE)
  esalandcover <- rast(esalandcover)
  attributes(shp)$years <- 2015:2020
  result <- .calc_landcover(shp, esalandcover)
  expect_equal(
    names(result),
    c("classes", "year", "area", "percentage")
  )
  expect_equal(
    unique(result$year),
    c("2015", "2016", "2017", "2018", "2019")
  )
  # expect_snapshot(
  #   result$area
  # )
})