<html><body><div style="color:#000; background-color:#fff; font-family:HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;font-size:16px"><div id="yui_3_16_0_1_1423868067048_4283" dir="ltr"><span id="yui_3_16_0_1_1423868067048_4454">Tom,</span></div><div id="yui_3_16_0_1_1423868067048_4283" dir="ltr"><span><br></span></div><div id="yui_3_16_0_1_1423868067048_4283" dir="ltr"><span id="yui_3_16_0_1_1423868067048_4288">Hopefully the higher spatial resolution data will also have higher vertical resolution. </span></div><div id="yui_3_16_0_1_1423868067048_4283" dir="ltr"><span><br></span></div><div id="yui_3_16_0_1_1423868067048_4283" dir="ltr"><span id="yui_3_16_0_1_1423868067048_4523">That said, if you have any idea where the stream channels are, for example a vector file, in the flat areas you can "burn" them into the raster using r.mapcalc. So you can rasterize vector streams (v.to.rast) and then thin the raster streams (r.thin) and use r.mapcalc to subtract a constant value where the raster stream map is. If you do not thin the raster you will get an error, as there will be ambiguous flow paths (endless loops).</span></div><div id="yui_3_16_0_1_1423868067048_4283" dir="ltr"><span><br></span></div><div id="yui_3_16_0_1_1423868067048_4283" dir="ltr"><span id="yui_3_16_0_1_1423868067048_4564">QGIS has a nice feature in its vector to raster tool. It will do the burn for you directly from the vector file into your DEM. You should make a backup of your DEM though first, since it edits the original file. You just add a field to the vector streams for a desired elevation. Then you select that elevation field in the vector to raster dialog and select your DEM as the output raster. After saying that yes it is OK to overwrite your DEM, it will burn the vector elevation. I have only used this to edit small areas on a raster, so I can not tell you if the burnt channel is sufficiently thinned or not. </span>One of the big problems with using flow algorithms is that bridges and roads will frequently act like dams, causing the flow to be rerouted outside of the actual channel. To prevent this you can digitize a vector map of places where you want to cut through these "dams" and set a field with desired elevations as described above.</div><div id="yui_3_16_0_1_1423868067048_4283" dir="ltr"><br></div><div id="yui_3_16_0_1_1423868067048_4283" dir="ltr"><span id="yui_3_16_0_1_1423868067048_4584">-Thayer</span></div><div id="yui_3_16_0_1_1423868067048_4283" dir="ltr"><span><br></span></div><br> <div style="font-family: HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif; font-size: 16px;" id="yui_3_16_0_1_1423868067048_4286"> <div style="font-family: HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif; font-size: 16px;" id="yui_3_16_0_1_1423868067048_4285"> <div dir="ltr" id="yui_3_16_0_1_1423868067048_4284"> <hr size="1" id="yui_3_16_0_1_1423868067048_4340"> <font size="2" face="Arial" id="yui_3_16_0_1_1423868067048_4287"> <b><span style="font-weight:bold;">From:</span></b> Thomas Adams <tea3rd@gmail.com><br> <b><span style="font-weight: bold;">To:</span></b> Thayer Young <thayeray@yahoo.com> <br><b><span style="font-weight: bold;">Cc:</span></b> "grass-user@lists.osgeo.org" <grass-user@lists.osgeo.org> <br> <b><span style="font-weight: bold;">Sent:</span></b> Friday, February 13, 2015 5:48 PM<br> <b><span style="font-weight: bold;">Subject:</span></b> Re: grass-user Digest, Vol 106, Issue 34<br> </font> </div> <div class="y_msg_container" id="yui_3_16_0_1_1423868067048_4293"><br><div id="yiv8685915947"><div id="yui_3_16_0_1_1423868067048_4292"><div dir="ltr" id="yui_3_16_0_1_1423868067048_4291">Thayer,<div id="yui_3_16_0_1_1423868067048_4396"><br clear="none"></div><div id="yui_3_16_0_1_1423868067048_4307">Thank you for the suggestion; I just tried it with the same result. I think the area just happens to be a difficult area, that is pretty flat…</div><div id="yui_3_16_0_1_1423868067048_4464"><br clear="none"></div><div id="yui_3_16_0_1_1423868067048_4308">Cheers!</div><div id="yui_3_16_0_1_1423868067048_4465">Tom</div></div><div class="yiv8685915947gmail_extra" id="yui_3_16_0_1_1423868067048_4300"><br clear="none"><div class="qtdSeparateBR"><br><br></div><div class="yiv8685915947yqt9771608455" id="yiv8685915947yqt86943"><div class="yiv8685915947gmail_quote" id="yui_3_16_0_1_1423868067048_4299">On Fri, Feb 13, 2015 at 1:54 PM, Thayer Young <span dir="ltr"><<a rel="nofollow" shape="rect" ymailto="mailto:thayeray@yahoo.com" target="_blank" href="mailto:thayeray@yahoo.com">thayeray@yahoo.com</a>></span> wrote:<br clear="none"><blockquote class="yiv8685915947gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;" id="yui_3_16_0_1_1423868067048_4298"><div id="yui_3_16_0_1_1423868067048_4297"><div style="color:#000;background-color:#fff;font-family:HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;font-size:16px;" id="yui_3_16_0_1_1423868067048_4296"><div dir="ltr" id="yui_3_16_0_1_1423868067048_4295"><span id="yui_3_16_0_1_1423868067048_4294">I don't know if you have tried this yet, but you may also want to look at r.terraflow. It does both D8 (single flow direction) and multiple flow direction (flow is partitioned, according to the steepness of slope, to all directions that are lower than the central cell). Supposedly you can also switch from MFD to SFD once flow exceeds a threshold, but I have not been successful at doing this. It would be a simple simulation of channelization though. <br clear="none"></span></div><div dir="ltr" id="yui_3_16_0_1_1423868067048_4466"><span><br clear="none"></span></div><div dir="ltr" id="yui_3_16_0_1_1423868067048_4468"><span id="yui_3_16_0_1_1423868067048_4467">Just set the D8 flag in the options tab, otherwise it will give you MFD by default. <br clear="none"></span></div><div dir="ltr" id="yui_3_16_0_1_1423868067048_4469"><br clear="none"><span></span></div><div dir="ltr" id="yui_3_16_0_1_1423868067048_4471"><span id="yui_3_16_0_1_1423868067048_4470">-Thayer</span></div><div dir="ltr"><span><br clear="none"></span></div><div dir="ltr"><span></span></div><br clear="none"> <div style="font-family:HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;font-size:16px;" id="yui_3_16_0_1_1423868067048_4474"> <div style="font-family:HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;font-size:16px;" id="yui_3_16_0_1_1423868067048_4473"> <div id="yui_3_16_0_1_1423868067048_4472"> <hr size="1"><font face="Arial"><b><span style="font-weight:bold;"></span></b></font>Date: Fri, 13 Feb 2015 11:12:30 -0700<br clear="none">From: Thomas Adams <<a rel="nofollow" shape="rect" ymailto="mailto:tea3rd@gmail.com" target="_blank" href="mailto:tea3rd@gmail.com" id="yui_3_16_0_1_1423868067048_4513">tea3rd@gmail.com</a>><br clear="none">To: "<a rel="nofollow" shape="rect" ymailto="mailto:grass-user@lists.osgeo.org" target="_blank" href="mailto:grass-user@lists.osgeo.org" id="yui_3_16_0_1_1423868067048_4475">grass-user@lists.osgeo.org</a>" <<a rel="nofollow" shape="rect" ymailto="mailto:grass-user@lists.osgeo.org" target="_blank" href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a>><br clear="none">Subject: [GRASS-user] Question about r.watershed and flow accumulation<br clear="none"> grid<br clear="none">Message-ID:<br clear="none"> <<a rel="nofollow" shape="rect" ymailto="mailto:CAGxgkWgefeTtcjZWKi3JP8wj79Eb0j3fZsf6hGzT8m9-dgce6w@mail.gmail.com" target="_blank" href="mailto:CAGxgkWgefeTtcjZWKi3JP8wj79Eb0j3fZsf6hGzT8m9-dgce6w@mail.gmail.com">CAGxgkWgefeTtcjZWKi3JP8wj79Eb0j3fZsf6hGzT8m9-dgce6w@mail.gmail.com</a>><br clear="none">Content-Type: text/plain; charset="utf-8"<br clear="none"><br clear="none">Hello all!<br clear="none"><br clear="none">I'm making use of the flow accumulation grid in GRASS 6.4.5 generated from<br clear="none">r.watershed using the SFD (D8) flow algorithm. The DEM has a 250m spatial<br clear="none">resolution. What I'm getting is a break in the flow accumulation in a few<br clear="none">locations which is causing me serious problems with subsequent processing<br clear="none">(with help from some here, I have put together some scripting to generate a<br clear="none">pixel connectivity file for a distributed hydrologic model).<br clear="none"><br clear="none">Besides going to a higher resolution DEM, are there any thoughts as to how<br clear="none">I can eliminate these flow accumulation breaks?<br clear="none"><br clear="none">Thank you,<br clear="none">Tom<br clear="none"><br clear="none">--<br clear="none">-------------- next part --------------<br clear="none">An HTML attachment was scrubbed...<br clear="none">URL: <<a rel="nofollow" shape="rect" target="_blank" href="http://lists.osgeo.org/pipermail/grass-user/attachments/20150213/9306343f/attachment-0001.html">http://lists.osgeo.org/pipermail/grass-user/attachments/20150213/9306343f/attachment-0001.html</a>><br clear="none"><br clear="none"></div> </div> </div> </div></div></blockquote></div></div><br clear="none"><br clear="all"><div><br clear="none"></div>-- <br clear="none"><div class="yiv8685915947gmail_signature">Thomas E Adams, III<div>718 McBurney Drive</div><div>Lebanon, OH 45036</div><div><br clear="none"></div><div>1 (513) 739-9512 (cell)</div><div><br clear="none"></div></div>
</div></div></div><br><br></div> </div> </div> </div></body></html>