如何将共享相同坐标的点数作为标记添加到 ppp 对象？

Question

在 spatstat 中，如果您有重复项或长度等于 ppp 对象中点数的整数（重数）的数字向量，则有许多函数返回逻辑向量。我找不到允许您将多重性结果作为标记直接添加到您的 ppp 对象的功能。你怎么能做到这一点？如何最终得到一个新的 ppp 对象，其中包含每对坐标中的事件计数作为标记？

另一种问法是，multiplicity.ppp 产生的结果向量的顺序是如何确定的？假设您已经在 ppp 对象 x 中测试了多重性，并基于此“multiple_x”创建了一个向量。

x <- ppp(x = mycoords[,1], y = mycoords[,2],
                   window = window, check = T)
multiple_x <- multiplicity.ppp(x)

您可以简单地使用原始坐标和新向量作为标记创建一个新的 ppp 对象吗？

x <- ppp(x = mycoords[,1], y = mycoords[,2],
                   window = window, check = T,
                   marks = multiple_x)

非常感谢。

Answer 1

您可以按照您的建议使用ppp 创建一个新的点模式。

但是使用成语marks(X) <- value 将标记分配给现有的点模式X 会容易得多。 （请参阅marks<-.ppp 的帮助）

例子：

Y <- X
marks(Y) <- multiplicity(X)

那么Y与X的模式相同，只是每个点都用一个整数标记，表示同一位置的点数。

如果您想折叠点模式X 中的重复点，只保留每个重复点的一个实例，但带有指示其多重性的标记，那么就这样做

Z <- Y[!duplicated(X)]

multiplicity.ppp 的帮助文件中给出了示例