1. Učit se
    2. /
  2. Kurzy
    3. /
  3. Vizualizace geoprostorových dat v R
Connected

cvičení

Převod mezi souřadnicovými referenčními systémy a projekcemi

Převod z jednoho CRS nebo projekce do jiného zajišťují metody spTransform() z balíčku rgdal. spTransform() funguje pro všechny objekty sp, včetně SpatialPolygonsDataFrame, ale nefunguje pro objekty raster. Transformace rastru je totiž o něco složitější – transformovaná obdélníková mřížka totiž po převodu obdélníkový tvar nezachová. Pokud tě transformace rastrů zajímá, podívej se na ?raster::projectRaster.

Samotná transformace je jednoduchá. Prvním argumentem funkce spTransform(), x, je prostorový objekt, který chceš transformovat, a druhým, CRS, je specifikace cílového CRS. CRS lze zadat pomocí PROJ4 řetězce, který bys mohl/a sestavit ručně, ale je mnohem pohodlnější ho převzít z existujícího objektu (například pomocí funkce proj4string()).

Je čas převést oba datové sety polygonů do stejného CRS.

Pokyny

100 XP
  • Transformuj neighborhoods tak, aby měl stejný CRS jako nyc_tracts, pomocí funkce spTransform() s argumentem CRS nastaveným na proj4string(nyc_tracts).
  • Ověř transformaci tím, že se podíváš na výstup head() funkce coordinates(neighborhoods).
  • Zkontroluj, zda jsou teď datové sety správně zarovnány: vykresli neighborhoods, pak nyc_tracts s add = TRUE a col = "red", a nakonec water s add = TRUE a col = "blue".