Shapefile 투영 변경

Question

NA에서 +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0으로 독일 모양 파일의 투영법을 변경하거나 할당하려하지만 어떻게 든 잘 작동하지 않습니다.

이 재현 예 : Shape 파일 및 기타 파일을 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) 

원하는 출력 :