在scipy中卷积不均匀间隔的向量

Question

我有一个测量光谱，即带有spec的1d数组spec.shape = (n, )。各个点对应于不均匀间隔的波长，存储在具有wl的1d阵列wl.shape = (n, )中。现在我需要将此频谱与预定义FWHM的高斯函数进行卷积。

我怎么能用numpy / scipy做到这一点？

目前，我有这个：

fwhm = .22
gaussian_window = scipy.signal.gaussian(spec.size, std=fwhm/2.3548201)
result = scipy.signal.fftconvolve(spec, gaussian_window)

但是，这并没有考虑输入数据的不均匀间距。

非常感谢任何帮助：）

Answer 1

听起来你可以interpolate将不均匀间隔的输入数据{{3}}放到统一网格上。如果光谱应该是光滑的，那么如果您选择与高斯窗口使用的网格相同的网格，您应该没问题。