我正在尝试更改或分配来自 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"))
所以,问题是我无法在 map 上绘制我的观察结果。见下文
使用
geocode
检测到 ponits函数来自 ggmap
包裹。任何想法如何改变形状文件的投影或谷歌坐标的投影将不胜感激!
最佳答案
我发现了我的错误。解决方案是:
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)
以及所需的输出:
关于r - 改变形状文件的投影,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/39979726/