Интерполяция на DataFrame в pandas

У меня есть DataFrame, скажем, поверхность волатильности с индексом как временем и столбцом как страйк. Как выполнить двумерную интерполяцию? Я могу reindex, но как мне быть с NaN? Я знаю, что можем fillna(method='pad'), но это даже не линейная интерполяция. Есть ли способ подключить наш собственный метод для интерполяции?