从纬度经度点

时间:2015-08-11 14:34:29

标签: r kml maptools

我正在尝试使用maptools从一组提取的数据中绘制一系列线条。 latlong数据位于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 = "")

我收到了一条消息

  

向量中的错误(“字符”,长度(spatialLines @ lines)):     对于“Lines”类的此对象,没有名称为“lines”的插槽

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

将数据推送到形状对象然后使用像rivers示例这样的现有代码会更容易吗?

1 个答案:

答案 0 :(得分:0)

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

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