> 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
_______________________________________________
grass-user mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/grass-user
