Changing the Projection of Shapefile

前端 未结 1 625
难免孤独
难免孤独 2021-02-13 09:50

I am trying to change or assign the projection of a Germany-Shapefile from NA to +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0, but s

1条回答
  •  遥遥无期
    2021-02-13 10:24

    I found my mistake. The solution would be:

    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)
    

    And the desired output:

    0 讨论(0)
提交回复
热议问题