我有两种不同仪器同时运行的气体浓度,因此它们的时间标记不同,但在相同的2小时块内(例如,O3有710个测量点,但CO2有2190):
我需要一种方法来平均或插入介于o3时间戳之间的二氧化碳浓度,因此阵列具有相同的点数,我可以生成散点图。我在这里使用python。有没有办法使用scipy或numpy插值功能?
我对编码和python很新(只有大约1个月的经验),所以感谢任何帮助。
答案 0 :(得分:1)
您可以使用scipy.interp.interp1d工具。根据您想要的插值类型,您可以将“kind ='cubic'”作为关键字参数添加到interp1d函数中。
import scipy.interp as interp
import numpy as np
f_o3 = interp.interp1d(o3_time, o3_conc)
f_co2 = interp.interp1d(co2_time, co2_conc)
new_time = np.linspace(start_time,stop_time,num_points)
new_o3 = f_o3(new_time)
new_co2 = f_co2(new_time)