tests/testthat/test-calculate_feature.R

In land-info-lab/viewscape: Viewscape Analysis

testthat::test_that("runs correctly", {

  #Load in DSM
  test_dsm <- terra::rast(system.file("test_dsm.tif",
                                      package ="viewscape"))

  # load canopy raster
  test_canopy <- terra::rast(system.file("test_canopy.tif",
                                         package ="viewscape"))

  #Load in the viewpoint
  test_viewpoint <- sf::read_sf(system.file("test_viewpoint.shp",
                                            package = "viewscape"))

  #Compute viewshed
  output <- viewscape::compute_viewshed(dsm = test_dsm,
                                        viewpoints = test_viewpoint,
                                        offset_viewpoint = 6,
                                        plot=FALSE)

  # calculate the percentage of canopy coverage
  test_canopy_proportion <- viewscape::calculate_feature(viewshed = output,
                                                         feature = test_canopy,
                                                         type = 2,
                                                         exclude_value=0)

  testthat::expect_type(test_canopy_proportion, "double")

})