使用Stienen集和扩张隔离组中的单例

Question

我想使用Steiner集（function count_transactions($array_transactions) { if(isset($array_transactions)){ $count = 0; foreach($array_transactions['data']['txs'] as $value){ echo "<p>Transactions id: " . $value['txid'] . "</p>"; $count++; } echo $count . "transactions found!"; } } ）形成点组。

有些单例需要删除，我正在考虑为此使用Stienen集（spatstat::dilation）：对于单例Stienen> Steiner。

这怎么办？

spatstat::stienenSet

Answer 1

您可以使用connected.ppp查找连接的组件并确定 这样的单例（抱歉，没有注释的快速代码）：

library(spatstat)
## Points within distance R are connected.
## This should happen when two discs of radius r=0.03 touch each other,
## i.e., R=r+r=0.06:
cc <- connected(redwood, R=0.06)
s <- split(cc)
np <- sapply(s, npoints)
X <- s[np==1]
X <- unmark(superimpose(X))
Y <- s[np>1]
Y <- unmark(superimpose(Y))
plot(dilation(X, 0.03), main = "", col = "orange")
plot(dilation(Y, 0.03), add = TRUE)