有没有更简单的方法来检测与边缘相交的Voronoi细胞?我将单位正方形的边界增加了0.1,并重新计算了单元格的面积。然后,我检查哪些单元格的面积增加了,因此自动成为边缘单元格。我也使用了tripack
但使用了函数voronoi.findrejectsites
,但是它需要一个三角剖分对象,并将边界作为点的凸包,而不是指定的边界(例如单位平方)。
我对这些点进行Voronoi细分,
library(spatstat)
x = runif(100, min = 0, max = 1)
y = runif(100, min = 0, max = 1)
dir = ppp(x=x, y=y, window = square(c(0,1)))
tess = dirichlet(dir)
plot(tess, main = "")
points(dir, pch=19, cex = 0.5)
然后得到 Plot of the Voronoi tessellation of 100 points
那么如何识别边缘(边界)的细胞?