tests/testthat/test_navigate.R

In nhdplusTools: NHDPlus Tools

test_that("navigation basics", {
  skip_on_cran()

  suppressWarnings(
  # 06287800
  net <-  navigate_network(list(featureSource = "nwissite", featureID = "USGS-11486500"),
                           "UM",
                           output = "flowlines",
                           trim_start = TRUE))

  expect_equal(nrow(net), 5)

  expect_warning(net2 <-  navigate_network(2562073,
                           "UM",
                           output = "flowlines",
                           trim_start = TRUE),
                 "trim_start ignored for comid start")


  start <- sf::st_sfc(sf::st_point(c(-121.6678, 42.14987)), crs = sf::st_crs(4269))

  suppressWarnings(
  net3 <-  navigate_network(start,
                           "UM",
                           output = "flowlines",
                           trim_start = TRUE))

  net_no_split <- navigate_network(start,
                                   "UM",
                                   output = "flowlines",
                                   trim_start = FALSE)

  expect_equal(names(net), names(net3))

  expect_true(length(sf::st_geometry(net_no_split)[[5]]) > length(sf::st_geometry(net)[[5]]))

  suppressWarnings(
  net4 <-  navigate_network(start,
                            "UM", network = net_no_split,
                            output = "flowlines",
                            trim_start = TRUE))

  expect_equal(names(net4), names(net))
  expect_equal(nrow(net4), nrow(net))

  expect_error(navigate_network(1234567,
                                "UM", network = net_no_split,
                                output = "flowlines",
                                trim_start = TRUE))

  expect_error(navigate_network("12345",
                                "UM",
                                output = "flowlines",
                                trim_start = TRUE))

  expect_error(navigate_network(12345.5,
                                "UM",
                                output = "flowlines",
                                trim_start = TRUE))

  suppressWarnings(
  net5 <-  navigate_network(start,
                            "UM", network = net_no_split,
                            output = "nwissite",
                            trim_start = TRUE))
})

test_that("walker", {
  skip_on_cran()
  source(system.file("extdata", "walker_data.R", package = "nhdplusTools"))
  hydro_location <- list(comid = 5329339,
                         reachcode = "18050005000078",
                         reach_meas = 30)

  hydro_location <- sf::st_sf(hydro_location,
                              geom = nhdplusTools::get_hydro_location(data.frame(hydro_location),
                                                                      walker_flowline))

  netD <- navigate_network(hydro_location,
                          mode = "DM", network = walker_flowline,
                          distance_km = 4, trim_start = TRUE)

  netM <- navigate_network(sf::st_geometry(hydro_location),
                          mode = "UM", network = walker_flowline,
                          distance_km = 4, trim_start = TRUE)

  expect_true(nrow(sf::st_coordinates(netD)) < nrow(sf::st_coordinates(netM)))

})