我有一个从gqig gis软件导出的MultiLineString Geogeson文件。一个小例子:
{
"type": "FeatureCollection",
"name": "route1",
"crs": { "type": "name", "properties": { "name": "urn:ogc:def:crs:EPSG::3857"
} },
"features": [
{ "type": "Feature", "properties": { "FID": 0 }, "geometry": { "type":
"MultiLineString", "coordinates": [ [ [ 1936131.287994222715497,
-4335318.772792792879045 ], [ -2633407.770391199737787,
1763382.609922708477825 ], [ -2922369.195528693497181,
4600947.908943663351238 ], [ -1640888.092745035886765,
5275789.498084637336433 ], [ -361201.781421858817339, 5970373.793290910311043
], [ -361201.781421858817339, 5970373.793290910311043 ] ] ]
} }
]
}
如何在长和纬度的数据帧绑定节点中转换它?预期结果 :
node long lat
1 1936131.287994222715497 -4335318.772792792879045
2 -2633407.770391199737787 1763382.609922708477825
我尝试了什么(创建列表):
route1 <- jsonlite::fromJSON(readr::read_file("routes/route1.geojson"))
library(sf)
可以读取GeoJSON。这会给你一个sf
对象。如果需要坐标,可以使用该st_coordinates()
功能。
library(sf)
sf <- sf::st_read( geo, quiet = T )
df <- as.data.frame( sf::st_coordinates( sf ) )
# X Y L1 L2
# 1 1936131.3 -4335319 1 1
# 2 -2633407.8 1763383 1 1
# 3 -2922369.2 4600948 1 1
# 4 -1640888.1 5275789 1 1
# 5 -361201.8 5970374 1 1
# 6 -361201.8 5970374 1 1
这个额外的L1
和L2
列告诉您每个坐标对在MULTILINESTRING中属于哪个线串。
本文收集自互联网,转载请注明来源。
如有侵权,请联系 [email protected] 删除。
我来说两句