Geostatistica con GIS open source: gli insediamenti neolitici del Tavoliere

Geostatistica con GIS open sourceGli insediamenti neolitici del Tavoliere

Tesi di laurea triennale in Geofisica applicata

Francesco de Virgilio

Relatore: Ch.mo Prof. Marcello Ciminale

Correlatrice: dott.ssa Mariangela Noviello

Dipartimento di Scienze della Terra e GeoambientaliUniversità degli Studi di Bari “Aldo Moro”

Geostatisticacon GIS open

source

Francescode Virgilio

Workflow

Compounds

Distinzione

Area, perim.

Accesso

Orientazione

Codice

WebGIS

Statistiche

Conclusioni

Insediamenti neolitici

compounds

source

Workflow

Compounds

Distinzione

Area, perim.

Accesso

Orientazione

Codice

WebGIS

Statistiche

Conclusioni

compounds

source

Workflow

Compounds

Distinzione

Area, perim.

Accesso

Orientazione

Codice

WebGIS

Statistiche

Conclusioni

compounds

source

2 Workflow

Compounds

Distinzione

Area, perim.

Accesso

Orientazione

Codice

WebGIS

Statistiche

Conclusioni

Il workflow dell’analisi geofisica

foto aeree database

SIT vettori dati derivati statistiche

magnetometria statistiche autom.

disegno

source

3 Workflow

Compounds

Distinzione

Area, perim.

Accesso

Orientazione

Codice

WebGIS

Statistiche

Conclusioni

foto aeree database

disegno

source

4 Workflow

Compounds

Distinzione

Area, perim.

Accesso

Orientazione

Codice

WebGIS

Statistiche

Conclusioni

foto aeree database

disegno

source

Workflow

5 Compounds

Distinzione

Area, perim.

Accesso

Orientazione

Codice

WebGIS

Statistiche

Conclusioni

Compounds: difficoltà e stato dell’arte

A. Niekamp — 2013Pietrele, Romania

source

Workflow

6 Compounds

Distinzione

Area, perim.

Accesso

Orientazione

Codice

WebGIS

Statistiche

Conclusioni

A. Laterza — 2013Tavoliere, Puglia

source

Workflow

6 Compounds

Distinzione

Area, perim.

Accesso

Orientazione

Codice

WebGIS

Statistiche

Conclusioni

perimetro

orientazione(accesso)

source

Workflow

6 Compounds

Distinzione

Area, perim.

Accesso

Orientazione

Codice

WebGIS

Statistiche

Conclusioni

perimetro

orientazione(accesso)

source

Workflow

6 Compounds

Distinzione

Area, perim.

Accesso

Orientazione

Codice

WebGIS

Statistiche

Conclusioni

source

Workflow

Compounds

7 Distinzione

Area, perim.

Accesso

Orientazione

Codice

WebGIS

Statistiche

Conclusioni

Distinzione di ditches e compounds: Jenks

Anglisano (Lucera): 47 geometrie

Dati geografici di partenza

source

Workflow

Compounds

8 Distinzione

Area, perim.

Accesso

Orientazione

Codice

WebGIS

Statistiche

Conclusioni

Anglisano (Lucera): 47 geometrie

Classificazione dei perimetri con intervalli naturali di Jenks, k = 5ditches: classi 1, 2 compounds: classi 3, 4, 5

source

Workflow

Compounds

9 Distinzione

Area, perim.

Accesso

Orientazione

Codice

WebGIS

Statistiche

Conclusioni

Anglisano (Lucera): 47 geometrie = 3 ditches + 44 compounds

ditches

compounds

source

Workflow

Compounds

9 Distinzione

Area, perim.

Accesso

Orientazione

Codice

WebGIS

Statistiche

Conclusioni

Anglisano (Lucera): 47 geometrie = 3 ditches + 44 compounds

ditches

compounds

source

Workflow

Compounds

Distinzione

10 Area, perim.

Accesso

Orientazione

Codice

WebGIS

Statistiche

Conclusioni

Derivazione aree e perimetro

Ccentroide

Fpunto piùdistante daC

Bbuffer da F

source

Workflow

Compounds

Distinzione

11 Area, perim.

Accesso

Orientazione

Codice

WebGIS

Statistiche

Conclusioni

� = 1�

xp = xC + r � cos�

yp = yC + r � cos�

source

Workflow

Compounds

Distinzione

12 Area, perim.

Accesso

Orientazione

Codice

WebGIS

Statistiche

Conclusioni

I = fi1(xi1 ; yi1)g

di2C < di1C

source

Workflow

Compounds

Distinzione

13 Area, perim.

Accesso

Orientazione

Codice

WebGIS

Statistiche

Conclusioni

I = fi1(xi1 ; yi1); i2(xi2 ; yi2)g

source

Workflow

Compounds

Distinzione

