Изменение проекции шейп-файла

Я пытаюсь изменить или назначить проекцию шейп-файла Германии с 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, будет высоко оценена!


r sp shapefile projection
person And_R    schedule 11.10.2016    source источник

Ответы (1)


arrow_upward
18
arrow_downward

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

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)

И желаемый результат: введите здесь описание изображения

person And_R    schedule 11.10.2016
comment
Ваша ссылка на данные больше не работает. Не могли бы вы предоставить образец с dput? - person Christoph; 14.06.2019