Как определить краевые ячейки тесселяции Вороного точек в единичном квадрате - PullRequest
0 голосов
/ 21 сентября 2018

Есть ли более простой способ обнаружить клетки Вороного, которые пересекаются с краем?Я увеличиваю границы единичного квадрата на 0,1 и пересчитываю площади ячеек.Затем я проверяю, какие области ячеек увеличились, так что автоматически становится краевойЯ также использовал tripack, но функцию voronoi.findrejectsites, но она требует объекта триангуляции и принимает границу как выпуклую оболочку точек, а не указанную границу, такую ​​как единичный квадрат.

Я делаю тесселяцию Вороного из точек следующим образом.

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)

Затем получаю Участок тесселяции Вороного из 100 баллов

Так какопределить ячейки на краю (границе)?

...