![](/img/trans.png)
[英]Plotting line shapefile on ggmap using fortify () or broom::tidy() producing a polygon-like output
[英]Wrong lat/long using fortify, broom::tidy
我在使用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位)
運行於:Windows> = 8 x64(內部版本9200)
其他附件包:
1掃帚_0.4.2 dplyr_0.7.4 rgeos_0.3-25 ggplot2_2.2.1。rgdal_1.2-13 sp_1.2-5
第一步之后
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")
然后按照其余步驟進行操作。 現在,您的軸應為十進制的經度和緯度
聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.