在三角形上使用 3D 等值线图绘制曲面

Question

rgl 的示例代码 几乎可以满足我的要求。然而，我的 (x,y) 不在矩形中，而是在矩形一半的三角形中（x=0..1，y=0..1，如果 x+y>1，则缺少 z）。

我的输入数据集采用以下格式：
x1 x2 x3 x4 x5 x6 z

，我想可视化（x4，x5，z）和（x2，x3，x6）等。如果如何在 R 中执行此操作，请提供任何帮助，我是初学者。

可以将最终绘图导出为 PDF、EPS 或 SVG 吗？

更新：好的，我看到我可以导出为 EPS。

更新2：多亏了卡尔，我几乎得到了我想要的东西。

唯一的问题是表面的边缘之一是锯齿状的。在齿之间，其中x+y≤4，应当保留表面，使得锯齿状边缘变得像其他边缘一样。我该怎么做？

R代码如下，输入数据在这里

library(rgl)
mat = matrix(scan("bpexamp.out"),ncol=9,byrow=T)
N <- 4
x <- c(0:N)
y <- c(0:N)
z <- mat[,9]
zlim <- range(y)
zlen <- zlim[2] - zlim[1] + 1
colorlut <- terrain.colors(zlen,alpha=0) # height color lookup table
col <- colorlut[ z-zlim[1]+1 ] # assign colors to heights for each point
open3d()
aspect3d(1,1,0.1)
surface3d(x, y, z, col)
axes3d()

Answer 1

这是一条评论 - 但我不知道如何将图像放入评论中，因此将其格式化为答案。我跑了你的东西并得到了这张图表：
我跑了你的东西并得到了这个图表：

那么，这里缺少或不准确的是什么？我们也许可以为您“修复”它......