REM 
REM
REM Inizializzazione
REM
g.mapset -c mapset=EcoreteLogica
g.region rast=uds_cod1 res=100
r.mask rast=uds_cod1

REM Aree inutilizzabili
REM Selezione artificiale senza strade
r.mapcalc e="fram1a=if(clc_cod1==1,1,null())" --o
REM Fondo in un unico layer le strade principali
r.patch i=StraPrin,SuperStra,AutoStra o=fram1b --o
REM Fonde i due dati
r.patch i=fram1a,fram1b o=fram1 --o

REM Zone vicine aree inutilizzabili 
r.buffer --o input=fram1 output=buff_fram1 dist=300
r.mapcalc e="fram2=if(buff_fram1==1,null(),if(buff_fram1==2,1))" --o

REM Ecotone interno ECO2 e aree core ECO1
r.mapcalc e="not_naturali=if(uds_cod1>2,null(),1)" --o
r.mapcalc e="naturali=if(uds_cod1>2,1,null())" --o
r.buffer --o input=not_naturali output=buffer_not_naturali dist=200
r.mapcalc e="eco2=if(buffer_not_naturali==2,1,null())" --o
r.mapcalc e="eco1=if(naturali==1 & isnull(eco2),1,null())" --o

REM ecotone esterno ECO3
r.buffer --o input=naturali output=buffer_naturali dist=200
r.mapcalc e="eco3=if(buffer_naturali==2 & isnull(fram1) & isnull(fram2),1,null())" --o

REM MATrice agricola
r.mapcalc e="mat=if(isnull(eco1)&&& isnull(eco2) &&& isnull(eco3) &&& isnull(fram1) &&& isnull(fram2) &&& uds_cod1==2,1,null())" --o

REM Mappa complessiva

REM ricodifico le strade in ordine di valenza ecologica
r.mapcalc e="eco1a=eco1*1" --o
r.mapcalc e="eco2a=eco2*2" --o
r.mapcalc e="eco3a=eco3*3" --o
r.mapcalc e="mata=mat*4" --o
r.mapcalc e="fram1a=fram1*6" --o
r.mapcalc e="fram2a=fram2*5" --o
REM Combino in una unica mappa
r.patch i=fram1a,fram2a,mata,eco3a,eco2a,eco1a o=ECO --o

REM POST-ELABORAZIONE
REM
REM Creare in GeoDati una cartella chiamata risultati

REM Ricodifico l'uso del suolo assegnando categorie alfanumeriche
REM
r.mapcalc e="uds_recod3=uds_cod3" --o
r.category map=uds_recod3 rules="D:\GeoDati\codici\cod3.txt" separator=":" 
REM Ricodifico la rete ecologica
r.category map=ECO rules="D:\GeoDati\codici\EcoRete.txt" separator=":" 

REM clacolo il report
r.report map=ECO,uds_recod3 units=h -n o="D:\GeoDati\risultati\reteEco.txt" --o
REM una forma alternativa di calcolo del report per importazione dentro excel
REM ATTENZIONE Excel deve essere settato con il separatore dei decimali sul punto
REM Altrimenti per avere il risultato in metri quadrati il dato va diviso per 1000000
REM per avere il risultato in ettari dividere per 1E+10
REM per fare grafici a torta potete lasciare il dato com'è
r.stats i=ECO,uds_recod3 separator=";" -nl o="D:\GeoDati\risultati\reteEco2.csv" --o