library("testthat")
library("mapi")
test_that("MAPI_Varicell", {
    keep <- c(2,3,6,9,10,11,16,18,19,20,21,23,26,27,29,31,33,34,38,41,46,50,54,58,61,63,65,71,
    72,73,76,78,79,81,85,92,93,94,97,98,99,101,103,113,115,119,120,121,124,127,130,134,142,143,
    151,152,159,160,176,185,189,191,195,196,197,198,199) # for reproductibility
    samples2 <- mapi::samples[keep,]
    grid <- MAPI_GridHexagonal(samples2, crs=3857, hw=250)
    expect_equal(nrow(grid), 1045)
    expect_equal(as.integer(sum(sf::st_area(grid))), 169686852)
    set.seed(1234) # for reproductibility
    grid.var <- MAPI_Varicell(grid, samples2, buf=250, var.coef=20)
    expect_equal(nrow(grid.var), 173)
    expect_equal(as.integer(sum(sf::st_area(grid.var))), 189126625)
    expect_equal(sd(sf::st_area(grid.var)), 963339.563502512)
})