Introduzione 

Le applicazioni della tecnologia LiDAR (Light Detection And Ranging, denominata anche Airborne Laser Scanner o ALS) nelle attività forestali sono giunte a un avanzato stadio di maturazione per ciò che riguarda la sperimentazione e la messa a punto di procedure di utilizzo pratico nel settore forestale. Il trasferimento di tali procedure alla pratica operativa è avvenuto rapidamente in Nord America e Nord Europa ([21], [11], [33], [22], [28], [34]), mentre procede più lentamente nel nostro Paese ([19]), soprattutto a causa della relativa scarsità di riprese LiDAR su estensioni sufficientemente ampie, legata al costo di acquisizione, in particolare per rilievi ad alta densità di punti ([23]). Laddove esistenti, tali riprese non sono inoltre quasi mai specificatamente orientate al monitoraggio del patrimonio forestale ([5]).

La disponibilità su scala nazionale dei dati LiDAR acquisiti nell’ambito del Piano Straordinario di Telerilevamento ([18]), pur se rilevati a risoluzioni di scansione non elevate in campagne finalizzate alla realizzazione di modelli digitali del terreno ad alta risoluzione, ha comunque consentito negli ultimi anni la messa a punto di metodi in grado di utilizzare proficuamente questa tipologia di dati ([4], [25]). Da essi è infatti possibile estrarre informazioni a livello di popolamento forestale che possono essere impiegate per ridurre in maniera consistente il numero di rilievi a terra negli inventari ([2]).

Obiettivi 

Uno degli obiettivi del progetto AlForLab è stato quindi l’impiego di tali dati LiDAR di pubblica disponibilità per la caratterizzazione dendrometrica dei boschi della Calabria, in particolare per la stima “automatica” del volume e delle fitomasse per unità di superficie ([1]) basata esclusivamente sulla disponibilità di dati LiDAR.

In particolare due sono state le finalità principali di questa attività di ricerca: (a) la messa a punto di un sistema di modelli di stima del volume e della biomassa legnosa, applicabili alle principali formazioni forestali della Calabria, utilizzando come variabili predittive alcune “metriche” elaborate da dati LiDAR; (b) l’implementazione di tali modelli in procedure informatiche facilmente utilizzabili dai tecnici forestali e dalle amministrazioni preposte alla pianificazione e gestione del patrimonio forestale regionale. Tali obiettivi sono stati costantemente coordinati con le altre attività AlForLab oggetto di altri contributi del presente numero speciale, in particolare quelle finalizzate alla realizzazione delle tariffe regionali calabresi di cubatura dei soggetti arborei ([26]), alla messa a punto di modelli di stima campionaria al suolo del volume legnoso tramite rilievi relascopici e anche alla realizzazione del sistema di supporto alle decisioni per la gestione forestale CFOR ([24]).

Materiali e metodi 

Area di studio

Nel progetto AlForLab ([1]) sono stati individuati quattro grandi scenari territoriali di studio, ciascuno designato come area forestale di manifestazione caratteristica di una delle quattro principali tipologie forestali regionali: le Faggete nell’area del Pollino, i cedui di Castagno nella Catena Costiera, le Pinete di Pino laricio nella Sila e le Abetine nelle Serre, per un totale di circa 140 000 ha.

In queste aree si è concentrato il rilievo di dati al suolo integrativi per la costruzione dei modelli. I dati di base derivano invece da una distribuzione estesa a tutta la regione essendo quelli derivanti dai rilievi al suolo della terza fase INFC 2005 ([6]), che hanno consentito di costruire i modelli anche per altre tipologie forestali, come ad esempio i querceti decidui, aventi una certa importanza nella regione, ma poco presenti nei quattro scenari-studio (Fig. 1). L’area di studio è pertanto spazialmente corrispondente a tutte le aree al suolo, per la fase di addestramento dei modelli, e a tutta la copertura regionale disponibile, pari a circa il 90% del territorio, per l’applicazione dei modelli e la produzione delle carte tematiche (Fig. S1 - Materiale supplementare).

