tests/testthat/test-join.R

In aarkilanian/dci: Calculate the Dendritic Connectivity in River Networks

test_that("Nodes within tolerance are joined", {

  # Create test river network
  rivers <- sf::st_as_sf(sf::st_sfc(
    sf::st_linestring(matrix(c(2,1,5,5), 2)),
    sf::st_linestring(matrix(c(3,2,5,5), 2)),
    sf::st_linestring(matrix(c(4,3,4,4), 2)),
    sf::st_linestring(matrix(c(3,3,3,4), 2)),
    sf::st_linestring(matrix(c(3,3,2,3), 2)),
    sf::st_linestring(matrix(c(3,3,1,2), 2)),
    sf::st_linestring(matrix(c(3,3,4,5), 2))
  )) %>%
    dplyr::rename("geometry" = "x") %>%
    sf::st_as_sf(wkt = "geometry")
  rivers <- sf::st_as_sf(rivers, wkt = "x")
  net <- sfnetworks::as_sfnetwork(rivers)

  # Create test nodes
  nodes <- sf::st_as_sf(sf::st_sfc(
    sf::st_point(c(2,2)),
    sf::st_point(c(2,7))
  )) %>%
    dplyr::rename("geometry" = "x") %>%
    sf::st_as_sf(wkt = "geometry")
  nodes <- sf::st_as_sf(nodes, wkt = "x") %>%
    dplyr::mutate(type = "Barrier") %>%
    dplyr::mutate(perm = 0)

  # Run test without tolerance
  expect_snapshot(join_attributes(net, nodes))

  # Run test with tolerance
  expect_snapshot(join_attributes(net, nodes, tolerance = 1))

})