python - python中的对数插值

Question

使用 numpy.interp 我能够计算在离散数据点处具有给定值的函数的一维分段线性插值。

返回我的对数插值是不是类似的函数？

Answer 1

在过去，我只是将普通插值包装在对数空间中，即

def log_interp(zz, xx, yy):
    logz = np.log10(zz)
    logx = np.log10(xx)
    logy = np.log10(yy)
    return np.power(10.0, np.interp(logz, logx, logy))

就我个人而言，我更喜欢 scipy interpolation functions (正如@mylesgallagher 提到的)，例如:

import scipy as sp
import scipy.interpolate

def log_interp1d(xx, yy, kind='linear'):
    logx = np.log10(xx)
    logy = np.log10(yy)
    lin_interp = sp.interpolate.interp1d(logx, logy, kind=kind)
    log_interp = lambda zz: np.power(10.0, lin_interp(np.log10(zz)))
    return log_interp

然后您可以将其作为任意值的函数调用。