Unire dati con CRS/proiezioni differenti
Ogni oggetto spaziale ha associato un sistema di riferimento delle coordinate (CRS). In genere viene impostato quando i dati vengono importati e viene letto direttamente dai file spaziali. È così che neighborhoods e nyc_tracts hanno ottenuto le informazioni sul loro sistema di coordinate.
Sia i pacchetti sp che raster hanno una funzione proj4string() che restituisce il CRS dell'oggetto su cui viene chiamata.
Cercare di lavorare con dati spaziali che usano CRS diversi è un po' come lavorare con un insieme di dati in miglia e un altro in chilometri. Misurano la stessa cosa, ma i numeri non sono direttamente confrontabili.
Diamo un'occhiata ai nostri due oggetti poligonali.
Questo esercizio fa parte del corso
Visualizzare dati geospaziali in R
Istruzioni dell'esercizio
- Chiama
proj4string()suneighborhoods, poi di nuovo sunyc_tracts. Verifica che le due stringhe siano diverse. - Guarda il
head()dellecoordinates()dineighborhoodse ripeti pernyc_tracts. Riesci a vedere il problema?nyc_tractsha coordinate x intorno a -70, mentreneighborhoodsè intorno a 1.000.000! - Fai il grafico di
neighborhoods, poi traccianyc_tractsconcol = "red"eadd = TRUEper aggiungerli sopra.
Esercizio pratico interattivo
Prova a risolvere questo esercizio completando il codice di esempio.
library(sp)
# proj4string() on nyc_tracts and neighborhoods
# coordinates() on nyc_tracts and neighborhoods
# plot() neighborhoods and nyc_tracts