Skip to content

Watershedding in WBW not producing the same output as whitebox_tools cli #480

@reasonman

Description

@reasonman

I'm trying to do what i think is a simple watershed exercise on a set of rasters from NHDPlus(precomputed fdir/fac). When I run it through python it seems to just be outputting the flow direction raster i gave as input. Running the same options from the cli produces what I expect, a small watershed encompassing a point.

Python 3.10, WBW 1.3.0

from whitebox_workflows import WbEnvironment, VectorGeometryType
wbe = WbEnvironment()
fdir = wbe.read_raster('12090205_fdir.tif')
pour_point = wbe.read_raster('pour_from_wbw.tif')
shed = wbe.watershed_from_raster_pour_points(fdir, pour_point, esri_pntr=True)
wbe.write_raster(shed, 'watershed_from_wbw_file.tif')

CLI(using the whitebox_tools via pip)

.venv/lib/python3.10/site-packages/whitebox/WBT/whitebox_tools -r=Watershed -v --d8_pntr=12090205_fdir.tif --pour_pts=pour_from_wbw.tif -o=watershed_from_wbw_file.tif --esri_pntr

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions