如何更改R中.shp文件的坐标参照系(CRS)?

时间:2018-05-04 20:24:45

标签: r spatial projection rgdal

我有一个.shp文件,我想更改其crs,我尝试使用spTransform,但在我的情况下它不起作用。可以在https://www.dropbox.com/s/8wfgf8207dbh79r/gpr_000b11a_e.zip?dl=0找到.shp文件。

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之后,shapefile的crs不会改变。谢谢你的帮助。

1 个答案:

答案 0 :(得分:1)

问题是您没有将变换后的形状归因于对象。 试试这个:

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

我建议您使用包sf来阅读和处理.shp文件,它易于使用且高效。

希望它有所帮助。