0

我在使用 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()

阿拉斯加河流地图

但是，如果您查看坐标轴，则纬度和经度的格式是我无法识别的。我也没有看到任何明显的转换因素。还有人经历过这个吗？

会话信息：
R 版本 3.4.2 (2017-09-28)
平台：x86_64-w64-mingw32/x64 (64-bit)
运行条件：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

4

1 回答 1

0

第一步之后

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

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

于 2017-10-10T18:46:37.467 回答