我是R的新手,无法绘制SpatialPointsDataFrame,最终希望创建最小的凸多边形。我已经尝试了几天但找不到任何解决这个问题的方法。
我将excel数据加载为TXT。文件有3列(纬度,经度,ID)和549行观察。然后我输入以下代码:
# Create coordinates variable
coords <- cbind(LAT = as.numeric(as.character(multi$LAT)), LONG = as.numeric(as.character(multi$LONG)))
# Create the SpatialPointsDataFrame
multi_SPDF <- SpatialPointsDataFrame(coords, data = data.frame(multi$ID), proj4string = CRS("+init=epsg:4326"))
#Plot the points
plot(multi_SPDF)
当我输入此内容时,会生成如下图:
我使用此链接中的类似代码制作了此代码:http://www.alex-singleton.com/R-Tutorial-Materials/7-converting-coordinates.pdf
如果有人能帮助我完成这项工作,我会非常感激。希望我提供了所有必要的信息。
修改
为了提供可重复的示例,我提取了数据的头部以复制到注释中,如下所示:
LAT LONG ID
1 -41.30853 174.7342 7
2 -41.30481 174.7353 6
3 -41.30681 174.7363 7
4 -41.30660 174.7360 10
5 -41.31400 174.7329 10
6 -41.31059 174.7350 6
当我单独在这6行上运行上面的代码时,它产生了一个垂直和水平分布点的图(正是我想要的!)
但是,相同的代码仍无法在我的整个数据集上运行。所以我认为问题可能出在我的excel文件而不是我的代码中。