如何在density.ppp()中使用eps?

时间:2015-07-21 03:23:55

标签: r resolution kernel-density spatstat

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

但是,eps仅指的是估算值的位置吗?或强度函数估计的区域单位?

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

1 个答案:

答案 0 :(得分:0)

参数eps通常可以是二维向量,给出x和y坐标中的像素大小。如果您只提供一个值(这可能是最常见的用法),则会在两个方向上使用此大小。因此,您的使用应该为您提供1000米×1000米像素和强度值的图像。该强度是每单位面积,在您的情况下,这意味着强度值是1米乘1米区域的预期点数。如果您想将整个分析重新调整为km,我建议您先使用rescale.ppp转换为km,然后在eps=1中使用density.ppp。如果您的数据位于ppp对象X中,则应执行以下操作:

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

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