tests/testthat/test-track_sensor.R

In Jean-Romain/lidR: Airborne LiDAR Data Manipulation and Visualization for Forestry Applications

context("track_sensor")

las = topography
ctg <- topography_ctg
opt_chunk_size(ctg) <- 300
opt_chunk_alignment(ctg) <- c(50,180)
opt_chunk_buffer(ctg) <- 40
opt_progress(ctg) <- FALSE

test_that("track_sensor works with Roussel2020", {

  flightlines <- track_sensor(las, Roussel2020(pmin = 10))

  expect_is(flightlines, "sf")
  expect_is(sf::st_geometry(flightlines), "sfc_POINT")
  expect_equal(dim(flightlines), c(8,4))
  expect_equal(names(flightlines)[1:3], c("gpstime", "PointSourceID", "SCORE"))
  expect_equal(mean(sf::st_coordinates(flightlines)[,3]), 3100, tol = 3)
  expect_equal(sf::st_crs(flightlines), st_crs(las))
})

test_that("track_sensor works with Gatziolis2019", {

  flightlines <- track_sensor(las, Gatziolis2019(deltaT = 0.5))

  expect_is(flightlines, "sf")
  expect_is(sf::st_geometry(flightlines), "sfc_POINT")
  expect_equal(dim(flightlines), c(9,4))
  expect_equal(names(flightlines)[1:3], c("gpstime", "PointSourceID", "SCORE"))
  expect_equal(mean(sf::st_coordinates(flightlines)[,3]), 3100, tol = 3)
  expect_equal(sf::st_crs(flightlines), st_crs(las))
})

test_that("track_sensor returns an empty sf", {

  flightlines <- track_sensor(las, Roussel2020(pmin = 50))

  expect_is(flightlines, "sf")
  expect_equal(dim(flightlines), c(0,4))
  expect_equal(names(flightlines)[1:3], c("gpstime", "PointSourceID", "SCORE"))
})

test_that("track_sensor works with a LAScatalog", {

  skip_on_cran()

  flightlines <- track_sensor(ctg, Roussel2020(pmin = 10))

  expect_is(flightlines, "sf")
  expect_is(sf::st_geometry(flightlines), "sfc_POINT")
  expect_equal(dim(flightlines), c(8,4))
  expect_equal(names(flightlines)[1:3], c("gpstime", "PointSourceID", "SCORE"))
  expect_equal(mean(sf::st_coordinates(flightlines)[,3]), 3100, tol = 3)
  expect_equal(sf::st_crs(flightlines), st_crs(las))
})

test_that("track_sensor handle errors", {

  las2 = las
  las2@data = las@data[, !"PointSourceID"]
  expect_error(track_sensor(las2, Roussel2020(pmin = 15)), "No 'PointSourceID' attribute found")

  las2@data$PointSourceID <- 0L
  expect_error(track_sensor(las2, Roussel2020(pmin = 15)), "'PointSourceID' attribute is not populated")

  las2$gpstime <- 0
  expect_error(track_sensor(las2, Roussel2020(pmin = 15)), "'gpstime' attribute is not populated")

  las2 = las
  las2$PointSourceID[15] <- 5L
  expect_warning(track_sensor(las2, Roussel2020(pmin = 15), thin_pulse_with_time = 0), "points with same gpstime")
})

test_that("track_sensor filter duplicates", {

  las = rbind(las, las[1:25])

  expect_warning(flightlines <- track_sensor(las, Roussel2020(pmin = 10)), "duplicated")

  expect_is(flightlines, "sf")
  expect_is(sf::st_geometry(flightlines), "sfc_POINT")
  expect_equal(dim(flightlines), c(8,4))
  expect_equal(names(flightlines)[1:3], c("gpstime", "PointSourceID", "SCORE"))
  expect_equal(mean(sf::st_coordinates(flightlines)[,3]), 3100, tol = 3)
  expect_equal(sf::st_crs(flightlines), st_crs(las))
})