test_that("raster_to_points works for all data type", {

    result_layer <- raster_to_points(landscape)
    result_stack <- raster_to_points(landscape_stack)
    result_list <- raster_to_points(landscape_list)

    expect_true(is.matrix(result_layer))
    expect_true(is.matrix(result_stack))
    expect_true(is.matrix(result_list))

    expect_equal(object = nrow(result_layer),
                 expected = terra::ncell(landscape))

    expect_equal(object = nrow(result_stack),
                 expected = terra::ncell(landscape_stack) *
                     terra::nlyr(landscape_stack))

    expect_equal(object = nrow(result_list),
                 expected = (terra::ncell(landscape_list[[1]]) +
                                 terra::ncell(landscape_list[[2]])) / 2 *
                     length(landscape_list))
})

test_that("raster_to_points considers NAs", {

    landscape_NA <- landscape
    landscape_NA[1, 1] <- NA

    result <- raster_to_points(landscape_NA, return_NA = FALSE)
    result_NA <- raster_to_points(landscape_NA, return_NA = TRUE)

    expect_true(any(!is.na(result)))
    expect_true(any(is.na(result_NA)))

    expect_equal(object = nrow(result),
                 expected = terra::ncell(landscape) - 1)
    expect_equal(object = nrow(result_NA),
                 expected = terra::ncell(landscape))
})