Daten aus unterschiedlichen CRS/Projektionen zusammenführen
Jedes räumliche Objekt hat ein zugehöriges Koordinatenreferenzsystem (CRS). In der Regel wird es beim Import der Daten gesetzt und direkt aus den räumlichen Dateien ausgelesen. So haben neighborhoods und nyc_tracts ihre Informationen zum Koordinatensystem erhalten.
Sowohl die Pakete sp als auch raster haben eine Funktion proj4string(), die das CRS des Objekts zurückgibt, auf dem sie aufgerufen wird.
Mit räumlichen Daten in unterschiedlichen CRS zu arbeiten, ist ein bisschen so, als würdest du mit einem Datensatz in Meilen und einem anderen in Kilometern arbeiten. Sie messen dasselbe, aber die Zahlen sind nicht direkt vergleichbar.
Schauen wir uns unsere beiden Polygonobjekte an.
Diese Übung ist Teil des Kurses
Geodaten in R visualisieren
Anleitung zur Übung
- Rufe
proj4string()aufneighborhoodsauf und dann erneut aufnyc_tracts. Überprüfe, dass die beiden Strings unterschiedlich sind. - Sieh dir das
head()dercoordinates()vonneighborhoodsan und wiederhole es fürnyc_tracts. Erkennst du das Problem?nyc_tractshat x‑Koordinaten um -70, aberneighborhoodsliegt bei etwa 1.000.000! - Plotte
neighborhoods, dann plotenyc_tractsmitcol = "red"undadd = TRUE, um sie darüberzulegen.
Interaktive Übung
Vervollständige den Beispielcode, um diese Übung erfolgreich abzuschließen.
library(sp)
# proj4string() on nyc_tracts and neighborhoods
# coordinates() on nyc_tracts and neighborhoods
# plot() neighborhoods and nyc_tracts