Markus, <br><br>Thank you. I&#39;m pretty sure this is the way. So I run:<br><br>r.to.vect input=mainchannel output=points feature=point --o<br>v.to.db map=points option=coor col=x,y (like the man page suggests)<br><br>but I get a number of errors like:<br>
<br>Column &#39;x&#39; not found<br>Error in db_execute_immediate()<br><br>where am I wrong now?<br><br><br><br><div class="gmail_quote">2010/6/2 Markus Neteler <span dir="ltr">&lt;<a href="mailto:neteler@osgeo.org">neteler@osgeo.org</a>&gt;</span><br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div class="im"><br>
</div>If you treat them as vectors, perhaps this could help?<br>
<br>
<br>
v.to.db<br>
...<br>
             slope: slope steepness of vector line or boundary<br></blockquote><div><br>I calculate slope as <br><br>S = (z2-z1)/(sqrt((x2-x1)^2+(y2-y1)^2))*100<br><br>where x1,y1 are the coords of point 1, x2,y2 are the coords of point 2, and z1, z2 are the elevations in the raster map associate at the points 1 and 2 respectively.<br>
<br>Margherita<br> </div><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
?<br>
<font color="#888888"><br>
Markus<br>
</font></blockquote></div><br>