[GRASS-user] elevation profile

Hamish hamish_b at yahoo.com
Fri Jul 16 01:42:25 EDT 2010


yet another two methods to try-

(spearfish dataset)

# v.out.ascii | r.profile method

g.region -d
v.in.region sp_box
v.type in=t_powerlines out=t_powerline type=boundary,line
v.overlay ain=t_powerline bin=sp_box out=sp_powerline atype=line op=and
v.build.polylines in=sp_powerline  out=sp_powerline_pl
g.remove v=sp_box,t_powerline,sp_powerline

v.out.ascii in=sp_powerline_pl format=standard | grep  '^ ' | \
  r.profile -g in=elevation.dem res=30



# v.to.rast MASK method

v.type in=t_powerlines out=t_powerline type=boundary,line
v.category in=t_powerline out=t_powerline_with_cat option=add
v.to.rast in=t_powerline_with_cat out=MASK use=val value=1
g.remove v=t_powerline,t_powerline_with_cat

r.out.xyz elevation.dem
g.remove MASK



      


More information about the grass-user mailing list