Python Интерполяция функции двух переменных

У менять имеется набор данных типа

x = [-3.1234, -3.1134, -3.1021, ...]
y = [-18.0123, -18.0223, -18.0236, ...]
z = [-1.432, -1.657, -1.876, ...]

Это как пример, точек около 800, со связью z = f(x,y). Функция не известна, промежутки неравномерны.

Моя цель добавить промежуточные точки, между заданными, то есть, как я понимаю, провести интерполяцию. Но надо сделать это так, чтобы функция получилась более-менее гладкой. То есть реализация через билинейную интерполяцию не пойдет.

Есть ли какие-то дополнительные возможности в python для реализации подобной задачи (scipy например)? Я видел много материала по поводу одномерного случая. Так же пробовал генерировать нейронкой, но там и вовсе нерабочий код получается или только билинейная.