def pitremove_example():
"""run function of TauDEM, take pitremove as an example.
Compare the max, min, and average of rawdem and filled DEM.
The result will be::
RawDEM: Max: 284.07, Min: 139.11, Mean: 203.92
FilledDEM: Max: 284.07, Min: 139.11, Mean: 203.93
"""
dem = '../tests/data/Jamaica_dem.tif'
wp = '../tests/data/tmp_results'
fel = 'dem_pitremoved.tif'
taudem_bin = None
mpi_bin = None
num_proc = 2
TauDEM.pitremove(num_proc,
dem,
fel,
wp,
mpiexedir=mpi_bin,
exedir=taudem_bin)
rawdem = RasterUtilClass.read_raster(dem)
feldem = RasterUtilClass.read_raster(wp + os.sep + fel)
print('RawDEM: Max: %.2f, Min: %.2f, Mean: %.2f' %
(rawdem.get_max(), rawdem.get_min(), rawdem.get_average()))
print('FilledDEM: Max: %.2f, Min: %.2f, Mean: %.2f' %
(feldem.get_max(), feldem.get_min(), feldem.get_average()))
def pitremove_simple_usage():
"""Simple usage of pitremove.
Workspace will be set as the base directory of input file.
"""
dem = '../tests/data/Jamaica_dem.tif'
fel = 'dem_pitremoved.tif'
num_proc = 2
TauDEM.pitremove(num_proc, dem, fel)
def pitremove_simple_usage():
"""Simple usage of pitremove.
Workspace will be set as the base directory of input file.
"""
dem = '../tests/data/Jamaica_dem.tif'
fel = '../tests/data/tmp_results/dem_pitremoved.tif'
num_proc = 2
TauDEM.pitremove(num_proc, dem, fel)
