tests/testthat/test-spReprojectVector.R

In USDAForestService/FIESTA: Forest Inventory Estimation and Analysis

test_that("Reproject of SF Vector Object works", {

  skip_on_cran()

  WYspplt <- spMakeSpatialPoints(xyplt = WYplt,
                                 xy.uniqueid = "CN",
                                 xvar = "LON_PUBLIC",
                                 yvar = "LAT_PUBLIC",
                                 xy.crs = 4269)

  new_prj <- "+proj=utm +zone=12 +ellps=GRS80 +datum=NAD83 +units=m +no_defs"

  # Reproject with new projection
  WYspplt.utm12 <- spReprojectVector(layer = WYspplt, crs.new = new_prj)

  # Retrieve new projection from SF object
  retrieved_prj <- sf::st_crs(WYspplt.utm12)

  # Check that retrieved projection contains input projection
  expect_contains(retrieved_prj, new_prj)
})