Scroll Top

Provando le nostre capacità

Introduzione

Dall’inizio del 2017, abbiamo raccolto tutti i tipi di dati spaziali riguardanti Paititi: immagini satellitari, mappe antiche e moderne, testimonianze scritte e orali, risultati di altre spedizioni, documenti autentici e leggende. È stato un compito difficile, che continuiamo a svolgere ancora oggi. Poiché lavoriamo con informazioni digitali, abbiamo anche bisogno di mappe in formato digitale: non immagini JPEG scansionate, ma in un formato adatto all’analisi geoinformativa. L’area del nostro interesse è una regione disabitata sui pendii orientali delle Ande, e i materiali cartografici relativi a essa sono difficili da reperire. Quando disponibili, sono spesso di piccola scala e poco informativi. Per capire meglio di cosa stiamo parlando, osservate queste due illustrazioni qui sotto:

Black Marble - Suomi NPP night image over Paititi research area
Un’immagine della Terra di notte scattata dal satellite meteorologico statunitense Suomi NPP. L’ovale rosso evidenzia l’area delle nostre ricerche su Paititi. I punti luminosi nella foto corrispondono agli insediamenti, mentre le aree scure indicano regioni disabitate. La parte superiore della foto mostra l’inizio di un bacino amazzonico quasi disabitato.

Insufficient data
La scritta “Datos Insuficientes” è piuttosto comune sulle mappe topografiche della nostra area di interesse, ammesso che si trovino mappe disponibili.

Solo il telerilevamento della Terra (di seguito ERS, Earth Remote Sensing) ci consente di ottenere informazioni affidabili e pertinenti sull’area di studio, mentre i sistemi di informazione geografica (di seguito GIS, Geographic Information System) raccolgono le conoscenze sul territorio in un unico luogo e ci aiutano nella ricerca di una risposta alla domanda principale: “Dove si trova Paititi?

Struttura del GIS “Alla ricerca di Paititi”

Per sfruttare al meglio i dati raccolti, abbiamo sviluppato il nostro GIS, che per bellezza chiameremo “Alla ricerca di Paititi”. Al centro della sua infrastruttura c’è un server privato virtuale che utilizza il sistema operativo Ubuntu. Sul server (che chiamiamo GIS server) abbiamo installato alcuni software di geoinformazione e mappatura open-source:

  • PostgreSQL—un sistema di gestione di database relazionale a oggetti con estensione PostGIS per supportare le funzionalità geografiche
  • GeoServer—software per la pubblicazione di geodati su un server (server di mappe)
Schema della struttura del sistema di informazione geografica "Alla ricerca di Paititi"
Schema della struttura del sistema di informazione geografica “Alla ricerca di Paititi”.

Gli oggetti geografici appartenenti a tipi diversi (ad esempio, “fiumi”, “strade inca”, “villaggi”, ecc.) sono memorizzati in tabelle separate, e le loro coordinate sono registrate in forma binaria non leggibile (per velocizzare l’elaborazione):

Screenshot della tabella in PostgreSQL
Un frammento di una tabella con una rete fluviale nel nostro database. Ogni record individuale corrisponde a un fiume o torrente. Naturalmente, è molto più facile visualizzare i fiumi su una mappa che interpretare coordinate sconosciute in forma binaria.

Una volta alla settimana, il database viene automaticamente sottoposto a backup e inviato a uno storage cloud per prevenire la perdita del nostro lavoro in caso di guasto del server GIS. Entro un giorno, possiamo ripristinare il progetto almeno allo stato di una settimana fa.

Memorizzare tutte le nostre informazioni spaziali in un unico database ci offre una grande opportunità: tutti i membri del team di Paititi Research possono connettersi da remoto ovunque nel mondo. Per fare ciò, è necessario installare il software open source QGIS (funziona su Windows, Mac OS e Linux) sul proprio computer. In QGIS è possibile visualizzare i nostri geodati nella forma consueta di una mappa composta da singoli Layer GIS, oltre a condurre varie analisi. Il modo di accesso è configurato individualmente per ogni utente, determinando cosa può essere visualizzato o modificato (ad esempio, aggiungere nuovi oggetti interpretando immagini satellitari).

