我的目标是根据地形高程数据(pp3的z轴)创建一个im对象，类似于bei.extra($elev和$grad).我首先从我的数据集创建了一个pp3对象:
x <- elev_test$x
y <- elev_test$y
z <- elev_test$z
X <- pp3(x, y, z, c(0, 102), c(0, 152), c(90, 102))
plot(X)
然后，我试图将pp3对象强制转换为类im，但as.im抛出了一个错误:
> Elev_im <- as.im(X)
Error in as.im.default(X) : Can't convert X to a pixel image
如何获得所需的im件物品？
像as.im这样的转换函数旨在将数据从一种格式转换为另一种表示相同信息的格式.
您正在try 从(x,y,z)个采样点的样本中推断曲面的空间高度.这不是格式之间的转换；这是某种内插.
我建议您使用位置的(x,y)个坐标创建一个二维点图案，并使用z个坐标作为该图案的标记.然后将这些印记打磨光滑，以获得表面.
类似于:
X <- ppp(elev_test$x, elev_test$y, c(0, 102), c(0, 152), marks=elev_test$z)
Z <- Smooth(X)
然后，您可以try 不同的参数到Smooth.ppp，以获得最佳的插值式.
另一种可能是使用nnmark而不是Smooth.ppp来获得由小平面组成的曲面(如果您的点来自稻田，这将更好).你得自己做这些决定，因为电脑不够聪明.