python 拉格朗日插值法处理异常数据

问题,run了以后错误提示如下图：

问题相关代码（拉格朗日插值法代码）

import pandas as pd
import numpy as np
# 拉格朗日插值
from scipy.interpolate import lagrange  #scipy.interpolate是内置工具包

def ploy (s,n,k=5):   #自定义插值函数
    y=s[list(range(n-k,n))+list(range(n+1,n+1+k))] #取出插值位置前后k个数据
    y=y[y.notnull()]  #剔除空值
    return lagrange(y.index,list(y))(n)

traj = pd.read_csv('DATASET-C-A.csv', header=None, usecols=[2,3,4]).iloc[:15]   #usecols表示选取那几列
traj.columns = ['timestamp', 'lon', 'lat']
traj['time_interval']=traj['timestamp']-traj['timestamp'].shift(1)
traj['time_interval'] = traj['timestamp'] - traj['timestamp'].shift(1)   #time_interval “时间间隔”   .shift(1)意思是把数据向下移动1位
index = traj[traj['time_interval'] >=6].index.to_list()
for i in index:
    timestamp = traj['timestamp'].loc[i-1] + 3
    insertRow = pd.DataFrame([[np.nan, np.nan, timestamp]], columns=['lon', 'lat', 'timestamp'])
    traj = pd.concat([traj[:i], insertRow, traj[i:]], ignore_index=True)
    traj['lon'][i]=ploy(traj['lon'],i)
    traj['lat'][i]=ploy(traj['lat'],i)
traj = traj.drop(['time_interval'], axis=1)

运行结果及报错内容

结果报错，求各位指导，哪里出了问题

错误提示不支持字符串数据减浮点数据

第11行
traj = pd.read_csv('DATASET-A.csv',header=None,usecols=[2,3,4],skiprows=[0]).iloc[:15]
加skiprows=[0]或者header=0，从数据第一行开始读即可，解决问题请联系我，我也有此类问题