test_that("get_tiles gets the same elevation tiles twice", {
skip_on_cran()
skip_if_offline()
dl_loc <- data.frame(
lat = c(44.04905, 44.04911),
lng = c(-74.01188, -74.01179)
)
dl_loc <- sf::st_as_sf(dl_loc, coords = c("lng", "lat"))
sf::st_crs(dl_loc) <- sf::st_crs(4326)
output_tif <- get_tiles(dl_loc,
services = c("elevation", "3DEPElevation"),
georeference = FALSE
)
expect_equal(names(output_tif), "elevation")
expect_equal(length(output_tif), 1)
expect_equal(length(output_tif[[1]]), 1)
stored_raster <- terra::rast("testdata/3DEP.tif")
test_raster <- terra::rast(output_tif[[1]])
expect_equal(as.vector(terra::crs(stored_raster)),
as.vector(terra::crs(test_raster)))
expect_equal(as.vector(terra::ext(stored_raster)),
as.vector(terra::ext(test_raster)))
expect_equal(
as.vector(terra::global(stored_raster, max))[[1]],
as.vector(terra::global(test_raster, max))[[1]],
tolerance = 0.01
)
})
test_that("get_tiles fails as expected", {
dl_loc <- data.frame(
lat = c(44.04905, 44.04911),
lng = c(-74.01188, -74.01179)
)
dl_loc <- sf::st_as_sf(dl_loc, coords = c("lng", "lat"))
sf::st_crs(dl_loc) <- sf::st_crs(4326)
expect_error(get_tiles(dl_loc,
services = c("elevation", "ortho", "USGSNAIPPlus"),
side_length = 4097
))
expect_error(get_tiles(dl_loc,
services = c("exactly_the_data_i_want")
))
expect_error(get_tiles(dl_loc,
services = c("elevation"),
side_length = 8001
))
})
test_that("specifying SR works", {
skip_on_cran()
skip_if_offline()
dat <- data.frame(
lat = c(44.04905, 44.04911),
lng = c(-74.01188, -74.01179)
)
dat <- sf::st_as_sf(dat, coords = c("lng", "lat"))
dat <- sf::st_set_crs(dat, 4326)
output_tif <- vector("list", 4)
output_tif[1] <- get_tiles(dat)
output_tif[2] <- get_tiles(dat, bboxSR = 4326)
output_tif[3] <- get_tiles(dat, imageSR = 4326)
output_tif[4] <- get_tiles(dat, bboxSR = 4326, imageSR = 4326)
expect_equal(
brio::read_file_raw(output_tif[[1]]),
brio::read_file_raw(output_tif[[2]])
)
expect_equal(
brio::read_file_raw(output_tif[[1]]),
brio::read_file_raw(output_tif[[3]])
)
expect_equal(
brio::read_file_raw(output_tif[[1]]),
brio::read_file_raw(output_tif[[4]])
)
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.