web-dev-qa-db-ja.com

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

Germany-ShapefileのプロジェクションをNAから+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0に変更または割り当てようとしていますが、どういうわけかうまくいきません。

再現可能な例:Shapefileおよびその他のファイルをダウンロードできます 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"))

したがって、問題はマップに観測をプロットできないことです。下記参照

enter image description heregeocodeパッケージのggmap関数を使用して、ポニーが検出されました。 enter image description here シェープファイルの投影またはGoogle座標の投影を変更する方法についてのアイデアは大歓迎です!

12
And_R

間違いを見つけました。解決策は次のとおりです。

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)

そして望ましい出力: enter image description here

12
And_R