r - eps 在密度.ppp() 中如何使用？

Question

我正在使用 density.ppp() 来获取一系列坐标的强度函数。数据是投影的，以米为单位。我试图估计点数/公里，所以在查看文档后，我使用 eps = 1000 将输出图像的分辨率设置为 1km2。

但是 eps 是否仅指估计值的位置？或者强度函数估计的面积单位？

我是否仍需要将值重新缩放 1000x1000 才能获得 km2？

Answer 1

参数eps通常可以是一个二维向量，给出 x 和 y 坐标中的像素大小。如果您只给出一个值（这可能是最常见的用法），则此大小将在两个方向上使用。因此，您的使用应该为您提供 1000 m x 1000 m 像素的图像，并具有强度值。此强度是每单位面积，在您的情况下，这意味着强度值是 1 m x 1 m 区域中的预期点数。如果您想将整个分析重新调整为以 km 为单位，我建议您首先使用rescale.pppto 转换为 km，然后使用eps=1in density.ppp。如果您的数据在一个ppp对象中X，您应该这样做：

Y <- rescale(X, 1000, unitname = "km")
density(Y, eps = 1)

我不记得调用了哪些内部函数density.ppp来传递eps给 as.mask，但如果你真的需要知道，我可以进一步调查。