Uno screenshot della finestra del GIS “Alla ricerca di Paititi”. La maggior parte è occupata dall’area della mappa, e sulla sinistra ci sono i layer che compongono questa mappa; in alto ci sono vari strumenti per modificare e analizzare i geodati. Tutti i layer sono caricati dal database tramite Internet.

 

Analisi Geoinformativa

Quindi inizia la fase più interessante—l’analisi geoinformativa: la creazione di nuove conoscenze a partire da quelle esistenti. Qui dimostriamo al meglio le nostre capacità. È importante notare il grande contributo che il modello digitale di elevazione (di seguito DEM, Digital Elevation Model) apporta alla nostra ricerca. Permette di estrarre le caratteristiche quantitative della superficie, simulare il movimento dei flussi d’acqua e calcolare l’illuminazione del territorio.

Modellazione della rete idrografica dettagliata

La rete fluviale è la base della nostra ricerca, con la quale abbiamo iniziato il nostro lavoro. I nomi dei fiumi si riferiscono a fonti letterarie; tutte le spedizioni descrivono i loro movimenti in relazione ai fiumi; con l’aiuto dei fiumi georeferenziamo mappe antiche; infine, lungo i fiumi le persone costruiscono e costruivano le loro abitazioni. Per collegare informazioni verbali occasionali a un riferimento di coordinate nel GIS, è fondamentale disporre di una mappa dettagliata della rete fluviale con i nomi effettivi dei corsi d’acqua. E una tale mappa pronta non si trova. Fortunatamente, simulare il movimento dei flussi d’acqua, avendo un DEM, non è un grosso problema. Questo è quello che abbiamo fatto. Naturalmente, dopo aver ricevuto una tale mappa, è necessario fornire nomi ai fiumi e torrenti, compito già molto più difficile: essendo nell’altro emisfero terrestre, non è sempre possibile scoprire come i locali chiamano un piccolo fiume (purtroppo, Internet non è onnipotente). Tuttavia, alla fine siamo riusciti a identificare tutti i fiumi che in qualche modo sono collegati alla ricerca di Paititi.

Rete fluviale costruita con il DEM che copre 44 mila kmq (area di ricerca di Paititi e dintorni)
La rete fluviale costruita con il DEM che copre 44 mila km² (area di ricerca di Paititi e dintorni).

Processo di georeferenziazione di una mappa antica (a destra) usando i meandri ben definiti della rete fluviale (a sinistra).
Il processo di georeferenziazione di una mappa antica (a destra) basandosi sui meandri distinti della rete fluviale (a sinistra).

 

Ottenere le caratteristiche morfometriche del rilievo

Il rilievo determina come scorrono i fiumi, come viene irrigata e illuminata la superficie e come crescono gli alberi. Per chi ha lavorato con mappe topografiche e DEM, dovrebbe essere evidente: il rilievo è la base su cui costruiamo i nostri insediamenti e conduciamo le nostre attività economiche.

Il DEM insieme alla rete idrografica fornisce la base di tutto lo studio. Senza una comprensione del rilievo è impossibile cercare Paititi.
Il DEM, insieme alla rete idrografica, fornisce la base di tutto lo studio. Senza una comprensione del rilievo, è impossibile cercare Paititi.

Pendenza della superficie

La pendenza determina il grado di variazione della superficie nelle direzioni orizzontale e verticale ed è misurata in gradi da 0 a 90. Abbiamo cercato di verificare quali valori di pendenza corrispondono agli insediamenti antichi e moderni nell’area di nostro interesse. I risultati sono interessanti: si scopre che gli insediamenti si trovano su superfici con una pendenza non superiore a 20°. Anche vivendo in alta montagna, le persone trovano superfici abbastanza pianeggianti per vivere. Pertanto, possiamo ridurre significativamente l’area di ricerca, escludendo con sicurezza le aree inabitabili con pendenze superiori a 60°.

Un frammento della mappa della pendenza della superficie. Le aree arancioni e rosse (pendenza superiore a 30°) possono essere escluse dall’analisi, concentrandosi sulle aree “verdi” (pendenza inferiore a 20°). Come si può vedere da questa illustrazione, l’area di studio sarà significativamente ridotta.

