多边形的行排序

Question

我的问题很简单。有没有一种自动的方式来为你的数据排序，以便它产生“干净”的多边形？我有生成环的函数（特别是 ahull 函数），我想要一种使用这些函数干净地生成多边形的方法。这是一个例子。

x <- c(1:3, 3:1, 1)
y <- c(1,1,1,3,3,2, 1)
xy <- cbind(x,y)

Sr1 <- Polygon(xy)
Srs1 = Polygons(list(Sr1), "s1")
SpP = SpatialPolygons(list(Srs1)) 
plot(SpP) 

z <- runif(7)
xyz <- cbind(x,y,z)
xyz <- xyz[order(z),]

xy <- xyz[,-3]
xy <- rbind(xy, xy[1,])

Sr1 <- Polygon(xy)
Srs1 = Polygons(list(Sr1), "s1")
SpP = SpatialPolygons(list(Srs1)) 
SpP = SpatialPolygons(list(Srs1)) 
plot(SpP)

这是我的真实数据：https://drive.google.com/file/d/0B8QG4cbDqH0UOUlobnlWaDgwOWs/edit?usp=sharing

Answer 1

从某种意义上说，您已经回答了自己的问题。

假设您有一组点，并且您使用alphahull 包中的ahull(...) 来生成凸包，您可以直接从ahull 对象中以正确的顺序提取边界上的点.这是一个例子：

library(sp)
library(alphahull)
set.seed(1)            # for reproducible example
X <- rnorm(100)
Y <- rnorm(100)
plot(X,Y)
XY <- cbind(X,Y)
hull <- ahull(XY,alpha=1)
plot(hull)

# extract the row numbers of the boundary points, in convex order.
indx=hull$arcs[,"end1"]  
points <- XY[indx,]                  # extract the boundary points from XY
points <- rbind(points,points[1,])   # add the closing point
# create the SpatialPolygonsDataFrame
SpP = SpatialPolygons(list(Polygons(list(Polygon(points)),ID="s1"))) 
plot(SpP)
points(XY)

编辑对 OP 提供数据集的回应。

ahull(...) 似乎在没有警告的情况下使用您的数据集失败 - 它不会产生任何凸包。经过一些实验，看起来问题与 x,y 值的大小有关。如果我将所有内容除以 1000，它会起作用。不知道这是怎么回事（也许其他人会提供见解？？）。无论如何，这是代码和结果：

library(sp)
library(alphahull)
df <- read.csv("ahull problem.csv")
hull <- ahull(df[2:3]/1000,alpha=2)
plot(hull)
# extract the row numbers of the boundary points, in convex order.
indx=hull$arcs[,"end1"]  
points <- df[indx,2:3]               # extract the boundary points from df
points <- rbind(points,points[1,])   # add the closing point
# create the SpatialPolygonsDataFrame
SpP = SpatialPolygons(list(Polygons(list(Polygon(points)),ID="s1"))) 
plot(SpP)
points(df[2:3])

还要注意alpha=2。使用此数据集设置 alpha=1 实际上会生成 2 个船体，一个具有 1 个点，一个具有所有其他点。设置 alpha=2 会创建 1 个船体。