Fig. 1 - I quattro scenari studio del progetto AlForLab (in verde chiaro) e i punti di campionamento al suolo (terza fase INFC-2005 più punti AlForLab) impiegati per la costruzione dei modelli LiDAR (in nero).

 Enlarge/Shrink  Open in Viewer

Dati LiDAR

I prodotti LiDAR rilevati nell’ambito del Piano Straordinario di Telerilevamento ambientale - PST e messi a disposizione del progetto AlForLab dal Ministero dell’ambiente e della tutela del territorio e del mare (MATTM) sono stati forniti in formato raster, come DTM (modello digitale del terreno) e DSM (modello digitale della superficie), ripartiti su ritagli (tile) aventi ciascuno estensione di 0°.01 × 0°.01 (latitudine × longitudine), non proiettati (in coordinate geografiche). I voli di ripresa sono stati effettuati in un arco temporale che è andato dal 2008 al 2014, pur se concentrati soprattutto nel 2011 (circa 81% del totale). La risoluzione geometrica del pixel, espressa in unità angolari, era di 0°.0001 e corrispondeva a “circa” 1 m. La risoluzione nativa della nuvola di impulsi LiDAR era stata in media di 1.6 pp m^-2, registrando fino a 4 echi di ritorno dal singolo impulso in uscita dal sensore. Per poter procedere alle successive elaborazioni è stato necessario riproiettare i tile su sistema di riferimento UTM-WGS84 33N (coordinate metriche) e ricampionarli a una risoluzione esatta del pixel di 1 × 1 m. Questa fase, dapprima realizzata con concatenazione semi-automatica di funzionalità in ambito QGIS, è stata successivamente ingegnerizzata e ottimizzata in ambiente R (per renderla usufruibile anche all’utenza), avendo cura di gestire, essenzialmente per sostituzione tramite interpolazione spaziale, i “no-data” già presenti nel dataset e quelli che si venivano a creare nelle zone di confine tra diversi tile a seguito del ricampionamento. Eventuali zone ampie di “no-data”, per le quali l’interpolazione spaziale avrebbe potuto dare luogo a valori non verosimili, sono state lasciate con i valori “no-data” immutati. Le principali caratteristiche tecniche dei voli e dei dati LiDAR sono riportate in Tab. 1.

Il Modello Digitale delle Chiome (CHM) necessario alle successive elaborazioni è stato calcolato come differenza algebrica tra DSM e DTM, ponendo a valore zero tutti i pixel aventi altezza inferiore a due metri, valore-soglia usuale per escludere la vegetazione non arborea ([20]). Porzioni del CHM attribuibili a elementi non arborei, quali pale eoliche, viadotti, costruzioni rurali e altre, sono state escluse dall’analisi, mentre valori eccessivamente elevati del CHM in corrispondenza dello strato arboreo, ritenuti outlier, sono stati modificati tramite interpolazione spaziale con i dati vicini (parametro mediana). Approfondimenti su metodi e procedure adottati per la produzione automatizzata di CHM il più possibile privi di dati aberranti sono riportati in Floris et al. ([9]). Un esempio di CHM è riportato in Fig. S2 (Materiale supplementare).

Rilievi dendrometrici sul terreno

Per mettere a punto i modelli oggetto di questo studio sono stati utilizzati sia dati al suolo INFC (anno di riferimento 2006, 311 osservazioni), sia dati al suolo rilevati appositamente nel progetto AlForLab (anno 2015, 143 osservazioni). Tutte le aree sono state accuratamente georeferenziate tramite posizionamento GPS con ricevitori di classe GIS (Trimble modelli PRO-XR, GeoXT 6000 e Juno SD), come media stazionaria di 100-170 posizioni e correzione differenziale in post-processing.

