如何确定单位正方形中Voronoi点细分的边缘像元

时间:2018-09-21 20:29:57

标签: r cell voronoi tessellation

有没有更简单的方法来检测与边缘相交的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

那么如何识别边缘(边界)的细胞?

0 个答案:

没有答案