тест на плитку для плитки Дирихле, используя R - PullRequest
0 голосов
/ 03 декабря 2018

Так что я могу взять очки и использовать библиотеки R deldir или spatstat::dirichlet, чтобы найти тесселяцию дирихле этих точек.

Теперь у меня нет точки в наборе, и я хочузнать индексы точек, образующих плитку дирихле, внутренняя часть которой не соответствует заданной точке.Я могу получить это, зная метку плитки (или индекс).

Существуют ли какие-либо библиотеки или методы для этого?Я думаю спатстат, но пока что-то там не нашел.

Ответы [ 2 ]

0 голосов
/ 13 декабря 2018

Если X является точечным шаблоном, а B является тесселяцией, тогда

M <- marks(cut(X, B))

возвращает коэффициент (вектор категориальных значений), идентифицирующий, какой тайл содержит каждую из точек X,В качестве альтернативы

M <- tileindex(X$x, X$y, B)

или

f <- as.function(B)
M <- f(X)
0 голосов
/ 04 декабря 2018

Функция cut.ppp() может взять точечный паттерн и найти, к какой мозаике мозаики относится каждая точка в паттерне.Ниже приведен код для простого примера точечного шаблона, который содержит только одну точку (0,5, 0,5).

library(spatstat)
dd <- dirichlet(cells)
plot.tess(dd, do.labels = TRUE)
xx <- ppp(.5, .5, window = Window(dd))
plot(xx, add = TRUE, col = "red", cex = 2, pch = 20)

yy <- cut(xx, dd)
yy
#> Marked planar point pattern: 1 point
#> Multitype, with levels = 
#>    1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 
#> 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42
#> window: rectangle = [0, 1] x [0, 1] units
marks(yy)
#> [1] 18
#> 42 Levels: 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 ... 42

Создано в2018-12-03 по представлению пакета (v0.2.1)

...