[GRASSLIST:4526] Re: roads cut slopes in right angle
Eric G. Miller
egm2 at jps.net
Wed Sep 18 03:10:35 EDT 2002
On Wed, Sep 18, 2002 at 09:20:45AM +0300, orkun wrote:
> Does anyone have an idea about how to query a road map to get a new map
> in which roads cut slopes in right angle ?
Convert vector roads to raster at resolution << than road width.
Create a raster elevation layer just along these road centerlines:
r.mapcalc 'elev_roads = if (!isnull(centerline), elevation)'
Now, grow the elev_roads to the width of the road (r.resample?) --
not sure if that will work well.
r.mapcalc 'cut_elev = if (!isnull(elev_roads2), elev_roads, elevation)'
Don't know if that'll work too well because the r.resample is liable
just create a NULL raster -- not sure how it's handling the NULL's in
begin 664 .signature
More information about the grass-user