在R中绘制样方图作为强度图

zu0ti5jz  于 2023-06-03  发布在  其他
关注(0)|答案(1)|浏览(201)

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

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

其中xhomicide是先前定义的 Dataframe 。
图1.每个样方的点数图:

图2示出了图1中所示的样方计数的相应强度图。图2是使用以下代码创建的:#添加每个样方图的强度(intensity(q,image=TRUE),main=NULL,las=1)
图2.每个样方的点数的强度图。

手工计算表明,图1和图2都显示了样方计数方法的正确结果。
请注意,图1和图2是在R v4.1.2下创建的。
我将R更新到v4.3.0并安装了spatstat包的当前版本。当在新版本的R下运行上述代码时,我获得了与图1相同的结果(R版本之间每个样方的计数值相同)。但是,当我运行intensity()函数时,结果如图3所示。
图3.使用R v4.3.0绘制的每个样方点数的强度图:

既然两个版本的R的代码是相同的,那么问题是出在新版本的R上,还是spatstat发生了什么变化?我很想知道这个问题是否与新版本的R有关,当前版本的spatstat,或者我的代码不再满足R或spatstat的要求。我找不到任何其他的评论是什么我是造成这个问题。

oaxa6hgo

oaxa6hgo1#

这是一个bug。谢谢你提醒我们注意。
此问题已在the github repository上提供的spatstat.geom开发版本(版本3.2-1.001)中得到修复

相关问题