使用加强的经度/经度错误，broom::tidy

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()

map of alaska rivers

但是，如果你看一下轴，就会发现经度和经度的格式我不认识。我也看不到它们的任何明显的转换因子。还有谁有过这样的经历？

会话信息：

R版本3.4.2 (2017-09-28)

平台: x86_64-w64-mingw32/x64 (64位)

运行环境: Windows >= 8 x64 (build 9200)

附加的其他软件包：

1 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")

然后按照剩下的步骤操作。现在，您的轴应该是以十进制度表示的纬度和经度