Una mappa della pendenza della superficie per le rovine scoperte nella valle del fiume Chunchusmayo. Le rovine si trovano in piccole aree verdi corrispondenti a superfici con una pendenza inferiore a 20°.

Frammento della mappa della pendenza della superficie per le rovine di Vilcabamba. Le rovine si trovano su un grande altopiano pianeggiante con una pendenza della superficie inferiore a 20°.
Un frammento della mappa della pendenza della superficie per le rovine di Vilcabamba. Tutte le strutture in pietra si trovano su un grande altopiano pianeggiante con una pendenza inferiore a 20°.

Esposizione della superficie

Ogni sito Inca conosciuto al centro del loro impero era situato per ricevere il sole del mattino (sia per ragioni religiose che pratiche). Sorge quindi la domanda: è possibile identificare aree del territorio in cui sia possibile incontrare l’alba? La risposta è : possiamo farlo calcolando l’esposizione del rilievo, che determina la direzione del pendio rispetto ai punti cardinali. Questo è un compito banale nei GIS.

Frammento della mappa dell'esposizione della superficie, in cui direzioni diverse sono codificate da colori differenti.
Un frammento della mappa dell’esposizione della superficie, dove le direzioni rispetto ai punti cardinali sono codificate da colori diversi.

Calcolo della radiazione solare incidente

Studiando le immagini spaziali e il DEM, abbiamo notato che alcune aree sono costantemente in ombra a causa delle peculiarità del rilievo. Poi ci è venuta un’altra idea: e se valutassimo la radiazione solare incidente sulla superficie (insolazione) per filtrare tali luoghi sfavorevoli alla vita? Per calcolare l’insolazione abbiamo utilizzato un software chiamato ArcGIS Pro, e il processo ha richiesto diversi giorni, poiché è molto dispendioso in termini di risorse. Il risultato è una mappa tematica raster della distribuzione della radiazione, dove ogni pixel contiene un numero espresso in watt per metro quadrato (W/m2). Ovviamente, i valori assoluti di insolazione non sono importanti per noi. Abbiamo confrontato i livelli di radiazione corrispondenti agli insediamenti moderni e alle rovine scoperte, eliminando alla fine le aree ovviamente sfavorevoli, quasi non illuminate dal sole.

Frammento della mappa rappresentante l'insolazione totale (diretta + diffusa), calcolata per l'intero anno 1500.
Un frammento della mappa che rappresenta la distribuzione della radiazione solare totale (diretta e diffusa) sulla superficie, calcolata per l’intero anno 1500. Il colore rosso corrisponde ai valori massimi, il giallo a quelli medi e il nero ai minimi.

Un frammento della mappa di insolazione per la valle del fiume Chunchusmayo. Le icone arancioni sulla mappa corrispondono agli insediamenti antichi trovati. Questi luoghi si trovano in aree con un’illuminazione superiore alla media (rosso). L’icona grigia corrisponde a un insediamento moderno situato nella zona “gialla”.

Modellazione 3D

L’aggiunta della terza dimensione ai dati spaziali aumenta significativamente la nostra capacità di interpretarli. Avendo un DEM, è possibile ricreare il rilievo in 3D e poi sovrapporre tutti gli altri dati.

Analisi visiva

Modalità di visualizzazione del DEM con immagini satellitari sovrapposte e foschia atmosferica realistica.
Modalità di visualizzazione del DEM con immagini satellitari sovrapposte e foschia atmosferica realistica.

Esempio di utilizzo della visualizzazione 3D per identificare luoghi specifici. A sinistra: fotografia della montagna sacra Apu Catinti. A destra: la stessa montagna trovata in Google Earth.
Un esempio di come l’uso della visualizzazione 3D del terreno ci aiuta a raccogliere le informazioni necessarie. A sinistra: fotografia della montagna sacra Apu Catinti che stavamo cercando di identificare (fonte: Quest for Paititi, sito della spedizione di Gregory Deyermenjian, copyright © 2004 Quest for Paititi). A destra: la stessa montagna trovata in Google Earth e successivamente aggiunta al nostro database.

