<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;background-color:#FFFFFF;font-family:Calibri,Arial,Helvetica,sans-serif;">
<p>I'm working with GRASS 7.0.3 on Ubuntu Linux and attempting to create and partition a longest flow path in a catchment comprised of 25km<sup>2</sup> basin, the original DEM is a 5 m resolution.</p>
<p><br>
</p>
<p>Most of the process appears to work well until I attempt to create segments of the vector map of the longest flow path this segment rules file.</p>
<p>====================</p>
<p></p>
<div>P 1 1 -10%</div>
<div>P 2 1 -20%</div>
<div>P 3 1 -30%</div>
<div>P 4 1 -40%</div>
<div>P 5 1 -50%</div>
<div>P 6 1 -60%</div>
<div>P 7 1 -70%</div>
<div>P 8 1 -80%</div>
<div>P 9 1 -90%</div>
<div>========================</div>
<div><br>
<span style="font-size: 12pt;">First I use r.watershed like this to create the various raster maps, </span></div>
<div><span style="font-size: 12pt;"></span><br>
<div>
<div>r.watershed elevation=topmod_demfill threshold=20000 drainage=topmod_drain \</div>
<div> stream=topmod_streams basin=topmod_basin convergence=5 --overwrite</div>
<div><br>
</div>
<div>Next I use r.stream.order to get the strahler numbers and then find the outlet coordinates</div>
<div><br>
</div>
<div>Next I use r.water.outlet as such to get the basin, </div>
<div><br>
</div>
<div>
<div>r.water.outlet input=topmod_drain output=topmod_basin \</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span> coordinates=$E_outflow,$N_outflow --overwrite</div>
<div><br>
</div>
<div>And find the longest flow path, in raster format as such</div>
<div><br>
</div>
<div>r.lfp input=topmod_drain output=topmod_lfp \</div>
<div> coordinates=$E_outflow,$N_outflow --overwrite</div>
<div><br>
</div>
<div>All seems to work well up to this point. I've found two ways to convert the raster map of the longest flow path to a vector map. The first is using v.lfp but it truncates the lfp map, almost in half, only shows the last half in the new vector lfp. Any
ideas why it might do this?</div>
<div><br>
</div>
<div>
<div>v.lfp input=topmod_lfp output=topmod_lfpv coordinates=$E_outflow,$N_outflow \</div>
<div> --verbose --overwrite</div>
<div><br>
</div>
<div>The other way is to use r.to.vect, which appears to work well, now the new vector lfp map (jamrlfp) is the same shape and length as the original raster lfp map. However the next job is to do the segmentation using the above rules file like this,</div>
<div><br>
</div>
<div>
<div>v.segment input=jamrlfp output=outlets rules=$segment_rules \</div>
<div><span class="Apple-tab-span" style="white-space:pre"></span> --overwrite --verbose</div>
<div><br>
</div>
<div>And for some reason it only creates the very first outlet (segment).</div>
<div><br>
</div>
<div>Sorry for the long missive, but I end up with two questions:</div>
<div><br>
</div>
<div>
<ol>
<li>Any suggestions as to why v.lfp is not doing the conversion correctly but r.to.vect appears to do it correctly?</li><li>Any suggestions as to why v.segment is not working correctly and creating the 9 segments correctly?</li></ol>
<div><br>
</div>
</div>
<div><br>
</div>
<div>Thanks,</div>
</div>
</div>
</div>
</div>
</div>
<p></p>
<p><br>
</p>
<div id="Signature">
<div id="divtagdefaultwrapper" style="font-size:12pt; color:#000000; background-color:#FFFFFF; font-family:Calibri,Arial,Helvetica,sans-serif">
<div class="BodyFragment"><font size="2"><span style="font-size:10pt">
<div class="PlainText">Dr. Jim Maas <br>
University of East Anglia <br>
<br>
</div>
<div class="PlainText"></div>
</span></font></div>
</div>
</div>
</div>
</body>
</html>