1. Nauka
    2. /
  2. Kursy
    3. /
  3. Wizualizacja danych geoprzestrzennych w R
Connected

ćwiczenie

Konwersja między układami CRS i odwzorowaniami

Konwersję między układami CRS lub odwzorowaniami obsługują metody spTransform() z pakietu rgdal. Funkcja spTransform() działa ze wszystkimi obiektami sp, w tym z SpatialPolygonsDataFrame, ale nie obsługuje obiektów raster. Wynika to z tego, że transformacja rastra jest nieco bardziej złożona – po przekształceniu prostokątna siatka przestaje być prostokątna. Jeśli chcesz dowiedzieć się więcej o transformowaniu rastrów, zajrzyj do ?raster::projectRaster.

Sama transformacja jest prosta. Pierwszy argument spTransform(), x, to obiekt przestrzenny do przekształcenia, a drugi, CRS, określa docelowy układ współrzędnych. Układ CRS można podać jako ciąg PROJ4 – można go skonstruować ręcznie, ale znacznie wygodniej jest pobrać go z istniejącego obiektu (np. za pomocą funkcji proj4string()).

Czas doprowadzić oba zbiory danych wielokątów do wspólnego układu CRS.

Instrukcje

100 XP
  • Przekształć neighborhoods tak, aby miał taki sam układ CRS jak nyc_tracts, używając spTransform() z argumentem CRS ustawionym na proj4string(nyc_tracts).
  • Sprawdź wynik transformacji, przeglądając head() z coordinates(neighborhoods).
  • Zweryfikuj, czy zbiory danych się pokrywają, rysując neighborhoods, następnie nyc_tracts z add = TRUE i col = "red", a na końcu water z add = TRUE i col = "blue".