Costruire un Digital Twin delle foreste del Québec a partire da 115.000 ortofoto

Le fondamenta prima di tutto

Un digital twin territoriale vive o muore in base ai suoi dati. Fondere migliaia di tessere di ortofoto eterogenee, assorbire sistemi di riferimento delle coordinate incompatibili, colmare le lacune di copertura, armonizzare i colori tra sensori e annate diverse. È il lavoro meno visibile, ma il più importante dal punto di vista strutturale. Finché queste fondamenta non sono solide, nulla di ciò che segue (rendering 3D, navigazione multiscala, analisi territoriale) può reggere.

Questo articolo è un resoconto sullo stato di avanzamento di questa fase, applicata al territorio del Québec.

A che punto siamo oggi

La pipeline di TerraLab: un'unica scena, tutte le fonti

TerraLab-Twin è un virtual globe planetario. Il suo compito è fondere in tempo reale fonti geospaziali eterogenee (ortofoto pubbliche, ortofoto drone ad alta risoluzione, modelli digitali del terreno, livelli vettoriali, immagini satellitari) su un unico terreno 3D navigabile, dalla scala planetaria fino al livello della singola particella.

Tre proprietà definiscono lo stato dell'arte a cui puntiamo:

• Continuità di scala: dall'intero Québec (10.000 km di diagonale) fino a una particella di 50 m, senza interruzioni, senza cambi di sistema di riferimento e senza cuciture ai confini amministrativi.
• Fusione dinamica multi-fonte: più fonti di qualità e risoluzione diverse coesistono sullo stesso terreno e vengono selezionate al volo in base alla zona, alla risoluzione richiesta e alla disponibilità dei dati.
• Astrazione completa a monte: l'utente finale non vede mai un formato di file, un codice EPSG o una zona MTM. Naviga, semplicemente.

Tutto il lavoro descritto di seguito consiste nel rispettare queste tre proprietà su un caso concreto: il territorio del Québec.

Tessere MRNF fuse da TerraLab WorldMapServer. Terreno continuo, senza cuciture visibili tra i fogli originali.

Cosa rende unici i dati del Québec

Il Québec offre uno dei set di dati pubblici più ricchi del Nord America per questo tipo di lavoro. Il Ministero delle Risorse Naturali e delle Foreste del Québec (MRNF, in precedenza MERN) pubblica oltre 115.000 tessere di ortofoto, per un totale di diverse centinaia di gigabyte di dati, a una risoluzione nativa di 20-30 cm per pixel. È circa trenta volte più fine delle immagini satellitari Sentinel-2 usate come complemento.

Ortofoto MRNF a risoluzione nativa (20-30 cm) su un settore forestale del Québec.

Questa ricchezza porta con sé le sue sfide: una famiglia eterogenea di proiezioni cartografiche, più formati di file, immagini acquisite su anni e stagioni diverse. È esattamente il tipo di set di dati per cui un virtual globe è pensato, laddove uno strumento GIS tradizionale chiederebbe all'utente di scegliere una proiezione, una stagione, un formato, e di convivere con le conseguenze.

E sappiamo gestire questa scala. Le 115.000 tessere, le centinaia di gigabyte, i formati misti, i dati multi-annata: l'architettura di TerraLab-Twin è progettata per assorbire nativamente questo volume, su una normale workstation da sviluppatore, senza bisogno di infrastrutture particolari.

Tessera MRNF originale a 20 cm. Alberi, strade e particelle chiaramente visibili prima di qualsiasi elaborazione.

Il punto dolente: le proiezioni cartografiche

Voglio fermarmi un momento su questo punto, perché per chiunque abbia lavorato nel geospaziale, si tratta di un incubo operativo.

I dati pubblici del Québec non arrivano in un'unica proiezione. Arrivano come una famiglia di varianti NAD83: una variante del Québec Lambert (EPSG:32198), varianti più moderne basate su CSRS, e dieci zone MTM (Modified Transverse Mercator) che suddividono il territorio per ragioni cartografiche storiche. In pratica, ogni tessera può arrivare in un sistema di riferimento delle coordinate diverso. Ogni datum ha le sue insidie. Ogni strumento della catena geospaziale ha il proprio comportamento ereditato attorno alla libreria PROJ. È estremamente comune, in produzione, vedere una pipeline fallire silenziosamente perché un'installazione obsoleta restituisce trasformazioni sbagliate di dieci metri o più.

