<div><div dir="auto">Ok great - thanks. I think I will use r.carve to carve the streams into the DTM first and then use this procedure you have described. Do you with the seed raster (for r.lake) if I convert the stream to raster, will ever cell be a seed point (i.e. Water will be placed into every cell at a height of x?)? If so, this is exactly what I need.</div><div dir="auto"><br></div><div dir="auto">Thanks</div><br><div class="gmail_quote"><div>On Luan 17 MFómh 2018 at 04:22, Vaclav Petras <<a href="mailto:wenzeslaus@gmail.com">wenzeslaus@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div><div>Hi Shane,<br></div><div><br><div class="gmail_quote"><div>On Fri, Sep 14, 2018 at 5:03 PM Shane Carey <<a href="mailto:careyshan@gmail.com" target="_blank">careyshan@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><br><div>Does anyone know is it possible to calculate the elevation above a river channel (actual river network that was digitised as opposed to being extracted from a DTM) from a DTM and create a polygon from it. <br></div><div><br></div><div>I need to calculate heights of 1m and 3m above a river channel on both sides of the channel and create a polygon from it.<br></div></div></blockquote><div><br></div><div>Maybe Height Above Nearest Drainage (HAND) is what you are looking for. It is basically done using:<br></div><div><br></div></div><div class="gmail_quote">r.stream.distance stream_rast=streams direction=drainage elevation=elevation method=downstream difference=above_stream<br></div><div class="gmail_quote"><br></div><div class="gmail_quote">See full workflow here:</div><div class="gmail_quote"><br></div><div class="gmail_quote"><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 class="gmail_quote"><br></div><div class="gmail_quote">You will need to do some post processing with r.lake or r.mapcalc and r.to.vect to get the polygons.<br></div><div class="gmail_quote"><br></div><div class="gmail_quote">Since you already have the streams and you want to use them, you need to convert them to raster and then also get stream direction for r.stream.distance by reclassifying, e.g. with r.mapcalc, output from something like this:<br></div><div class="gmail_quote"><br></div><div class="gmail_quote">v.to.rast in=streams output=streams_dir use=dir</div><div class="gmail_quote"><br></div><div class="gmail_quote">which you can see in context here (different hydrology tool):<br></div><div class="gmail_quote"><br><a href="http://ncsu-geoforall-lab.github.io/geospatial-modeling-course/grass/simwe.html" target="_blank">http://ncsu-geoforall-lab.github.io/geospatial-modeling-course/grass/simwe.html</a><br><br></div><div class="gmail_quote">However, I would be careful about how the streams fit with the terrain you are using to be sure results make sense.</div></div></div></div><div><div><div><div class="gmail_quote"><br></div><div class="gmail_quote">Vaclav<br></div><div class="gmail_quote"><div><br></div></div></div></div></div>
</blockquote></div></div><div dir="ltr">-- <br></div><div class="gmail_signature" data-smartmail="gmail_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><div><span style="color:rgb(61,133,198)"><b><span style="background-color:rgb(255,255,255)"><i>GIS and Data Solutions Consultant</i></span></b></span><br></div></div></div></div></div></div></div>