使用R

时间:2018-12-03 16:43:11

标签: r delaunay spatstat dirichlet

因此,我可以取点并使用Rdeldirspatstat::dirichlet来查找这些点的dirichlet镶嵌。

现在我没有一个点在集合中,我想知道构成dirichlet磁贴的点的索引,我的非设定点位于该点内部。我可以通过了解磁贴标签(或索引)到达那里。

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

2 个答案:

答案 0 :(得分: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创建

答案 1 :(得分:0)

如果X是点图案,而B是棋盘格,则

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

返回一个因子(分类值的向量),该因子标识哪个图块包含X的每个点。或者,

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

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