IniziaInizia gratis

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

Visualizza il corso

Istruzioni dell'esercizio

  • Chiama proj4string() su neighborhoods, poi di nuovo su nyc_tracts. Verifica che le due stringhe siano diverse.
  • Guarda il head() delle coordinates() di neighborhoods e ripeti per nyc_tracts. Riesci a vedere il problema? nyc_tracts ha coordinate x intorno a -70, mentre neighborhoods è intorno a 1.000.000!
  • Fai il grafico di neighborhoods, poi traccia nyc_tracts con col = "red" e add = TRUE per 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

Modifica ed esegui il codice