У мене є ряд форм-файлів у різних CRS (в основному WGS84 lat / lon), які я хотів би перетворити на загальну проекцію (ймовірно, Albers Equal Area Conic, але я можу попросити допомоги щодо вибору в іншому питанні, як тільки моя проблема стане кращою -визначений).
Я провів кілька місяців, займаючись просторовою статистикою в R, але це було 5 років тому. За все життя я не можу пригадати, як перетворити sp
об'єкт (наприклад SpatialPolygonsDataFrame
) з однієї проекції на іншу.
Приклад коду:
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
Тепер у мене є SpatialPolygonsDataFrame
відповідна інформація про проекцію, але я хотів би перетворити її на потрібну проекцію. Я пригадую, що для цього була якась неінтуїтивно названа функція, але я не можу пригадати, що це таке.
Зауважте, що я не хочу просто змінювати CRS, а змінювати координати на відповідність ("перепроектувати", "перетворити" тощо).
Редагувати
За винятком AK / HI, які дратують розміщення в Мексиці для цього формату файлів:
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