使用 R 外推空间点模式

Question

我想知道是否有一个内置函数可以在 R 中的“父”窗口之外推断点模式。例如，让我们生成一个空间点模式“X”：

require(spatstat)
X <- ppp(runif(200), runif(200), 
c(0,1), c(0,1), unitname=c("metre","metre"))

让我们重新采样数据：

a <- quadratresample(X, nx=25, ny=5, replace=F, nsamples = 1)

但新点是在同一区域/空间窗口内生成的

> a
 planar point pattern: 200 points 
window: rectangle = [0, 1] x [0, 1] metres  

我的问题是：如何在比原始窗口大的新窗口中重新采样 200 个点（1 x 1 m）；换句话说，我将如何在保持相同的重采样密度的同时将这小组 200 个空间点外推到更大的范围；假设我想在 5 x 5 m 范围内查看总共 1,000 个数据点？

Answer 1

这很容易，因为spatstat 为我们提供了所有正确的工具。您目前有一个 1x1 网格。您需要一个由 25 个 1x1 网格构成的 5x5 网格。我们可以使用nsamples 参数对这些网格的点进行采样：

a <- quadratresample(X, nx = 25, ny = 5, replace = F, nsamples = 25)

现在我们有一个包含 25 个ppps 的列表。正如您所指出的，所有这些都将在同一个 1x1 窗口中。为了将它们变成一个网格，我们将它们适当地移动，从 x 偏移 0 到 4 个单位，在 y 偏移中从 0 到 4 个单位：

for (i in seq_along(a)) {
    a[[i]] = shift(a[[i]], vec = c((i - 1) %% 5, (i - 1) %/% 5))
}

要将它们组合起来，请使用superimpose：

b = superimpose(a)

这会在 5x5 窗口中提供一个具有 200 * 25 = 5000 个点的 ppp 对象，这会保留原始的每单位平方 200 个点。