Burned flow direction

[tcalegario]

Dear dtarb!

This package is awesome. Thanks for developing it.

I`m trying to obtain burned flow direction file using the .exe procedure described on https://hydrology.usu.edu/taudem/taudem5/TauDEMRScript.txt

But I'm not being well succeeded on it. Can you give me some lights?

1. Filled Dem (work perfectly)

system("mpiexec -n 8 pitremove -z ./intermediarios/NASAdem_S2000_30m_raster.tif -fel ./intermediarios/burned/demFill.tif")

2. D8 (work perfectly)

system("mpiexec -n 8 D8Flowdir -p ./intermediarios/burned/demFillDir8.tif -sd8 ./intermediarios/burned/demFillSlopeDir8.tif -fel ./intermediarios/burned/demFill.tif")

3. D8 Burned based on streamnet (do not work)

system("mpiexec -n 8 D8Flowdir -p ./intermediarios/burned/demFillDir8.tif -sd8 ./intermediarios/burned/demFillSlopeDir8.tif -fel ./intermediarios/burned/demFill.tif -sfdr ./intermediarios/streamnetRaster30m.tif")

Thanks for your attention!

[dtarb]

@tcalegario what is the error that you get. The most likely problem with this third step is that streamnetRaster30m.tif is not constructed correctly for use by D8Flowdir. How did you obtain or calculate it? More broadly what are you trying to do?

[tcalegario]

@dtarb, thanks for your reply!

to obtain stremnetRaster30m.tif, I converted a polyline stream file to raster on ArcGis respecting the grid of "dem".

I'm trying to make my accumulated flow coincide with the mapped rivers that are in a shapefile.

[tcalegario]

@dtarb, I tried again moving "mapped streams" input to "pitremove" tool. But the result of flow Accumulation was not what I'm expecting (Figure below).

system("mpiexec -n 8 pitremove -z ./demProjected.tif -fel ./demPFillBurn.tif -sfdr ./stream30m.tif")

system("mpiexec -n 8 D8Flowdir -p ./Dir8.tif -sd8 ./SlopeDir8.tif -fel ./demPFillBurn.tif" ,show.output.on.console=F,invisible=F)

system("mpiexec -n 8 AreaD8 -p ./Dir8.tif -ad8 ./accumul_nc.tif -nc")

When I try to use -sfrd optional input on "D8Flowdir" I got error 65533.

system("mpiexec -n 8 D8Flowdir -p ./Dir8.tif -sd8 ./SlopeDir8.tif -fel ./demPFillBurn.tif -sfdr ./stream30m.tif" ,show.output.on.console=F,invisible=F)

Cause of my difficulty, I obtained flow accumulation using "whiteBox" package on R, and apparently worked fine (Figure Below). But I would like to compare both procedures (Taudem x WhiteBox). If there is a way to obtain "similar" result as "WhiteBox" please let me know. My goal is to automating watershed delineation on R.

Thanks a lot for your time.

[lucarraro]

Hello,

The R package traudem was just published on CRAN, and it might help you solve this issue.

What traudem does:

• it guides installation of TauDEM and its dependencies GDAL and MPI for different OSs
• 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