The example code here for rgl does almost what I want. However my (x,y) are not in a rectangle but in a triangle that is half of the rectangle (x=0..1, y=0..1 and z is missing if x+y>1). My input
computational geometry - Efficiently determining if 3D points are within a surface composed of polygons - Mathematica Stack Exchange
plotting - Creating Optimal Smoothed Contour-Type Plots (X,Y,Z) Using Mathematica 8.04 - Mathematica Stack Exchange
Free variable mass nonlocal systems, jerks, and snaps, and their implications in rotating fluids in rockets
Making Contour Maps in R
r - Adding Contour Lines to 3D Plots - Stack Overflow
swiftui - How to draw surface contour plot (Swift language) - Stack Overflow
Plot 3D data in R - Stack Overflow
plotting - Create smooth colour map from given array for ListPlot3D - Mathematica Stack Exchange
Use Python to plot Surface graphs of irregular Datasets
Add contour plot to surface plot in R - Stack Overflow