Matt Bartos

- You can `del startnodes` after line 1233, then set line 1251 to something like `startnodes = np.flatnonzero(indegree == 0)`. - Moving some other code around to reduce the amount...

Greetings, Glad this has been helpful for you. I believe you can just do: ```python dem_coords = dem.coords dem_coords_reshaped = [dem_coords[:,0].reshape(dem.shape), dem_coords[:,1].reshape(dem.shape)] ``` For verification, see: ```python import numpy as...

Greetings. Because z1 has the same coordinates as the current grid, you can just do: ```python # Add z1 to grid grid.add_gridded_data(z1, data_name='z1', affine=grid.affine, shape=z1.shape, crs=grid.crs, nodata=np.nan) # Write to...

Hmmm, it's hard to say without a traceback, or without looking at the data source you are using. One issue is that there could be cycles (loops) in the flow...

I don't have enough info to give an answer. Feel free to post the dataset.

Yes, external upload would be preferred.

Greetings, The bbox format should be (xmin, ymin, xmax, ymax), so I believe your bbox needs to be (0, 0, 1291, 1291). Let me know if that fixes the issue....

Greetings, For the first image it looks like you are trying to delineate at the x-coordinate 145.29, which is less than than the minimum x-coordinate of your dataset view (145.290066...)....

Basically, if you know what the river network is supposed to look like, you'll just want to find an approximate threshold accumulation value that captures that structure. Using a percentile...

Interesting. What version of numpy are you using? I'm on 1.18.3 and the following still seem to work: ```python import numpy as np x = np.random.randint(100, size=1000, dtype=np.int64) _ =...