使用plotKML函数绘制Lat Long坐标的矢量图

时间:2016-11-27 23:47:16

标签: r plot google-earth

我有一个数据框,在22个变量中有3000多个观测值。

第5:22列包含每个观察的唯一lat长坐标,它们没有配对,即lat long lat long lat long lat long ......等等。下面是我的数据帧的样本

> head(sites,5)
          name    cellName      lat       lon Longitude Latitude        X1       X2        X3       X4        X5
1 Parteen Silo CE0001U21A1 52.69373 -8.609117 -8.609117 52.69373 -8.609428 52.71170 -8.604171 52.71145 -8.599072
2 Parteen Silo CE0001U21B1 52.69373 -8.609117 -8.609117 52.69373 -8.594033 52.67827 -8.598788 52.67689 -8.603870
3 Parteen Silo CE0001U21C1 52.69373 -8.609117 -8.609117 52.69373 -8.637082 52.69970 -8.634897 52.70260 -8.631898
4 Parteen Silo CE0001U21A2 52.69373 -8.609117 -8.609117 52.69373 -8.609381 52.70900 -8.604914 52.70879 -8.600579
5 Parteen Silo CE0001U21B2 52.69373 -8.609117 -8.609117 52.69373 -8.596295 52.68059 -8.600337 52.67942 -8.604657
        X6        X7       X8        X9      X10       X11      X12       X13      X14 Longitude.1 Latitude.1
1 52.71064 -8.594290 52.70929 -8.589977 52.70745 -8.586270 52.70518 -8.583286 52.70255   -8.609117   52.69373
2 52.67605 -8.609117 52.67577 -8.614364 52.67605 -8.619446 52.67689 -8.624202 52.67827   -8.609117   52.69373
3 52.70523 -8.628178 52.70749 -8.623855 52.70932 -8.619065 52.71066 -8.613961 52.71146   -8.609117   52.69373
4 52.70810 -8.596515 52.70696 -8.592849 52.70540 -8.589698 52.70347 -8.587161 52.70123   -8.609117   52.69373
5 52.67870 -8.609117 52.67846 -8.613578 52.67870 -8.617897 52.67942 -8.621940 52.68059   -8.609117   52.69373 

可以找到以dput格式输入数据的链接here

为了在Google地球中绘制这些点,我相信您需要将它们投影到使用以下代码完成的坐标系:

coordinates(sites) <- ~lon+lat
proj4string(sites) <- CRS("+proj=longlat +datum=WGS84")

当我的文件只包含1:4列时,我可以使用

轻松绘制网站详细信息
plotKML(sites)

但现在我希望为每次观察绘制第5至22列,但我无法弄明白。

当每排多个纬线长时,投影和投影片是否不起作用?

当空间点数据框中有多个坐标时,有人可以评论plotKML如何用于绘制多边形吗?

0 个答案:

没有答案