具有混乱形状的空间多边形

时间:2017-05-30 11:48:07

标签: r polygon openstreetmap spatial sp

我是空间数据的新手。我的目标是从openstreetmap获取gpx个文件并绘制匈牙利边界的多边形。我可以在Járások下成功绘制边界,但Kistérségek下有一些边界,这些边界混乱,包含点之间的循环。

图像

Plot of kistérségek

特定边界的可重现示例

library(XML)
library(magrittr)
library(sp)

parsed <- xmlParse("http://osmrm.openstreetmap.de/gpx.jsp?relation=1368104") %>% xmlToList()
coord <- do.call(rbind, parsed$rte)
name <- coord[1, 1]
coord <- coord[-(1:2), ]
coord <- apply(coord, 2, as.numeric)

poly <- Polygons(list(Polygon(coord)), name)
sp <- SpatialPolygons(list(poly), proj4string = CRS("+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0"))
plot(sp)

结果

Result

错误是因为我错过了什么,还是仅仅是openstreetmap数据的问题?任何帮助表示赞赏。

0 个答案:

没有答案