tests/testthat/test-get_elev_point.R

In USEPA/elevatr: Access Elevation Data from Various APIs

context("get_elev_point")
library(sf)
library(elevatr)
library(terra)
data("pt_df")
data("sf_big")
skip_on_cran()
#skip_on_os("solaris")
ll_prj  <- 4326
aea_prj <- 5072

sf_sm <- st_as_sf(pt_df, coords = c("x", "y"), crs = ll_prj)
sf_sm_prj <- st_transform(sf_sm, crs = aea_prj) 
bad_sf <- st_as_sf(data.frame(x = 1000, y = 1000), coords = c("x", "y"), 
                   crs = ll_prj)
blank_raster <- rast(sf_sm, nrow = 5, ncol = 5, vals = 1)
sf_sm_raster <- rasterize(sf_sm, rast(sf_sm, nrow = 10, ncol = 10))

test_that("get_elev_point returns correctly", {
  
  epqs_df <- get_elev_point(locations = pt_df, prj = ll_prj, src = "epqs", serial = TRUE)
  epqs_sf <- get_elev_point(locations = sf_sm, src = "epqs", serial = TRUE)
  epqs_sf_prj <- get_elev_point(locations = sf_sm_prj, src = "epqs", serial = TRUE)
  epqs_ft <- get_elev_point(locations = sf_sm, src = "epqs", units = "feet", serial = TRUE)
  epqs_m <- get_elev_point(locations = sf_sm, src = "epqs", units = "meters", serial = TRUE)
  epqs_df_aws <- get_elev_point(locations = pt_df, prj = ll_prj, src = "aws", serial = TRUE)
  epqs_sf_aws <- get_elev_point(locations = sf_sm, src = "aws", serial = TRUE)
  epqs_sf_aws_z <- get_elev_point(locations = sf_sm, src = "aws", z = 4, serial = TRUE)
  epqs_sf_aws <- get_elev_point(locations = sf_sm, src = "aws", serial = TRUE)
  epqs_ft_aws <- get_elev_point(locations = sf_sm, src = "aws", units = "feet", serial = TRUE)
  epqs_rast <- get_elev_point(locations = blank_raster, serial = TRUE)
  epqs_sf_raster <- get_elev_point(locations = sf_sm_raster, serial = TRUE)
  
  
  
  # class
  expect_is(epqs_df, "sf")
  expect_is(epqs_sf, "sf")
  expect_is(epqs_sf_prj, "sf")
  expect_is(epqs_rast, "sf")
  expect_is(epqs_sf_raster, "sf")
  
  # crs
  expect_equal(st_crs(sf_sm)$wkt,st_crs(epqs_sf)$wkt)
  expect_equal(st_crs(sf_sm_prj)$wkt,st_crs(epqs_sf_prj)$wkt)
  expect_equal(st_crs(sf_sm)$wkt,st_crs(epqs_sf_aws)$wkt)
  
  
  # units
  expect_equal(epqs_ft$elev_units[1],"feet")
  expect_equal(epqs_m$elev_units[1],"meters")
  expect_equal(epqs_ft_aws$elev_units[1],"feet")
  expect_equal(epqs_sf_aws$elev_units[1],"meters")
  
  # num points from SpatRaster
  expect_equal(nrow(epqs_rast), 25)
  expect_equal(nrow(epqs_sf_raster), 5)
})