Sammanfoga data från olika CRS/projektioner
Varje rumsligt objekt har ett associerat koordinatreferenssystem (CRS). Det sätts vanligtvis när data importeras och läses direkt från de rumsliga filerna. Det är på det sättet som neighborhoods och nyc_tracts fick sin koordinatsysteminformation.
Både sp- och raster-paketen har en funktion proj4string() som returnerar CRS för det objekt den anropas på.
Att arbeta med rumsliga data som använder olika CRS liknar att jämföra ett dataset i miles med ett annat i kilometer. De mäter samma sak, men talen går inte att jämföra direkt.
Låt oss titta närmare på våra två polygonobjekt.
Den här övningen är en del av kursen
Visualisering av geospatial data i R
Övningsinstruktioner
- Anropa
proj4string()påneighborhoods, och sedan igen pånyc_tracts. Kontrollera att de två strängarna är olika. - Ta en titt på
head()avcoordinates()förneighborhoodsoch upprepa förnyc_tracts. Ser du problemet?nyc_tractshar x-koordinater runt -70, medanneighborhoodsligger runt 1 000 000! - Plotta
neighborhoods, och plotta sedannyc_tractsmedcol = "red"ochadd = TRUEför att lägga dem ovanpå.
Interaktiv övning med praktiskt arbete
Testa den här övningen genom att slutföra den här exempelkoden.
library(sp)
# proj4string() on nyc_tracts and neighborhoods
# coordinates() on nyc_tracts and neighborhoods
# plot() neighborhoods and nyc_tracts