从纬度经度点画线

Question

我正在尝试使用 maptools 从一组提取的数据中绘制一系列线条。 lat 和 long 数据位于 long 数据集 72000 点中，另一列包含每个单独行的标签。每个点都重复此标签。

# from the sp vignette:
M2 = cbind(c(28.55841,25.2479,46.24100,25.24810,25.24570,25.24670,25.24840,28.55820,25.2460),      
       c(77.11510,55.3591,6.11329,55.35931,55.36831,55.38071,55.35801,77.11601,55.3684)
       )

L2 = Line(M2)

# I want to create a loop here to add a new line each time there is a break ie when the line label changes.  
# Perhaps there is a better way to organise this step? 

S2 = Lines(list(L2), ID="b")

library(maptools)

kmlLines(S2, kmlfile = "zrivers.kml", name = "R Lines",
     description = "Hello!", col = "blue", visibility = 1, lwd = 1,
     kmlname = "", kmldescription = "")

我收到消息

Error in vector("character", length(spatialLines@lines)) : no slot of name "lines" for this object of class "Lines"

我没有找到创建 Line 对象的简单方法。

将数据推入形状对象然后使用现有代码（如 rivers 示例）是否更容易？

Answer 1

我尝试了创建形状文件的想法，然后使用 kmlLines 函数进行处理。我发现由于我缺乏知识，我只是遇到了新问题。

最后，我编写了一些 R 代码来处理整个数据集并创建所需的 KML 输出。我选择使用 KML 文件中的格式，最终结果非常令人满意！

从纬度经度点画线

Draw line from latitude Longitude points

r

kml

maptools