library("testthat") library("mapi") test_that("MAPI_GridHexagonal", { library("data.table") library("sf") data("samples") grid <- MAPI_GridHexagonal(samples, crs=3857, hw=250) expect_equal(nrow(grid), 1071) samplesPlus <- data.table::data.table(ind=samples$ind, x=samples$x*1000, y=samples$y*1000, errRad=samples$errRad*10) samples.ll <- sf::st_transform(sf::st_as_sf(samplesPlus, crs=3857, coords=c("x","y")), crs=4326) samples4326 <- cbind(sf::st_drop_geometry(samples.ll), sf::st_coordinates(samples.ll)) data.table::setnames(samples4326, "X", "x") data.table::setnames(samples4326, "Y", "y") grid.ll <- MAPI_GridHexagonal(samples4326, crs=4326, hw=250000) expect_equal(nrow(grid.ll), 2432) })