Îmbinarea datelor din CRS-uri/proiecții diferite
Fiecare obiect spațial are asociat un sistem de referință de coordonate (CRS). De obicei, acesta este setat la importul datelor și este citit direct din fișierele spațiale. Exact așa au obținut neighborhoods și nyc_tracts informațiile despre sistemul lor de coordonate.
Atât pachetul sp, cât și raster oferă funcția proj4string(), care returnează CRS-ul obiectului pe care îl primește ca argument.
A lucra cu date spațiale care folosesc CRS-uri diferite seamănă puțin cu a compara distanțe — unele în mile, altele în kilometri. Măsoară același lucru, dar valorile nu sunt direct comparabile.
Hai să aruncăm o privire asupra celor două obiecte de tip poligon.
Acest exercițiu face parte din cursul
Vizualizarea datelor geospațiale în R
Instrucțiuni pentru exercițiu
- Apelează
proj4string()peneighborhoods, apoi penyc_tracts. Verifică dacă cele două șiruri sunt diferite. - Uită-te la
head()alcoordinates()pentruneighborhoodsși repetă pentrunyc_tracts. Observi problema?nyc_tractsare coordonate x în jurul valorii -70, în timp ceneighborhoodsare valori în jurul lui 1.000.000! - Reprezintă grafic
neighborhoods, apoi adaugănyc_tractscucol = "red"șiadd = TRUE.
Exercițiu interactiv practic
Încearcă acest exercițiu completând acest cod de exemplu.
library(sp)
# proj4string() on nyc_tracts and neighborhoods
# coordinates() on nyc_tracts and neighborhoods
# plot() neighborhoods and nyc_tracts