tests/testthat/test_get_hydro_location.R

In hydroloom: Utilities to Weave Hydrologic Fabrics

test_that("get location", {
  if (!requireNamespace("hydrogeofetch", quietly = TRUE)) skip("Missing hydrogeofetch")

  source(system.file("extdata", "sample_flines.R", package = "hydrogeofetch"))

  points <- sf::st_sfc(sf::st_sfc(list(sf::st_point(c(-76.86934, 39.49328)),
    sf::st_point(c(-76.91711, 39.40884)),
    sf::st_point(c(-76.88081, 39.36354))),
  crs = 4326))

  sample_flines <- sf::st_cast(sample_flines, "LINESTRING", warn = FALSE)

  indexes <- index_points_to_lines(sample_flines, points)

  locations <- get_hydro_location(indexes, sample_flines)

  expect_equal(as.numeric(sf::st_coordinates(locations)[, 1:2]),
    c(-76.8693957911233, -76.9176139963277, -76.8810037292214, 39.4932572053652,
      39.4090934738461, 39.3632976014239), tolerance = 1e-6)

  points <- sf::st_sfc(sf::st_sfc(list(sf::st_point(c(-76.86934, 39.49328))),
    crs = 4326))

  indexes <- index_points_to_lines(sample_flines, points)

  locations <- get_hydro_location(indexes, sample_flines)

  expect_equal(sf::st_coordinates(locations),
    sf::st_coordinates(sf::st_sfc(sf::st_point(c(-76.8694, 39.49326)), crs = 4326)),
    tolerance = 0.001)

})