Per ciascuna delle aree di saggio al suolo è stato calcolato il volume riferito all’ettaro dei soggetti arborei con diametro ≥ 4.5 cm. Nelle aree INFC, soggette a cavallettamento totale e a rilievo ipsometrico secondo le modalità riportate da Di Cosmo et al. ([7]), si è pervenuti a tale volume applicando individualmente a ciascun albero i modelli a doppia entrata INFC ([29]), assegnando a ciascun albero misurato per il diametro l’altezza della curva guida tariffaria corrispondente all’indicatore di tariffa TaCAL ([26]), calcolato sugli alberi modello dell’altezza rilevati nella specifica area di saggio nella campagna INFC (in media 8 soggetti per unità campionaria). Pertanto i valori di altezza ottenuti possono essere leggermente diversi da quelli utilizzati per le stime INFC. Si ricorda che il volume cui si fa riferimento nei modelli INFC e nel sistema tariffario TaCAL è quello del fusto e dei rami, in entrambi i casi fino alla sezione diametrica minima di 5 cm. Il coefficiente di espansione dei dati unitari all’ettaro è dato dal rapporto 1/a, con a pari alla superficie in ha dell’area di saggio (raggio 4 m per gli alberi con d ≥ 4.5 cm, raggio 13 m per gli alberi con d ≥ 9.5 cm).

Nelle unità campionarie integrative AlForLab, che erano di tipo relascopico con misurazione del diametro di tutti gli alberi inclusi (alberi “IN”) e con rilievo dell’altezza di 8 alberi, è stato calcolato individualmente il volume di tutti gli alberi IN di ciascuna prova di numerazione angolare (PNA), previa assegnazione a ciascun albero dell’altezza della curva guida tariffaria TaCAL corrispondente all’indicatore di tariffa individuato sugli alberi modello misurati in altezza, con procedura analoga a quella adottata per le aree INFC. I dati di volume riferiti all’ettaro sono stati ottenuti mediante “ricostruzione del popolamento virtuale” secondo il metodo di Bitterlich ([3]). L’idoneità delle prove di numerazione angolare per la calibrazione di modelli che usino dati telerilevati, soprattutto quando la variabile in stima sia il volume, è stata dimostrata in diversi studi ([12], [13], [16], [15], [31], [27], [30]).

Le principali statistiche descrittive del volume a ettaro calcolato dai rilievi al suolo, riferite alle aree aventi copertura LiDAR, sono riportate in Tab. 2.

Estrazione delle metriche CHM per ciascuna unità campionaria

In corrispondenza delle coordinate di ogni unità campionaria al suolo (sia INFC che AlForLab), sul supporto LiDAR CHM è stata individuata un’area circolare corrispondente alla superficie dell’unità campionaria individuata al suolo. Nel caso delle aree di saggio INFC è stato individuato un cerchio con raggio fisso di 13 m, mentre, per le aree relascopiche AlForLab, che sono da ritenere “a raggio variabile”, è stata dimensionata e individuata sul CHM un’area circolare di esplorazione relascopica locale con raggio di volta in volta variabile calcolato con la procedura messa a punto da Matérn ([17]) (vedi [27]).

Una volta prodotti i CHM per i singoli tile interessati da punti di campionamento al suolo, sono state estratte con funzionalità GIS di Map Algebra (statistiche zonali), per l’estensione di ciascuna area di saggio (INFC) o PNA (AlForLab), le principali metriche presumibilmente correlate al volume legnoso: media, mediana, minimo, massimo, deviazione standard dei valori di altezza (in metri). Si sottolinea che tali calcoli delle metriche includono tutti i pixel del CHM locale estratto, compresi quelli che evidenziano altezza del CHM pari a zero (echi del suolo non coperto da strato arboreo o con altezza della vegetazione inferiore a 2 m). In tal modo, come da esperienze realizzate in altri studi, le metriche estratte interpretano implicitamente anche la densità locale, parametro che influisce, con le altezze dei soggetti presenti, a “spiegare” in termini previsionali il volume a ettaro riscontrato al suolo.

