REM
REM ATTENZIONE! cambiare nei nomi file "18" con il vostro numero ambito paesaggistico
REM
REM ATTNZIONE! verificate che il disco esetrno sui cui lavorate abbia assegnato la lettera D:\
REM Atrimenti modifcare nei nomi file
REM
REM Importazione geodato raster DEM
REM
r.import -e input="D:\GeoDati\P18\D18.tif" output="DEM" --o
REM
REM Importazione ed elaborazione Geodato Uso del Suolo Regione Toscana
v.in.ogr -e input="D:\GeoDati\P18\UdsPae18.shp" output="uds" --o
REM Aggiungo un nuovo campo al database
v.db.addcolumn map=uds@PERMANENT columns="uso_cod3 INT"
REM converto il campo lv3_13 in intero
v.db.update map=uds@PERMANENT column=uso_cod3 query_column="CAST(lv3_13 AS int)"
REM Aggiungo un nuovo campo al database
v.db.addcolumn map=uds@PERMANENT columns="uso_cod2 INT"
REM converto il campo lv2_13 in intero
v.db.update map=uds@PERMANENT column=uso_cod2 query_column="CAST(lv2_13 AS int)"
REM Aggiungo un nuovo campo al database
v.db.addcolumn map=uds@PERMANENT columns="uso_cod1 INT"
REM converto il campo lv1_13 in intero
v.db.update map=uds@PERMANENT column=uso_cod1 query_column="CAST(lv1_13 AS int)"
REM Definisco la regione cartografica e imposto una risoluzione di 5 metri
g.region vect=uds res=5
REM Conversione da vettore a raster
v.to.rast input=uds@PERMANENT output=uds_cod1 use=attr attribute_column=uso_cod1 --o
v.to.rast input=uds@PERMANENT output=uds_cod2 use=attr attribute_column=uso_cod2 --o
v.to.rast input=uds@PERMANENT output=uds_cod3 use=attr attribute_column=uso_cod3 --o
REM
REM
REM
REM Elaborazione del geodato CorineLandCover
v.in.ogr -e input="D:\GeoDati\P18\clc2012ivPae18.shp" output="clc" --o
REM Aggiungo un nuovo campo al database
v.db.addcolumn map=clc@PERMANENT columns="uso_cod INT"
REM converto il campo clc06 in intero
v.db.update map=clc@PERMANENT column=uso_cod query_column="CAST(clc06 AS int)"
REM Aggiungo un nuovo campo al database
v.db.addcolumn map=clc@PERMANENT columns="uso_cod3 INT"
REM estraggo le prime tre cifre di uso_cod
v.db.update map=clc@PERMANENT column=uso_cod3 query_column="CASE WHEN uso_cod>9999 THEN uso_cod/100 WHEN uso_cod>999 then uso_cod/10 ELSE uso_cod end"
REM Aggiungo un nuovo campo al database
v.db.addcolumn map=clc@PERMANENT columns="uso_cod2 INT"
REM estraggo le prime 2 cifre di uso_cod3
v.db.update map=clc@PERMANENT column=uso_cod2 query_column="uso_cod3/10"
REM Aggiungo un nuovo campo al database
v.db.addcolumn map=clc@PERMANENT columns="uso_cod1 INT"
REM estraggo la prima cifre di uso_cod3
v.db.update map=clc@PERMANENT column=uso_cod1 query_column="uso_cod3/100"
REM Definisco la regione cartografica e imposto una risoluzione di 10 metri
g.region vect=uds res=10
REM Conversione da vettore a raster
v.to.rast input=clc@PERMANENT output=clc_cod1 use=attr attribute_column=uso_cod1 --o
v.to.rast input=clc@PERMANENT output=clc_cod2 use=attr attribute_column=uso_cod2 --o
v.to.rast input=clc@PERMANENT output=clc_cod3 use=attr attribute_column=uso_cod3 --o
REM
REM Importazione strade
REM Definisco la regione cartografica e imposto una risoluzione di 10 metri
g.region vect=uds res=10
REM Strade Principali
v.in.ogr -r input="D:\GeoDati\StradePrincipali\StradePrincipali.shp" output="StraPrin" --o
v.to.rast input=StraPrin output=StraPrin use=val --o
REM Strade Principali Importanti
v.in.ogr -r input="D:\GeoDati\StradePrincipali\StradePrincipaliImportanti.shp" output="SuperStra" --o
v.to.rast input=SuperStra output=SuperStra use=val --o
REM AutoStrade
v.in.ogr -r input="D:\GeoDati\StradePrincipali\AutoSuper.shp" output="AutoStra" --o
v.to.rast input=AutoStra output=AutoStra use=val --o
REM StradeLocaliConnessione
v.in.ogr -r input="D:\GeoDati\StradeSecondarie\StradeLocaliConnessione.shp" output="StraLoc" --o
v.to.rast input=StraLoc output=StraLoc use=val --o
REM Strade secondarie
v.in.ogr -r input="D:\GeoDati\StradeSecondarie\StradeSecondarie.shp" output="StraSec" --o
v.to.rast input=StraSec output=StraSec use=val --o