The goal of traudem is to wrap the TauDEM CLI from R.
- It provides a guide to installation of TauDEM and its dependencies GDAL and MPI, depending on your operating system.
- It checks that TauDEM and its dependencies are correctly installed and included to the PATH.
- It provides wrapper commands for calling TauDEM methods from R.
- It does not provide an automatic installation of TauDEM and its
dependencies. Please refer to the
vignette("taudem-installation")for instructions on this. - It does not provide extensive documentation on the functioning of TauDEM commands. Please refer to the original TauDEM documentation for that.
Yes! Some TauDEM methods have dedicated wrappers like
taudem_pitremove(), with argument names that we strove to make
informative. For other methods, you can use taudem_exec() and provide
the arguments as TauDEM CLI would expect them. You could make a PR to
this repository to add more dedicated wrappers.
Refer to
vignette("taudem-installation", package = "traudem").
In particular, after installing the TauDEM CLI and dependencies as well
as the R package, please run traudem::taudem_sitrep().
traudem::taudem_sitrep()
#> ✔ Found GDAL version 3.0.4.
#> ✔ Found mpiexec (OpenRTE) 4.0.3 (MPI).
#> ✔ Found TauDEM path (/usr/local/taudem).
#> ✔ Found TauDEM executables directory (/usr/local/taudem).
#> ✔ Found all TauDEM executables.
#> ℹ Testing TauDEM on an example file (please wait a bit)...
#> ── TauDEM output ───────────────────────────────────────────────────────────────
#> PitRemove version 5.3.9
#> Input file DEM.tif has projected coordinate system.
#> Nodata value input to create partition from file: nan
#> Nodata value recast to float used in partition raster: nan
#> Processes: 1
#> Header read time: 0.008763
#> Data read time: 0.003737
#> Compute time: 0.033597
#> Write time: 0.002236
#> Total time: 0.048333
#> This run may take on the order of 1 minutes to complete.
#> This estimate is very approximate.
#> Run time is highly uncertain as it depends on the complexity of the input data
#> and speed and memory of the computer. This estimate is based on our testing on
#> a dual quad core Dell Xeon E5405 2.0GHz PC with 16GB RAM.
#> ── End of TauDEM output ────────────────────────────────────────────────────────
#> ✔ Was able to launch a TauDEM example!
#> ! Double-check above output for serious error messages.You can install traudem from CRAN:
install.packages("traudem")You can install the development version of traudem from R-universe:
# Enable repository from lucarraro
options(
repos = c(
lucarraro = 'https://lucarraro.r-universe.dev',
CRAN = 'https://cloud.r-project.org'
)
)
# Download and install traudem in R
install.packages('traudem')Or from GitHub with:
# install.packages("devtools")
devtools::install_github("lucarraro/traudem")This is how you would use traudem to launch TauDEM’s PitRemove:
library(traudem)
test_dir <- withr::local_tempdir()
fs::file_copy(
system.file("test-data", "DEM.tif", package = "traudem"),
file.path(test_dir, "DEM.tif")
)
output <- taudem_pitremove(file.path(test_dir, "DEM.tif"))
#> PitRemove version 5.3.9
#> Input file /tmp/RtmpoEf2qy/file110784aada156/DEM.tif has projected coordinate system.
#> Nodata value input to create partition from file: nan
#> Nodata value recast to float used in partition raster: nan
#> Processes: 1
#> Header read time: 0.008310
#> Data read time: 0.002408
#> Compute time: 0.023698
#> Write time: 0.001950
#> Total time: 0.036366
#> This run may take on the order of 1 minutes to complete.
#> This estimate is very approximate.
#> Run time is highly uncertain as it depends on the complexity of the input data
#> and speed and memory of the computer. This estimate is based on our testing on
#> a dual quad core Dell Xeon E5405 2.0GHz PC with 16GB RAM.
output
#> [1] "/tmp/RtmpoEf2qy/file110784aada156/DEMfel.tif"We ran the example above in a temporary directory that withr
automatically deletes. If you want to automatically get rid of some of
the intermediary files created by TauDEM in one of your pipelines, you
might be interested in withr::local_tempfile().
If you wanted to run this same code without seeing the messages from
TauDEM, you can either use the quiet argument:
test_dir <- withr::local_tempdir()
fs::file_copy(
system.file("test-data", "DEM.tif", package = "traudem"),
file.path(test_dir, "DEM.tif")
)
output <- taudem_pitremove(file.path(test_dir, "DEM.tif"), quiet = TRUE)
output
#> [1] "/tmp/RtmpoEf2qy/file110783a666c67/DEMfel.tif"Or set the traudem.quiet option (options(traudem.quiet = TRUE) or
for just the session, withr::local_options(traudem.quiet = TRUE)):
withr::local_options(traudem.quiet = TRUE)
test_dir <- withr::local_tempdir()
fs::file_copy(
system.file("test-data", "DEM.tif", package = "traudem"),
file.path(test_dir, "DEM.tif")
)
output <- taudem_pitremove(file.path(test_dir, "DEM.tif"))
output
#> [1] "/tmp/RtmpoEf2qy/file110784599c1a1/DEMfel.tif"