owin - 如何绘制等高线图并显示每个等高线级别的标记平均值（与 x，y 相关）

Question

我在 spatstat 中有一个由点模式 X 的强度函数生成的等高线图（如“树木的位置”）。该点模式中的每个 x,y 坐标都标有相应的第三个向量（如“树的直径”）。-->cf image（当然代表树的竖线可以省略）

我想用不同的颜色显示轮廓每个级别的标记（直径）的平均值。建议？谢谢！

Answer 1

您实际上是在要求一种非参数回归。

这是使用该函数进行的快速而简单的计算，rhohat并在数据集上进行了演示longleaf。

首先计算强度估计：Z <- density(longleaf)产生图像Z。接下来在对命令Z的调用中将其视为协变量：rhohat

f <- rhohat(unmark(longleaf), Z)

和

g <- rhohat(unmark(longleaf), Z, weights=marks(longleaf)).

现在取比率，h <- eval.fv(g/f)并绘制它，plot(h)。这显示了估计的平均树木直径作为森林密度的函数。要将这个函数 h 应用于 Z 的原始轮廓，您首先要转换h为一个真正的函数，H <- as.function(h)然后再评估hZ <- eval.im(H(Z))和最后plot(hZ)。