REM inizializzazione
g.mapset mapset="idoneita" -c
g.region rast=uds_cod3 res=10
r.mapcalc e="uds_idoneita=float(uds_cod3)" --o
r.mask rast=uds_idoneita
REM Mappa idoneità Boitani
r.reclass i=uds_idoneita o=IdBoitani rules="D:\GeoDati\script\idoneita\mini_id.txt" --o
REM riscalo mappa Boitani nell'intervallo 0-1
REM la mappa IdBoitani è in formatointeger. Per avere un output decimale debbo utilizzare la funzione float(...)
r.mapcalc e="idoneita=float(IdBoitani)/3" --o
REM Valutazione percolazione
r.neighbors i=idoneita o=permeabilita size=35 -c --o
REM Valutazione frammentazione
r.reclass i=uds_idoneita o=framm rules="D:\GeoDati\script\idoneita\mini_fram.txt" --o
r.neighbors i=framm o=artificiale size=35 -c --o
REM Classi di permabilità
r.mapcalc e="c_permeabilita=if(permeabilita<0.30,10,if(permeabilita<0.60,20,30))" --o
REM Classi di frammentazione
r.mapcalc e="c_artificiale=if(artificiale<0.15,3,if(artificiale<0.60,2,1))" --o
r.cross i=c_permeabilita,c_artificiale o=rete --o
REM
REM |0|category 10; category 1 Barriere critiche in aree a bassa idoneita
REM |1|category 10; category 2 Aree a bassa idoneità con barriere
REM |2|category 10; category 3 Aree a bassa idoneità
REM |3|category 20; category 1 Barriere critiche in aree a media idoneità
REM |4|category 20; category 2 Aree a media idoneità con barriere
REM |5|category 20; category 3 Aree a media idoneità
REM |6|category 30; category 2 Aree ad alta idoneità con barriere
REM |7|category 30; category 3 Aree ad alta idoneità
REM

REM
REM Dalla rete specie specifica alla mappa di vocazione
REM
REM Barriere
REM Selezione artificiale senza strade
r.mapcalc e="edif=if(uds_cod2==11,1,null())" --o
REM Fondo in un unico layer con le strade principali
r.patch i=edif,StraPrin,SuperStra,AutoStra o=fram1 --o
REM calcola la distanza
r.grow.distance i=fram1 d=distArt --o
REM
REM Valutazioni sfocate (Fuzzy)
REM
REM Valutazione fuzzy della permeabilità
r.mapcalc e="fuzzyPerm=graph(permeabilita,0,0,0.6,1,1,1)" --o
REM
REM valutazione fuzzy disturboda strade
REM
r.mapcalc e="fuzzyArt=graph(distArt,0,0,100,0,500,1)" --o
REM
REM aggregazione
REM
r.mapcalc e="vocazione=(fuzzyArt+fuzzyPerm)/2" --o