Lettura di uno shapefile
Gli shapefile sono uno dei modi più comuni per condividere dati spaziali e si leggono facilmente in R usando readOGR() del pacchetto rgdal. readOGR() ha due argomenti importanti: dsn e layer. Esattamente cosa passare a questi argomenti dipende dal tipo di dati che stai leggendo. Nel video hai visto che, per gli shapefile, dsn deve essere il percorso della cartella che contiene i file che compongono lo shapefile e layer è il nome del file dello specifico shapefile (senza estensione).
Per la tua mappa, ti servono i confini dei quartieri. Abbiamo scaricato le Neighborhood Tabulation Areas, come definite dalla City of New York, dalla Open Data Platform del Department of City Planning. Il download era un archivio zip e abbiamo messo il risultato dell’estrazione del file scaricato nella tua directory di lavoro.
Userai la funzione dir() di base R per esaminare il contenuto della directory di lavoro, quindi leggerai lo shapefile in R.
Questo esercizio fa parte del corso
Visualizzare dati geospaziali in R
Istruzioni dell'esercizio
- Usa
dir()senza argomenti per trovare il nome della directory dello shapefile. - Usa
dir(), passando il percorso alla directory dello shapefile, per vedere i file al suo interno. - Ora conosci la directory e il nome del file. Usa
readOGR()per leggere lo shapefile dei quartieri in un oggetto chiamatoneighborhoods. - Controlla il contenuto chiamando
summary()suneighborhoods. - Controlla il contenuto tracciando il grafico di
neighborhoods.
Esercizio pratico interattivo
Prova a risolvere questo esercizio completando il codice di esempio.
library(sp)
library(rgdal)
# Use dir() to find directory name
# Call dir() with directory name
# Read in shapefile with readOGR(): neighborhoods
# summary() of neighborhoods
# Plot neighborhoods