【发布时间】:2018-11-18 09:40:58
【问题描述】:
我尝试了几天来找到一种在 R 中使用 sf 包的方法,但没有成功。我想绘制与this example 中的示例 2 类似的区域,但没有任何点。我将我的区域加载为SpatialpolygonsDataFrame,然后我使用设防来获取纬度和长度,如下所示:
area<-readOGR(dsn="/home/ubuntu/..",layer="area")
f_area<-fortify(area)
head(f_area)
long lat order hole piece id group
1 116.1045 57.23717 1 FALSE 1 0 0.1
2 116.1551 57.21548 2 FALSE 1 0 0.1
3 116.2420 57.14505 3 FALSE 1 0 0.1
4 116.1706 57.12011 4 FALSE 1 0 0.1
5 116.1222 57.12006 5 FALSE 1 0 0.1
6 116.0756 57.09926 6 FALSE 1 0 0.1
从这一点来看,我很困惑我必须做什么才能得到我想要的结果。
感谢您的帮助。
【问题讨论】:
-
如果你使用
sf,使用read_sf会更容易 -
您想要使用 sf 的解决方案,还是说您不想使用 sf?
-
@sebdalgarno 我想要一个 sf 的解决方案,因为我从未使用过它。
-
@geo_dd 该解决方案对您有用吗?
标签: r dictionary grid area sf