REM 
REM
REM Inizializzazione: creo un  nuovo mapset di lavoro, imposto la regione cartografica e la risoluzione
REM
g.mapset -c mapset=Esempio01
g.region rast=uds_cod1 res=5
r.mask rast=uds_cod1
REM
REM Selezionare aree agricole a rischio danno fauna selvatica -> entro un buffer di 200 metri da aree naturali
REM
REM Creo mediante selezione le mappe delle aree naturali e agricole
REM
r.mapcalc e="nat=if(uds_cod1==3,1,null())" --o
r.mapcalc e="agr=if(uds_cod1==2,1,null())" --o
REM
REM Creo un buffer di 200 metri dalle aree naturali
REM
r.buffer --o input=nat output=buff_nat dist=200
REM
REM Creo mediante selezione logica mappa delle aree agricole AND entro buffer 200 metri aree naturali
REM
r.mapcalc e="rischio=if(agr==1 &&& buff_nat==2,1,null())" --o
REM
REM Creo un report delle colture a rischio in ettari
REM
r.report -n map=rischio,uds_cod3 units=h
REM End