Québec: Digital Twin forestale, 115.000 ortofoto
Pubblicato il
Le fondamenta prima di tutto
Il Québec pubblica oltre 115.000 tile ortofotografiche da 20 a 30 cm per pixel, diverse centinaia di gigabyte, e noi le abbiamo inserite tutte in un unico terreno 3D continuo. Quasi nessuna di queste attività era 3D. Si trattava di unire tile che non condividono la stessa proiezione, colmare gli spazi tra loro e uniformare i colori tra sensori e anni di acquisizione.
Questa è la parte che nessuno mostra, ed è quella che decide se tutto ciò che viene dopo funzionerà. Finché i dati sottostanti non sono continui, la vista 3D, lo zoom dalla provincia alla particella e ogni calcolo territoriale poggiano sul nulla. Questo post è un rapporto sullo stato di avanzamento di questa fase, applicata al Québec.
La pipeline TerraLab: una scena, ogni fonte
TerraLab-Twin è un globo virtuale planetario. Il suo compito è fondere in tempo reale fonti geospaziali eterogenee (ortofoto pubbliche, ortofoto da drone ad alta risoluzione, modelli digitali del terreno, layer vettoriali, immagini satellitari) su un unico terreno 3D navigabile, dalla scala planetaria fino a quella della particella.
Tre proprietà definiscono lo stato dell'arte a cui puntiamo:
- Continuità della scala: dall'intero Québec (10.000 km di diagonale) fino a una particella di 50 m, senza interruzioni, cambi di CRS o giunzioni ai confini amministrativi.
- Fusione dinamica di più fonti: più fonti di qualità e risoluzione variabili coesistono sullo stesso terreno e vengono selezionate al volo in base all'area, 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. Deve solo navigare.
Tutto il lavoro descritto qui sotto serve a rendere concrete queste tre proprietà in un caso specifico: il territorio del Québec.
Tile MRNF fuse da TerraLab WorldMapServer. Terreno continuo, senza giunzioni visibili tra i fogli originali.
Cosa rende unici i dati del Québec
Il Québec offre uno dei dataset pubblici più ricchi del Nord America per questo tipo di lavoro. Il Ministero delle Risorse Naturali e delle Foreste del Québec (MRNF, precedentemente MERN) pubblica oltre 115.000 tile ortofotografiche, per diverse centinaia di gigabyte di dati complessivi, a una risoluzione nativa di 20-30 cm per pixel. È circa trenta volte più fine delle immagini satellitari Sentinel-2 utilizzate come complemento.
Ortofoto MRNF alla risoluzione nativa di 20-30 cm su un settore forestale del Québec.
Questa ricchezza comporta sfide specifiche: una famiglia eterogenea di proiezioni cartografiche, diversi formati di file, immagini acquisite in più anni e stagioni. È esattamente il tipo di dataset per cui è progettato un globo virtuale, mentre uno strumento GIS tradizionale chiederebbe all'utente di scegliere una proiezione, una stagione e un formato, accettandone le conseguenze.
E sappiamo gestire questa scala. Le 115.000 tile, le centinaia di gigabyte, i formati misti e i dati di più annate: l'architettura TerraLab-Twin è progettata per assorbire nativamente questo volume, su una workstation da sviluppatore, senza richiedere infrastrutture speciali.
Tile MRNF originale a 20 cm. Alberi, strade e particelle sono chiaramente visibili prima di qualsiasi elaborazione.
La parte più difficile: le proiezioni cartografiche
Voglio soffermarmi su questo punto, perché per chiunque abbia lavorato nel geospaziale, è 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 tile può arrivare con un diverso sistema di riferimento delle coordinate.
Gli strumenti sono la seconda metà del problema. Ogni datum ha le sue insidie e ogni strumento della catena geospaziale porta con sé il proprio comportamento storico rispetto alla libreria PROJ. In produzione è estremamente comune vedere una pipeline fallire silenziosamente perché un'installazione obsoleta restituisce trasformazioni con errori di dieci metri o più.
È esattamente ciò che un globo virtuale astrae. TerraLab-Twin conosce una sola proiezione finale: WGS84 sull'ellissoide planetario. Tutta la complessità delle proiezioni 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 interruzioni ai confini delle zone.
Una volta risolto questo problema difficile, tutto il resto diventa possibile: fusione di più fonti, navigazione continua, 3D in tempo reale.
Strumento interno di tiling: ogni fonte viene riproiettata in WGS84, quindi composta nella scena 3D.
Perché WGS84 in un globo virtuale e perché è importante
La scelta di WGS84 sull'ellissoide non è casuale. Una proiezione cartografica classica è, per costruzione, un compromesso: appiattisce un oggetto 3D (l'ellissoide terrestre) su un piano 2D, introducendo sempre distorsioni (negli angoli, nelle aree, nelle distanze o in tutte e tre). Accettabile alle piccole scale, intollerabile a scala continentale.
Un globo virtuale aggira 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 è che si può zoomare liberamente dall'intero Québec fino a una particella di 50 m, senza cambiare CRS, senza giunzioni ai confini delle zone MTM, senza allungamenti polari e senza artefatti di allineamento ai bordi. È questo che rende la navigazione continua, qualcosa che nessuno strumento GIS 2D tradizionale può offrire.
Calibrazione e fusione dei colori: 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:
- Ortofoto aerea MRNF: priorità 1, 20-30 cm, il layer di dettaglio
- Immagini satellitari Sentinel-2: priorità 2, 10 m, coprono le aree senza copertura aerea (in genere le regioni costiere)
- Base cartografica planetaria Blue Marble: fallback finale, garantisce che nessuna area resti mai vuota
Ognuna di queste tre fonti ha la propria firma colorimetrica: sensori diversi, condizioni di acquisizione diverse, date di cattura diverse e processi upstream diversi. La sfida non è solo visualizzarle insieme, ma assicurarsi che i confini restino invisibili all'occhio e che le tile adiacenti condividano lo stesso intervallo cromatico.
Alcune tile MRNF sono inoltre incomplete: ai margini delle campagne di volo, la copertura aerea si interrompe al centro di una tile, lasciando zone nodata. Senza elaborazione, questi vuoti appaiono come bande nere nel terreno 3D. È un altro motivo per cui la pipeline deve fondere più fonti: colmare automaticamente questi vuoti con Sentinel-2 dove l'ortofoto aerea non è sufficiente.
Tile MRNF con copertura parziale. La pipeline colma automaticamente i vuoti con Sentinel-2.
Per Sentinel-2, creiamo una composizione usando una mediana su più date su un intervallo di riflettanza fisso: questo rimuove gli outlier (nuvole, ombre, giornate atipiche) e garantisce che tutte le tile condividano la stessa calibrazione.
La dominante blu delle campagne del 2015
In Québec, alcune campagne di volo più datate (in particolare il 2015) presentano una dominante blu sistematica. L'origine esatta del difetto non è documentata dal fornitore dei dati, ma l'anomalo rapporto blu-rosso (B/R ~1.30 contro ~0.88 per le campagne normali) interessa una larga parte delle tile di quell'anno. La pipeline filtra per stagione (acquisizioni estive per ottenere tonalità verdi coerenti), ma il problema persiste nelle aree in cui esistono solo tile del 2015.
Per queste aree sono possibili diversi approcci: correzione automatica del colore (matching degli istogrammi tra tile adiacenti) 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 costi significativi). Nessuna di queste opzioni è ideale da sola. La combinazione più probabile consiste nel filtrare per annata, associando un fallback Sentinel-2 alle aree prive di alternative.
Il tempo gioca a nostro favore. Il MRNF continua a pubblicare regolarmente nuove campagne aeree e le ortofoto del 2023 e 2024 sono già disponibili come dati aperti per alcune regioni. Le aree coperte solo da tile difettose del 2015 saranno probabilmente sostituite da acquisizioni più recenti e meglio calibrate.
In alto: tile MRNF del 2015 con dominante blu (B/R ~1.30). In basso: 2016 con colori normali (B/R ~0.88).
Tile del 2015 con dominante blu accanto a tile normalmente calibrate: un mosaico visibile che interrompe l'omogeneità del terreno.
Dove l'ortofoto del Québec incontra Sentinel-2
Al confine tra le due fonti, la transizione geometrica è impercettibile. La transizione colorimetrica, invece, resta visibile in alcune aree costiere.
Quella giunzione cromatica è un compromesso deliberato. La fusione a livello di pixel lungo centinaia di chilometri di costa costerebbe più elaborazione di quanto valga il guadagno visivo, soprattutto alle scale di navigazione tipiche del Digital Twin. La pipeline privilegia copertura e continuità rispetto alla perfezione estetica ai margini.
In alto: solo Sentinel-2. In basso: composizione Québec + Sentinel-2. Zona di fusione visibile lungo la costa.
La pipeline include un passaggio di matching dei colori che attenua queste transizioni. L'esempio seguente, dall'area di Rivière-du-Loup, mostra la differenza: senza matching dei colori, il confine tra l'ortofoto del Québec e Sentinel-2 è chiaramente visibile. Con il matching dei colori, la transizione diventa molto più sottile.
Rivière-du-Loup. In alto: senza matching dei colori. In basso: con matching dei colori.
Penisola della Gaspésie, vista ampia. Zona ciano (a sinistra) = tile del 2015 con dominante blu, terraferma = ortofoto del Québec, oceano = Sentinel-2.
Stato attuale
Al 7 aprile 2026, siamo in linea con il programma. La famiglia NAD83 viene assorbita a monte e la pipeline di ingestione sta elaborando progressivamente le circa 115.000 tile 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 vengono temporaneamente riempite con Sentinel-2. La composizione a tre layer è stata convalidata sulla costa della Gaspésie e il volume di dati viene gestito agevolmente da una workstation standard da sviluppatore.
Risoluzione attuale
Il motore produce attualmente 10 m per pixel. È una scelta deliberata per questa prima iterazione: lavorare prima alla scala globale ci consente di convalidare l'intera pipeline (proiezioni, fusione, calibrazione dei colori, volume dei dati) prima di aumentare la risoluzione zona per zona. I dati MRNF di origine a 20-30 cm sono disponibili. Il raffinamento fino a 30 cm avverrà caso per caso, in base alle necessità.
Elevazione
Per quanto riguarda l'elevazione, il terreno 3D si basa su un SRTM (Shuttle Radar Topography Mission) Digital Elevation Model. Non è la fonte più precisa disponibile, ma copre l'intero territorio e ci permette di convalidare la catena completa. Per le aree prioritarie è possibile passare a DTM del Québec più precisi.
Ecco il risultato attuale: l'ortofoto del Québec a 10 m drappeggiata sul rilievo SRTM, senza vegetazione 3D e senza dati vettoriali. Questa è la base grezza del terreno su cui verrà costruito tutto il resto.
Prime viste 3D del Québec. Ortofoto a 10 m + rilievo SRTM. Dall'alto verso il basso: scattering atmosferico, illuminazione diretta, nessuna illuminazione.
Una conclusione importante: la qualità del rendering non può superare quella dei dati di origine. Dove i dati MRNF sono puliti e correttamente calibrati (campagne successive al 2015), il risultato è eccellente. Dove presentano difetti (campagne del 2015 con dominante blu), la pipeline non può inventare una calibrazione corretta dei colori.
Disponibile oggi:
- Ortofoto aerea del Québec (MRNF) fusa con Sentinel-2
- Rilievo 3D tramite SRTM
- Pipeline completa convalidata sulla costa della Gaspésie
Prossimamente:
- Risoluzione di 20-30 cm allo zoom ravvicinato
- Dati vettoriali OSM del Québec
- Alberi 3D su tutto il territorio del Québec
Hai un progetto geospaziale in mente?
Sviluppiamo strumenti IA per QGIS e gemelli digitali 3D su misura. Parlaci del tuo progetto.