14 Area, perim.

Accesso

Orientazione

Codice

WebGIS

Statistiche

Conclusioni

I = fi1(xi1 ; yi1); i2(xi2 ; yi2); i3(xi3 ; yi3)g

source

Workflow

Compounds

Distinzione

15 Area, perim.

Accesso

Orientazione

Codice

WebGIS

Statistiche

Conclusioni

I = fi1(xi1 ; yi1); i2(xi2 ; yi2); i3(xi3 ; yi3); i4(xi4 ; yi4); : : :g

source

Workflow

Compounds

Distinzione

16 Area, perim.

Accesso

Orientazione

Codice

WebGIS

Statistiche

Conclusioni

source

Workflow

Compounds

Distinzione

17 Area, perim.

Accesso

Orientazione

Codice

WebGIS

Statistiche

Conclusioni

source

Workflow

Compounds

Distinzione

Area, perim.

18 Accesso

Orientazione

Codice

WebGIS

Statistiche

Conclusioni

Derivazione dell’accesso

Segmento rettilineo più lungo tra quelli che compongono ilperimetro.

d = 0d < dp3;p4

d = dp3;p4

source

Workflow

Compounds

Distinzione

Area, perim.

19 Accesso

Orientazione

Codice

WebGIS

Statistiche

Conclusioni

d = dp3;p4

d > dp6;p7

d = dp3;p4

source

Workflow

Compounds

Distinzione

Area, perim.

20 Accesso

Orientazione

Codice

WebGIS

Statistiche

Conclusioni

d = dp3;p4

d < dp9;p10

d = dp9;p10

source

Workflow

Compounds

Distinzione

Area, perim.

21 Accesso

Orientazione

Codice

WebGIS

Statistiche

Conclusioni

d = dp9;p10

d > dp12;p13

d = dp9;p10

p12p13

source

Workflow

Compounds

Distinzione

Area, perim.

22 Accesso

Orientazione

Codice

WebGIS

Statistiche

Conclusioni

d = dp9;p10

d < dp14;p1

d = dp14;p1

source

Workflow

Compounds

Distinzione

Area, perim.

23 Accesso

Orientazione

Codice

WebGIS

Statistiche

Conclusioni

d = din ;in+1

in = in+1

get in+1

fin (xin ; yin ) j in 2 I and n = 1g

din ;in+1 > d

register d , in

source

Workflow

Compounds

Distinzione

Area, perim.

Accesso

24 Orientazione

Codice

WebGIS

Statistiche

Conclusioni

Calcolo dell’orientazione

(xe ; ye) = ( 12 (x1 + x2) ; 12 (y1 + y2))

source

Workflow

Compounds

Distinzione

Area, perim.

Accesso

25 Orientazione

Codice

WebGIS

Statistiche

Conclusioni

� = 45�

xp = xC + r � cos�

yp = yC + r � cos�

source

Workflow

Compounds

Distinzione

Area, perim.

Accesso

26 Orientazione

Codice

WebGIS

Statistiche

Conclusioni

source

Workflow

Compounds

Distinzione

Area, perim.

Accesso

Orientazione

27 Codice

WebGIS

Statistiche

Conclusioni

Dalla logica al codice

d = din;in+1

in = in+1

get in+1

fin (xin ; yin) j in 2 I and n = 1g

din;in+1 > d

register d , in

1 # i t e r a t e on a l l open compounds2 f o r compound i n cur_shp . he l p e r compoundsa r ea_se t \3 . f i l t e r ( t ype= ’ compound ’ , open=True ) :4 # get s i d e s and r e l a t i v e l e n g t h s as d i c t i o n a r y5 s i d e s = ge t_s i d e_d i c t ( compound , 3857)6 # get l o n g e s t s i d e i n a r e a po l y gon as a L i n e S t r i n g7 a c c e s s_ l i n e s t r = max ( s i d e s , key= s i d e s . ge t )89 # get a c c e s s l e n g h t as p r o j e c t e d v a l u e10 p r o j_ a c c e s s_ l i n e s t r = a c c e s s_ l i n e s t r11 p r o j_ a c c e s s_ l i n e s t r . t r a n s f o rm ( cur_shp . p r o j )1213 # get the c e n t r o i d o f the a c c e s s s i d e14 f e a t u r e_ c e n t r o i d = compound . p o l y . c e n t r o i d1516 # get compound ’ s f a r t h e s t p o i n t from c e n t r o i d17 max_point = Po i n t ( compound . p o l y . c on v e x_hu l l . e x t e n t [ 2 ] ,18 compound . p o l y . c on v e x_hu l l . e x t e n t [ 3 ] ,19 s r i d =3857)20 r a d i u s = max_point . d i s t a n c e ( f e a t u r e_ c e n t r o i d )2122 # draw c a r d i n a l p o i n t s a round the compound e v e r y23 # 45 degree , and r o t a t e them by 12 deg r e e24 # to a l i g n p e r p e n d i c u l a r l y to N25 c a r d i n a l_ p t s = get_round_ver tex (26 45 ,27 r a d i u s ,28 f e a t u r e_ c e n t r o i d . x ,29 f e a t u r e_ c e n t r o i d . y ,30 3857 ,31 12)3233 # c r e a t e " cake s l i c e s " u s i n g c a r d i n a l p o i n t s34 p o l y g o n_ l i s t = [ ]35 f o r i , i t em i n enumerate ( c a r d i n a l_ p t s ) :36 p o i n t s = ( f e a t u r e_ c e n t r o i d . coords ,37 i t em . coords ,38 c a r d i n a l_ p t s [ i � 1 ] . coords ,39 f e a t u r e_ c e n t r o i d . c oo r d s )40 p o l y g o n_ l i s t . append ( Po lygon ( p o i n t s , s r i d =3857))41 s e c t o r s = Mu l t i Po l y gon ( p o l y g o n_ l i s t , s r i d =3857)4243 # get a c c e s s s i d e c e n t r o i d44 a c c e s s_c e n t r o i d = a c c e s s_ l i n e s t r . c e n t r o i d45 a c c e s s_c e n t r o i d . t r a n s f o rm (3857)