Esempio di successo nell'uso della visualizzazione 3D per trovare una montagna con cinque vette.
Un altro esempio di successo nell’uso della visualizzazione 3D per cercare una montagna con cinque vette. A sinistra: foto di Carlos Neuenschwander Landa dal suo libro “Paititi en la bruma de la historia”. A destra: modello 3D del terreno.

Analisi computerizzata della visibilità

Studiando il territorio in modalità tridimensionale, abbiamo scoperto un fatto interessante: tutte le rovine nella valle del fiume Chunchusmayo erano visibili l’una all’altra.

Visibilità reciproca tra le rovine nella valle del fiume Chunchusmayo
Visibilità reciproca tra le rovine nella valle del fiume Chunchusmayo.

Successivamente, abbiamo trovato una conferma della nostra teoria. Gli Inca utilizzavano regolarmente metodi di comunicazione semplici in linea diretta di vista: fumo e luce solare riflessa durante il giorno e fuochi di guardia di notte. Armati di questo nuovo indizio per la ricerca di Paititi, abbiamo iniziato a testare la teoria nella pratica, evidenziando le aree che possono essere osservate dalle rovine già conosciute.

Visibilità dalle aree evidenziate in giallo dalla cima della montagna sacra Apu Catinti.
Le aree evidenziate in giallo corrispondono ai territori visibili dalla cima della montagna sacra Apu Catinti.

Creazione di mappe tematiche

Le cosiddette mappe tematiche rappresentano con grande completezza elementi, proprietà o fenomeni che non sono mostrati sulle mappe geografiche generali. Siamo interessati a due proprietà del paesaggio per le quali non si trovano mappe disponibili: la vivibilità del territorio e la sua percorribilità a piedi.

Mappa della percorribilità

Abbiamo bisogno di questa mappa per pianificare una spedizione a terra, per tracciare il percorso nel modo più ottimale (in termini di durata e facilità di movimento). I principali fattori che influenzano la percorribilità sono:

  • densità della vegetazione
  • pendenza della superficie

Per valutare la densità della foresta, abbiamo utilizzato le mappe VCF (Vegetation Continuous Fields) create dall’Università del Maryland basate sui dati dei satelliti Landsat. Per la pendenza abbiamo utilizzato il DEM. Per combinare questi due set di dati, sono stati riclassificati agli stessi valori. La mappa risultante mostra l’indice di percorribilità con un valore da 1 a 5, dove 5 rappresenta aree facilmente percorribili (pianeggianti e senza foresta) e 1 quelle assolutamente impraticabili (pendii ripidi con foresta densa).

Un frammento della mappa della percorribilità del terreno.

Mappa di vivibilità del territorio

Abbiamo creato diverse mappe per fattori individuali che riteniamo influenzino la vivibilità delle aree per le persone. Questi fattori includono mappe già note come la pendenza, l’insolazione e l’esposizione della superficie. Abbiamo anche cercato di valutare l’umidità dell’area di studio attraverso il calcolo del cosiddetto indice di umidità topografica. Tuttavia, dopo numerosi esperimenti, ci siamo concentrati su due fattori principali, dai quali abbiamo creato una mappa di vivibilità del territorio:

  • pendenza della superficie
  • insolazione

Il principio di combinazione di queste due mappe è lo stesso usato per la mappa della percorribilità del terreno. Il risultato è un indice di vivibilità che assume valori da 1 a 5, dove 5 corrisponde ad aree pianeggianti e ben illuminate, adatte all’attività umana.

Frammento della mappa di vivibilità per l'area intorno a Machu Picchu.
A destra—un frammento della nostra mappa di vivibilità per l’area intorno a Machu Picchu. Come si può vedere, il famoso sito Inca si trova nella zona “verde” con un alto indice di vivibilità.

L’indice di vivibilità è codificato sulla nostra mappa in diversi colori come segue:

  1. Rosso—poco adatto
  2. Arancione—non adatto
  3. Giallo—leggermente non adatto
  4. Verde chiaro—adatto
  5. Verde scuro—molto adatto

