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
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
CLI(using the whitebox_tools via pip)