2D-интерполяция большой нерегулярной сетки в регулярную сетку

У меня есть сетка нерегулярных данных размером 2048x2048 zi = f (xi, yi) , которые, по сути, являются тремя независимыми наборами 2048 вещественных значения. Мне нужно плавно интерполировать (возможно, бикубическим сплайном) это в регулярную сетку wi = f (ui, vi) , где ui и vi - целые значения из От 0 до 2047.

Я пробовал данные griddata, которые, кажется, хорошо работают с изображениями меньше 1000x1000, но взрываются, когда вы доходите до 1500x1500 (очевидно, ошибки qhull памяти для сетки Делоне). Я просмотрел некоторые из функций ndimage , а именно geometry_transform , RectBivariateSpline и map_coordinates , но все они, похоже, принимают на вход регуляризованные данные . Я мог что-то упустить и просто реализовать это неправильно!

Я пытаюсь использовать Python / SciPy для выполнения того же сценария Matlab, который я делал, используя tformarray и makeresampler . Есть предложения относительно того, какую функцию я могу использовать для обработки этого большого набора данных? Спасибо!