Analisi dei dati e studio dei modelli di stima del volume

Considerato che le epoche dei rilievi a terra erano sensibilmente diverse tra le due fonti di dati, ed entrambe diverse dal riferimento temporale delle acquisizioni LiDAR, si è attuata una sincronizzazione dei dati di volume rilevati al suolo rispetto all’anno di riferimento 2011, applicando ai due dataset (INFC e AlForLab) valori di incremento corrente di volume (positivo o negativo) ricavati, rispettivamente, dalle stime di incremento corrente INFC (per unità campionaria) per il primo gruppo ([6]) e dagli incrementi rilevati in campo con metodo diretto (prelievo di carotine incrementali con succhiello di Pressler su 8 alberi modello per PNA) per il gruppo delle aree AlForLab.

A tale sincronizzazione è seguita una necessaria selezione riguardante: (a) punti INFC privi di copertura LiDAR PST-MATTM; (b) punti rivelatisi outlier (posizionamento GPS palesemente incongruente, interventi di taglio o schianti o cambi di uso del suolo intervenuti tra l’epoca del rilievo a terra e il volo LiDAR).

È stato di conseguenza possibile utilizzare 240 punti INFC. Dei punti AlForLab rilevati al suolo (PNA relascopiche) è stata fatta una analoga selezione, che ha portato a un set effettivamente utilizzabile di 137 punti sui 143 rilevati. Nel complesso il dataset di calibrazione è risultato quindi costituito da 377 osservazioni.

Un primo passo di analisi è consistito in una esplorazione delle relazioni tra le variabili esplicative potenziali di cui sopra e la variabile dipendente osservata, producendo gli scatterplot tra ciascuna delle metriche CHM prese in considerazione e il volume totale a ettaro (V). Al tempo stesso è stata eseguita una analisi preliminare di regressione lineare multipla su tutto il dataset, senza ripartizione delle osservazioni in base al tipo di formazione forestale, utilizzando come potenziali variabili esplicative tutte le metriche estratte in precedenza dal CHM. La metrica Hmean (media delle altezze CHM nell’area di saggio o PNA) si è confermata essere la variabile di gran lunga più significativa per spiegare in termini previsionali il volume al suolo attraverso dati ALS LiDAR, come peraltro riscontrato in numerose esperienze precedenti ([8], [4], [32]) e dalla letteratura esistente ([14], [31]), soprattutto quando i dati CHM LiDAR derivino da scansioni ALS a bassa risoluzione.

Per tale motivo e tenuto conto che la numerosità non elevata del dataset suggeriva di contenere il numero di variabili predittive da includere nei modelli, tutte le successive analisi riguardano esclusivamente detta variabile Hmean, pur non escludendo che in futuri eventuali affinamenti dei modelli possa essere recuperato il contributo (seppur presumibilmente limitato) di altre variabili esplicative.

Si è quindi provveduto alla ripartizione dell’intero dataset di unità campionarie in subset riconducibili a formazioni di una determinata specie principale prevalente (livello dei modelli di massimo dettaglio), di gruppi di specie consimili anche in riferimento all’area altitudinale di vegetazione (livello di dettaglio intermedio), fino a giungere a una semplice ripartizione conifere/latifoglie/formazioni miste (livello dei modelli di dettaglio minimo). In Tab. 3 sono riportate le statistiche descrittive, riferite alla metrica Hmean, per i diversi subset prodotti.

La stima delle fitomasse

