-> 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)