It should be explicitly specified that the returned graph is a nx.DiGraph. This could be done by adding the create_using=nx.DiGraph to the list of arguments of the function above.

Adding a list of outlet IDs can be helpful for multiple basin domain extraction

There are cases where non-contributing areas are part of the catchment geometry, with the difference that no river segments are defined for them. This can be easily detected as part of the prepare_cat() function.

The following argument needs to get out of hard-coded state, and should only be used if provided as an input. Meaning, the area is only calculated if this variable is entered as an input (for now, using **kwargs argument):

The unit of the area needs to be defined by the user. Otherwise, the default unit of CRS with EPSG value of 6933 will be considered. Furthermore, the fid of the area variable needs to be given by the user.

If the dtype of input columns of the rivs object is not int or float, the following lines will cause an issue:

It is better to get the data types of each column as a separate argument (optional), and do necessary checks given the object's data type.

The package misses relevant readthedocs documentation.

An aggregation criteria applicable to both river DAGs and sub-basins Polygons is needed to target sub-basins below a certain threshold

Request by: @zktesemma