tests/testthat/testthat_changeCoords.R

In PhilipMostert/inlabruSDMs: Fit Models Derived from Point Processes to Species Distributions using 'inlabru'

testthat::test_that('changeCoords changes the coordinate variables of a list of data.', {
  
  ##Generate arbitrary data:
  #Arbitrary projection
  projection <- '+proj=tmerc'
  x <- c(16.48438,  17.49512,  24.74609, 22.59277, 16.48438)
  y <- c(59.736328125, 55.1220703125, 55.0341796875, 61.142578125, 59.736328125)
  xy <- cbind(x, y)
  xy <- cbind(x, y)
  SpatialPoly <- st_sfc(st_polygon(list(xy)), crs = projection)
  
  ##Old coordinate names
  oldcoords <- c('x','y')
  #Make random points
  #Random presence only dataset
  PO <- st_as_sf(st_sample(SpatialPoly, 100, crs = projection))
  st_geometry(PO) <- 'geometry'
  
  #Random presence absence dataset
  PA <- st_as_sf(st_sample(SpatialPoly, 100, crs = projection))
  st_geometry(PA) <- 'geometry'
  PA$PAresp <- sample(x = c(0,1), size = nrow(PA), replace = TRUE)
  #Add trial name
  PA$trial <- sample(x = c(1,2,3), size = nrow(PA), replace = TRUE)
  PA$x <- st_coordinates(PA)[,1]
  PA$y <- st_coordinates(PA)[,2]
  st_geometry(PA) <- NULL
  
  ##Make PA a data.frame object
  PA <- data.frame(PA)
  
  spData <- list(PO,PA)  
  
  #Set new coordinate names
  newcoords <- c('long', 'lat')
  
  mod <- changeCoords(data = spData, oldcoords = oldcoords, newcoords = newcoords)
  
  expect_equal(colnames(mod[[1]]), 'geometry')
  expect_true(all(newcoords%in%names(mod[[2]])))
  expect_false(any(oldcoords%in%names(mod[[2]])))
  
})