所以我可以拿分并使用R
图书馆deldir
或spatstat::dirichlet
找到这些点的狄利克雷镶嵌。
现在我有一个点不在集合中,我想知道形成狄利克雷图块的点的索引,我的非设置点在内部。我可以通过了解磁贴标签(或索引)到达那里。
是否有任何库或方法可以做到这一点?我在想spatstat,但还没有找到东西。
该函数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 包(v0.2.1)于 2018 年 12 月 3 日创建
如果X
是一个点模式并且B
是一个镶嵌,那么
M <- marks(cut(X, B))
返回一个因子(分类值的向量),标识哪个图块包含 的每个点X
。或者,
M <- tileindex(X$x, X$y, B)
或者
f <- as.function(B)
M <- f(X)