R:多边形组合算法

Question

我有一个包含不同多边形的文件，定义为 data.frame用id , lng和 lat .使用以下代码，我可以在 map 上绘制这些多边形。

    library(ggmap)
    library(ggplot2)

    map <- get_googlemap(center = c(lon = 9.26, lat = 47.3), zoom=10)
    xy <- data.frame(id = c(rep(1,4), rep(2,4)), lng = c(9,9.5,9.5,9,9.25,9.7,9.7,9.24), lat= c(47.1,47.1,47.4,47.4,47.2,47.2,47.5,47.5))
    p <- ggmap(map) + geom_polygon(data=xy, aes(x=xy$lng, y=xy$lat,group=xy$id),fill='red',alpha= 0.2)
    print(p)

在绘制多边形之前是否有合并多边形的功能？ 我想绘制一个多边形来覆盖至少一个多边形下的所有内容。因此，我需要在交点处创建新点，凸包不会被切割。

Answer 1

如果你的意思是:

I want to plot one polygon which covers everything that's under AT LEAST one of the polygons.

then chull(凸包)可能就是您想要的。

举个例子(虽然不是 ggplot):

# create data: (x,y)-coords
x<-c(1:9,9:1)
y<-runif(length(x))*10
# plot `em
plot(x,y,type="l")
# these are the indexes of points that are the convex hull
print(chull(x,y))
these.idx<-chull(x,y)
# plot the convex hull using the indexes
lines(x[these.idx],y[these.idx],type="l",col="red",lwd=2)

对于您的 ggplot 场景，如果您的代码像下面这样修改，它会绘制一个多边形(尽管很难知道这是否是您所追求的):

library(ggmap)
library(ggplot2)

map <- get_googlemap(center = c(lon = 9.26, lat = 47.3), zoom=10)
xy <- data.frame(id = c(rep(1,4), rep(2,4))
    , lng = c(9,9.5,9.5,9,9.25,9.7,9.7,9.24)
    , lat= c(47.1,47.1,47.4,47.4,47.2,47.2,47.5,47.5))
ch.idx<-chull(xy$lng,xy$lat)
xy<-xy[ch.idx,]
# removing:: ,group=xy$id
p <- ggmap(map) 
p <- p + geom_polygon(data=xy, aes(x=xy$lng, y=xy$lat),fill='red',alpha= 0.2)
print(p)

(这更像是评论而不是答案，但我没有“代表”来做那件事抱歉)