waveファイルを波形にしたい。

Question

waveファイルを、x軸を時間にした波形で表したいのですが、録音した時間以上の数値が表示されます。7秒録音したwaveファイルを以下のコードで波形表示したところ、波のx軸（時間）は７前後のはずが10以上になってしまいます。

書いたコードは以下の通りです。

import matplotlib.pyplot as plt
import wave
import numpy as np

file = "filename" 

wf = wave.open(file, "rb")
data = np.frombuffer(wf.readframes(wf.getnframes()), dtype='int16')
RATE = wf.getframerate()#サンプリングレート
wf.close()

x = np.arange(data.shape[0]) / RATE
plt.plot(x, data)
plt.show()

何故録音時間と同じ時間が表示されないのでしょうか。
ご教授お願いいたします。

追記　.wavファイルをアップロードしました。

la.wav,https://www.axfc.net/u/3982348

Comments

[can110]

当方環境にてhttp://www.non-fiction.jp/2015/08/17/sin_wave/で生成した.wavファイルを用いて質問コードを実行しましたが、時間の範囲は0～5と正しくplot表示されました。
.wavファイルを提示すると回答得られやすくなるかと思います。

[yuuki_yui_]

指摘ありがとうございます。.wavファイルの提示ですが、質問内容の場所には音声ファイルをアップロードできないようで、どのように提示すれば良いでしょうか。

[can110]

githubが無難かと。アカウントなければどこか適当なアップローダでしょうか。

Answer 1

提示.wavのサンプル幅は一般的な16bitではなく24bitになっています。
その場合は以下のように読込処理を行う必要があります。
詳細はPythonでのwavファイル操作を参照ください。

import matplotlib.pyplot as plt
import wave
import numpy as np
from struct import unpack

file = "la.wav" 
wf = wave.open(file, "rb")

samplewidth = wf.getsampwidth()
buf = wf.readframes(wf.getnframes())

# サンプル幅により必要な処理が異なる！
if samplewidth == 2:  # 16bit
    data = np.frombuffer(buf, dtype='int16')
elif samplewidth == 3:# 24bit
    read_frames = wf.getnframes()
    nbyte = samplewidth
    data = [unpack("<i",
            bytearray([0]) + buf[nbyte * idx: nbyte * (idx + 1)])[0]
            for idx in range(read_frames)]
    data = np.array(data, dtype='int16')
elif samplewidth == 4:# 32bit
    data = np.frombuffer(buf, dtype='int32')

RATE = wf.getframerate()#サンプリングレート

wf.close()

x = np.arange(data.shape[0]) / RATE
plt.plot(x, data)
plt.show()

Comments

[yuuki_yui_]

サンプル幅によって処理が異なるんですね。詳細のページも紹介してくださりありがとうございました。音声処理についてもっと勉強します。