Подготовить нужно только растр аккумуляции (с помощью стандартного <r.watershed>). Остальные сгенерит <r.stream.extract>.
Вот последовательность команд:
Код: Выделить всё
r.hydrodem in=dem out=dem_h
r.watershed ele=dem_h acc=acc
r.stream.extract elevation=dem_h accumulation=acc threshold=100 stream_rast=streams stream_vect=streams dir=dir stream_length=200
r.stream.basins -l dir=dir stream_rast=streams basins=basins_main
r.stream.basins dir=dir stream_rast=streams basins=basins_all
r.stream.order stream_rast=streams direction=dir elevation=dem_h accumulation=acc stream_vect=streams_order strahler=streams_strahler