tests/testthat/test_get_dem.R

In FluvialGeomorph/fluvgeo: A Package for Fluvial Geomorphic Analysis

dem_plot <- function(dem, sf_lines) {
  terra::plot(dem)
  terra::lines(vect(sf_lines))
}

test_that("WGS84 input - dem is a terra::SpatRaster", {
  xs_mapedit <- sf::st_read(system.file("extdata", "shiny", "xs_mapedit.shp",
                                package = "fluvgeodata"), quiet = TRUE)
  xs <- sf::st_transform(xs_mapedit, crs = 3857) # Web Mercator
  dem <- get_dem(xs)
  #dem_plot(dem, xs)
  expect_true("SpatRaster" %in% class(dem))
})

test_that("web mercator input - dem is a terra::SpatRaster", {
  xs <- sf::st_read(system.file("extdata", "shiny", "xs.shp",
                       package = "fluvgeodata"), quiet = TRUE)
  dem <- get_dem(xs)
  #terra::plot(dem)
  expect_true("SpatRaster" %in% class(dem))
})

test_that("expect error - WGS84 input", {
  xs_mapedit <- sf::st_read(system.file("extdata", "shiny", "xs_mapedit.shp",
                                package = "fluvgeodata"), quiet = TRUE)
  expect_error(get_dem(xs_mapedit))
})