更改 Shapefile 的投影
Changing the Projection of Shapefile
我正在尝试将 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"))
所以,问题是我无法在地图上绘制我的观察结果。见下文
这些点是使用 ggmap
包中的 geocode
函数检测到的。
任何想法如何更改 shapefile 的投影或 google 坐标的投影将不胜感激!
我发现了我的错误。解决方案是:
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)
以及所需的输出:
我正在尝试将 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"))
所以,问题是我无法在地图上绘制我的观察结果。见下文
ggmap
包中的 geocode
函数检测到的。
我发现了我的错误。解决方案是:
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)
以及所需的输出: