tests/testthat/test-aggregate-polys.R

In mapscanner: Print Maps, Draw on Them, Scan Them Back in

context("ms_aggregate_polys")

test_that ("aggregate", {
               geometry <- list (sf::st_point (cbind (0, 0)),
                                 sf::st_point (cbind (0, 1)),
                                 sf::st_point (cbind (1, 0)))
               pts <- sf::st_sf (a = 1:3, geometry = geometry)
               overlapping_polys <- sf::st_buffer (pts, 0.75)

               requireNamespace ("dplyr") # produces output
               expect_silent (x <- ms_aggregate_polys (overlapping_polys))
               expect_is (x, "sf")
               expect_is (x$geometry, "sfc_POLYGON")
               expect_equal (nrow (x), 3)
               expect_equivalent(1:3, x$n)
})