source

Workflow

Compounds

Distinzione

Area, perim.

Accesso

Orientazione

Codice

28 WebGIS

Statistiche

Conclusioni

Interfaccia alle funzioni: webGIS

source

Workflow

Compounds

Distinzione

Area, perim.

Accesso

Orientazione

Codice

29 WebGIS

Statistiche

Conclusioni

source

Workflow

Compounds

Distinzione

Area, perim.

Accesso

Orientazione

Codice

30 WebGIS

Statistiche

Conclusioni

source

Workflow

Compounds

Distinzione

Area, perim.

Accesso

Orientazione

Codice

31 WebGIS

Statistiche

Conclusioni

source

Workflow

Compounds

Distinzione

Area, perim.

Accesso

Orientazione

Codice

32 WebGIS

Statistiche

Conclusioni

source

Workflow

Compounds

Distinzione

Area, perim.

Accesso

Orientazione

Codice

33 WebGIS

Statistiche

Conclusioni

source

Workflow

Compounds

Distinzione

Area, perim.

Accesso

Orientazione

Codice

34 WebGIS

Statistiche

Conclusioni

source

Workflow

Compounds

Distinzione

Area, perim.

Accesso

Orientazione

Codice

WebGIS

35 Statistiche

Conclusioni

Affidabilità: numero di ditches e compounds

ang2 ipp

mlr2 mt

plm11 pn

compounds

(n) Lat13

current

ang2 ipp

mlr2 mt

plm11 pn

insediamenti

ditches(n)

source

Workflow

Compounds

Distinzione

Area, perim.

Accesso

Orientazione

Codice

WebGIS

36 Statistiche

Conclusioni

Affidabilità: aree e perimetri

ang2 ipp

mlr2 mt

plm11 pn

�105

current

ang2 ipp

mlr2 mt

plm11 pn

insediamenti

source

Workflow

Compounds

Distinzione

Area, perim.

Accesso

Orientazione

Codice

WebGIS

37 Statistiche

Conclusioni

N NE E SE SSW W NW

orientazione dei compounds

frequenza(n)

current

H0 : �21 = �22 s2 =

1N � 1

N∑i=1

(xi � x)2 Fcal =s22s21

= 1:07

Ha : �21 6= �22 s21 = 438:79 Fcal < Ftab

s22 = 470:84 H0 verificata

source

Workflow

Compounds

Distinzione

Area, perim.

Accesso

Orientazione

Codice

WebGIS

37 Statistiche

Conclusioni

N NE E SE SSW W NW

orientazione dei compounds

frequenza(n)

current

H0 : �21 = �22 s2 =

1N � 1

N∑i=1

(xi � x)2 Fcal =s22s21

= 1:07

Ha : �21 6= �22 s21 = 438:79 Fcal < Ftab

s22 = 470:84 H0 verificata

source

Workflow

Compounds

Distinzione

Area, perim.

Accesso

Orientazione

Codice

WebGIS

Statistiche

38 Conclusioni

Conclusioni

Dati: 11 insediamenti – 155 compounds

Metodo standard Metodo proposto

tempi � 2 mesi � 30min

affidabilità 100% 100%riproducibilità parziale completa

software installazione locale internet (webGIS)costi licenza GIS (e 2850) open source: e 0

source

Workflow

Compounds

Distinzione

Area, perim.

Accesso

Orientazione

Codice

WebGIS

Statistiche

38 Conclusioni

Conclusioni

tempi � 2 mesi � 30minaffidabilità 100% 100%

