r - 在 R 中将样方计数绘制为强度图

Question

我有一个城市犯罪点数据集。通过 spatstat 包使用 quadratcount()，我绘制了该数据集的样方计数。图 1 显示了每个样方的数据点计数图像，是通过以下代码创建的:

#Quadrat Count Map 
q <- quadratcount(xhomicide, 4, 8) 
plot(q) 

其中 xhomicide 是先前定义的数据框。

图 1. 每样方点数图:

图 2 显示了图 1 中所示样方计数的相应强度图。图 2 是使用以下代码创建的: #添加每个样方的强度 绘图(强度(q，图像= TRUE)，主要= NULL，las = 1)

图2.每样方点数的强度图:

手工计算表明，图1和图2均显示了样方计数法的正确结果。

请注意，图 1 和图 2 是在 R v4.1.2 下创建的。

我将 R 更新到 v4.3.0 并安装了当前版本的 spatstat 软件包。在新版本的 R 下运行上述代码时，我获得了与图 1 相同的结果(每个样方的计数值在 R 版本之间是相同的)。然而，当我运行强度()函数时，结果如图 3 所示。

图 3. 使用 R v4.3.0 的每个样方点数的强度图:

由于两个版本的 R 之间的代码相同，是新版本 R 的问题还是 spatstat 发生了一些变化？我很想知道问题是否在于新版本的 R、当前版本的 spatstat，或者我的代码不再满足 R 或 spatstat 的要求。我找不到任何其他评论来解释我造成这个问题的原因。

Answer 1

这是一个错误。 感谢您引起我们的注意。

它已在 spatstat.geom 的开发版本(版本3.2-1.001)中修复，可在 the github repository 获取。