LoslegenKostenlos loslegen

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

Kurs anzeigen

Anleitung zur Übung

  • Rufe proj4string() auf neighborhoods auf und dann erneut auf nyc_tracts. Überprüfe, dass die beiden Strings unterschiedlich sind.
  • Sieh dir das head() der coordinates() von neighborhoods an und wiederhole es für nyc_tracts. Erkennst du das Problem? nyc_tracts hat x‑Koordinaten um -70, aber neighborhoods liegt bei etwa 1.000.000!
  • Plotte neighborhoods, dann plote nyc_tracts mit col = "red" und add = 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

Code bearbeiten und ausführen