Łą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
Instrukcje do ćwiczenia
- Wywołaj
proj4string()na obiekcieneighborhoods, a następnie nanyc_tracts. Sprawdź, czy obie zwrócone wartości różnią się od siebie. - Przyjrzyj się wynikom
head()funkcjicoordinates()dlaneighborhoods, a potem powtórz to dlanyc_tracts. Widzisz problem? Współrzędne x obiektunyc_tractswynoszą około -70, podczas gdy wneighborhoods– około 1 000 000! - Narysuj wykres dla
neighborhoods, a następnie dodajnyc_tractsz argumentamicol = "red"iadd = 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