Nothing
test_that("Creating WORLD-COVER-2021 cubes from TERRASCOPE", {
class_cube <- .try(
{
sits_cube(
source = "TERRASCOPE",
collection = "WORLD-COVER-2021",
roi = c(
lon_min = 137.15991,
lon_max = 138.18467,
lat_min = -33.85777,
lat_max = -32.56690
),
progress = FALSE
)
},
.default = NULL
)
testthat::skip_if(purrr::is_null(class_cube),
message = "TERRASCOPE is not accessible"
)
expect_true(all(sits_bands(class_cube) %in% c("class")))
expect_equal(nrow(class_cube), 4)
bbox_cube <- sits_bbox(class_cube, as_crs = "EPSG:4326")
bbox_cube_1 <- sits_bbox(.tile(class_cube), as_crs = "EPSG:4326")
expect_true(bbox_cube["xmax"] >= bbox_cube_1["xmax"])
expect_true(bbox_cube["ymax"] >= bbox_cube_1["ymax"])
rast <- .raster_open_rast(class_cube$file_info[[1]]$path[1])
cube_nrows <- .tile_nrows(class_cube)
expect_true(.raster_nrows(rast) == cube_nrows)
})
test_that("Creating WORLD-CEREAL-2021 cubes from TERRASCOPE",{
# get roi for an MGRS tile
bbox_22LBL <- sits_mgrs_to_roi("22LBL")
# retrieve the world cereal map for the chosen roi
world_cereal_2021 <- sits_cube(
source = "TERRASCOPE",
collection = "WORLD-CEREAL-2021",
roi = bbox_22LBL
)
# cut the 3 x 3 degree grid to match the MGRS tile 22LBL
world_cereal_2021_20LBL <- sits_cube_copy(
cube = world_cereal_2021,
roi = bbox_22LBL,
multicores = 6,
output_dir = tempdir()
)
roi_wc <- sits_bbox(world_cereal_2021_20LBL)[,1:4]
roi_20LBL <- .bbox(bbox_22LBL, as_crs = "EPSG:4326")[,1:4]
expect_equal(roi_wc[["xmin"]], roi_20LBL[["xmin"]], tolerance = 0.001)
sumwc <- summary(world_cereal_2021_20LBL)
expect_true(all(sumwc[["class"]] %in% c("Non_Cropland", "Cropland")))
expect_true(all(sumwc[["value"]] %in% c("0", "100")))
})
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.