シェイプファイルの投影を変更する

Question

ドイツのシェイプファイルの投影をNAから+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0に変更しようとしましたが、何とかうまく機能しません。

再現例： シェープファイルやその他のファイルはhereをダウンロードすることができます。

私が試したことは以下の通りです：

だから、

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ファンクションggmapパッケージを使用してポニットが検出されました。 シェイプファイルの投影方法やGoogle座標の投影方法を変更する方法があれば、非常に感謝しています。

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) 

所望の出力：