Modification de la projection du fichier de formes
J'essaie de modifier ou d'assigner la projection d'un fichier de formes Allemagne de NA à +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0, mais ça ne marche pas bien.
Exemple reproductible: Shapefile et d'autres fichiers peuvent être téléchargés ici :
Ce que j'ai essayé est le suivant:
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"))
Donc, le problème est que je ne peux pas tracer mes observations sur la carte. Voir ci-dessous
Les ponits ont été détectés en utilisant la fonction geocode du package ggmap. 
Toute idée sur la façon de modifier la projection du shapefile ou la projection des coordonnées google serait très appréciée!
J'ai trouvé mon erreur. La solution serait:
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)
Et la sortie souhaitée: 