Sotto-selezionare i neighborhood
Non ti servono tutti quei neighborhood extra a New York, quindi farai una sotto-selezione mantenendo solo quelli nel New York County. Sai già come fare!
neighborhoods è uno SpatialPolygonsDataFrame e nel Capitolo 2 hai imparato a sotto-selezionare in base alla colonna nello slot data. La chiave era creare un vettore logico, poi sotto-selezionare lo SpatialPolygonsDataFrame come un data frame.
Come puoi identificare i neighborhood giusti? Prova con:
head(neighborhoods@data)
CountyFIPS è un codice numerico che identifica la contea. Se riesci a trovare il codice per il New York County, puoi mantenere solo le righe con quel valore.
Questo esercizio fa parte del corso
Visualizzare dati geospaziali in R
Istruzioni dell'esercizio
- L'oggetto
nyc_tracts_mergeha anch'esso i codici di contea nella colonnaCOUNTYFP. Trova i valori conunique()per individuare il codice del New York County. - Sotto-seleziona
neighborhoodsaggiungendo una condizione logica che verifichi seneighborhoods$CountyFIPSha il valore corretto. - Modifica il grafico per usare
manhat_hoodsal posto di neighborhoods. - Aggiungi un layer
tm_text(), mappandotexta"NTAName".
Esercizio pratico interattivo
Prova a risolvere questo esercizio completando il codice di esempio.
library(tmap)
# Find unique() nyc_tracts_merge$COUNTYFP
# Add logical expression to pull out New York County
manhat_hoods <- neighborhoods[___, ]
tm_shape(nyc_tracts_merge) +
tm_fill(col = "estimate") +
tm_shape(water) +
tm_fill(col = "grey90") +
# Edit to use manhat_hoods instead
tm_shape(neighborhoods) +
tm_borders() +
# Add a tm_text() layer