[GRASS-user] slope along defined section (e.g. road, stream,)

[Hamish]

> Tim Michelsen wrote:
> > 3) add the slope to the attribute table of the stream/road vector

Frank Broniewski wrote:
....
> Why don't you just calculate a slope map for your required region and
> use  d.profile to see your slope values as a profile. Or you could
use
> v.report to see whats under your stream or road.

that gives me an idea, using v.out.ascii + r.profile:


#spearfish dataset
VMAP=railroads
QMAP=slope
OUT=rr_slope_points

v.out.ascii "$VMAP" format=standard | grep '^ [0-9]' | \
   grep -v ^' 1 ' | awk '{print $1 "," $2}' | \
   r.profile -g -c in="$QMAP" --quiet 2> /dev/null | \
   grep -v ' \* ' | cut -f1,2,4,5 -d' ' | \
   v.in.ascii x=1 y=2 out="$OUT" fs=space \
     columns='x double, y double, slope double, GRASSRGB varchar(11)'

d.vect -a "$OUT" size=0 width=2


(line vertices only; bug: includes profiles between line jumps...)


Hamish



      ____________________________________________________________________________________
Be a better friend, newshound, and 
know-it-all with Yahoo! Mobile.  Try it now.  http://mobile.yahoo.com/;_ylt=Ahu06i62sR8HDtDypao8Wcj9tAcJ