使用 wavio Python 3.5 读取 24 位 WAV

Question

我正在使用WarrenWeckesser的 wavio，因为我需要在 python 中读取 24 位 wav 文件。 我拥有的 wav 文件是由一些仪器生成的，我试图在没有任何标准化或缩放的情况下获取原始值。

在 wavio 模块中，完成工作的代码是这样的：

    if sampwidth == 3:
        a = _np.empty((num_samples, nchannels, 4), dtype=_np.uint8)
        raw_bytes = _np.fromstring(data, dtype=_np.uint8)
        a[:, :, :sampwidth] = raw_bytes.reshape(-1, nchannels, sampwidth)
        a[:, :, sampwidth:] = (a[:, :, sampwidth - 1:sampwidth] >> 7) * 255
        result = a.view('<i4').reshape(a.shape[:-1]

有人能解释一下它实际上在做什么吗（我对 numpy 和数组切片比较陌生）。 我明白其中的大部分，但我不明白这里发生了什么：

    a[:, :, sampwidth:] = (a[:, :, sampwidth - 1:sampwidth] >> 7) * 255

在我的例子中，它进行了从 24 位到 32 位的转换，但我无法确定它是在缩放数据，还是只是在不更改任何原始值的情况下填充它。

Answer 1

的形状a是(num_samples, nchannels, 4)和sampwidth == 3 ，从而使线相同

a[:, :, 3:] = (a[:, :, 2:3] >> 7) * 255

这与

a[:, :, 3] = (a[:, :, 2] >> 7) * 255

我们可以对外部的两个循环进行去向量化：

for i in range(num_samples):
    for j in range(nchannels):
        a[i, j, 3] = (a[i, j, 2] >> 7) * 255

a的_np.uint8是_np.uint8 ，所以a[...] >> 7只能在值 <128 时给出 0，或者在 ≥128 时给出 1，所以上面变成：

for i in range(num_samples):
    for j in range(nchannels):
        v = a[i, j, 2]
        a[i, j, 3] = 255 if v >= 128 else 0

如果数据是 24 位小端整数，这相当于对 32 位进行符号扩展。