Per la stima delle fitomasse epigee riferite all’ettaro si sarebbe potuto procedere con approccio analogo a quanto fatto per la costruzione dei modelli LiDAR di stima del volume, mettendo in relazione i valori di fitomassa (e non più di volume) a ettaro riscontrati nelle unità campionarie col medesimo set di metriche LiDAR locali del CHM. Infatti nota la specie e le dimensioni (d, h) dei soggetti presenti, al pari del volume individuale, specifiche equazioni INFC 2005 di Tabacchi et al. ([29]) per le fitomasse avrebbero potuto fornire i dati individuali e per sommatoria i dati globali riferiti a ettaro, come nel caso del volume legnoso. In tal caso però, con perequazioni prive di vincoli espliciti di interdipendenza reciproca tra volume e (rispettive) fitomasse per unità di superficie, si sarebbe potuti giungere a determinare modelli potenzialmente incongruenti tra loro.

Gli Autori citati in precedenza riportano peraltro anche i coefficienti medi, diversi per ciascuna specie, di trasformazione del volume individuale nelle rispettive fitomasse anidre (legnosa e totale) e anche rispetto al peso fresco. Nel presente studio sono stati adottati tali coefficienti medi e il passaggio da volume a fitomasse è avvenuto calcolando un coefficiente di trasformazione ponderato per ciascun modello (gruppo specifico), sulla base delle composizioni medie in specie osservate nei punti di campionamento al suolo afferenti al modello/gruppo stesso. Si è ritenuto potesse essere utile fornire, oltre ai due coefficienti relativi alle componenti di peso secco, anche il coefficiente di trasformazione in peso fresco del volume legnoso, facendo ricorso anche in questo caso ai valori riportati in Tabacchi et al. ([29]). Dal punto di vista operativo, per ottenere una descrizione digitale (“carta”) della fitomassa presente in un determinato ambito, è sufficiente moltiplicare i valori della carta del volume a ettaro per il coefficiente idoneo riportato in Tab. S1 (Materiale supplementare).

La produzione di carte tematiche

Le funzioni di stima messe a punto, dirette per i volumi e basate su coefficienti medi per le fitomasse, sono state implementate in ambiente GIS open source (QGIS), per la produzione di carte tematiche del volume e delle fitomasse a ettaro. Tutte le elaborazioni necessarie a tal fine sono state automatizzate tramite script di R eseguibili in ambiente QGIS. Le informazioni indispensabili richieste all’utente per la produzione di una carta del volume (e delle fitomasse) sono l’area di interesse (in formato shapefile) e il livello di dettaglio (alto, medio, basso) per il quale si desidera avere le stime. Altre informazioni facoltative sono uno shapefile dei comparti assestamentali (particelle, comprese, eventuali lotti di utilizzazione) che, se presente, consente l’evidenziamento delle stime anche a livello di comparto, con calcolo di valori medi e totali per ciascuno di essi, e uno shapefile delle formazioni forestali classificate secondo codifica AlForLab, in assenza del quale viene comunque impiegata la classificazione Corine Land Cover IV livello. Un esempio di interfaccia utente del modulo R che guida alla produzione delle carte del volume e delle fitomasse è riportato in Fig. S3 (Materiale supplementare). Si sottolinea che tale possibilitàè conseguibile in qualsiasi ambito territoriale boschivo della Calabria, purché dotato di copertura LiDAR MATTM-PST (oltre il 90% del territorio).

Risultati e discussione 

Fig. 2 - Andamenti delle funzioni di regressione di V sulla variabile esplicativa LiDAR Hmean, per i modelli a dettaglio meno elevato: formazioni di Conifere e di Latifoglie.

 Enlarge/Shrink  Open in Viewer

Nello sviluppo dei modelli in studio la variabile in stima (dipendente) era ovviamente il volume a ettaro (V, m³ ha^-1) riscontrato al suolo nelle unità campionarie. Osservando gli andamenti di V in funzione di Hmean (Fig. 2), tenuto conto di numerose precedenti esperienze di modellizzazione di questa relazione, si è scelto di adottare una struttura funzionale dei modelli del tipo (eqn. 1):

dove i coefficienti a e b sono diversi per ciascun tipo di formazione forestale considerata.

