Filtrowanie dzielnic
Nie potrzebujesz wszystkich dzielnic Nowego Jorku – wystarczy przefiltrować te należące do hrabstwa New York County. Już wiesz, jak to zrobić!
neighborhoods to obiekt SpatialPolygonsDataFrame. W rozdziale 2 nauczyłeś się filtrować tego typu obiekty na podstawie kolumny w slocie data. Kluczem było stworzenie wektora logicznego, a następnie przefiltrowanie SpatialPolygonsDataFrame jak ramki danych.
Jak zidentyfikować właściwe dzielnice? Sprawdź:
head(neighborhoods@data)
Kolumna CountyFIPS zawiera numeryczny kod identyfikujący hrabstwo. Gdy znajdziesz kod dla New York County, zachowasz tylko wiersze z tą wartością.
To ćwiczenie jest częścią kursu
Wizualizacja danych geoprzestrzennych w R
Instrukcje do ćwiczenia
- Obiekt
nyc_tracts_mergerównież zawiera kody hrabstw – w kolumnieCOUNTYFP. Znajdź wartościunique(), aby ustalić kod dla New York County. - Przefiltruj
neighborhoods, dodając warunek logiczny sprawdzający, czyneighborhoods$CountyFIPSma odpowiednią wartość. - Zmodyfikuj wykres tak, aby używał
manhat_hoodszamiastneighborhoods. - Dodaj warstwę
tm_text(), mapująctextna"NTAName".
Interaktywne ćwiczenie praktyczne
Spróbuj tego ćwiczenia, uzupełniając ten przykładowy kod.
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