更改Shapefile的投影

Question

我试图更改或指定一个德国-shape文件的投影NA来+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0,但不知何故,它不能很好地工作.

可重复的示例:Shapefile和其他文件可以在这里下载:

我尝试的是以下内容:

library(maptools)
library(sp)
library(rgeos)
library(rgdal)
projection.x <- CRS("+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +towgs84=0,0,0")
mapG <- readShapePoly("vg2500_lan.shp", verbose=TRUE, proj4string=projection.x)
summary(mapG)
mapG <- spTransform(mapG, CRS("+proj=longlat +ellps=WGS84 +datum=WGS84"))

所以,问题是我无法在地图上绘制我的观察结果.见下文

使用包装中的geocode功能检测到ponits ggmap. 任何想法如何改变shapefile的投影或谷歌坐标的投影将非常感激!

Answer 1

我发现了自己的错误.解决方案是:

mapG <- readOGR("vg2500_lan.shp", layer="vg2500_lan")
summary(mapG)
germG <- spTransform(mapG, CRS("+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0"))
plot(germG, axes=T)

和期望的输出: