context("split_emis")
data(net)
net <- sf::st_transform(net, 31983)
g <- make_grid(net, 500) #500m in degrees
netsf <- sf::st_as_sf(net)[, "ldv"]
x <- split_emis(netsf, g)
dim(x)
test_that("split_emis works", {
expect_equal(round(as.numeric(sum(x$ldv))), 1946955)
})
g$A <- rep(letters, length = 20)[1:nrow(g)]
g$B <- rev(g$A)
netsf <- sf::st_as_sf(net)[, c("ldv", "hdv")]
test_that("split_emis works", {
expect_equal(round(as.numeric(sum(split_emis(netsf,
g,
add_column = c("A", "B"))$ldv
))), 1946955)
})
test_that("split_emis message", {
expect_equal(round(as.numeric(split_emis(net = netsf, distance = 1000, verbose = TRUE)$ldv[1])),
1603)
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.