R - 如何从插值等高线图中获取数据

Question

我使用“R”已经有一段时间了。我的实际问题位于 3d 绘图和插值区域。

我有一个包含 3d 数据的数据框。假设它是一个区域和相应的高程。我只有几个兴趣点并使用插值功能。 3d 图为我提供了热图和等高线。

现在我对恒定高度的轮廓线感兴趣，我想得到它的 X 和 Y 坐标。例如，实际的插值数据并不包含所有正好为零的点。

test = data.frame(x = rep(c(1,2,3), 3), y = c(rep(1,3),rep(2,3),rep(3,3)), z = rnorm(9, 0, 1))
mat = as.matrix(interp(test$x, test$y, test$z, nx = 10, ny = 10, duplicate = "strip"))[[3]]

image(mat)
contour(mat, add = TRUE)
contour(mat, add = TRUE, levels = 0, lwd = 3)

我尝试了几个拟合函数，如 lm、loess、lowess、nls 等。但它们都不适合我。要么是二维数据，要么是给我与等高线完全不同的数据点。

有什么想法吗？

最好的问候，

拉勒

Answer 1

好的，找了几天，终于找到了解决办法。

“grDevices”包中的“contourLines”函数返回所有内插线的列表。