<br><br><div class="gmail_quote">On Thu, Apr 19, 2012 at 5:58 PM, Rich Shepard <span dir="ltr"><<a href="mailto:rshepard@appl-ecosys.com">rshepard@appl-ecosys.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div class="im">On Thu, 19 Apr 2012, Margherita Di Leo wrote:<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Please post information about the region, <a href="http://r.info" target="_blank">r.info</a> of the map and the<br>
command line of r.watershed that you use. Also, do you use a mask? Pankaj<br>
Kr Sharma, please avoid to reply attaching the whole digest.<br>
</blockquote>
<br></div>
madi, Micha, Pankaj:<br>
<br>
The Level 6 HUC basins (a.k.a. HUC-12 or sub-watersheds) are defined by<br>
the USGS to enclose areas of 10k-40k acres (15.625-62.5 square miles;<br>
40.5-161.9 square kilometers). These are a bit too large, but vector<br>
boundaries for 14- or 16-digit hydrologic units are not available for this<br>
area (or for my other project area which is approximately 125 square miles,<br>
323.75 square kilometers) so I tried different threshold values for<br>
r.watershed and 150000 seems to provide additional resolution without<br>
generating a separate drainage for each individual stream segment. Yes, it's<br>
arbitrary. But for the purposes of these two projects it is a practical<br>
balance between too coarse and too fine.<br>
<br>
The region was set by the base, 10m DEM:<br>
<br>
projection: 99 (Lambert Conformal Conic; units in International Feet)<br>
zone: 0<br>
datum: nad83<br>
ellipsoid: grs80<br>
north: 735780<br>
south: 655650<br>
west: 2042100<br>
east: 2138580<br>
nsres: 30<br>
ewres: 30<br>
rows: 2671<br>
cols: 3216<br>
cells: 8589936<br>
<br>
The r.watershed command line (all on one line but wrapped in the MUA):<br>
<br>
r.watershed -f elev=dem10m.avg acc=dem.acc drain=dem.drain bas=dem.basin<br>
stream=dem_str.seg thresh=150000 --o<br>
<br>
And <a href="http://r.info" target="_blank">r.info</a> for dem.basin:<br>
<br>
Type of Map: raster Number of Categories: 94<br>
Data Type: CELL<br>
Rows: 2671<br>
Columns: 3216<br>
Total Cells: 8589936<br>
Projection: Lambert Conformal Conic<br>
N: 735780 S: 655650 Res: 30<br>
E: 2138580 W: 2042100 Res: 30<br>
Range of data: min = 2 max = 94<br>
Data Source:<br>
dem10m.avg<br>
Data Description:<br>
generated by r.watershed<br>
Comments:<br>
Processing mode: MFD<br>
Memory mode: All in RAM<br>
r.watershed -f elevation="dem10m.avg" accumulation="dem.acc" drainag\<br>
e="dem.drain" basin="dem.basin" stream="dem_str.seg" threshold=15000\<br>
0 convergence=5 memory=300<br>
<br>
If the white areas with vector streams (the blue lines) represent basins,<br>
then that's OK. What puzzles me is why a couple of basins would be white<br>
while the others are different colors and the small white areas without<br>
mapped streams are appently external to the basins of interest.<br></blockquote><div><br></div><div>You can query the white cells and discover that they are null values. This is because r.watershed delineates only the basins that fall entirely into the computational region. You didn't provide information about any mask, but you probably either set a too small mask or should enlarge the computational region in order to delineate also the other basins. In addiction, consider to provide at least a small buffer along the ridge line in order to allow the algorithm to work properly. </div>
<div><br></div><div>Regards,</div><div>madi</div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>
As far as I know, there's no attribute table for raster maps that hold<br>
descriptions of the generated basins, only information on the contents of<br>
each cell.<span class="HOEnZb"><font color="#888888"><br>
<br>
Rich</font></span><div class="HOEnZb"><div class="h5"><br>
<br>
______________________________<u></u>_________________<br>
grass-user mailing list<br>
<a href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/grass-user" target="_blank">http://lists.osgeo.org/<u></u>mailman/listinfo/grass-user</a><br>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br>Ing. Margherita Di Leo, Ph.D.<br>