Subset pe cartiere
Nu ai nevoie de toate cartierele din New York, așa că vei selecta doar cartierele din New York County. Știi deja cum să faci asta!
neighborhoods este un SpatialPolygonsDataFrame, iar în Capitolul 2 ai învățat cum să faci subset pe baza unei coloane din slotul data. Cheia era să creezi un vector logic, apoi să aplici subsetul pe SpatialPolygonsDataFrame ca pe un data frame.
Cum poți identifica cartierele potrivite? Verifică:
head(neighborhoods@data)
CountyFIPS este un cod numeric care identifică județul. Dacă afli codul pentru New York County, poți păstra doar rândurile cu acea valoare.
Acest exercițiu face parte din cursul
Vizualizarea datelor geospațiale în R
Instrucțiuni pentru exercițiu
- Obiectul
nyc_tracts_mergeconține și coduri de județ în coloanaCOUNTYFP. Găsește valorileunique()pentru a identifica codul pentru New York County. - Aplică subset pe
neighborhoodsadăugând o condiție logică care verifică dacăneighborhoods$CountyFIPSare valoarea corectă. - Editează graficul pentru a folosi
manhat_hoodsîn loc deneighborhoods. - Adaugă un strat
tm_text(), care mapeazătextla"NTAName".
Exercițiu interactiv practic
Încearcă acest exercițiu completând acest cod de exemplu.
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