使用 R 测试 Dirichlet tile 的 in-tile

Question

所以我可以取点并使用 R 库 deldir 或 spatstat::dirichlet 来找到这些点的狄利克雷镶嵌。

现在我有一个点不在集合中，我想知道形成 dirichlet tile 的点的索引，我的 not-in-set-point 在其内部。我可以通过了解磁贴标签（或索引）到达那里。

是否有任何库或方法可以做到这一点？我在想 spatstat，但还没有找到任何东西。

Answer 1

函数 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

^{由 reprex package (v0.2.1)}

创建于 2018-12-03

Answer 2

如果X是点图案，B是镶嵌，那么

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

returns 一个因子（分类值的向量）标识哪个图块包含 X 的每个点。或者，

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

或

f <- as.function(B)
M <- f(X)

使用 R 测试 Dirichlet tile 的 in-tile

test for in-tile for Dirichlet tile, using R

r

delaunay

dirichlet

spatstat