Zacznij terazZacznij za darmo

Łączenie danych z różnych układów CRS/rzutowań

Każdy obiekt przestrzenny ma przypisany układ współrzędnych (CRS). Zazwyczaj jest on ustawiany podczas importu danych i odczytywany bezpośrednio z plików przestrzennych. W ten właśnie sposób obiekty neighborhoods i nyc_tracts uzyskały informacje o swoich układach współrzędnych.

Zarówno pakiet sp, jak i raster udostępniają funkcję proj4string(), która zwraca układ CRS przekazanego obiektu.

Próba pracy z danymi przestrzennymi opartymi na różnych układach CRS przypomina porównywanie odległości – jednych podanych w milach, drugich w kilometrach. Obie miary opisują to samo, ale liczby nie są ze sobą wprost porównywalne.

Przyjrzyjmy się naszym dwóm obiektom poligonowym.

To ćwiczenie jest częścią kursu

Wizualizacja danych geoprzestrzennych w R

Zobacz kurs

Instrukcje do ćwiczenia

  • Wywołaj proj4string() na obiekcie neighborhoods, a następnie na nyc_tracts. Sprawdź, czy obie zwrócone wartości różnią się od siebie.
  • Przyjrzyj się wynikom head() funkcji coordinates() dla neighborhoods, a potem powtórz to dla nyc_tracts. Widzisz problem? Współrzędne x obiektu nyc_tracts wynoszą około -70, podczas gdy w neighborhoods – około 1 000 000!
  • Narysuj wykres dla neighborhoods, a następnie dodaj nyc_tracts z argumentami col = "red" i add = TRUE.

Interaktywne ćwiczenie praktyczne

Spróbuj tego ćwiczenia, uzupełniając ten przykładowy kod.

library(sp)

# proj4string() on nyc_tracts and neighborhoods



# coordinates() on nyc_tracts and neighborhoods



# plot() neighborhoods and nyc_tracts

Edytuj i uruchom kod