[GRASS-user] running r.watershed taking in to account water bodies

Hi,
I would like to run r.watershed and get the
basin=name Output map: unique label for each watershed basin

considering the water bodies (lakes and rivers) as an ancillary layer.

If i set water bodies as NULL value in the elevation data-set

I get nodata computation along the water-bodies borders (the truncate basin are not compute)

If i set - water bodies as MASK

I just receive a mask-out result of the basin calculation (resulting in splitting of a basin in case of river)

Is there a way to get basin raster taking in to account the actual water-bodies?

Thanks

Best

Giuseppe Amatulli, Ph.D.

Center for Earth Observation
ESC Room 119A
Yale University
P.O. Box 208109
New Haven, CT, 06520-8109

Teaching: spatial-ecology.net
Work: http://sbsc.yale.edu/training-and-support