-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathindex.html
711 lines (645 loc) · 28.6 KB
/
index.html
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
<!DOCTYPE html>
<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8"/>
<meta name="viewport" content="width=1024, user-scalable=no">
<meta name="author" content="Enrico Ferreguti - enricofer@gmail.com">
<meta name="venue" content="FOSS4G Italia 2021">
<meta name="date" content="20/09/2021">
<title>Geoprocessing in SQL</title>
<!-- Required stylesheet -->
<link rel="stylesheet" media="screen" href="core/deck.core.css">
<!-- Extension CSS files go here. Remove or add as needed. -->
<link rel="stylesheet" media="screen" href="extensions/goto/deck.goto.css">
<link rel="stylesheet" media="screen" href="extensions/menu/deck.menu.css">
<link rel="stylesheet" media="screen" href="extensions/navigation/deck.navigation.css">
<link rel="stylesheet" media="screen" href="extensions/status/deck.status.css">
<link rel="stylesheet" media="screen" href="extensions/scale/deck.scale.css">
<!-- Style theme. More available in /themes/style/ or create your own. -->
<link rel="stylesheet" media="screen" href="themes/style/swiss.css">
<!-- Transition theme. More available in /themes/transition/ or create your own. -->
<link rel="stylesheet" media="screen" href="themes/transition/horizontal-slide.css">
<!-- Basic black and white print styles -->
<link rel="stylesheet" media="print" href="core/print.css">
<link href='https://cdnjs.cloudflare.com/ajax/libs/codemirror/5.26.0/codemirror.min.css' rel='stylesheet' />
<link href='https://cdnjs.cloudflare.com/ajax/libs/codemirror/5.33.0/theme/duotone-light.min.css' rel='stylesheet' />
<link href='https://cdnjs.cloudflare.com/ajax/libs/codemirror/5.33.0/theme/monokai.min.css' rel='stylesheet' />
<link href='https://cdnjs.cloudflare.com/ajax/libs/codemirror/5.33.0/theme/base16-light.min.css' rel='stylesheet' />
<link href='https://cdnjs.cloudflare.com/ajax/libs/codemirror/5.33.0/theme/neat.min.css' rel='stylesheet' />
<style>
.editor {
width: 100%;
margin-bottom: 15px !important;
}
.wrapper {
display:inline-block;
position:relative;
}
.wrapper .hidden {
display:none;
}
.wrapper button{
position:absolute;
top:4px;
right:4px;
z-index: 100;
height: 25px;
width: 25px;
background-image: url(img/play-icon.png);
background-repeat: no-repeat;
background-position: center;
background-size: cover;
}
.CodeMirror {
height: auto;
font-family: Courier, monospace !important;
font-size: 20px !important;
font-weight: 700 !important;
background-color: #eee
}
</style>
<!-- Required Modernizr file -->
<script src="modernizr.custom.js"></script>
</head>
<body>
<div class="deck-container">
<!-- Begin slides. Just make elements with a class of slide. -->
<section class="slide">
<h2>INTRODUZIONE AI DATABASE GEOSPAZIALI</h2>
<h2>https://enricofer.github.io/spatial_sql_workshop/</h2>
<h3>FOSS4G Italia 2021</h3>
<p><a href="https://creativecommons.org/licenses/by-sa/3.0/it/">(C)2021 CC-BY-SA</a></br>
Enrico Ferreguti enricofer@gmail.com</br>
Luca Lanteri lklanteri@gmail.com</br>
Rocco Pispico r.pispico@gmail.com</br></p>
</section>
<section class="slide">
<h1>Esercitazione guidata</h1>
</section>
<section class="slide">
<h2>SINTASSI SQL</h2>
<p>In linguaggio SQL sono possibili parecchi tipi di Query:</p>
<p>query di creazione tabella:</p>
<div class ='wrapper editor'>
<textarea class="code" id='code1' name='code'>
CREATE TABLE tab(id INT PRIMARY KEY NOT NULL, campo1 CHAR(50), campo2 INT);
</textarea>
<button class="hidden"></button>
</div>
<p>query di inserimento dati:</p>
<div class ='wrapper editor'>
<textarea class="code" id='code2' name='code'>
INSERT INTO tab (id, campo1, campo2) VALUES (1, 'valore 1', 120);
</textarea>
<button class="hidden"></button>
</div>
<p>query di cancellazione tabella:</p>
<div class ='wrapper editor'>
<textarea class="code" id='code3' name='code'>
DROP TABLE tab;
</textarea>
<button class="hidden"></button>
</div>
<p>query di aggiornamento:</p>
<div class ='wrapper editor'>
<textarea class="code" id='code4' name='code'>
UPDATE tab SET campo1 = 'valore2' WHERE id = 1;
</textarea>
<button class="hidden"></button>
</div>
</section>
<section class="slide">
<h2>SQL SELECT</h2>
<p>ma quello che ci interessa oggi è la query di ricerca</p>
<div class ='wrapper editor'>
<textarea class="code" id='code4' name='code'>
SELECT id, campo1 -- prendi id e campo1
FROM tabella -- da tabella
WHERE campo2 > 100; -- dalle righe che soddisfano una condizione</textarea>
<button class="hidden"></button>
</div>
<div class ='wrapper editor'>
<textarea class="code" id='code4' name='code'>
SELECT *
FROM tabella -- elenca tutti i campi e tutte le righe di una tabella</textarea>
<button class="hidden"></button>
</div>
<p><br/>I dialetti SQL sono molti e ci sono grandi differenze di dettaglio.
Ma la struttura della sintassi principale per l'estrazione dei dati
rimane praticamente la stessa per tutti i RDBMS.<br/>
Oggi lavoreremo con SQlite, nella sua versione con estensioni spaziali
<a href="http://www.gaia-gis.it/gaia-sins/">Spatialite</a>
il cui sviluppo è open source ed è curato da
<a href="mailto:a.furieri@lqt.it">Alessandro Furieri</a></p>
</section>
<section class="slide">
<h2>L'ambiente di lavoro del Workshop</h2>
<p>Per facilitare e velocizzare il workshop le esercitazioni sono completamente basate su un una versione speciale di Spatialite compilata per l'utilizzo da browser. Il database è residente nella memoria del browser,
ma è facilmente scaricabile per l'utilizzo con QGIS per chi lo desidera:<a href="data/piem.sqlite">dati esercitazione</a></p>
<p>E' possibile una visualizzazione tabellare:</p>
<div class ='wrapper editor'>
<textarea class="code" id='code' name='code'>
SELECT * FROM comuni LIMIT 10</textarea>
<button></button>
</div>
<p>o una visualizzazione su mappa: Le tabelle contengono un campo *geom* con le informazioni spaziali memorizzate in latitudine e longitudine (WGS84 epsg:4326) e visualizzabili
producendo una geometria <a href="https://it.wikipedia.org/wiki/GeoJSON">geojson</a> per mezzo della funzione AsGeoJSON()<br/></p>
<div class ='wrapper editor'>
<textarea class="code" id='code' name='code'>
SELECT *, AsGeoJSON(geom) AS geojson
FROM comuni</textarea>
<button></button>
</div>
<p>Il formato <a href="https://it.wikipedia.org/wiki/GeoJSON">geojson</a> è un formato testuale facilmente interpretabile da javascript come un'oggetto (dati serializzati) contentente dati geografici e relativi attributi.
Il formato è uno standard codificato di internet e costituisce una specie di lingua universale per l'interscambio di dati geografici nel web</p>
</section>
<section class="slide">
<h2>I dati dell’esercitazione</h2>
<p>
Dati ISTAT (http://www.istat.it/it/archivio/104317)</br>
</br>
Limiti provinciali del Piemonte (province poligonale)</br>
Località abitate (localita poligonale)</br>
Tabella TIPO_LOC Tipologia di località 2001/2011. (tipo_localita.csv)</br>
Il campo può assumere i seguenti valori:</br>
1. centro abitato</br>
2. nucleo abitato</br>
3. località produttiva</br>
4. case sparse</br>
</br>
Da Geoportale Regione Piemonte (http://www.geoportale.piemonte.it/cms/)</br>
Limiti amministrativi comunali (comuni poligonale)</br>
Una selezione della viabilità principale (autostrade lineare)</br>
</p>
</section>
<section class="slide">
<img src="img/dbmanager1.png" alt="DB manager">
</section>
<section class="slide">
<h1>SELECT ... FROM ... WHERE ...</h1>
</section>
<section class="slide">
<h2>query di ricerca SELECT</h2>
<p>Disponiamo delle seguenti tabelle, visualizzabili eseguendo le query:</p>
<div class ='wrapper editor'>
<textarea class="code" id='code' name='code'>
SELECT * FROM comuni -- Confini comunali del Piemonte (geometrie poligonali)</textarea>
<button></button>
</div>
<div class ='wrapper editor'>
<textarea class="code" id='code' name='code'>
SELECT * FROM viab_princ -- Viabilità principale (geometrie lineari)</textarea>
<button></button>
</div>
<div class ='wrapper editor'>
<textarea class="code" id='code' name='code'>
SELECT * FROM localita -- Località ISTAT (geometrie poligonali)</textarea>
<button></button>
</div>
<div class ='wrapper editor'>
<textarea class="code" id='code' name='code'>
SELECT * FROM province -- Province (geometrie poligonali)</textarea>
<button></button>
</div>
<div class ='wrapper editor'>
<textarea class="code" id='code' name='code'>
SELECT * FROM tipo_localita -- Tabella di decodifica tipo di località (senza geometria)</textarea>
<button></button>
</div>
<p>L'elenco delle tabelle e dei relativi campi disponibili può essere visualizzato in ogni momento premendo il bottone "elenco tabelle"</p>
</section>
<section class="slide">
<h2>Esercitazione SELECT</h2>
<div class ='wrapper editor'>
<textarea class="code" id='code' name='code'>
--ES 3
SELECT
t1.descrizione, localita.id,
localita.geom, localita.loc2011,
localita.pro_com, localita.tipo_loc,
localita.popres as residenti, localita.denominazi as nome
FROM
localita,
tipo_localita as t1
WHERE
localita.tipo_loc = t1.tipo_loc
--and popres > 20000
LIMIT 10
</textarea>
<button></button>
</div>
<p>Indicare i campi delle tavole che si desidera selezionare, è possibile indicare anche un ALIAS.
Elenco delle tavole, anche in questo caso è possibile indicare un ALIAS.
Indicare la clausola di JOIN tra le due tabelle.
</p>
<p>Durante la fase di test è consigliato impostare un limite massimo di risultati con la clausola: Limit <numero di valori desiderati>
</br> ad esempio: Limit 10</b>
</p>
</section>
<section class="slide">
<h1>FUNZIONI GEOSPAZIALI</h1>
</section>
<section class="slide">
<h2>Funzioni geometriche di Spatialite</h2>
<p>L'estensione spaziale di SQlite è Spatialite, e dispone di una serie molto completa di funzioni per la manipolazione di informazioni geografiche: <a href="http://www.gaia-gis.it/gaia-sins/spatialite-sql-4.2.0.html">Reference List</a></p>
<p>Le funzioni basilari che utilizzeremo nel workshop sono le seguenti</p>
<ul>
<li>condizioni topologiche (condizioni spaziali): dentro, fuori, contenuto da, coincidente etc...</li>
<li>manipolazione geometrica (aggregazione): unione, differenza, intersezione, buffer, generalizzazione</li>
<li>misurazione geometrica: lunghezza, area, distanza</li>
<li>estrazione geometrica: centroidi, linee di unione</li>
<li>trasformazioni: riproiezioni, conversioni di formati geometrici</li>
</ul>
</section>
<section class="slide">
<h2>Il campo geometria</h2>
<p>I <a href="https://www.gaia-gis.it/gaia-sins/spatialite-cookbook/html/wkt-wkb.html">tipi geometrici</a> possono essere i seguenti:</p>
<ul>
<li>geometria puntuale (POINT e MULTIPOINT)</li>
<li>geometria lineare (LINESTRING e MULTILINESTRING)</li>
<li>geometria poligonale (POLYGON e MULTIPOLYGON)</li>
</ul>
<p>L'informazione spaziale è contenuta in almeno un campo geometrico. Nelle tabelle del workshop il campo geografico è per convenzione il campo "geom", ma nulla vieta che possa avere un altro nome.<br/>Il formato di memorizzazione è binario, e quindi non visualizzabile "per gli umani" ma disponiamo di funzioni che permettono di trasformarlo in un formato più amichevole:</p>
<div class ='wrapper editor'>
<textarea class="code" id='code' name='code'>
SELECT
geom,
ST_AsText(geom) AS wkt,
Hex(ST_AsBinary(geom)) AS wkb
FROM localita</textarea>
<button></button>
</div>
</section>
<section class="slide">
<h2>Esercitazione 4</h2>
<div class ='wrapper editor'>
<textarea class="code" id='code' name='code'>
--ES 4
SELECT
l.id,
ST_Centroid(l.geom) as geom, asgeojson(ST_Centroid(l.geom)) as geojson,
l.loc2011 as codice_localita, l.pro_com as codice_comune, l.tipo_loc,
tl.descrizione as desc_tipo_localita, l.denominazi as nome_localita,
l.popres as popolazione_residente, l.maschi as maschi,
l.popres - l.maschi as femmine
FROM
localita l, tipo_localita tl
WHERE
--localita.tipo_loc = t1.tipo_loc;
cast(l.tipo_loc as integer) = cast(tl.tipo_loc as integer)
LIMIT 10;
</textarea>
<button></button>
</div>
In questa query utilizziamo la funzione ST_Centroid per trasformare un elemento poligonale in puntuale.
Esiste anche la funzione ST_Pointonsurface che restituisce un punto interno alla geometria.
</section>
<section class="slide">
<h2>VISUALIZZAZIONE</h2>
<p>campi calcolati ed alias</p>
<div class ='wrapper editor'>
<textarea class="code" id='code' name='code'>
select
id,
geom,
st_astext(geom) as wkt,
asgeojson(geom) as geojson
from
localita</textarea>
<button></button>
</div>
</section>
<section class="slide">
<h2>ESERCITAZIONE 6 - Le subquery</h2>
<div class ='wrapper editor'>
<textarea class="code" id='code' name='code'>
select
c.id, c.comune, l.nome_localita, l.numero_famiglie, l.numero_abitazioni
from
comuni c,
(
SELECT
l.id,
ST_Centroid(l.geom) as geom, asgeojson(ST_Centroid(l.geom)) as geojson,
l.loc2011 as codice_localita, l.pro_com as codice_comune, l.tipo_loc, tl.descrizione as desc_tipo_localita,
l.denominazi as nome_localita,
l.popres as popolazione_residente, l.maschi as maschi, l.popres - l.maschi as femmine,
l.famiglie as numero_famiglie, l.abitazioni as numero_abitazioni, l.edifici as numero_edifici, l.shape_area as superficie_in_mq
FROM
localita l, tipo_localita tl
WHERE
--localita.tipo_loc = t1.tipo_loc;
cast(l.tipo_loc as integer) = cast(tl.tipo_loc as integer)
) as l
where
st_intersects(l.geom,c.geom)
and c.comune like 'Torino';</textarea>
<button></button>
</div>
In questa query facciamo dei calcoli: deriviamo il numero di femmine come sottrazione tra la popolazione totale e i maschi.
</section>
<section class="slide">
<h2>ESERCITAZIONE 7 - test geospaziali</h2>
<div class ='wrapper editor'>
<textarea class="code" id='code' name='code'>
--ES 7
SELECT
c.comune, count(*), AsGeoJSON(c.geom) as geojson
FROM
comuni AS c, localita AS l
WHERE
c.cod_pro = 5 and
st_intersects(st_envelope(c.geom),st_envelope(l.geom)) and
st_intersects(c.geom,st_centroid(l.geom))
GROUP BY c.comune;</textarea>
<button></button>
</div>
<p>In questa query viene calcolato il numero di località presente in ogni comune.
Dato che l'operazione viene fatta sulla combinazione dei dati di due tabelle (record tabella1 x record tabella2),
il numero di dati su cui si va ad operare è enorme, quindi è opportuno restringere il campo di ricerca inserendo delle condizioni in WHERE</p>
<p>In tutti i database la valutazione delle condizioni WHERE si dice che è "lazy" ovvero alla prima condizione false si tralascia la valutazione delle restanti condizioni</p>
<p>Quindi prima restringiamo la ricerca ad un codica provincia, poi valutiamo se le estensioni delle geometrie coincidono con <strong>st_envelope</strong> e solo allora andremo ad effettuare il test sull'intersezione tra geometrie con <strong>st_intersects</strong> che è una funzione molto pesante dal punto di vista computazionale</p>
</section>
<section class="slide">
<h2>Tipi geometrici</h2>
<img src="http://vos.openlinksw.com/owiki/wiki/VOS/VirtGeoSPARQLEnhancementDocs/VirtWktGeomerticObjects.png"/>
</section>
<section class="slide">
<h2>Visualizzazione in mappa</h2>
<p>Il risultato della query può essere visualizzato in mappa quando la tabella di output contiene un campo "geojson" contente le geometrie da visualizzare in formato <a href="http://geojson.org/">geojson:</a></p>
<p>Il formato geojson è l'espressione testuale di un'oggetto javascript prodotta secondo lo standard "json" ed è facilmente "digeribile" dalle applicazioni web</p>
<p>L'altra condizione perchè la geometria possa essere visualizzata correttamente nel visualizzatore è che il dato geografico deve essere espresso in gradi di latitudine e longitudine. Le tabelle predefinite sono nel formato ETRS89-epsg4258, sono quindi compatibili con il visualizzatore </p>
<div class ='wrapper editor'>
<textarea class="code" id='code' name='code'>
SELECT
AsGeoJSON(
ST_GeomFromText('LINESTRING(7 45,8 45,8 44,7 44, 7 45)')
) AS geojson
</textarea>
<button></button>
</div>
</section>
<section class="slide">
<h2>Condizioni topologiche</h2>
<p></p>
<img src="https://www.gaia-gis.it/spatialite-2.1/resource/mbrs.gif">
<p></p>
<div class ='wrapper editor'>
<textarea class="code" id='code' name='code'>
SELECT
comuni.comune,
AsGeoJSON(geom) as geojson
FROM comuni
WHERE ST_Intersects(geom,ST_GeomFromText('POLYGON((7 45,8 45,8 44,7 44, 7 45))'))</textarea>
<button></button>
</div>
</section>
<section class="slide">
<h2>AGGREGAZIONE SPAZIALE</h2>
<p>Le funzioni spaziali possono essere usate per collezionare le geometrie aggregate da una query di raggruppamento</p>
<div class ='wrapper editor'>
<textarea class="code" id='code' name='code'>
SELECT
cod_pro,
COUNT(comune) as num_comuni,
SUM(shape_area)/10000 as area_ha,
AVG(shape_area)/10000 as area_media,
AsGeoJSON(ST_union(geom)) as geojson
FROM comuni
GROUP BY cod_pro</textarea>
<button></button>
</div>
<p>L'intersezione tra geometrie lineari genera geometrie puntuali:</p>
<div class ='wrapper editor'>
<textarea class="code" id='code' name='code'>
SELECT
asgeojson(st_intersection(viab_princ.geom,padanasuperiore.geom)) as geojson,
viab_princ.id,
padanasuperiore.id
FROM viab_princ, padanasuperiore
WHERE st_intersects(viab_princ.geom,padanasuperiore.geom)</textarea>
<button></button>
</div>
</section>
<section class="slide">
<h2>MISURE, TRASFORMAZIONI E BUFFER</h2>
<p>E' possibile acquisire dati dimensionali delle geometrie inserite lunghezze, aree, distanze</p>
<div class ='wrapper editor'>
<textarea class="code" id='code' name='code'>
SELECT comune,St_area(geom) AS superficie, St_length(geom) AS perimetro FROM comuni</textarea>
<button></button>
</div>
<div class ='wrapper editor'>
<textarea class="code" id='code' name='code'>
SELECT id,St_length(geom) AS lunghezza FROM padanasuperiore</textarea>
<button></button>
</div>
<div class ='wrapper editor'>
<textarea class="code" id='code' name='code'>
SELECT id,St_distance(geom,ST_GeomFromText('POINT(7 45)')) AS dist_pd FROM centri</textarea>
<button></button>
</div>
<p>Essendo le geometrie espresse in latitudine e longitudine, I dati geometrici sono espressi in gradi. E' quindi necessario trasformare le coordinate in un sistema di riferimento proiettato per ottenere le corrette dimensioni in metri</p>
<div class ='wrapper editor'>
<textarea class="code" id='code' name='code'>
SELECT id,name,St_length(st_transform(geom,3003)) AS lunghezza FROM viab_princ
limit 10</textarea>
<button></button>
</div>
<p>La funzione di buffer genera un poligono i cui vertici distano ad una distanza costante da una geometria generatrice</p>
<div class ='wrapper editor'>
<textarea class="code" id='code' name='code'>
SELECT asgeojson(st_buffer(viab_princ.geom,0.03)) as geojson FROM viab_princ
limit 10</textarea>
<button></button>
</div>
</section>
<section class="slide">
<h1>ELABORAZIONI</h1>
</section>
<section class="slide">
<h2>Quali sono i centri serviti dalla ferrovia?</h2>
<p>Individuazione di un buffer di 10 km dalla ferrovia</p>
<div class ='wrapper editor'>
<textarea class="code" id='code' name='code'>
SELECT AsGeoJson(st_buffer(st_union(geom), 0.03)) AS geojson
FROM padanasuperiore
</textarea>
<button></button>
</div>
<p>individuazione dei comuni ricadenti nel buffer appena individuato</p>
<div class ='wrapper editor'>
<textarea class="code" id='code' name='code'>
SELECT *,ASgeoJSON(geom) as geojson
FROM comuni
WHERE st_intersects(geom, (
SELECT st_buffer(st_union(geom), 0.03)
FROM padanasuperiore
))
</textarea>
<button></button>
</div>
<p>metodo alternativo meno performante ma che recupera i dati </p>
<div class ='wrapper editor'>
<textarea class="code" id='code' name='code'>
SELECT comuni.comune,st_distance(comuni.geom,padanasuperiore.geom),ASgeoJSON(comuni.geom) as geojson
FROM comuni JOIN padanasuperiore
ON st_intersects(comuni.geom, st_buffer(padanasuperiore.geom, 0.03))
GROUP BY comuni.comune
</textarea>
<button></button>
<p>NB: i dati geografici sono espressi in latitudine e longitudine quindi la distanza di buffer va espressa in gradi (0.01 gradi corrispondono circa a 950 m)</p>
</div>
</section>
<section class="slide">
<h2>Qual'è la popolazione dell'area metropolitana di Torino?</h2>
<p>Individuazione dei comuni limitrofi</p>
<div class ='wrapper editor'>
<textarea class="code" id='code' name='code'>
SELECT comune, asgeojson(geom) as geojson
FROM comuni
WHERE st_touches(geom, (
SELECT geom
FROM comuni
WHERE comune = 'Torino'
))
</textarea>
<button></button>
<p>Usando st_intersects al posto di st_touches si include anche il comune stesso di Torino</p>
</div>
<div class ='wrapper editor'>
<textarea class="code" id='code' name='code'>
SELECT comune, asgeojson(geom) as geojson
FROM comuni
WHERE st_intersects(geom, (
SELECT geom
FROM comuni
WHERE comune = 'Torino'
))
</textarea>
<button></button>
</div>
</section>
<section class="slide">
<h2>Quanto è lungo Corso Francia?</h2>
<p>Individuazione del fiume: Trattandosi di un grafo idrologico, il fiume è individuato da una serie di archi:</p>
<div class ='wrapper editor'>
<textarea class="code" id='code' name='code'>
SELECT name,asgeojson(geom) as geojson
FROM viab_princ
WHERE name = 'Corso Francia'
</textarea>
<button></button>
</div>
<p>Metodo 1: quindi si può raggruppare sommando la lunghezza di ugni arco</p>
<div class ='wrapper editor'>
<textarea class="code" id='code' name='code'>
SELECT
name,
st_length(st_transform(st_union(geom),3003)),
asgeojson(st_union(geom)) as geojson
FROM viab_princ
WHERE name = 'Corso Francia'
GROUP BY name
</textarea>
<button></button>
</div>
<p>Nota bene: Questa query non funzionerebbe in postgis poichè nel caso di raggruppamento Postgresql richiede di specificare una funzione di aggregazione per ogni campo</p>
</section>
<section class="slide">
<h1>QGIS</h1>
</section>
<section class="slide">
<h2>Spatialite in QGIS</h2>
<p>E' possibile replicare in ambiente QGIS le esercitazioni del workshop. In tal caso è necessario: </p>
<ul>
<li>Scaricare in locale il database dell'esercitazione: <a href="https://github.com/enricofer/spatiasql.js/blob/gh-pages/data/veneto.sqlite?raw=true">veneto.sqlite</a></li>
<li>Stabilire una connessione con il database </li>
<li>Aprire DB MANAGER</li>
<li>Incollare le query nella finestra SQL ed aggiungerle in mappa</li>
</ul>
<img src="img/QGIS_spatialite.png"/>
</section>
<section class="slide">
<h2>Virtual layers in QGIS</h2>
<p>Un'altro modo per usare in modo produttivo le query geospaziali in QGIS sono i <a href="https://docs.qgis.org/2.14/it/docs/user_manual/working_with_vector/virtual_layers.html">Virtual Layers:</a></p>
<p>I virtual Layers sono una caratteristica poco nota ma molto potente di QGIS, e permettono di creare relazioni SQL tra i layer correntemente caricati da qualunque datasource. Il virtual layer viene rigenerato in modo automatico e trasparente ogni qualvolta vengono aggiornati i layers collegati al layer virtuale: </p>
<p>La funzionalità è resa possibile replicando i layers in un database Spatialite locale gestito direttamente da QGIS, dunque la sintassi da utilizzare nelle query è quella di Spatialite. Su questa funzionalità è disponibile un breve ma chiaro <a href="https://pigrecoinfinito.wordpress.com/2016/04/06/qgis-virtual-layer/">tutorial di Totò Fiandaca</a></p>
<img src="img/QGIS_virtual_layers.png"/>
</section>
<!-- End slides. -->
<!-- Begin extension snippets. Add or remove as needed. -->
<!-- deck.navigation snippet -->
<div aria-role="navigation">
<a href="#" class="deck-prev-link" title="Previous">←</a>
<a href="#" class="deck-next-link" title="Next">→</a>
</div>
<!-- deck.goto snippet -->
<form action="." method="get" class="goto-form">
<label for="goto-slide">Go to slide:</label>
<input type="text" name="slidenum" id="goto-slide" list="goto-datalist">
<datalist id="goto-datalist"></datalist>
<input type="submit" value="Go">
</form>
<p class="deck-status deck-progress-10">
<span class="deck-status-current"></span> / <span class="deck-status-total"></span> − <span class="var-venue">GFOSS 2021</span> − <span class="var-title">Introduzione ai database geospaziali</span>
</p>
<!-- End extension snippets. -->
</div>
<!-- Required JS files. -->
<script src="jquery.min.js"></script>
<script src="core/deck.core.js"></script>
<!-- Extension JS files. Add or remove as needed. -->
<script src="extensions/menu/deck.menu.js"></script>
<script src="extensions/goto/deck.goto.js"></script>
<script src="extensions/status/deck.status.js"></script>
<script src="extensions/navigation/deck.navigation.js"></script>
<script src="extensions/scale/deck.scale.js"></script>
<script src='https://cdnjs.cloudflare.com/ajax/libs/codemirror/5.26.0/codemirror.min.js'></script>
<script src='https://cdnjs.cloudflare.com/ajax/libs/codemirror/5.33.0/keymap/sublime.min.js'></script>
<script src='https://cdnjs.cloudflare.com/ajax/libs/codemirror/5.33.0/mode/sql/sql.min.js'></script>
<script src='https://cdnjs.cloudflare.com/ajax/libs/codemirror/5.33.0/addon/scroll/simplescrollbars.min.js'></script>
<!-- Initialize the deck. You can put this in an external file if desired. -->
<script>
$(function() {
$.deck('.slide');
});
window.onload = function () {
var editors = {};
var sqlareas = document.getElementsByClassName('code');
for (var elem of sqlareas) {
var editor = CodeMirror.fromTextArea(elem, {
mode: 'text/x-sql',
keyMap: 'sublime',
//theme: 'neat',
viewportMargin: 0,
lineWrapping: true,
scrollbarStyle: "null",
addModeClass: true,
//readOnly: "nocursor"
});
editor.target = elem.getAttribute('target')
editors[elem.id] = editor
var button = elem.nextSibling.nextSibling.nextSibling;
console.log(button)
button.editor = editor;
button.addEventListener('click', function(event) {
var url = 'browse.html?wasm';
console.log(event.target.editor.getValue());
url += '&query='+event.target.editor.getValue();
if (event.target.editor.getValue().toLowerCase().includes('geojson')){
url += '&tab=m'
} else {
url += '&tab=t'
}
url += '&title=INTRODUZIONE AI DATABASE GEOSPAZIALI - WORKSHOP GFOSS 2021'
console.log(encodeURI(url));
var win = window.open(encodeURI(url), '_blank');
win.focus();
});
}
console.log(editors)
}
</script>
</body>
</html>