我是空间数据的新手。我的目标是从openstreetmap获取gpx
个文件并绘制匈牙利边界的多边形。我可以在Járások
下成功绘制边界,但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)
错误是因为我错过了什么,还是仅仅是openstreetmap数据的问题?任何帮助表示赞赏。