#!/bin/sh
# geoprocessing_spearfish_aree_linee.sh
# geoprocessing sullo spearfish: combinazione di aree e linee
# per GRASS 6.0/6.1
# C Paolo Zatelli 20/3/2006
# distribuibile sotto licenza GNU

eval `g.gisenv`
: ${GISDBASE?} ${LOCATION_NAME?} ${MAPSET?}
LOCATION="$GISDBASE/$LOCATION_NAME/$MAPSET"

# metto la regione a quella di default
g.region -d

# apro il monitor x0
d.mon start=x0

# visualizzo
d.erase

# fields
d.vect map=fields type=area,face display=shape icon=basic/x size=8 layer=1 color=black fcolor=gray

# fields D.Portillo#1
d.vect map=fields type=point,line,boundary,centroid,area,face display=shape,attr attrcol=label where="label='D. Portillo#1'" icon=basic/x size=8 layer=1 color=yellow fcolor=gray llayer=1 lcolor=red bgcolor=none bcolor=none lsize=8 font=romans xref=left yref=center 

# roads
d.vect map=roads color=red

echo "Premere invio per continuare"
read IN

# crea una mappa vettoriale con il campo di D. Portillo#1
v.extract input=fields@PERMANENT output=fields_portillo1 type=point,line,boundary,centroid,area,face layer=1 new=-1 'where=label='\''D. Portillo#1'\''' --overwrite

# fields_portillo1 in giallo
d.vect map=fields_portillo1 type=point,line,boundary,centroid,area,face display=shape,attr attrcol=label icon=basic/x size=8 layer=1 color=black fcolor=yellow llayer=1 lcolor=red bgcolor=none bcolor=none lsize=8 font=romans xref=left yref=center 

# roads
d.vect map=roads color=red

echo "Premere invio per continuare"
read IN

######################
# select: seleziona gli elementi di un tematismo vettoriale che intersecano un'area data senza ritagliarla
# crea il file roads_portillo1 che contiene la strada che passa nel campo di D. Portillo #1
v.select ainput=roads@PERMANENT atype=point,line,boundary,centroid,area alayer=1 binput=fields_portillo1 btype=point,line,boundary,centroid,area blayer=1 output=roads_portillo_select operator=overlap --overwrite

# visualizzazione
d.erase
# fields
d.vect map=fields type=area,face display=shape icon=basic/x size=8 layer=1 color=black fcolor=gray
# fields_portillo1 in giallo
d.vect map=fields_portillo1 type=point,line,boundary,centroid,area,face display=shape,attr attrcol=label icon=basic/x size=8 layer=1 color=black fcolor=yellow llayer=1 lcolor=black bgcolor=none bcolor=none lsize=8 font=romans xref=left yref=center 

# roads
d.vect map=roads
# la strada selezionata in giallo
d.vect map=roads_portillo_select color=red

# DB associato
echo "select * from roads_portillo_select"|db.select

echo "Premere invio per continuare"
read IN

######################
# overlay: ritaglia un tematismo vettoriale su un area data
# ritaglia la strada usando il campo di Portillo #1
v.overlay ainput=roads@PERMANENT atype=line alayer=1 binput=fields_portillo1 btype=area blayer=1 output=roads_portillo_overlay operator=and --overwrite

# visualizzazione
d.erase
# fields
d.vect map=fields type=area,face display=shape icon=basic/x size=8 layer=1 color=black fcolor=gray
# fields_portillo1 in giallo
d.vect map=fields_portillo1 type=point,line,boundary,centroid,area,face display=shape,attr attrcol=label icon=basic/x size=8 layer=1 color=black fcolor=yellow llayer=1 lcolor=black bgcolor=none bcolor=none lsize=8 font=romans xref=left yref=center 

# roads
d.vect map=roads
# la strada selezionata in giallo
d.vect map=roads_portillo_overlay color=red

# DB associato
echo "select * from roads_portillo_overlay"|db.select


echo " "
echo "Fine geoprocessing linee ed aree"
echo " "