Questa scelta presenta il vantaggio di adottare un’unica struttura funzionale comune a tutti i modelli di stima, con una funzione dotata di elasticità (assume andamento rettilineo per b=1) e con positive ripercussioni anche sulla facilità di implementazione in procedure di calcolo automatico e produzione di layer cartografici in continuo.

Una preliminare analisi statistica, effettuata imponendo alla struttura funzionale di cui sopra anche la presenza di un parametro di intercetta, ha messo in evidenza (per tutti i tipi di composizione) che tale parametro non è significativamente diverso da zero (se la metrica locale CHM Hmean è pari a zero, il volume a ettaro stimato è a sua volta pari a zero), confermando definitivamente la possibilità di adottare una configurazione del modello come indicato nella eqn. 1.

L’elenco dei tipi di formazione e i relativi coefficienti dei modelli messi a punto sono riportati in Tab. 4. Per i modelli a dettaglio elevato, i grafici con le funzioni di perequazione sono riportati in Fig. 3.

Fig. 3 - Andamenti delle funzioni di regressione di V sulla variabile esplicativa LiDAR Hmean, per i modelli di stima LiDAR del volume a ettaro a dettaglio elevato.

 Enlarge/Shrink  Open in Viewer

Si osservi che i modelli, in particolare quelli a livello di dettaglio elevato, sono costruiti su numerosità di osservazioni a terra in alcuni casi esigue (ad es., Eucalitti e altri gruppi), ragione per cui le stime prodotte avranno in alcuni casi margini di incertezza ampi; al tempo stesso una certa prudenza andrà posta nel determinare l’intervallo di validità dei modelli stessi, che non dovrà in linea di principio estendersi troppo dall’estremo superiore dei valori di Hmean realmente osservati e desumibili dalla Tab. 3 e dalla Fig. 3.

Le incertezze di stima dei modelli realizzati, espresse dai parametri root mean square error (RMSE) ed errore standard di stima (ESS%), sono riportate in Tab. 5. I valori di ESS% sono abbastanza in linea con quelli ottenuti in altri studi nei quali siano state adottate analoghe modalità di stima ([12], [13], [31], [27]). Si noti il fisiologico aumento di dette incertezze nel passare dai modelli completi di tutte le osservazioni usati in addestramento a quelli della cross-validazione. Tale aumento di errore è più marcato per le specie con pochi dati di calibrazione al suolo (Eucalitteti, Querceti sempreverdi), nei quali una singola osservazione può incidere maggiormente sulla predizione, mentre è minore e spesso trascurabile nei modelli delle specie più rappresentate da osservazioni al suolo: nelle Faggete, ad esempio, ESS% passa da 29.58% a 30.25%, restando di fatto invariato ai fini applicativi.

Una considerazione va fatta al riguardo di tali valori di incertezza: essi sono riferiti a stime su superfici analoghe a quelle usate nella costruzione dei modelli, cioè le aree di saggio al suolo, che nel caso in studio hanno estensione di alcune centinaia di metri quadrati. I principi di propagazione degli errori di stima per aggregazione ([10]), confermati da alcuni studi recenti in analogo ambito ([27]), mostrano che, per modelli non distorti da errori sistematici, all’aumento delle dimensioni dell’area stimata corrisponde una rapida diminuzione dell’incertezza di stima che già alle superfici tipiche delle unità gestionali sub-particellari ad elevato dettaglio spaziale (1-4 ha) si attesta su valori ampiamente compatibili con l’uso pratico delle stime stesse (4-8%).

