python - 在 python scipy 版本 1.1 中计算信噪比

Question

我在网上查了一下，似乎scipy.stats中的signaltonoise比率函数已被弃用，并且在1.1版本中不可用。 scipy 包中是否还有其他等效方法，因为我无法在网上找到它。

如果不是 scipy 那么是否有其他推荐的库用于此类计算？

Answer 1

如scipy issue #609 on github中所示，signaltonoise 函数

[...] is not useful except for backwards compatibility. The reason is that there's a Matlab signal-to-noise function http://www.mathworks.com/help/signal/ref/snr.html which means something different. This is not good, because scipy clones the Matlab interface of other signal-related functions, and this incompatibility apparently has no offsetting benefit.

如果您确实需要此函数以实现向后兼容性，可以在 history of scipy repository 中找到简短的实现。 as(在此处转载，没有文档注释， license ):

def signaltonoise(a, axis=0, ddof=0):
    a = np.asanyarray(a)
    m = a.mean(axis)
    sd = a.std(axis=axis, ddof=ddof)
    return np.where(sd == 0, 0, m/sd)