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

Question

我有一个测量光谱，即一维数组spec 和spec.shape = (n, )。各个点对应于不均匀间隔的波长，存储在一维数组wl 和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 放到一个统一的网格上。如果光谱应该是平滑的，那么只要您选择与高斯窗口使用的网格相同的网格，就可以了。