Dem module tests

.Rbuildignore

@@ -13,3 +13,4 @@
 ^CITATION.cff$
 ^doc$
 ^Meta$
+^tests/testthat/fixtures/

data-raw/bucharest_test.R

@@ -0,0 +1,11 @@
+dem <- terra::unwrap(bucharest_dem)
+river <- bucharest_osm$river_surface
+crs <- "epsg:32635"
+
+valley <- get_valley(dem, river, crs)
+
+if (!dir.exists("tests/testthat/fixtures")) {
+  dir.create("tests/testthat/fixtures", recursive = TRUE)
+}
+
+saveRDS(valley, "tests/testthat/fixtures/expected_valley.rds")

tests/testthat/test-valley.R

@@ -0,0 +1,54 @@
+test_that("STAC asset urls are correctly retrieved", {
+  skip_on_ci()
+  bb <- bucharest_osm$bb
+  ep <- "https://earth-search.aws.element84.com/v1"
+  col <- "cop-dem-glo-30"
+
+  asset_urls_retrieved <- get_stac_asset_urls(bb, endpoint = ep,
+                                              collection = col)
+  asset_urls_retrieved_default <- get_stac_asset_urls(bb)
+
+  expasset1 <- paste0("s3://copernicus-dem-30m/",
+                      "Copernicus_DSM_COG_10_N44_00_E026_00_DEM/",
+                      "Copernicus_DSM_COG_10_N44_00_E026_00_DEM.tif")
+  expasset2 <- paste0("s3://copernicus-dem-30m/",
+                      "Copernicus_DSM_COG_10_N44_00_E025_00_DEM/",
+                      "Copernicus_DSM_COG_10_N44_00_E025_00_DEM.tif")
+  expected_asset_urls <- c(expasset1, expasset2)
+
+  expect_equal(expected_asset_urls, asset_urls_retrieved)
+  expect_equal(expected_asset_urls, asset_urls_retrieved_default)
+})
+
+
+test_that("raster data are correctly retrieved and merged", {
+  skip_on_ci()
+  bb <- bucharest_osm$bb
+  rp1 <- paste0("s3://copernicus-dem-30m/",
+                "Copernicus_DSM_COG_10_N44_00_E026_00_DEM/",
+                "Copernicus_DSM_COG_10_N44_00_E026_00_DEM.tif")
+  rp2 <- paste0("s3://copernicus-dem-30m/",
+                "Copernicus_DSM_COG_10_N44_00_E025_00_DEM/",
+                "Copernicus_DSM_COG_10_N44_00_E025_00_DEM.tif")
+  raster_paths <- c(rp1, rp2)
+
+  dem <- load_raster(bb, raster_paths) |> CRiSp::reproject(32635)
+  expected_dem <- terra::unwrap(bucharest_dem)
+  expect_equal(terra::values(dem), terra::values(expected_dem))
+  expect_true(all.equal(terra::ext(dem), terra::ext(expected_dem),
+                        tolerance = 1e-4))
+  expect_equal(terra::crs(dem), terra::crs(expected_dem))
+})
+
+test_that("valley polygon is correctly constructed", {
+  skip_on_ci()
+  dem <- terra::unwrap(bucharest_dem)
+  river <- bucharest_osm$river_surface
+  crs <- "epsg:32635"
+
+  valley <- get_valley(dem, river, crs)
+  expected_valley <- readRDS("fixtures/expected_valley.rds")
+
+  expect_equal(valley, expected_valley, tolerance = 1e-4)
+
+})