Transcript Introduction to spatial data analysis (using ArcGIS)

-> construction, analysis, modelling and
mapping of spatial data in a GIS environment
Introduction to spatial data
analysis (using ArcGIS)
Scuola di Dottorato in Economia, Università di Roma La Sapienza,
2013/2014
Instructors: Filippo Celata, Federico Martellozzo e Luca Salvati
([email protected])
Dati spaziali
- Confini (municipal, county, district,
provinces, urban areas, neighbourhoods,
census tracks, land parcels)
Punti
Linee
Vettoriali
Poligoni
Superfici
Geodati: i layer di base (basemap)
- Rete di trasporto (roads, railroads, trails,
canals, etc.).
Vector
- Singoli siti, “points of interest” (firms,
airports, schools, etc. etc. etc.)
- Idrografia (surface water, rivers, lakes,
canals, etc.)
Raster
- Copertura del suolo (forest, wetlands,
etc,) o uso del suolo (agriculture, urban,
recreational, etc.).
Raster
- Altimetria
Tipologie di geodati
Risoluzione spaziale dei dati (vs. scala di analisi)
Managing geodata(base)
in a GIS
Aggiungere layer all’area di lavoro
Shapefile
Table
Coordinate system
Spatial index
Spatial index
Vettoriali:
Geometry
Add data - c/lezgis14/start/zoneurbanistiche.shp
Layer file
Index of geometry
Add data - c/lezgis14/basemap/…(tutti)
Raster:
File - Add data - c/lezgis14/start/raster.tif
[ File – Save as – c/lezgis14/ - (nome) – Salva ]
Tool bar
L’acquisizione di dati da fonti esterne
Raster
Georiferimento tramite coordinate (latitudine-longitudine, x/y)
Tools - Add XY data: c/lezgis14/AddXY/BBreakfastXY.dbf
- Geocoding services : georiferimento di dati
puntuali tramite informazioni nominali con indirizzo
Formati numerici delle coordinate
Geografiche:
Decimali:
75 59 32.483 W
-7599.2356
ID
HOT EL
INDIRIZZO
CAMERE
LETT I
addetti
MediaDoppia
…
40390353547738
Regina Hotel Baglioni
Via V. Veneto, 72
151
310
108,72
495,0
75 59.32483’ West -75,992356
40390353547732
Marriott Grand Hotel Flora
Via V. Veneto, 191
156
360
57,72
415,0
…
40390353547731
Excelsior
Via V. Veneto, 125
327
616
235,44
812,3
…
75° 59’ 32.483’’ W -75.992356
40390353547723
Ambasciatori Palace
Via V. Veneto, 62
154
366
110,88
425,0
…
40390353526305
De Russie
Via del Babbuino, 9
130
370
93,6
819,0
…
40390353525546
Sheraton Roma
V. le del Pattinaggio, 100
680
1360
251,6
570,5
…
40390353525545
Sheraton Golf
V. le Parco de' Medici, 167
272
543
100,64
388,0
…
40390353525522
Plaza
Via del Corso, 126
200
547
144
598,0
…
40390353525516
Parco Dei Principi
Via G. Frescobaldi, 5
198
364
142,56
580,0
…
…
…
…
…
…
…
…
…
Converti coordinate geografiche in decimali:
DD = (Seconds/3600) + (Minutes/60) + Degrees
[Negativi = - (Seconds/3600) - (Minutes/60) + Degrees]
Per visualizzare esportare i valori delle COORDINATE:
- Calculate geometry / X/Y coordinates of a point (decimali)
- Write feature to text file (crea un file di testo con ID e
coordinate geografiche dei diversi punti di input)
…
-Geocoding services :
http://batchgeo.com/it/
- Apri geocode.xls
- Copia e incolla su
batchgeo.com (solo le prime
righe..)
- Imposta le opzioni e esegui
- Salva mappa(+dati)
- Open in google earth
-Export da google earth (save as .kmz or .kml)
-Da google a GIS (kml to shapefile)
Zonums’ Kml2Shp
online:
-Upload / Process /
Export
-Importa in ArcGIS
su layer già
georiferiti
-File / data / export
data (same coord.
syst. as the
dataframe)
-Da Gis a google (shapefile to .kml)
Zonums’ “shp2kml”
Usa layer in geographic coordinate system WGS1984 (!)
-Gmap editing e collaborazione
Table join: associare tabelle
esterne a layer georiferiti
Editare la
tabella
La tabella degli attributi
(tasto destro sul layer / open
attribute table)
Esportare la tabella
Importare tabelle nel formato corretto (.dbf)
- Table join: associare
gli attributi di una
tabella esterna alla
tabella del layer,
utilizzando una colonna
corrispondente nelle
due tabelle
Esercizio: carta della densità della popolazione nelle
zone urbanistiche di Roma
- Scarica tabella popolazione per zona urbanistica dal sito del Comune di
Roma: link / tasto dx mouse su tavola 3
- Apri il file su excel, sistema (vedi slide 28), salva e chiudi
- Aggiungi il file excel all’area di lavoro di Arcmap
Associazione di dati tabellari (Table join) 1
Associa a zoneurbanistiche.shp la
tabella excel: tasto destro sullo
shapefile / Join and relates / Join / Join
attributes from a table
- Specifica il field per il join nello
shapefile (Zone_urban) e nella tabella
/Esegui / Verifica
Verifica che il join sia riuscito (problemi
e indizi)
-Utilizza la simbologia per creare il
tematismo: tasto destro sul layer zone
urbanistiche / properties / symbology /
quantities / graduated colors, Value:
“POP”, normalization: “areaKM”
(- Sistema il layout e esporta la carta
come immagine a 300dpi)
Fonti di geodati: confini amministrativi
-Rete trasportistica
OPENSTREETMAP (www.openstreetmap.org)
OPENSTREETMAP (www.openstreetmap.org)
2005
Oggi
Openstreetmap editor(s)
OSM geodata
download: OSM,
cloudmade.com,
geofabrik.de
-Land cover e land use
-Corine land cover codes
-Corine land cover codes
Georeferencing e
vettorializzazione
Selezioni e queries
- Selection queries:
- Selezioni manuali
Select by attributes
(=, <, >…)
Sulla tabella
Sulla carta
Select by
location
Selezioni multiple: sulla tabella tieni premuto CTRL / sulla carta
tieni premuto MAIUSCOLO
Select by location:
Intersect
Are within a distance of
Are within
Are completely within
Contain
Completely contain
Have their centroid in
Share a line segment with
Touch the boundary of
Are identical to
Are crossed by the outline of
-Per creare un layer o uno shapefile contenente solo
gli elementi selezionati
1. Create new layer from selected features: tasto destro
sul layer che contiene la selezione / selection / create layer from selected
features
2. Crea uno shapefile
con gli elementi
selezionati tramite:
tasto destro sul layer
creato al punto 1. /
data / export data
(come shapefile)
3. Esporta come layer:
tasto destro / save as
layer file
Selezioni
Fai: Crea selezione contenente solo il bene culturale “colosseo” = seleziona
il bene culturale “colosseo” (sulla carta o sulla tabella) / tasto destro sul
layer / selection / create layer from selected features / modifica nome layer
in “colosseo”
Fai: Crea layer contenente solo la zona urbanistica “centro storico”
Fai: Seleziona tutti i beni culturali classificati come “MUSEI” (Selezione da
tabella / select by attributes: “TIPO” = ‘MUSEI’)
Fai: Seleziona tutti i beni culturali classificati con visitatori paganti (Select
by attributes – “paganti” > 0)
Fai: Seleziona tutti i beni culturali localizzati nella zona urbanistica “centro
storico” (Select / Select by location/ beniculturale che “are completely
within” il source layer “centro storico”)
Fai: Seleziona tutti i beni culturali in un raggio di 5 km dal Colosseo (Select
by location tutti i beniculturale che “are within a distance” da “colosseo” <
5000 mt)