Как изменить систему координат (CRS) для файла .shp в R? - PullRequest
Новая
       61

Как изменить систему координат (CRS) для файла .shp в R?

0 голосов
/ 04 мая 2018

У меня есть файл .shp, и я хочу изменить его crs, я пытался использовать spTransform, но он не работает в моем случае. Файл .shp можно найти по адресу https://www.dropbox.com/s/8wfgf8207dbh79r/gpr_000b11a_e.zip?dl=0. 

library(rgdal)
shpfile <- readOGR(dsn="D:/Map",layer = "gpr_000b11a_e")
crs(shpfile)

CRS arguments:
 +proj=longlat +datum=NAD83 +no_defs +ellps=GRS80 +towgs84=0,0,0

spTransform(shpfile, CRS("+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0"))
crs(shpfile)

CRS arguments:
 +proj=longlat +datum=NAD83 +no_defs +ellps=GRS80 +towgs84=0,0,0

Проблема: после spTransform crs для шейп-файла не меняется. Спасибо за любую помощь.

1 Ответ

0 голосов
/ 05 мая 2018

Проблема в том, что вы не приписали преобразованную форму объекту. Попробуйте это: 

shpfile <- spTransform(shpfile, 
CRS("+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0"))

Я рекомендую вам пакет sf, для чтения и обработки .shp файлов, он прост в использовании и эффективен.

Надеюсь, это поможет.

...