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
Instrucțiuni pentru exercițiu
- Transformă
neighborhoodspentru a avea același CRS canyc_tracts, folosindspTransform()cu argumentulCRSsetat laproj4string(nyc_tracts). - Verifică transformarea uitându-te la
head()alcoordinates(neighborhoods). - Verifică dacă seturile de date se suprapun corect reprezentând grafic
neighborhoods, apoinyc_tractscuadd = TRUEșicol = "red", și în finalwatercuadd = TRUEșicol = "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