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

Question

我有一个测量光谱，即一个一维数组spec。spec.shape = (n, )各个点对应于不均匀间隔的波长，存储在 1d 数组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

听起来您可以将不均匀间隔的输入数据插入到统一的网格中。如果光谱应该是平滑的，那么只要您选择与高斯窗口使用的网格相同的网格，就可以了。