ÎncepețiÎncepe gratuit

Conversia dintr-un CRS/proiecție în altul

Procesul de conversie dintr-un CRS sau o proiecție în alta este gestionat de metodele spTransform() din pachetul rgdal. spTransform() are metode pentru toate obiectele sp, inclusiv SpatialPolygonsDataFrame, însă nu funcționează pe obiecte de tip raster. Acest lucru se datorează faptului că transformarea unui raster este puțin mai complexă: grila rectangulară transformată nu va mai fi rectangulară. Poți consulta ?raster::projectRaster dacă ești curios(oasă) despre transformarea rasterelor.

Transformarea este simplă. Primul argument al funcției spTransform(), x, este obiectul spațial de transformat, iar al doilea, CRS, specifică CRS-ul dorit. CRS poate fi specificat printr-un șir PROJ4, pe care l-ai putea construi manual, dar este mult mai simplu să îl preiei dintr-un obiect existent (de exemplu, cu funcția proj4string()).

A venit momentul să aduci cele două seturi de date cu poligoane în același CRS.

Acest exercițiu face parte din cursul

Vizualizarea datelor geospațiale în R

Vezi cursul

Instrucțiuni pentru exercițiu

  • Transformă neighborhoods pentru a avea același CRS ca nyc_tracts, folosind spTransform() cu argumentul CRS setat la proj4string(nyc_tracts).
  • Verifică transformarea uitându-te la head() al coordinates(neighborhoods).
  • Verifică dacă seturile de date se suprapun corect reprezentând grafic neighborhoods, apoi nyc_tracts cu add = TRUE și col = "red", și în final water cu add = TRUE și col = "blue".

Exercițiu interactiv practic

Încearcă acest exercițiu completând acest cod de exemplu.

library(sp)
library(raster)

# Use spTransform on neighborhoods: neighborhoods



# head() on coordinates() of neighborhoods


# Plot neighborhoods, nyc_tracts and water


Editează și rulează codul