Z2 = np.array([0, 1, 0.5, 0])
p2d = np.vstack([X2, Y2]).T
d2d = Delaunay(p2d)
fig = mlab.figure(1, bgcolor=(1, 0.7, 1), fgcolor=(0.5, 0.5, 0.5))
tmesh = mlab.triangular_mesh(X2, Y2, Z2, d2d.vertices, scalars=Y2, colormap
    ='jet')
mlab.outline(extent=(0, 1, 0, 1, 0, 1))
mlab.axes(extent=(0, 1, 0, 1, 0, 1))
mlab.show()
Mayavi How to do Delaunay Triangulation with xz data instead of xy data