Ho un certo numero di shapefile in diversi CRS (principalmente WGS84 lat / lon) che vorrei trasformare in una proiezione comune (probabilmente Albers Equal Area Conic, ma potrei chiedere aiuto per scegliere un'altra domanda una volta che il mio problema è migliorato -defined).
Ho trascorso alcuni mesi a fare cose sulle statistiche spaziali in R, ma è stato 5 anni fa. Per la mia vita, non riesco a ricordare come trasformare un sp
oggetto (ad esempio SpatialPolygonsDataFrame
) da una proiezione all'altra.
Codice di esempio:
P4S.latlon <- CRS("+proj=longlat +datum=WGS84")
hrr.shp <- readShapePoly("HRR_Bdry"), verbose=TRUE, proj4string=P4S.latlon)
# Shapefile available at
# http://www.dartmouthatlas.org/downloads/geography/hrr_bdry.zip
# but you must rename all the filenames to have the same
# capitalization for it to work in R
Ora ho un SpatialPolygonsDataFrame
con le informazioni di proiezione appropriate, ma mi piacerebbe trasformarle nella proiezione desiderata. Ricordo che per questo esiste una funzione in qualche modo nominata in modo non intuitivo, ma non ricordo di cosa si tratti.
Nota che non voglio solo cambiare il CRS ma cambiare le coordinate in modo che corrispondano ("riproiettare", "trasformare", ecc.).
modificare
Escludendo AK / HI che sono fastidiosamente posizionati in Messico per questo shapefile:
library(taRifx.geo)
hrr.shp <-
subset(hrr.shp, !(grepl( "AK-" , hrr.shp@data$HRRCITY ) |
grepl( "HI-" , hrr.shp@data$HRRCITY )) )
proj4string(hrr.shp) <- P4S.latlon