Question

我在使用broom :: tidy功能时遇到了一些麻烦。这可能是我的shapefile特有的（这是阿拉斯加所有河流的地图）。

我使用readOGR函数加载文件

shp <- readOGR(dsn = ".","mv_navigable_water_ln")

然后使用broom :: tidy函数将空间数据强化为数据框。

shp_points <- tidy(shp, region="group")

这一切似乎都有效。当我使用ggplot绘制点时，它给出了一个空间有意义的地图。

ggplot(data=shape_points, aes(x=long, y=lat, group=group)) +
  geom_path() +
  coord_fixed()

https://github.com/LearnBoost/websocket.io/issues/55

但是如果你看一下轴，那么lat和long的格式是我无法识别的。我也没有看到任何明显的转换因素。还有其他人经历过这个吗？

会话信息：
R版本3.4.2（2017-09-28）
平台：x86_64-w64-mingw32 / x64（64位）
运行于：Windows＆gt; = 8 x64（build 9200）
其他附件包：
map of alaska rivers broom_0.4.2 dplyr_0.7.4 rgeos_0.3-25 ggplot2_2.2.1 rgdal_1.2-13 sp_1.2-5

Answer 1

第一步之后

shp <- readOGR(dsn = ".","mv_navigable_water_ln")

尝试将您的CRS转换为EPSG：4326

library(sp)

shp <- spTransform(shp,"+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0")

然后按照其余步骤进行操作。现在你的轴应该是十进制度的纬度和经度

错了纬度/长度使用强化，扫帚::整洁

1 个答案: