Python 双向插值的替代方法

Python 双向插值的替代方法,python,scipy,interpolation,Python,Scipy,Interpolation,我编写了一些代码,根据两个标准执行插值,保险金额和可扣除金额%。我当时正努力一次完成插值,所以将滤波分离。hf表包含我用来作为插值结果基础的已知数据。df表包含需要基于hf插值的开发因子的新数据 现在,我的工作是首先根据ded_数量百分比过滤每个表,然后执行插值到空数据帧中,并在每个循环后追加 我觉得这样做效率很低,有更好的方法来实现这一点,希望听到一些关于我可以做的改进的反馈。谢谢 测试数据如下 import pandas as pd from scipy import interpolate

我编写了一些代码,根据两个标准执行插值,保险金额和可扣除金额%。我当时正努力一次完成插值,所以将滤波分离。hf表包含我用来作为插值结果基础的已知数据。df表包含需要基于hf插值的开发因子的新数据

现在,我的工作是首先根据ded_数量百分比过滤每个表,然后执行插值到空数据帧中,并在每个循环后追加

我觉得这样做效率很低,有更好的方法来实现这一点,希望听到一些关于我可以做的改进的反馈。谢谢

测试数据如下

import pandas as pd
from scipy import interpolate

known_data={'AOI':[80000,100000,150000,200000,300000,80000,100000,150000,200000,300000],'Ded_amount':['2%','2%','2%','2%','2%','3%','3%','3%','3%','3%'],'factor':[0.797,0.774,0.739,0.733,0.719,0.745,0.737,0.715,0.711,0.709]}
new_data={'AOI':[85000,120000,130000,250000,310000,85000,120000,130000,250000,310000],'Ded_amount':['2%','2%','2%','2%','2%','3%','3%','3%','3%','3%']}

hf=pd.DataFrame(known_data)
df=pd.DataFrame(new_data)

deduct_fact=pd.DataFrame()
for deduct in hf['Ded_amount'].unique():
    deduct_table=hf[hf['Ded_amount']==deduct]
    aoi_table=df[df['Ded_amount']==deduct]
    x=deduct_table['AOI']
    y=deduct_table['factor']
    f=interpolate.interp1d(x,y,fill_value="extrapolate")
    xnew=aoi_table[['AOI']]
    ynew=f(xnew)
    append_frame=aoi_table
    append_frame['Factor']=ynew
    deduct_fact=deduct_fact.append(append_frame)

是的,有一种方法可以更有效地做到这一点,而不必制作一堆中间数据帧并附加它们。请查看以下代码:

from scipy import interpolate
known_data={'AOI':[80000,100000,150000,200000,300000,80000,100000,150000,200000,300000],'Ded_amount':['2%','2%','2%','2%','2%','3%','3%','3%','3%','3%'],'factor':[0.797,0.774,0.739,0.733,0.719,0.745,0.737,0.715,0.711,0.709]}
new_data={'AOI':[85000,120000,130000,250000,310000,85000,120000,130000,250000,310000],'Ded_amount':['2%','2%','2%','2%','2%','3%','3%','3%','3%','3%']}

hf=pd.DataFrame(known_data)
df=pd.DataFrame(new_data)

# Create this column now
df['Factor'] = None

# I like specifying this explicitly; easier to debug
deduction_amounts = list(hf.Ded_amount.unique())
for deduction_amount in deduction_amounts:
    # You can index a dataframe and call a column in one line
    x, y = hf[hf['Ded_amount']==deduction_amount]['AOI'], hf[hf['Ded_amount']==deduction_amount]['factor']

    f = interpolate.interp1d(x, y, fill_value="extrapolate")

    # This is the most important bit. Lambda function on the dataframe
    df['Factor'] = df.apply(lambda x: f(x['AOI']) if x['Ded_amount']==deduction_amount else x['Factor'], axis=1)

lambda函数的工作方式是: 它逐行遍历“Factor”列,并根据其他列的条件为其提供一个值

如果扣减金额匹配,它将返回df的AOI列的插值(这就是您所谓的xnew),否则它将返回相同的内容