Блог

Я пытаюсь изменить или назначить проекцию шейп-файла Germany с NA на proj=longlat datum=WGS84 no_defs ellps=WGS84 towgs84=0,0,0 , но почему-то это работает плохо.

Воспроизводимый пример: шейп-файл и другие файлы можно загрузить здесь:

Я попробовал следующее:

 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, будет высоко оценена!

Я нашел свою ошибку. Решение было бы:

 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)
  

И желаемый результат: