<div dir="ltr"><div>Hi,<br>I would like to run r.watershed and get the<br>basin=name  Output map: unique label for each watershed basin<br><br>considering the water bodies (lakes and rivers) as an  ancillary layer.<br><br>
If i set water bodies as NULL value in the elevation data-set<br></div><div>  I get nodata computation along the water-bodies borders  (the truncate basin are not compute)<br><br>If i set - water bodies as MASK<br></div><div>
  I just receive a mask-out result of the basin calculation (resulting in splitting of a basin in case of river)<br><br></div><div>Is there a way to get basin raster taking in to account the actual water-bodies?<br><br></div>
<div>Thanks<br></div><div>Best<br></div><div><div><div><br>-- <br><div dir="ltr">Giuseppe Amatulli, Ph.D.
   <br><br>Center for Earth Observation                         <br>ESC Room 119A                                           <br>Yale University                                            <span>      </span><br>P.O. Box 208109 <span>                                       </span><br>
New Haven, CT, 06520-8109  <br><div>
   Teaching: <a href="http://spatial-ecology.net" target="_blank">spatial-ecology.net</a>
  </div> 
  
   Work:  <a href="http://www.exeter.ac.uk/esi/people/casalegno/" target="_blank"></a><a href="http://sbsc.yale.edu/training-and-support" target="_blank">http://sbsc.yale.edu/training-and-support</a><br><a href="http://www.spatial-ecology.net" target="_blank"></a></div>

</div></div></div></div>