Per quanto attiene alla realizzazione delle carte tematiche, è opportuno che il livello di dettaglio al quale vengono prodotte sia congruente con la risoluzione e l’accuratezza della compartimentazione tipologica: l’applicazione di modelli di singola specie a una classificazione delle specie poco precisa potrebbe dar luogo a stime distorte. Di conseguenza, in assenza di una carta dettagliata e accurata delle formazioni forestali (congruente con i tipi di modelli impiegati) è preferibile scegliere un grado medio o basso di dettaglio. Un esempio di carte in continuo (raster) del volume e della fitomassa legnosa, riferite a una proprietà comunale dell’estensione di circa 4300 ha, è riportato in Fig. 4, mentre in Fig. 5 è presentata una cartografia della stessa area ma con i valori medi a ettaro e totali aggregati sui poligoni delle diverse formazioni forestali.

Fig. 4 - Piani informativi raster GIS del volume a ettaro e delle fitomasse a ettaro nel comune di Stilo (Serre calabresi): volume a ettaro (in alto a sinistra), fitomassa anidra di fusto e rami grossi (in alto a destra), fitomassa anidra epigea totale (in basso a sinistra), fitomassa fresca di fusto e rami grossi (in basso a destra). Risoluzione 30 m.

 Enlarge/Shrink  Open in Viewer

Fig. 5 - Piano informativo vettoriale GIS del volume a ettaro per le diverse formazioni forestali (comune di Stilo). Il valore all’interno di ogni poligono è il volume legnoso totale stimato per il comparto (m³).

 Enlarge/Shrink  Open in Viewer

Conclusioni 

I modelli di stima presentati in questo studio rappresentano un primo esempio di utilizzo di dati pubblicamente disponibili (INFC-CFS, PST-MATTM) per la realizzazione di prodotti di stima automatica su base LiDAR della consistenza georeferenziata delle risorse forestali a supporto delle attività forestali, validi su ampia scala (inventario forestale regionale), ma configurabili, senza sostanziale variazione della precisione, a dettaglio territoriale più elevato, ad esempio nel caso di inventari forestali di comprensorio o Parco, di piano di assestamento forestale a scala aziendale o sovraziendale.

Una naturale evoluzione dello studio consiste nell’implementare le procedure automatizzate di produzione di piani tematici GIS dei volumi e delle fitomasse, già messe a disposizione in questo caso come script di R (Box S1 in Materiale supplementare), in codice nativo incorporato in QGIS, facilitandone l’uso e ampliando in questo modo la platea di utenti potenziali. Tale attivitàè già in corso e costituirà una componente del citato sistema di supporto alle decisioni CFOR.

Ringraziamenti 

Questo studio è stato realizzato nell’ambito del Progetto “ALForLab” (PON03PE_00024_1) cofinanziato dal Programma Operativo Nazionale Ricerca e Competitività (PON R&C) 2007-2013, attraverso il Fondo Europeo di Sviluppo Regionale (FESR) e risorse nazionali (Fondo di Rotazione (FDR) - Piano di Azione e Coesione (PAC) MIUR).

Gli autori sono particolarmente grati al Corpo Forestale dello Stato per aver messo a disposizione dello studio i dati nativi dei rilievi inventariali INFC 2005 della Calabria, ed al Ministero dell’Ambiente e della Tutela del Territorio e del Mare, per aver messo a disposizione dello studio i dati LiDAR rilevati nell’ambito del Piano Straordinario di Telerilevamento Ambientale.

References

Supplementary Material

Tab. S1 - Coefficienti moltiplicatori di trasformazione del volume a ettaro in fitomassa a ettaro, calcolati sulle composizioni medie ponderate dei gruppi modellistici AlForLab.

Fig. S1 - Copertura voli LiDAR della Regione Calabria (Piano Straordinario di Telerilevamento ambientale - MATTM 2014).

Fig. S2 - Ortofotocarta (a sinistra) e CHM LiDAR (a destra) del comune di Stilo (Serre calabresi).

Fig. S3 - Interfaccia utente del modulo di QGIS (implementato come script R) per la stima del volume e delle fitomasse e la produzione dei relativi piani informativi digitali.

Box S1 - Script di R per la produzione di piani tematici GIS dei volumi e delle fitomasse.