<div dir="ltr"><div dir="ltr"><div>Hi Shane,<br></div><div dir="ltr"><br><div class="gmail_quote"><div dir="ltr">On Fri, Sep 14, 2018 at 5:03 PM Shane Carey <<a href="mailto:careyshan@gmail.com">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 dir="ltr"><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">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">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 class="gmail_quote"><br></div><div class="gmail_quote">Vaclav<br></div><div class="gmail_quote"><div><br></div></div></div></div></div>