Questo è esattamente ciò che un virtual globe astrae. TerraLab-Twin conosce solo un'unica proiezione finale: WGS84 sull'ellissoide planetario. Tutta la complessità proiettiva viene assorbita a monte, nella pipeline di ingestione. L'utente finale non vede mai un codice EPSG. Vede un territorio 3D continuo e navigabile, dall'intero Québec fino a una singola particella forestale, senza che nulla si rompa ai confini di zona.

Una volta risolto questo problema difficile, tutto il resto diventa possibile: fusione multi-fonte, navigazione continua, 3D in tempo reale.

Strumento interno di tiling: ogni fonte viene riproiettata in WGS84, poi composita nella scena 3D.

Perché WGS84 in un virtual globe, e perché è importante

La scelta di WGS84 sull'ellissoide non è arbitraria. Una proiezione cartografica classica è, per costruzione, un compromesso: appiattisce un oggetto 3D (l'ellissoide terrestre) su un piano 2D, il che introduce sempre una distorsione (di angoli, aree, distanze, o di tutti e tre). Accettabile a piccola scala, intollerabile a scala continentale.

Un virtual globe elude il problema mantenendo tutto in 3D fino in fondo. Il motore lavora sull'ellissoide reale, non su una proiezione piana. La conseguenza diretta per l'utente: si può zoomare liberamente dall'intero Québec fino a una particella di 50 m, senza cambiare sistema di riferimento, senza cuciture ai confini delle zone MTM, senza stiramenti polari, senza artefatti di allineamento ai margini. È ciò che rende la navigazione continua, qualcosa che nessuno strumento GIS 2D tradizionale può offrire.

Calibrazione del colore e fusione: la sfida visiva

Una volta gestite le proiezioni, la sfida successiva è visiva. Un digital twin 3D del Québec non fonde una, ma tre fonti di immagini complementari:

1. Ortofoto aerea MRNF: priorità 1, 20-30 cm, il livello di dettaglio
2. Immagini satellitari Sentinel-2: priorità 2, 10 m, colma le aree senza copertura aerea (tipicamente le zone costiere)
3. Basemap planetaria Blue Marble: fallback finale, garantisce che nessuna area sia mai vuota

Ciascuna di queste tre fonti ha la propria firma colorimetrica: sensori diversi, condizioni di acquisizione diverse, date di ripresa diverse, elaborazioni a monte diverse. La sfida non è solo mostrarle insieme, ma garantire che i confini restino invisibili all'occhio e che le tessere vicine condividano la stessa gamma cromatica.

Alcune tessere MRNF sono anche incomplete: ai margini delle campagne di volo, la copertura aerea si interrompe a metà di una tessera, lasciando zone nodata. Senza elaborazione, questi vuoti appaiono come bande nere nel terreno 3D. Questo è un altro motivo per cui la pipeline deve fondere più fonti: colmare queste lacune con Sentinel-2 ovunque l'ortofoto aerea sia insufficiente.

Tessera MRNF con copertura parziale. La pipeline colma automaticamente le lacune con Sentinel-2.

Sul lato Sentinel-2, componiamo usando una mediana multi-data su un intervallo di riflettanza fisso: questo elimina i valori anomali (nuvole, ombre, giornate atipiche) e garantisce che tutte le tessere condividano la stessa calibrazione.

Sul lato Québec, alcune campagne di volo più vecchie (in particolare quella del 2015) mostrano una dominante blu sistematica. L'origine esatta del difetto non è documentata dal fornitore dei dati, ma il rapporto blu/rosso anormalmente alto (B/R ~1,30 contro ~0,88 per le campagne normali) colpisce gran parte delle tessere di quell'anno. La pipeline filtra per stagione (acquisizioni estive per toni di verde coerenti), ma il problema persiste per le aree in cui esistono solo tessere del 2015.

Per queste aree sono possibili diversi approcci: correzione cromatica automatica (histogram matching tra tessere vicine), oppure sostituzione con immagini satellitari (Sentinel-2 come prima opzione, gratuita a 10 m, o fonti commerciali come Airbus Pleiades e Maxar per una risoluzione più fine, ma a un costo significativo). Nessuna di queste opzioni è ideale da sola. La combinazione più probabile è un filtraggio per annata abbinato a un fallback su Sentinel-2 per le aree senza alternativa. Sul lato positivo, il MRNF continua a pubblicare regolarmente nuove campagne aeree (le ortofoto 2023 e 2024 sono già disponibili come open data per alcune regioni). Con il tempo, le aree coperte solo da tessere difettose del 2015 saranno probabilmente sostituite da acquisizioni più recenti e meglio calibrate.

