GeoArchaeo - Analisi Geostatistica per l'Archeologia
GeoArchaeo è un plugin QGIS per l'analisi geostatistica di dati archeologici puntuali:
statistiche, variogramma, kriging con mappa dell'incertezza, pattern recognition e pianificazione
di nuovi campionamenti. Per la teoria e gli esercizi usa il pulsante Tutorial nel pannello.
Come funziona il pannello
In alto scegli una volta sola il layer di punti e il campo numerico: tutte le schede lavorano
su questa selezione. Sotto, un riepilogo mostra quanti punti validi ci sono, il CRS e l'estensione.
Le schede seguono l'ordine di lavoro consigliato; ogni scheda conserva i propri risultati e il grafico.
1. Dati
Cosa fa: statistiche descrittive, istogramma e note di interpretazione (asimmetria, valori anomali)
Perché: capire la distribuzione prima di interpolare (dati molto asimmetrici → valutare una trasformazione)
2. Variogramma
Cosa fa: misura come la differenza tra i valori cresce con la distanza e adatta un modello (sferico, esponenziale, gaussiano, Matérn)
Parametri: distanza massima (pulsante Auto = 1/3 della diagonale), modello, eventuale verifica dell'anisotropia; i parametri manuali sono opzionali
Risultato: Nugget, Sill, Range e qualità del fit, con una lettura in parole
3. Kriging
Cosa fa: interpola i valori su una griglia e produce due raster: la stima e la sua incertezza (deviazione standard)
Prima di partire: la scheda mostra quante celle avrà la griglia (max 500 per lato) e se il variogramma calcolato è utilizzabile
Dopo: cross-validation leave-one-out (RMSE, MAE, ME, R², MSDR) con interpretazione
4. Pattern (ML)
Metodi: K-Means (gruppi, numero scelto o automatico), DBSCAN (gruppi per densità + punti isolati), Isolation Forest (anomalie)
Risultato: layer di punti con etichetta e legenda per gruppo; Random Forest supervisionato è in Processing
5. Campionamento
Cosa fa: propone dove campionare/scavare in futuro, dentro un'area di studio opzionale
Criteri: copertura (maximin), riduzione dell'incertezza (varianza di kriging), griglia stratificata, adattivo
Risultato: layer con i punti proposti (stelle gialle) e tabella delle coordinate
6. Report
Cosa fa: riunisce in un HTML (o PDF) statistiche, parametri, cross-validation e grafici interattivi delle analisi eseguite
Nota: ogni sezione indica se è disponibile o non ancora calcolata
Processing Toolbox
Nel gruppo GeoArchaeo della Toolbox trovi anche: Co-Kriging (regressione sui dati secondari +
kriging dei residui), kriging spazio-temporale, analisi composizionale (CLR/ILR/ALR), batch kriging per dataset
grandi, pattern recognition con Random Forest supervisionato e il design di campionamento. Tutti gli algoritmi
sono utilizzabili nei modelli grafici e da Python.
Tips
- Usa sempre un CRS proiettato in metri (es. UTM): con coordinate geografiche il pannello ti avvisa
- Per aree grandi aumenta la dimensione della cella: il numero di celle è indicato prima di avviare
- Il pulsante "Carica dati di esempio" aggiunge i dataset inclusi nel plugin (cartella
test_layers)
- Senza QtWebEngine il grafico è un'anteprima statica: "Apri grafico interattivo nel browser" apre la versione Plotly
Risoluzione Problemi
Nessun campo nella tendina? Il campo è di tipo testo: convertilo in numero con il Calcolatore di campi
Il kriging è lento? Aumenta la dimensione della cella o riduci il numero di punti
R² basso? I dati hanno poca correlazione spaziale: la mappa è indicativa, guarda il raster dell'incertezza
DBSCAN mette tutto in "isolati"? Prova K-Means