#### first part from chapter 5

g.region rast=elev_ned_30m -p

# visibility from a new 32 story tower being built in downtown
echo "642212,224767,165"|v.in.ascii -z out=tower_165m fs=, z=3
r.los elev_ned_30m out=tower_165_los coord=642212,224767 \
      obs=165 max=50000

#### second part from chapter 6
g.region rast=tower_165_los -p
# create a binary viewshed map and convert it to vector map
r.mapcalc "tower_165_losbin=if(tower_165_los)"
r.to.vect -s tower_165_losbin out=tower_165_los feat=area
v.select ain=poi_names_wake bin=tower_165_los out=tower_poi
d.erase
d.his h=tower_165_los i=elevation_shade
d.vect streets_wake
d.vect tower_165m siz=10 col=orange icon=basic/marker
d.vect tower_poi icon=basic/circle fcol=yellow
d.vect tower_poi disp=attr attrcol=class