<div dir="ltr"><div dir="ltr"><div dir="ltr"><div>These are steps based on:</div><div><br></div><div><a href="https://grasswiki.osgeo.org/wiki/From_GRASS_GIS_novice_to_power_user_(workshop_at_FOSS4G_Boston_2017)#Hydrology:_Estimating_inundation_extent_using_HAND_methodology" target="_blank">https://grasswiki.osgeo.org/wiki/From_GRASS_GIS_novice_to_power_user_(workshop_at_FOSS4G_Boston_2017)#Hydrology:_Estimating_inundation_extent_using_HAND_methodology</a></div><div><br></div><div>You need r.stream.distance module from Addons:</div><br>g.extension r.stream.distance<br><div><br></div><div>Get drainage and streams from your DEM (your carved DEM):<br></div><div><br></div>r.watershed elevation=dem accumulation=flowacc drainage=drainage stream=streams threshold=100000</div><div dir="ltr"><br></div><div>Compute height above nearest drainage/stream (HAND):<br></div><div dir="ltr"><br></div><div dir="ltr">r.stream.distance stream_rast=streams direction=drainage elevation=elevation method=downstream difference=hand</div><div dir="ltr"><br></div><div>Use r.lake not on the original DEM, but on the HAND and start flooding ("lake") from the streams:<br></div><div dir="ltr"><br></div><div dir="ltr">r.lake elevation=hand water_level=3 lake=flood_3m seed=streams<br><div><br></div><div>Convert to vector if desired:<br></div><div><br></div><div>r.to.vect -s input=flood_3m output=flood_3m type=area<br></div><div><br></div><div>The difference to the r.grow+r.mapcalc method [1] is that this uses an addon module (there should be no problem installing it) and that r.grow uses euclidean distance for what is later used for height difference while r.steam.distance follows drainage and further that r.lake floods only the cells accessible to water unlike the r.mapcalc expression which just looks at height. The two methodological differences can be summarized as "not respecting the surrounding terrain enough." Anyway, the r.grow+r.mapcalc method can get you quite far and I would be interested in the comparison (will differ for different terrains).<br></div><div><br></div><div>Best,</div><div>Vaclav<br></div><div><br></div><div>[1] <a href="https://lists.osgeo.org/pipermail/grass-user/2018-September/079134.html" target="_blank">https://lists.osgeo.org/pipermail/grass-user/2018-September/079134.html</a><br></div><div><br></div>On Fri, Sep 21, 2018 at 11:39 AM Shane Carey <<a href="mailto:careyshan@gmail.com" target="_blank">careyshan@gmail.com</a>> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Hi All,</div><div><br></div><div>I have used r.carve to carve out the rivers of a DTM - a really super job. I now need to pour 3meters of water into every cell in the river and see how for this water extends out - onto the floodplain. <br></div><div><br></div><div>I was trying to use r.lake to do this, but unsure as to how r.lake will work to pour 3 meters of water in every cell along the river network. <br></div><div><br></div><div>Any advice on this would be great. It is for the creation of a floodplain.</div><div><br></div><div>Thanks</div><div><br></div><div><br></div><div><div><div class="m_110382330510375694gmail-m_-3297281078989513694gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>Le gach dea ghui,<br></div><div><span style="color:rgb(61,133,198)"><b><span style="background-color:rgb(255,255,255)"><i>Shane Carey</i></span></b></span></div><br></div></div></div></div></div></div></div>
</div></div>
_______________________________________________<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="https://lists.osgeo.org/mailman/listinfo/grass-user" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/grass-user</a></blockquote></div></div></div></div>