riproducibilità parziale completasoftware installazione locale internet (webGIS)costi licenza GIS (e 2850) open source: e 0

source

Workflow

Compounds

Distinzione

Area, perim.

Accesso

Orientazione

Codice

WebGIS

Statistiche

38 Conclusioni

Conclusioni

riproducibilità parziale completa

software installazione locale internet (webGIS)costi licenza GIS (e 2850) open source: e 0

source

Workflow

Compounds

Distinzione

Area, perim.

Accesso

Orientazione

Codice

WebGIS

Statistiche

38 Conclusioni

Conclusioni

riproducibilità parziale completasoftware installazione locale internet (webGIS)

costi licenza GIS (e 2850) open source: e 0

source

Workflow

Compounds

Distinzione

Area, perim.

Accesso

Orientazione

Codice

WebGIS

Statistiche

38 Conclusioni

Conclusioni

riproducibilità parziale completasoftware installazione locale internet (webGIS)costi licenza GIS (e 2850) open source: e 0

source

Workflow

Compounds

Distinzione

Area, perim.

Accesso

Orientazione

Codice

WebGIS

Statistiche

39 Conclusioni

Grazie

So Long, and Thanks for All the Fish(D. Adams)

Geostatistica con GIS open source: gli insediamenti neolitici del Tavoliere

Technology

Transcript of Geostatistica con GIS open source: gli insediamenti neolitici del Tavoliere

Il Tavoliere news dicembre 2014

Geostatistica e interpolazione spaziale Practical Guide to Geostatistical Mapping of Environmental Variables Experimental Variogram Variogram Model Modelli di semivariogramma Vieneadaatounmodello

INDICE CAPITOLO 1 - I GIOCHI DI TAVOLIERE: TRA STORIA E ...scacchicecina.altervista.org/assets/tesi-gioco-scacchi-nella... · Si distinguono dai giochi da tavolo, quali le carte,

Il Tavoliere news marzo 2015

Il Tavoliere delle Puglie e il Basso Volturno tra le Due Guerre · 2017. 9. 20. · In copertina: Opera Nazionale per i Combattenti (a cura di), O.N.C. 36 Anni dell’Opera Nazionale

Il Tavoliere News gennaio 2015

Caso di studio: l’applicazione della geostatistica nell ...

La città moderna e ben dotata di strade e ferrovie, si trova al centro del Tavoliere, fra il Gargano e i Monti della Daunia. La provincia è bagnata.

LUNGO IL 40° PARALLELO NORD · nomadi nordici pre-camuni, primi neolitici...) e poterle, così, retrodatare anche di alcuni millenni. (Marisa Grande, Dai simboli universali alla

Archeologia e Calcolatori 28, 2017, 165-177 · siti preistorici: più precisamente, procedendo da N a S, si tratta degli abitati neolitici di Masainas (zona acquedotto), Is Solinas

Insegnare il mito nella scuola media inferioretesi.supsi.ch/911/1/ANZALONE_Lavinia_2012.pdfmiti agrari. Come dichiara Ries (2005) “i millenni neolitici costituiscono l’humus in

Ambito 10/ TAVOLIERE SALENTINO · 2014. 10. 17. · INTERPRETAZIONE IDENTITARIA E STATUTARIA A partire dalle descrizioni tematiche sviluppate nella sezione A, insieme alla rappresentazione

Creazione di un database ed elaborazione geostatistica dei ... · Strumenti e tecniche di analisi di dati ambientali 2. Acquisizione dati e problematiche 3. Implementazione dati:

Catalogo tavoliere dei sapori

Creazione di un database ed elaborazione geostatistica dei dati...5.2 Geostatistica Siti contaminati: Caratterizzazione, Bonifica e Analisi di Rischio Maria Elena Piccione, Nicoletta

LUNGO LE STRADE DEL VINO N. 6 e N. 8 TRA ANTICO E … · la Premurgia, la Murgia (alta e dei Trulli), l’Anfiteatro Tarantino, il Salento delle Murge e il Tavoliere Salentino. La

Itinerario e curiosità attraverso i sapori mantovani · 2015. 12. 25. · I “mantovani” neolitici dovevano essere abili nel procurarsi il cibo: lo dimostrano gli abbondanti ritrovamenti

Dalla Statistica alla Geostatistica - Aracne editrice · Dalla Statistica alla Geostatistica: Giancarlo Ciotoli Dipartimento di Scienze della Terra, Università di Roma “La Sapienza”

L’uso della geostatistica ambientale per la caratterizzazione geochimica dei suoli del sito inquinato di interesse nazionale di Tito, di Pietro Lucia, Achille Palma, Adriano Sofo,

· PRESENTAZIONE Poesia delle radici: un sostantivo, ... «Tu cerchi la Puglia negli ... Già ondeggia il grano nel Tavoliere