Delaunay triangulation in best fit plane

Hi I have a set of points scattered in 3D x,y, z, with an associated scalar property f. But these points are very nearly planar, and I want the triangulation of f in a “best fit” plane, which I can then import as a VTK file into ParaView using WriteVTK for viewing. Does anyone know of a Julia package / have examples on how to do this?