De buurten subselecteren
Je hebt al die overbodige buurten in New York niet nodig, dus je selecteert alleen de buurten in New York County. Je weet al hoe dat moet!
neighborhoods is een SpatialPolygonsDataFrame en in Hoofdstuk 2 heb je geleerd hoe je kunt subselecteren op basis van de kolom in de data-slot. De sleutel was een logische vector maken en vervolgens het SpatialPolygonsDataFrame subselecteren zoals een data frame.
Hoe identificeer je de juiste buurten? Bekijk:
head(neighborhoods@data)
CountyFIPS is een numerieke code die de county aangeeft. Als je de code voor New York County weet te vinden, kun je alleen de rijen met die waarde behouden.
Deze oefening maakt deel uit van de cursus
Georuimtelijke data visualiseren in R
Oefeninstructies
- Het object
nyc_tracts_mergeheeft ook county-codes in de kolomCOUNTYFP. Zoek deunique()-waarden om de code voor New York County te vinden. - Subselecteer
neighborhoodsdoor een logische test toe te voegen die controleert ofneighborhoods$CountyFIPSde juiste waarde heeft. - Pas je plot aan om
manhat_hoodste gebruiken in plaats van neighborhoods. - Voeg een
tm_text()-laag toe en maptextnaar"NTAName".
Praktische interactieve oefening
Probeer deze oefening eens door deze voorbeeldcode in te vullen.
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