In alto: tessera MRNF 2015 con dominante blu (B/R ~1,30). In basso: 2016 con colori normali (B/R ~0,88).

Tessere del 2015 con dominante blu accanto a tessere calibrate normalmente: un effetto patchwork visibile che rompe l'omogeneità del terreno.

Al confine tra le fonti Québec e Sentinel-2, la transizione geometrica è impercettibile. La transizione colorimetrica, invece, resta visibile in alcune aree costiere. È un compromesso deliberato: una fusione a livello di pixel su centinaia di chilometri di costa rappresenterebbe un costo di elaborazione sproporzionato rispetto al guadagno visivo, specialmente alle scale di navigazione tipiche del digital twin. La pipeline privilegia la copertura e la continuità rispetto alla perfezione estetica ai margini.

In alto: solo Sentinel-2. In basso: composito Québec + Sentinel-2. Zona di fusione visibile lungo la linea di costa.

La pipeline include un passaggio di color matching che attenua queste transizioni. L'esempio qui sotto, dalla zona di Rivière-du-Loup, mostra la differenza: senza color matching, il confine tra l'ortofoto del Québec e Sentinel-2 è chiaramente visibile. Con il color matching, la transizione diventa molto più sottile.

Rivière-du-Loup. In alto: senza color matching. In basso: con color matching.

Penisola della Gaspésie, vista ampia. Zona ciano (a sinistra) = tessere 2015 con dominante blu, terraferma = ortofoto Québec, oceano = Sentinel-2.

A che punto siamo

Al 7 aprile 2026, siamo in linea con i tempi. La famiglia NAD83 è assorbita a monte, e la pipeline di ingestione elabora progressivamente le ~115.000 tessere MRNF. La copertura del territorio è parziale in questa fase: le aree già elaborate (per esempio la penisola della Gaspésie) sono navigabili con l'ortofoto del Québec, mentre le altre aree sono riempite nel frattempo con Sentinel-2. Il compositing a tre livelli è stato validato sulla costa della Gaspésie, e il volume di dati funziona senza problemi su una normale workstation da sviluppatore.

Il motore produce attualmente output a 10 m per pixel. È una scelta deliberata per questa prima iterazione: lavorare prima alla scala globale permette di validare l'intera pipeline (proiezioni, fusione, calibrazione del colore, volume di dati) prima di aumentare la risoluzione zona per zona. I dati sorgente MRNF a 20-30 cm ci sono già. L'affinamento fino a 30 cm avverrà caso per caso, in base alle necessità.

Sul lato altimetrico, il terreno 3D si basa su un modello digitale di elevazione SRTM (Shuttle Radar Topography Mission). Non è la fonte più fine disponibile, ma copre l'intero territorio e permette di validare l'intera catena. Un passaggio a modelli digitali del terreno (DTM) del Québec più precisi è realizzabile per le aree prioritarie.

Ecco il risultato attuale: l'ortofoto del Québec a 10 m adagiata sul rilievo SRTM, senza vegetazione 3D e senza dati vettoriali. Questa è la base grezza del terreno su cui verrà costruito tutto il resto.

Prime vedute 3D del Québec. Ortofoto 10 m + rilievo SRTM. Dall'alto in basso: diffusione atmosferica, illuminazione diretta, senza illuminazione.

Una conclusione importante: la qualità del rendering non può superare la qualità dei dati sorgente. Dove i dati MRNF sono puliti e correttamente calibrati (campagne successive al 2015), il risultato è eccellente. Dove presentano dei difetti (campagne 2015 con dominante blu), la pipeline non può inventare una calibrazione cromatica corretta.

Cosa è già attivo oggi:

• Ortofoto aerea del Québec (MRNF) fusa con Sentinel-2
• Rilievo 3D tramite SRTM
• Pipeline completa validata sulla costa della Gaspésie

In arrivo:

• Risoluzione 20-30 cm in zoom ravvicinato
• Dati vettoriali OSM del Québec
• Alberi 3D su tutto il territorio del Québec