LoslegenKostenlos loslegen

Die Stadtteile subsetten

Du brauchst nicht all die überflüssigen Stadtteile in New York, also subsettest du nur die Stadtteile im New York County. Du weißt schon, wie das geht!

neighborhoods ist ein SpatialPolygonsDataFrame, und in Kapitel 2 hast du gelernt, wie man basierend auf einer Spalte im data-Slot subsettet. Der Trick war, zuerst einen logischen Vektor zu erstellen und dann den SpatialPolygonsDataFrame wie einen Data Frame zu subsetten.

Wie findest du die richtigen Stadtteile? Schau dir an:

head(neighborhoods@data)

CountyFIPS ist ein numerischer Code, der den County identifiziert. Wenn du den Code für New York County herausfindest, kannst du nur die Zeilen mit diesem Wert behalten.

Diese Übung ist Teil des Kurses

Geodaten in R visualisieren

Kurs anzeigen

Anleitung zur Übung

  • Das Objekt nyc_tracts_merge hat ebenfalls County-Codes in der Spalte COUNTYFP. Ermittle die unique()-Werte, um den Code für New York County zu finden.
  • Subsette neighborhoods, indem du eine Logik hinzufügst, die prüft, ob neighborhoods$CountyFIPS den richtigen Wert hat.
  • Bearbeite deinen Plot, sodass manhat_hoods statt neighborhoods verwendet wird.
  • Füge eine tm_text()-Ebene hinzu und mappe text auf "NTAName".

Interaktive Übung

Vervollständige den Beispielcode, um diese Übung erfolgreich abzuschließen.

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
Code bearbeiten und ausführen