Sloučení dat z různých CRS/projekcí
Každý prostorový objekt má přiřazený souřadnicový referenční systém (CRS). Ten se obvykle nastaví při importu dat a načte se přímo z prostorových souborů. Přesně takto získaly objekty neighborhoods a nyc_tracts informace o svém souřadnicovém systému.
Balíčky sp i raster obsahují funkci proj4string(), která vrátí CRS objektu, na který ji zavoláš.
Pracovat s prostorovými daty v různých CRS je trochu jako porovnávat dataset v mílích s datasetem v kilometrech. Obojí měří totéž, ale čísla přímo porovnat nelze.
Podívejme se na naše dva polygonové objekty.
Toto cvičení je součástí kurzu
Vizualizace geoprostorových dat v R
Pokyny k cvičení
- Zavolej
proj4string()naneighborhoods, a pak znovu nanyc_tracts. Ověř, že se oba řetězce liší. - Prohlédni si
head()hodnotcoordinates()objektuneighborhoodsa totéž zopakuj pronyc_tracts. Vidíš problém? Objektnyc_tractsmá souřadnice x kolem -70, zatímconeighborhoodskolem 1 000 000! - Vykresli
neighborhoods, a pak vykreslinyc_tractsscol = "red"aadd = TRUE, aby se zobrazily přes ně.
Interaktivní cvičení na vyzkoušení si v praxi
Vyzkoušejte si toto cvičení dokončením tohoto ukázkového kódu.
library(sp)
# proj4string() on nyc_tracts and neighborhoods
# coordinates() on nyc_tracts and neighborhoods
# plot() neighborhoods and nyc_tracts