# chapter 6 - "Shortest path maps"

g.region n=229590 s=221790 w=635010 e=648330 nsres=30 ewres=30 -pa

# fetch 30m DEM
r.in.wms output=elevation_meters mapserver=http://wms.jpl.nasa.gov/wms.cgi layers=us_ned styles=real -o
r.shaded.relief elevation_meters shad=elevation_meters_shaded

d.erase
d.rast elevation_meters_shaded
d.vect streets_hnet
d.vect streets_hnet layer=2
d.vect lakes fcol=blue type=area
d.vect streets_hnet width=4 where="SPEED > 25 AND SPEED < 40" col=176:164:44
d.vect streets_hnet width=2 where="SPEED > 25 AND SPEED < 40" col=255:250:150
d.vect streets_hnet width=6 where="SPEED >=40" col=176:164:44
d.vect streets_hnet width=4 where="SPEED >=40" col=242:191:36
d.vect streets_hnet layer=2 disp=cat
# Rex: 8; Dorothea Dix: 5

# ID as first number, then cat1 and cat2
#echo "1 8 5" | v.net.path streets_hnet out=spath \
#     afcolumn=navcost

d.vect hospitals disp=attr attrcol=NAME bg=white lsize=10
d.vect hospitals fcol=red icon=basic/diamond
d.vect spath col=red width=2