L’area iniziale della nostra ricerca era di circa 1300 km². Dopo aver mappato la vivibilità del territorio, abbiamo ridotto l’area di ricerca di molte volte, risparmiandoci molte ore di lavoro infruttuoso.

Interpretazione delle immagini satellitari

Questo è l’ultimo stadio dei nostri studi d’ufficio. Riducendo significativamente l’area di ricerca di Paititi attraverso le analisi spaziali precedenti, siamo stati in grado di procedere all’interpretazione dei dati di telerilevamento da vari satelliti e aerei (leggi di più in “Geodati utilizzati da noi”). Ora possiamo studiare in dettaglio il territorio della selva peruviana alta alla ricerca di qualcosa di “sospetto”, utilizzando diversi range spettrali: dal visibile (lunghezza d’onda media di 500 nm) al microonde (lunghezza d’onda di circa 20 cm).

Processo di interpretazione delle immagini in QGIS. Le linee rosse indicano strutture lineari identificate.
Il processo di interpretazione delle immagini in QGIS. Le linee rosse mostrano le strutture lineari identificate (possibili strade).

Abbiamo effettuato l’interpretazione visivamente, utilizzando il miglior strumento per questo: i nostri occhi e il cervello. Abbiamo esperienza nel rilevamento tramite Machine Learning, e abbiamo cercato di automatizzare questo processo, ma ci siamo scontrati con diversi problemi:

  • Non siamo in grado di formalizzare chiaramente ciò che stiamo cercando in una forma comprensibile per il computer. Questo può includere, ad esempio, variazioni nel colore e nella luminosità della copertura forestale, un certo gioco di ombra e luce… Abbiamo un’idea di come applicare tecniche di apprendimento automatico per cercare siti archeologici in aree disabitate, ma richiede uno studio scientifico approfondito e molta esperienza nella scrittura di nuovi algoritmi di apprendimento automatico.
  • I dati di telerilevamento che abbiamo, e che sono adatti per ricerche dettagliate, hanno una risoluzione spettrale insufficiente. Abbiamo cercato di collegare le proprietà spettrali della vegetazione con le caratteristiche della superficie sottostante. Ad esempio, se alcune strutture in pietra sono semi-sepolte sotto la copertura forestale, certamente cambiano il regime di umidità e nutrizione del suolo in un determinato luogo. Questo dovrebbe influire sul contenuto di clorofilla nelle foglie (misurabile dallo spazio), sulla quantità di biomassa, che alla fine porta a un cambiamento nella riflettanza rispetto agli alberi circostanti. Purtroppo, le immagini satellitari pubblicamente disponibili più adatte a risolvere questo problema hanno una risoluzione spaziale bassa (circa 10 m). È necessario ordinare immagini satellitari da satelliti commerciali ad altissima risoluzione (30 cm).

Puoi vedere alcune delle potenziali rovine trovate nella pagina “Risultati attuali”.

Prospettive future

L’acquisizione di dati di telerilevamento con la migliore risoluzione spaziale e spettrale (e tali dati sono disponibili sul mercato commerciale) ci consentirà di accelerare e semplificare significativamente il compito di trovare nuovi siti archeologici. Tuttavia, solo una spedizione sul campo potrà fornire risposte definitive. L’opzione migliore, a cui aspiriamo, è l’acquisizione di un veicolo aereo senza pilota (UAV), dotato di uno scanner laser (lidar), che fornirà un’immagine dettagliata della superficie sotto la copertura di una fitta foresta.

Phoenix LiDAR - separa il suolo dalla foresta
I droni dotati di lidar consentono la creazione di cosiddetti cloud di punti, dai quali è possibile ottenere informazioni sulla superficie sotto la foresta. Fonte: screenshot da “Webinar: Popular UAV LiDAR Use Cases”, Phoenix LiDAR, sotto licenza YouTube standard / Modificato dall’originale.

Riferimenti

  1. LEE, Vincent R. Forgotten Vilcabamba: Final Stronghold of the Incas, 2000, USA, Empire Publishing, ISBN: 0-9677109-1, p.113
  2. LEE, Vincent R. Forgotten Vilcabamba: Final Stronghold of the Incas, 2000, USA, Empire Publishing, ISBN: 0-9677109-1, p.165