0

我是 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 行上运行上面的代码时,它产生了一个点垂直和水平分布的图(正是我想要的!)

6 次观察图

但是,相同的代码仍然不适用于我的整个数据集。所以我认为问题可能出在我的 excel 文件中,而不是我的代码中。

4

0 回答 0