inst/tinytest/_test-runGdal.R

In MODIS: Acquisition and Processing of MODIS Products

lap = file.path(tempdir(), "MODIS_ARC")
odp = file.path(lap, "PROCESSED")
jnk = capture.output(
  suppressMessages(
    MODISoptions(localArcPath = lap, outDirPath = odp
                 , save = FALSE, quiet = TRUE, checkTools = FALSE)
  )
)

jnk = capture.output(
  orgStruc(system.file("external", package = "MODIS")
           , pattern = "^MOD13A2.A2016145.h18v04.006.2016166145124.hdf$"
           , quiet = TRUE)
)


### 0 default settings ----

jnk = capture.output(
  tfs0 <- runGdal("MOD13A2", collection = "006"
                 , begin = "2016145", end = "2016145"
                 , tileH = 18, tileV = 4
                 , SDSstring = "111", job = "test-runGdal"
                 , checkIntegrity = FALSE, forceDownload = FALSE
                 , overwrite = TRUE)
)

rst0 = raster::stack(unlist(tfs0))

expect_equivalent(
  dim(rst0)
  , target = c(1200, 1410, 3)
  , info = "default output has expected dimensions (# rows, cols, layers)"
)


### 1 custom settings ----

jnk = capture.output(
  tfs1 <- runGdal(
    "MOD13A2", collection = "006"
    , begin = "2016145", end = "2016145"
    , tileH = 18, tileV = 4
    , SDSstring = "1", job = "test-runGdal"
    , outProj = 32632, pixelSize = 1000
    , checkIntegrity = FALSE, forceDownload = FALSE
    , overwrite = TRUE
  )
)

rst1 = raster::stack(unlist(tfs1))

expect_true(
  raster::nlayers(rst1) == 1
  , info = "customized output has expected # layers"
)

expect_equivalent(
  raster::res(rst1)
  , target = c(1000, 1000)
  , info = "customized output has expected resolution"
)

expect_true(
  sf::st_crs(rst1) == sf::st_crs(32632)
  , info = "customized output inherits specified crs"
)