Многовариантная сплайн-интерполяция в python / scipy?

Есть ли модуль библиотеки или другой простой способ реализовать многомерную сплайн-интерполяцию в Python?

В частности, у меня есть набор скалярных данных на регулярной трехмерной сетке, которые мне нужно интерполировать в небольшом количестве точек, разбросанных по всей области. Для двух измерений я использовал scipy.interpolate.RectBivariateSpline , и я, по сути, ищу расширение этого для трехмерных данных.

Я нашел процедуры N-мерной интерполяции: недостаточно хорош: я бы предпочел сплайны LinearNDInterpolator для гладкости, и у меня слишком много точек данных (часто более миллиона), например, для работы радиальной базовой функции.

Если кто-нибудь знает библиотеку Python, которая может это сделать, или, возможно, одну библиотеку на другом языке, которую я мог бы вызвать или перенести, я был бы очень признателен.