如何绘制带孔的曲面 z=f(x,y)？

Question

我想画一个类似于左边的图（目标是说明插值：右边的图）。

我用受限的 Delaunay 镶嵌（使用 delaunay 包）来绘制这个图，应用于在带孔的网格上给定的点

(x,y)

，然后我构建了一个 3D rgl 网格。我想知道是否存在更直接的方法。我试了

rgl::perps3d.deldir

但是填坑了，即使可以去掉这部分，这种方式也实在是太慢了

这是带有孔和

z

值的网格：

x <- y <- seq(-3, 3, by = 0.1)
Grid <- expand.grid(X = x, Y = y)
toremove <- with(Grid, X >= -0.5 & X <= 1 & Y >= -0.5 & Y <= 1)
Grid <- Grid[!toremove, ]
z <- with(Grid, dt(X, df = 1) * dt(Y, df = 1))

我还希望能够在第一个图（黑色插值部分）中添加另一个图。