在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)

然而，这完全没有考虑输入数据的间隔不均匀。任何帮助都将不胜感激 :)

- andreas-h

Answer 1

看起来你可以将不均匀间隔的输入数据插值到一个均匀的网格上。如果频谱应该是平滑的，那么只要选择与高斯窗口使用的相同网格，就应该没问题。

- Greg